Source author record

Eric B. Isaacs

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

8works
4topics
4close collaborators

Published work

8 publications

preprint2022arXiv

Assessing the accuracy of compound formation energies with quantum Monte Carlo

Accurately predicting the formation energy of a compound, which describes its thermodynamic stability, is a key challenge in materials physics. Here, we employ many-body quantum Monte Carlo (QMC) with single-reference trial functions to compute the formation energy of two electronically disparate compounds, the intermetallic VPt2_2 and the semiconductor CuI, for which standard density functional theory (DFT) predictions using both the Perdew-Burke Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) density functional approximations deviate markedly from available experimental values. For VPt2_2, we find an agreement between QMC, SCAN, and PBE0 estimates, which therefore remain in disagreement with the much less exothermic experimental value. For CuI, the QMC result agrees with neither SCAN nor PBE pointing towards DFT exchange-correlation biases, likely related to the localized Cu 3d3d electrons. Compared to the behavior of some density functional approximations within DFT, spin-averaged QMC exhibits a smaller but still appreciable deviation when compared to experiment. The QMC result is slightly improved by incorporating spin-orbit corrections for CuI and solid I2_2, so that experiment and theory are brought into imperfect but reasonable agreement within about 120~meV/atom.

preprint2022arXiv

Testing the r2^2SCAN density functional for the thermodynamic stability of solids with and without a van der Waals correction

A central aim of materials discovery is an accurate and numerically reliable description of thermodynamic properties, such as the enthalpies of formation and decomposition. The r2^2SCAN revision of the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (meta-GGA) balances numerical stability with high general accuracy. To assess the r2^2SCAN description of solid-state thermodynamics, we evaluate the formation and decomposition enthalpies, equilibrium volumes, and fundamental bandgaps of more than 1,000 solids using r2^2SCAN, SCAN, and PBE, as well as two dispersion-corrected variants, SCAN+rVV10 and r2^2SCAN+rVV10. We show that r2^2SCAN achieves accuracy comparable to SCAN and often improves upon SCAN's already excellent accuracy. Whereas SCAN+rVV10 is often observed to worsen the formation enthalpies of SCAN, and makes no substantial correction to SCAN's cell volume predictions, r2^2SCAN+rVV10 predicts marginally less-accurate formation enthalpies than r2^2SCAN, and slightly more-accurate cell volumes than r2^2SCAN. The average absolute errors in predicted formation enthalpies are found to decrease by a factor of 1.5 to 2.5 from the GGA level to the meta-GGA level. Smaller decreases in error are observed for decomposition enthalpies. For formation enthalpies r2^2SCAN improves over SCAN for intermetallic systems. For a few classes of systems -- transition metals, intermetallics, weakly-bound solids, and enthalpies of decomposition into compounds -- GGAs are comparable to meta-GGAs. In total, r2^2SCAN and r2^2SCAN+rVV10 can be recommended as stable, general-purpose meta-GGAs for materials discovery.

preprint2020arXiv

Compositional phase stability of correlated electron materials within DFT+DMFT

Predicting the compositional phase stability of strongly correlated electron materials is an outstanding challenge in condensed matter physics, requiring precise computations of total energies. In this work, we employ the density functional theory plus dynamical mean-field theory (DFT+DMFT) formalism to address local correlations due to transition metal d electrons on compositional phase stability in the prototype rechargeable battery cathode material LixCoO2, and detailed comparisons are made with the simpler DFT+U approach (i.e. the Hartree-Fock solution of the DMFT impurity problem). Local interactions are found to strongly impact the energetics of the band insulator LiCoO2, most significantly via the Eg orbitals, which are partially occupied via hybridization with O p states. We find CoO2 and Li1/2CoO2 to be moderately correlated Fermi liquids with quasiparticle weights of 0.6-0.8 for the T2g states, which are most impacted by the interactions. As compared to DFT+U, DFT+DMFT considerably dampens the increase in total energy as U is increased, which indicates that dynamical correlations are important to describe this class of materials despite the relatively modest quasiparticle weights. Unlike DFT+U, which can incorrectly drive LixCoO2 towards spurious phase separating or charge ordered states, DFT+DMFT correctly captures the system's phase stability and does not exhibit a strong charge ordering tendency. Most importantly, the error within DFT+U varies strongly as the composition changes, challenging the common practice of artificially tuning U within DFT+U to compensate the errors of Hartree-Fock. DFT+DMFT predicts the average intercalation voltage decreases relative to DFT, opposite to the result of DFT+U, which would yield favorable agreement with experiment in conjunction with the overprediction of the voltage by the strongly constrained and appropriately normed DFT functional.

preprint2020arXiv

Inverse Design of Ultralow Lattice Thermal Conductivity Materials Via Lone Pair Cation Coordination Environment

The presence of lone pair (LP) electrons is strongly associated with the disruption of lattice heat transport, which is a critical component of strategies to achieve efficient thermoelectric energy conversion. By exploiting an empirical relationship between lattice thermal conductivity κLκ_L and the bond angles of pnictogen group LP cation coordination environments, we develop an inverse design strategy based on a materials database screening to identify chalcogenide materials with ultralow κLκ_L for thermoelectrics. Screening the ∼\sim 635,000 real and hypothetical inorganic crystals of the Open Quantum Materials Database based on the constituent elements, nominal electron counting, LP cation coordination environment, and synthesizability, we identify 189 compounds expected to exhibit ultralow κLκ_L. As a validation, we explicitly compute the lattice dynamical properties of two of the compounds (Cu2_2AgBiPbS4_4 and MnTl2_2As2_2S5_5) using first-principles calculations and successfully find both achieve ultralow κLκ_L values at room temperature of ∼\sim 0.3--0.4 W/(m⋅\cdotK) corresponding to the amorphous limit. Our data-driven approach provides promising candidates for thermoelectric materials and opens new avenues for the design of phononic properties of materials.

preprint2020arXiv

Prediction of Li intercalation voltages in rechargeable battery cathode materials: effects of exchange-correlation functional, van der Waals interactions, and Hubbard UU

Quantitative predictions of the Li intercalation voltage and of the electronic properties of rechargeable battery cathode materials are a substantial challenge for first-principles theory due to the possibility of (1) strong correlations associated with localized transition metal dd electrons and (2) significant van der Waals (vdW) interactions in layered systems, both of which are not accurately captured by standard approximations to density functional theory (DFT). Here, we perform a systematic benchmark of electronic structure methods based on the widely-used generalized gradient approximation of Perdew, Burke, and Ernzerhof (PBE) and the new strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation for battery cathode materials. Studying layered Lix_xTiS2_2, Lix_xNiO2_2, and Lix_xCoO2_2, olivine Lix_xFePO4_4, and spinel Lix_xMn2_2O4_4, we compute the voltage, crystal structure, and electronic structure with and without extensions to incorporate on-site Hubbard interactions and vdW interactions. Within pure DFT (i.e., without corrections for on-site Hubbard interactions), SCAN is a significant improvement over PBE for describing cathode materials, decreasing the mean absolute voltage error by more than 50%. Although explicit vdW interactions are not critical and in cases even detrimental when applied in conjunction with SCAN, Hubbard UU corrections are still in general necessary to achieve reasonable agreement with experiment. We show that no single method considered here can accurately describe the voltage and overall structural, electronic, and magnetic properties (i.e., errors no more than 5% for voltage, volume, band gap, and magnetic moments) of battery cathode materials, motivating a strong need for improved electronic structure approaches for such systems.

preprint2019arXiv

Electronic Structure and Phase Stability of Yb-filled CoSb3_3 Skutterudite Thermoelectrics from First Principles

Filling the large voids in the crystal structure of the skutterudite CoSb3_3 with rattler atoms RR provides an avenue for both increasing carrier concentration and disrupting lattice heat transport, leading to impressive thermoelectric performance. While the influence of RR on the lattice dynamics of skutterudite materials has been well studied, the phase stability of RR-filled skutterudite materials and the influence of the presence and ordering of RR on the electronic structure remain unclear. Here, focusing on the Yb-filled skutterudite Ybx_xCo4_4Sb12_{12}, we employ first-principles methods to compute the phase stability and electronic structure. Yb-filled CoSb3_3 exhibits (1) a mild tendency for phase separation into Yb-rich and Yb-poor regions and (2) a strong tendency for chemical decomposition into Co--Sb and Yb--Sb binaries (i.e., CoSb3_3, CoSb2_2, and YbSb2_2). We find that, at reasonable synthesis temperatures, configurational entropy stabilizes single-phase solid solutions with limited Yb solubility, in agreement with experiments. Filling CoSb3_3 with Yb increases the band gap, enhances the carrier effective masses, and generates new low-energy ``emergent'' conduction band minima, which is distinct from the traditional band convergence picture of aligning the energies of existing band extrema. The explicit presence of RR is necessary to achieve the emergent conduction band minima, though the rattler ordering does not strongly influence the electronic structure. The emergent conduction bands are spatially localized in the Yb-rich regions, unlike the delocalized electronic states at the Brillouin zone center that form the unfilled skutterudite band edges.

preprint2016arXiv

Electronic correlations in monolayer VS2_2

The layered transition metal dichalcogenide vanadium disulfide (VS2_2), which nominally has one electron in the 3d3d shell, is potent for strong correlation physics and is possibly another realization of an effective one-band model beyond the cuprates. Here monolayer VS2_2 in both the trigonal prismatic and octahedral phases is investigated using density functional theory plus Hubbard UU (DFT+UU) calculations. Trigonal prismatic VS2_2 has an isolated low-energy band that emerges from a confluence of crystal field splitting and direct V--V hopping. Within spin density functional theory, ferromagnetism splits the isolated band of the trigonal prismatic structure, leading to a low-band-gap S=1/2S=1/2 ferromagnetic Stoner insulator; the octahedral phase is higher in energy. Including the on-site interaction UU increases the band gap, leads to Mott insulating behavior, and for sufficiently high values stabilizes the ferromagnetic octahedral phase. The validity of DFT and DFT+UU for these two-dimensional materials with potential for strong electronic correlations is discussed. A clear benchmark is given by examining the experimentally observed charge density wave (CDW) in octahedral VS2_2, for which DFT grossly overestimates the bond length differences compared to known experiments; the presence of CDWs is also probed for the trigonal prismatic phase. Finally, we investigate why only the octahedral phase has been observed in experiments and discuss the possibility of realizing the trigonal prismatic phase. Our work suggests trigonal prismatic VS2_2 is a promising candidate for strongly correlated electron physics that, if realized, could be experimentally probed in an unprecedented fashion due to its monolayer nature.

preprint2013arXiv

Ideal strength and phonon instability of strained monolayer materials

The ideal strength of monolayer materials possessing semimetallic, semiconducting, and insulating ground states is computed using density functional theory. Here we show that, as in graphene, a soft mode occurs at the K-point in BN, graphane, and MoS2_2, while not in silicene. The transition is first-order in all cases except graphene. In BN and graphane the soft mode corresponds to a Kekul{é}-like distortion similar to that of graphene, while MoS2_2 has a distinct distortion. The phase transitions for BN, graphane, and MoS2_2 are not associated with the opening of a band gap, which indicates that Fermi surface nesting is not the driving force. We perform an energy decomposition that demonstrates why the soft modes at the K-point are unique and how strain drives the phonon instability.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.