Source author record

Chris Wolverton

Chris Wolverton appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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.

Published work

15 published item(s)

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 VPt$_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 VPt$_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 $3d$ 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 I$_2$, so that experiment and theory are brought into imperfect but reasonable agreement within about 120~meV/atom.

preprint2022arXiv

Identification of High-Dielectric Constant Compounds from Statistical Design

The discovery of high-dielectric materials is crucial to increasing the efficiency of electronic devices and batteries. Here, we report three previously unexplored materials with very high dielectric constants (69 $<$ $ε$ $<$ 101) and large band gaps (2.9$<$ $E_{\text{g}}$(eV) $<$ 5.5) obtained by screening materials databases using statistical optimization algorithms aided by artificial neural networks (ANN). Two of these new dielectrics are mixed-anion compounds (Eu$_5$SiCl$_6$O$_4$ and HoClO), and are shown to be thermodynamically stable against common semiconductors via phase-diagram analysis. We also uncovered four other materials with relatively large dielectric constants (20$<$$ε$$<$40) and band gaps (2.3$<$$E_{\text{g}}$(eV)$<$2.7). While the ANN training data is obtained from Materials Project, the search-space consists of materials from Open Quantum Materials Database (OQMD) - demonstrating a successful implementation of cross-database materials design. Overall, we report dielectric properties of 17 materials calculated using ab-initio calculations, that were selected in our design workflow. The dielectric materials with high dielectric properties predicted in this work open up further experimental research opportunities.

preprint2022arXiv

Testing the r$^2$SCAN 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 r$^2$SCAN 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 r$^2$SCAN description of solid-state thermodynamics, we evaluate the formation and decomposition enthalpies, equilibrium volumes, and fundamental bandgaps of more than 1,000 solids using r$^2$SCAN, SCAN, and PBE, as well as two dispersion-corrected variants, SCAN+rVV10 and r$^2$SCAN+rVV10. We show that r$^2$SCAN 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, r$^2$SCAN+rVV10 predicts marginally less-accurate formation enthalpies than r$^2$SCAN, and slightly more-accurate cell volumes than r$^2$SCAN. 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 r$^2$SCAN 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, r$^2$SCAN and r$^2$SCAN+rVV10 can be recommended as stable, general-purpose meta-GGAs for materials discovery.

preprint2021arXiv

Recent Advances and Applications of Deep Learning Methods in Materials Science

Deep learning (DL) is one of the fastest growing topics in materials data science, with rapidly emerging applications spanning atomistic, image-based, spectral, and textual data modalities. DL allows analysis of unstructured data and automated identification of features. Recent development of large materials databases has fueled the application of DL methods in atomistic prediction in particular. In contrast, advances in image and spectral data have largely leveraged synthetic data enabled by high quality forward models as well as by generative unsupervised DL methods. In this article, we present a high-level overview of deep-learning methods followed by a detailed discussion of recent developments of deep learning in atomistic simulation, materials imaging, spectral analysis, and natural language processing. For each modality we discuss applications involving both theoretical and experimental data, typical modeling approaches with their strengths and limitations, and relevant publicly available software and datasets. We conclude the review with a discussion of recent cross-cutting work related to uncertainty quantification in this field and a brief perspective on limitations, challenges, and potential growth areas for DL methods in materials science. The application of DL methods in materials science presents an exciting avenue for future materials discovery and design.

preprint2020arXiv

Accurate and scalable multi-element graph neural network force field and molecular dynamics with direct force architecture

Recently, machine learning (ML) has been used to address the computational cost that has been limiting ab initio molecular dynamics (AIMD). Here, we present GNNFF, a graph neural network framework to directly predict atomic forces from automatically extracted features of the local atomic environment that are translationally-invariant, but rotationally-covariant to the coordinate of the atoms. We demonstrate that GNNFF not only achieves high performance in terms of force prediction accuracy and computational speed on various materials systems, but also accurately predicts the forces of a large MD system after being trained on forces obtained from a smaller system. Finally, we use our framework to perform an MD simulation of Li7P3S11, a superionic conductor, and show that resulting Li diffusion coefficient is within 14% of that obtained directly from AIMD. The high performance exhibited by GNNFF can be easily generalized to study atomistic level dynamics of other material systems.

preprint2020arXiv

Direct Visualization of Electric Field induced Structural Dynamics in Monolayer Transition Metal Dichalcogenides

Layered transition metal dichalcogenides (TMDs) offer many attractive features for next-generation low-dimensional device geometries. Due to the practical and fabrication challenges related to in situ methods, the atomistic dynamics that give rise to realizable macroscopic device properties are often unclear. In this study, in situ transmission electron microscopy techniques are utilized in order to understand the structural dynamics at play, especially at interfaces and defects, in the prototypical film of monolayer MoS2 under electrical bias. Through our sample fabrication process, we clearly identify the presence of mass transport in the presence of a lateral electric field. In particular, we observe that the voids present at grain boundaries combine to induce structural deformation. The electric field mediates a net vacancy flux from the grain boundary interior to the exposed surface edge sites that leaves molybdenum clusters in its wake. Following the initial biasing cycles, however, the mass flow is largely diminished, and the resultant structure remains stable over repeated biasing. We believe insights from this work can help explain observations of non-uniform heating and preferential oxidation at grain boundary sites in these materials.

preprint2020arXiv

Microscopic Mechanisms of Glass-Like Lattice Thermal Transport in Cubic Cu$_{12}$Sb$_{4}$S$_{13}$ Tetrahedrites

Materials based on cubic tetrahedrites (Cu$_{12}$Sb$_{4}$S$_{13}$) are useful thermoelectrics with unusual thermal and electrical transport properties, such as very low and nearly temperature-independent lattice thermal conductivity ($κ_{L}$). We explain the microscopic origin of the glass-like $κ_{L}$ in Cu$_{12}$Sb$_{4}$S$_{13}$ by explicitly treating anharmonicity up to quartic terms for both phonon energies and phonon scattering rates. We show that the strongly unstable phonon modes associated with trigonally coordinated Cu atoms are anharmonically stabilized above approximately $100$ K and continue hardening with increasing temperature, in accord with experimental data. This temperature induced hardening effect reduces scattering of heat carrying acoustic modes by reducing the available phase space for three-phonon processes, thereby balancing the conventional $\propto T$ increase in scattering due to phonon population and yielding nearly temperature-independent $κ_{L}$. Furthermore, we find that very strong phonon broadening lead to a qualitative breakdown of the conventional phonon-gas model and modify the dominant heat transport mechanism from the particle-like phonon wave packet propagation to incoherent tunneling described by off-diagonal terms in the heat-flux operator, which are typically prevailing in glasses and disordered crystals. Our work paves the way to a deeper understanding of glass-like thermal conductivity in complex crystals with strong anharmonicity.

preprint2020arXiv

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

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 $d$ 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 Li$_x$TiS$_2$, Li$_x$NiO$_2$, and Li$_x$CoO$_2$, olivine Li$_x$FePO$_4$, and spinel Li$_x$Mn$_2$O$_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 $U$ 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

Developing an improved Crystal Graph Convolutional Neural Network framework for accelerated materials discovery

The recently proposed crystal graph convolutional neural network (CGCNN) offers a highly versatile and accurate machine learning (ML) framework by learning material properties directly from graph-like representations of crystal structures ("crystal graphs"). Here, we develop an improved variant of the CGCNN model (iCGCNN) that outperforms the original by incorporating information of the Voronoi tessellated crystal structure, explicit 3-body correlations of neighboring constituent atoms, and an optimized chemical representation of interatomic bonds in the crystal graphs. We demonstrate the accuracy of the improved framework in two distinct illustrations: First, when trained/validated on 180,000/20,000 density functional theory (DFT) calculated thermodynamic stability entries taken from the Open Quantum Materials Database (OQMD) and evaluated on a separate test set of 230,000 entries, iCGCNN achieves a predictive accuracy that is significantly improved, i.e., 20% higher than that of the original CGCNN. Second, when used to assist high-throughput search for materials in the ThCr2Si2 structure-type, iCGCNN exhibited a success rate of 31% which is 310 times higher than an undirected high-throughput search and 2.4 times higher than that of the original CGCNN. Using both CGCNN and iCGCNN, we screened 132,600 compounds with elemental decorations of the ThCr2Si2 prototype crystal structure and identified a total of 97 new unique stable compounds by performing 757 DFT calculations, accelerating the computational time of the high-throughput search by a factor of 130. Our results suggest that the iCGCNN can be used to accelerate high-throughput discoveries of new materials by quickly and accurately identifying crystalline compounds with properties of interest.

preprint2019arXiv

Electronic Structure and Phase Stability of Yb-filled CoSb$_3$ Skutterudite Thermoelectrics from First Principles

Filling the large voids in the crystal structure of the skutterudite CoSb$_3$ with rattler atoms $R$ provides an avenue for both increasing carrier concentration and disrupting lattice heat transport, leading to impressive thermoelectric performance. While the influence of $R$ on the lattice dynamics of skutterudite materials has been well studied, the phase stability of $R$-filled skutterudite materials and the influence of the presence and ordering of $R$ on the electronic structure remain unclear. Here, focusing on the Yb-filled skutterudite Yb$_x$Co$_4$Sb$_{12}$, we employ first-principles methods to compute the phase stability and electronic structure. Yb-filled CoSb$_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., CoSb$_3$, CoSb$_2$, and YbSb$_2$). We find that, at reasonable synthesis temperatures, configurational entropy stabilizes single-phase solid solutions with limited Yb solubility, in agreement with experiments. Filling CoSb$_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 $R$ 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

Lithium transport through Lithium-ion battery cathode coatings

The surface coating of cathodes using insulator films has proven to be a promising method for high-voltage cathode stabilization in Li-ion batteries. However, there is still substantial uncertainty about how these films function, specifically with regard to important coating design principles such as lithium solubility and transport through the films. This study uses Density Functional Theory to examine the diffusivity of interstitial lithium in crystalline α-$AlF_3$, α-$Al_2O_3$, m-$ZrO_2$, c-MgO, and α-quartz $SiO_2$, which provide benchmark cases for further understanding of insulator coatings in general. In addition, we propose an Ohmic electrolyte model to predict resistivities and overpotential contributions under battery operating conditions. For the crystalline materials considered we predict that Li+ diffuses quite slowly, with a migration barrier larger than 0.9 eV in all crystalline materials except α-quartz $SiO_2$, which is predicted to have a migration barrier of 0.276 eV along <001>. These results suggest that the stable crystalline forms of these insulator materials, except for oriented α-quartz $SiO_2$, are not practical for conformal cathode coatings. Amorphous $Al_2O_3$ and $AlF_3$ have higher Li+ diffusivities than their crystalline counterparts. Our predicted amorphous $Al_2O_3$ resistivity (1789 MΩm) is near the top of the range of fitted resistivities extracted from previous experiments on nominal $Al_2O_3$ coatings (7.8 to 913 MΩm) while our predicted amorphous $AlF_3$ resistivity (114 MΩm) is close to the middle of the range. These comparisons support our framework for modeling and understanding the impact on overpotential of conformal coatings in terms of their fundamental thermodynamic and kinetic properties, and support that these materials can provide practical conformal coatings in their amorphous form.

preprint2016arXiv

Prediction of superconducting iron-bismuth intermetallic compounds at high pressure

The synthesis of materials in high-pressure experiments has recently attracted increasing attention, especially since the discovery of record breaking superconducting temperatures in the sulfur-hydrogen and other hydrogen-rich systems. Commonly, the initial precursor in a high pressure experiment contains constituent elements that are known to form compounds at ambient conditions, however the discovery of high-pressure phases in systems immiscible under ambient conditions poses an additional materials design challenge. We performed an extensive multi component $ab\,initio$ structural search in the immiscible Fe--Bi system at high pressure and report on the surprising discovery of two stable compounds at pressures above $\approx36$ GPa, FeBi$_2$ and FeBi$_3$. According to our predictions, FeBi$_2$ is a metal at the border of magnetism with a conventional electron-phonon mediated superconducting transition temperature of $T_{\rm c}=1.3$ K at 40 GPa. In analogy to other iron-based materials, FeBi$_2$ is possibly a non-conventional superconductor with a real $T_{\rm c}$ significantly exceeding the values obtained within Bardeen-Cooper-Schrieffer (BCS) theory.

preprint2016arXiv

Ultralow Thermal Conductivity in Full-Heusler Semiconductors

Semiconducting half- and, to a lesser extent, full-Heusler compounds are promising thermoelectric materials due to their compelling electronic properties with large power factors. However, intrinsically high thermal conductivity resulting in a limited thermoelectric efficiency has so far impeded their widespread use in practical applications. Here, we report the computational discovery of a class of hitherto unknown stable semiconducting full-Heusler compounds with ten valence electrons ($X_2YZ$, $X$=Ca, Sr, and Ba; $Y$= Au and Hg; $Z$=Sn, Pb, As, Sb, and Bi) through high-throughput $ab-initio$ screening. These new compounds exhibit ultralow lattice thermal conductivity $κ_{\text{L}}$ close to the theoretical minimum due to strong anharmonic rattling of the heavy noble metals, while preserving high power factors, thus resulting in excellent phonon-glass electron-crystal materials.

preprint2015arXiv

Superconductivity in metastable phases of phosphorus-hydride compounds under high pressure

Hydrogen-rich compounds have been extensively studied both theoretically and experimentally in the quest for novel high-temperature superconductors. Reports on sulfur-hydride attaining metallicity under pressure and exhibiting superconductivity at temperatures as high as 200 K have spurred an intense search for room-temperature superconductors in hydride materials. Recently, compressed phosphine was reported to metallize at pressures above 45 GPa, reaching a superconducting transition temperature (T$_{c}$) of 100 K at 200 GPa. However, neither the exact composition nor the crystal structure of the superconducting phase have been conclusively determined. In this work the phase diagram of PH$_n$ ($n=1,2,3,4,5,6$) was extensively explored by means of {\it ab initio} crystal structure predictions using the Minima Hopping Method (MHM). The results do not support the existence of thermodynamically stable PH$_n$ compounds, which exhibit a tendency for elemental decomposition at high pressure even when vibrational contributions to the free energies are taken into account. Although the lowest energy phases of PH$_{1,2,3}$ display T$_{c}$'s comparable to experiments, it remains uncertain if the measured values of T$_{c}$ can be fully attributed to a phase-pure compound of PH$_n$.

preprint2012arXiv

Thermodynamic Stability of Co-Al-W L12 γ'

Co-based superalloys in the Co-Al-W system exhibit coherent L12 Co3(Al,W) γ' precipitates in an fcc Co γmatrix, analogous to Ni3Al in Ni-based systems. Unlike Ni3Al however, experimental observations of Co3(Al,W) suggest that it is not a stable phase at 1173K. Here, we perform an extensive series of density functional theory (DFT) calculations of the γ' Co3(Al,W) phase stability, including point defect energetics and finite-temperature contributions. We first confirm and extend previous DFT calculations of the metastability of L12 Co3(Al0.5W0.5) γ' at 0K with respect to HCP Co, B2 CoAl, and D019 Co3W using the special quasi-random structure (SQS) approach to describe the Al/W solid solution, employing several exchange/correlation functionals, structures with varying degrees of disorder, and newly developed larger SQS. We expand the validity of this conclusion by considering the formation of antisite and vacancy point defects, predicting defect formation energies similar in magnitude to Ni3Al. However, in contrast to the Ni3Al system, we find that substituting Co on Al sites is thermodynamically favorable at 0K, consistent with experimental observation of Co excess and Al deficiency in γ' with respect to the Co3(Al0.5W0.5) composition. Lastly, we consider vibrational, electronic, and magnetic contributions to the free energy, finding that they promote the stability of γ', making the phase thermodynamically competitive with the convex hull at elevated temperature. Surprisingly, this is due to the relatively small finite-temperature contributions of one of the γ' decomposition products, B2 CoAl, effectively destabilizing the Co, CoAl, and Co3W three phase mixture, thus stabilizing the γ' phase.