BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Paul Erhart

Paul Erhart contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-scicond-mat.mes-hallphysics.comp-phphysics.chem-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

18 published item(s)

preprint2022arXiv

The moiré potential in twisted transition metal dichalcogenide bilayers

Moiré superlattices serve as a playground for emerging phenomena, such as localization of band states, superconductivity, and localization of excitons. These superlattices are large and are often modeled in the zero angle limit, which obscures the effect of finite twist angles. Here, by means of first-principles calculations we quantify the twist-angle dependence of the moiré potential in the MoS$_2$ homobilayer and identify the contributions from the constituent elements of the moiré potential. Furthermore, by considering the zero-angle limit configurations, we show that the moiré potential is rather homogeneous across the transition metal dichalcogenides (TMDs) and briefly discuss the separate effects of potential shifts and hybridization on the bilayer hybrid excitons. We find that the moiré potential in TMDs exhibits both an electrostatic component and a hybridization component, which are intertwined and have different relative strengths in different parts of the Brillouin zone. The electrostatic component of the moiré potential is a varying dipole field, which has a strong twist angle dependence. In some cases, the hybridization component can be interpreted as a tunneling rate but the interpretation is not generally applicable over the full Brillouin zone.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantitative predictions of thermodynamic hysteresis: Temperature-dependent character of the phase transition in Pd-H

The thermodynamics of phase transitions between phases that are size-mismatched but coherent differs from conventional stress-free thermodynamics. Most notably, in open systems such phase transitions are always associated with hysteresis. In spite of experimental evidence for the relevance of these effects in technologically important materials such as Pd hydride, a recipe for first-principles-based atomic-scale modeling of coherent, open systems has been lacking. Here, we develop a methodology for quantifying phase boundaries, hysteresis, and coherent interface free energies using density-functional theory, alloy cluster expansions, and Monte Carlo simulations in a constrained ensemble. We apply this approach to Pd--H and show that the phase transition changes character above approximately 400 K, occurring with an at all times spatially homogeneous hydrogen concentration, i.e., without coexistence between the two phases. Our results are consistent with experimental observations but reveal aspects of hydride formation in Pd nanoparticles that have not yet been accessible in experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient construction of linear models in materials modeling and applications to force constant expansions

Linear models, such as force constant (FC) and cluster expansions, play a key role in physics and materials science. While they can in principle be parametrized using regression and feature selection approaches, the convergence behavior of these techniques, in particular with respect to thermodynamic properties is not well understood. Here, we therefore analyze the efficacy and efficiency of several state-of-the-art regression and feature selection methods, in particular in the context of FC extraction and the prediction of different thermodynamic properties. Generic feature selection algorithms such as recursive feature elimination with ordinary least-squares (OLS), automatic relevance determination regression, and the adaptive least absolute shrinkage and selection operator can yield physically sound models for systems with a modest number of degrees of freedom. For large unit cells with low symmetry and/or high-order expansions they come, however, with a non-negligible computational cost that can be more than two orders of magnitude higher than that of OLS. In such cases, OLS with cutoff selection provides a viable route as demonstrated here for both second-order FCs in large low-symmetry unit cells and high-order FCs in low-symmetry systems. While regression techniques are thus very powerful, they require well-tuned protocols. Here, the present work establishes guidelines for the design of protocols that are readily usable, e.g., in high-throughput and materials discovery schemes. Since the underlying algorithms are not specific to FC construction, the general conclusions drawn here also have a bearing on the construction of other linear models in physics and materials science.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hot-Carrier Generation in Plasmonic Nanoparticles: The Importance of Atomic Structure

Metal nanoparticles are attractive for plasmon-enhanced generation of hot carriers, which may be harnessed in photochemical reactions. In this work, we analyze the coherent femtosecond dynamics of photon absorption, plasmon formation, and subsequent hot-carrier generation through plasmon dephasing using first-principles simulations. We predict the energetic and spatial hot-carrier distributions in small metal nanoparticles and show that the distribution of hot electrons is very sensitive to the local structure. Our results show that surface sites exhibit enhanced hot-electron generation in comparison to the bulk of the nanoparticle. While the details of the distribution depend on particle size and shape, as a general trend lower-coordinated surface sites such as corners, edges, and {100} facets exhibit a higher proportion of hot electrons than higher-coordinated surface sites such as {111} facets or the core sites. The present results thereby demonstrate how hot carriers could be tailored by careful design of atomic-scale structures in nanoscale systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Luminescence quenching via deep defect states: A recombination pathway via oxygen vacancies in Ce-doped YAG

Luminescence quenching via non-radiative recombination channels limits the efficiency of optical materials such as phosphors and scintillators and therefore has implications for conversion efficiency and device lifetimes. In materials such as Ce-doped yttrium aluminum garnet (YAG:Ce), quenching shows a strong dependence on both temperature and activator concentration, limiting the ability to fabricate high-intensity white-light emitting diodes with high operating temperatures. Here, we reveal by means of first-principles calculations an efficient recombination mechanism in YAG:Ce that involves oxygen vacancies and gives rise to thermally activated concentration quenching. We demonstrate that the key requirements for this mechanism to be active are localized states with strong electron-phonon coupling. These conditions are commonly found for intrinsic defects such as anion vacancies in wide band-gap materials. The present findings are therefore relevant to a broad class of optical materials and shine light on thermal quenching mechanisms in general.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Defects from phonons: Atomic transport by concerted motion in simple crystalline metals

Point defects play a crucial role in crystalline materials as they do not only impact the thermodynamic properties but are also central to kinetic processes. While they are necessary in thermodynamic equilibrium spontaneous defect formation in the bulk is normally considered highly improbable except for temperatures close to the melting point. Here, we demonstrate by means of atomistic simulations that processes involving concerted atomic motion that give rise to defect formation are in fact frequent in body-centered cubic metals even down to about 50% of the melting temperature. It is shown that this behavior is intimately related to the anharmonicity of the lattice vibrations and a flat energy landscape along certain crystallographic directions, a feature that is absent in, e.g., face-centered cubic lattice structures. This insight has implications for our general understanding of these materials and furthermore provides a complementary explanation for the so-called anomalous diffusion in group 4 transition metals.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

icet - A Python library for constructing and sampling alloy cluster expansions

Alloy cluster expansions (CEs) provide an accurate and computationally efficient mapping of the potential energy surface of multi-component systems that enables comprehensive sampling of the many-dimensional configuration space. Here, we introduce \textsc{icet}, a flexible, extensible, and computationally efficient software package for the construction and sampling of CEs. \textsc{icet} is largely written in Python for easy integration in comprehensive workflows, including first-principles calculations for the generation of reference data and machine learning libraries for training and validation. The package enables training using a variety of linear regression algorithms with and without regularization, Bayesian regression, feature selection, and cross-validation. It also provides complementary functionality for structure enumeration and mapping as well as data management and analysis. Potential applications are illustrated by two examples, including the computation of the phase diagram of a prototypical metallic alloy and the analysis of chemical ordering in an inorganic semiconductor.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Structurally driven asymmetric miscibility in the phase diagram of W-Ti

Phase diagrams for multi-component systems represent crucial information for understanding and designing materials but are very time consuming to assess experimentally. Computational modeling plays an increasingly important role in this regard but has been largely focused on systems with matching lattice structures and/or stable boundary phases. Here, using a combination of density functional theory calculations, alloy cluster expansions, free energy integration, and Monte Carlo simulations, we obtain the phase diagram of W-Ti, a system that features metastable boundary phases on both sides of the phase diagram. We find that the mixing energy on the body-centered cubic (BCC) lattice is asymmetric and negative with a minimum of about -120 meV/atom, whereas for the hexagonal closed packed (HCP) lattice the mixing energy is positive. By combining these data with a model for the vibrational free energy, we propose a revision of the W-rich end of the phase diagram with a much larger solubility of Ti in BCC-W than previous assessments. Finally, by comparison with the W-V and W-Re system we demonstrate how strongly asymmetric phase diagrams can arise from a subtle energy balance of stable and metastable lattice structures.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The hiphive package for the extraction of high-order force constants by machine learning

The efficient extraction of force constants (FCs) is crucial for the analysis of many thermodynamic materials properties. Approaches based on the systematic enumeration of finite differences scale poorly with system size and can rarely extend beyond third order when input data is obtained from first-principles calculations. Methods based on parameter fitting in the spirit of interatomic potentials, on the other hand, can extract FC parameters from semi-random configurations of high information density and advanced regularized regression methods can recover physical solutions from a limited amount of data. Here, we present the hiPhive Python package, that enables the construction of force constant models up to arbitrary order. hiPhive exploits crystal symmetries to reduce the number of free parameters and then employs advanced machine learning algorithms to extract the force constants. Depending on the problem at hand both over and underdetermined systems are handled efficiently. The FCs can be subsequently analyzed directly and or be used to carry out e.g., molecular dynamics simulations. The utility of this approach is demonstrated via several examples including ideal and defective monolayers of MoS$_2$ as well as bulk nickel.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Thermal conductivity in intermetallic clathrates: A first principles perspective

Inorganic clathrates such as Ba$_8$Ga$_{x}$Ge$_{46-x}$ and Ba$_8$Al$_{x}$Si$_{46-x}$ commonly exhibit very low thermal conductivities. A quantitative computational description of this important property has proven difficult, in part due to the large unit cell, the role of disorder, and the fact that both electronic carriers and phonons contribute to transport. Here, we conduct a systematic analysis of the temperature and composition dependence of low-frequency modes associated with guest species in Ba$_8$Ga$_{x}$Ge$_{46-x}$ and Ba$_8$Al$_{x}$Si$_{46-x}$ ("rattler modes"), as well as of thermal transport in stoichiometric Ba$_8$Ga$_{16}$Ge$_{30}$. To this end, we account for phonon-phonon interactions by means of temperature dependent effective interatomic force constants (TDIFCs), which we find to be crucial in order to achieve an accurate description of the lattice part of the thermal conductivity. While the analysis of the thermal conductivity is often largely focused on the rattler modes, here, it is shown that at room temperatures modes with $\hbarω\gtrsim\,10\,\text{meV}$ account for 50\%\ of lattice heat transport. Finally, the electronic contribution to the thermal conductivity is computed, which shows the Wiedemann-Franz law to be only approximately fulfilled. As a result, it is crucial to employ the correct prefactor when separating electronic and lattice contributions for experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

atomicrex - A general purpose tool for the construction of atomic interaction models

We introduce atomicrex, an open-source code for constructing interatomic potentials as well as more general types of atomic-scale models. Such effective models are required to simulate extended materials structures comprising many thousands of atoms or more, because electronic structure methods become computationally too expensive at this scale. atomicrex covers a wide range of interatomic potential types and fulfills many needs in atomistic model development. As inputs, it supports experimental property values as well as \textit{ab initio} energies and forces, to which models can be fitted using various optimization algorithms. The open architecture of atomicrex allows it to be used in custom model development scenarios beyond classical interatomic potentials while thanks to its Python interface it can be readily integrated e.g., with electronic structure calculations or machine learning algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

libvdwxc: A library for exchange-correlation functionals in the vdW-DF family

We present libvdwxc, a general library for evaluating the energy and potential for the family of vdW-DF exchange--correlation functionals. libvdwxc provides an efficient implementation of the vdW-DF method and can be interfaced with various general-purpose DFT codes. Currently, the GPAW and Octopus codes implement interfaces to libvdwxc. The present implementation emphasizes scalability and parallel performance, and thereby enables \textit{ab initio} calculations of nanometer-scale complexes. The numerical accuracy is benchmarked on the S22 test set whereas parallel performance is benchmarked on ligand-protected gold nanoparticles ($\text{Au}_{144}(\text{SC}_{11}\text{NH}_{25})_{60}$) up to 9696 atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Mechanism of Re precipitation in irradiated W-Re alloys from kinetic Monte Carlo simulations

High-temperature, high-dose, neutron irradiation of W results in the formation of Re-rich clusters at concentrations one order of magnitude lower than the thermodynamic solubility limit. These clusters may eventually transform into brittle W-Re intermetallic phases, which can lead to high levels of hardening and thermal conductivity losses. Standard theories of radiation enhanced diffusion and precipitation cannot explain the formation of these precipitates and so understanding the mechanism by which nonequilibrium clusters form under irradiation is crucial to predict materials degradation and devise mitigation strategies. Here we carry out a thermodynamic study of W-Re alloys and conduct kinetic Monte Carlo simulations of Re cluster formation in irradiated W-2Re alloys using a generalized Hamiltonian for crystals containing point defects parameterized entirely with electronic structure calculations. Our model incorporates recently-gained mechanistic information of mixed-interstitial solute transport, which is seen to control cluster nucleation and growth by forming quasi-spherical nuclei after an average incubation time of 20 s at 1800 K. These nuclei are seen to grow by attracting more mixed interstitials bringing solute atoms, which in turns attracts vacancies leading to recombination and solute agglomeration. The clusters grow to a maximum size of approximately 4-nm radius, and are not fully dense with Re, containing 50% or less near the center. Our simulations are in reasonable agreement with recent atom probe examinations of ion irradiated W-2Re systems at 773 K.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Finite-temperature properties of non-magnetic transition metals: Comparison of the performance of constraint-based semi and nonlocal functionals

We assess the performance of nonempirical, truly nonlocal and semi-local functionals with regard to structural and thermal properties of $3d$, $4d$, and $5d$ non-magnetic transition metals. We focus on constraint-based functionals and consider the new consistent-exchange van der Waals density functional version vdW-DF-cx [Phys. Rev. B 89, 035412 (2014)], the semi-local PBE [Phys. Rev. Lett. 77, 3865 (1996)] and PBEsol functionals [Phys. Rev. Lett. 100, 136406 (2008)] as well as the AM05 meta-functional [Phys. Rev. B 72, 085108 (2005)]. Using the quasi-harmonic approximation structural parameters, elastic response, and thermal expansion at finite temperatures are computed and compared to experimental data. We also compute cohesive energies explicitly including zero-point vibrations. It is shown that overall vdW-DF-cx provides an accurate description of thermal properties and retains a level of transferability and accuracy that is comparable to or better than some of the best constraint-based semi-local functionals. Especially, with regard to the cohesive energies the consistent inclusion of spin polarization effects in the atoms turns out to be crucial and it is important to use the rigorous spin-vdW-DF-cx formulation [Phys. Rev. Lett. 115, 136402 (2015)]. This demonstrates that vdW-DF-cx has general-purpose character and can be used to study systems that have both sparse and dense electron distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermal expansion and transport in van der Waals solids from first-principles calculations

The lattice thermal expansion and conductivity in bulk Mo and W-based transition metal dichalcogenides are investigated by means of density functional and Boltzmann transport theory calculations. To this end, a recent van der Waals density functional (vdW-DF-CX) is employed, which is shown to yield excellent agreement with reference data for the structural parameters. The calculated in-plane thermal conductivity compares well with experimental room temperature values, when phonon-phonon and isotopic scattering are included. To explain the behavior over the entire available temperature range one must, however, include additional (temperature independent) scattering mechanisms that limit the mean free path. Generally, the primary heat carrying modes have mean free paths of $1\,μ\text{m}$ or more, which makes these materials very susceptible to structural defects. The conductivity of Mo and W-based TMDs is primarily determined by the chalcogenide species and increases in the order Te-Se-S. While for the tellurides and selenides the transition metal element has a negligible effect, the conductivity of WS$_2$ is notably higher than for MoS$_2$, which can be traced to the much larger phonon band gap of the former. Overall the present study provides a consistent set of thermal conductivities that reveal chemical trends and constitute the basis for future investigations of van der Waals solids.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A first-principles study of co-doping in lanthanum bromide

Co-doping of Ce-doped LaBr$_3$ with Ba, Ca, or Sr improves the energy resolution that can be achieved by radiation detectors based on these materials. Here, we present a mechanism that rationalizes of this enhancement that on the basis of first principles electronic structure calculations and point defect thermodynamics. It is shown that incorporation of Sr creates neutral $V_\text{Br}$-Sr$_\text{La}$ complexes that can temporarily trap electrons. As a result, Auger quenching of free carriers is reduced, allowing for a more linear, albeit slower, scintillation light yield response. Experimental Stokes shifts can be related to different Ce$_\text{La}$-Sr$_\text{La}$-$V_\text{Br}$ triple complex configurations. Co-doping with other alkaline as well as alkaline earth metals is considered as well. Alkaline elements are found to have extremely small solubilities on the order of 0.1 ppm and below at 1000 K. Among the alkaline earth metals the lighter dopant atoms prefer interstitial-like positions and create strong scattering centers, which has a detrimental impact on carrier mobilities. Only the heavier alkaline earth elements combine matching ionic radii with sufficiently high solubilities. This provides a rationale for the experimental finding that improved scintillator performance is exclusively achieved using Sr, Ca, or Ba. The present mechanism demonstrates that co-doping of wide gap materials can provide an efficient means for managing charge carrier populations under out-of-equilibrium conditions. In the present case dopants are introduced that manipulate not only the concentrations but the electronic properties of intrinsic defects without introducing additional gap levels. This leads to the availability of shallow electron traps that can temporarily localize charge carriers, effectively deactivating carrier-carrier recombination channels. The principles of this ... [continued]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A variational polaron self-interaction corrected total-energy functional for charge excitations in wide-band gap insulators

We conduct a detailed investigation of the polaron self-interaction (pSI) error in standard approximations to the exchange-correlation (XC) functional within density-functional theory (DFT). The pSI leads to delocalization error in the polaron wave function and energy, as calculated from the Kohn-Sham (KS) potential in the native charge state of the polaron. This constitutes the origin of the systematic failure of DFT to describe polaron formation in band insulators. It is shown that the delocalization error in these systems is, however, largely absent in the KS potential of the closed-shell neutral charge state. This leads to a modification of the DFT total-energy functional that corrects the pSI in the XC functional. The resulting pSIC-DFT method constitutes an accurate parameter-free {\it ab initio} methodology for calculating polaron properties in insulators at a computational cost that is orders of magnitude smaller than hybrid XC functionals. Unlike approaches that rely on parametrized localized potentials such as DFT+$U$, the pSIC-DFT method properly captures both site and bond-centered polaron configurations. This is demonstrated by studying formation and migration of self-trapped holes in alkali halides (bond-centered) as well as self-trapped electrons in an elpasolite compound (site-centered). The pSIC-DFT approach consistently reproduces the results obtained by hybrid XC functionals parametrized by DFT+$G_0W_0$ calculations. Finally, we generalize the pSIC approach to hybrid functionals, and show that in stark contrast to conventional hybrid calculations of polaron energies, the pSIC-hybrid method is insensitive to the parametrization of the hybrid XC functional. On this basis, we further rationalize the success of the pSIC-DFT approach.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Surface potential at a ferroelectric grain due to asymmetric screening of depolarization fields

Nonlinear screening of electric depolarization fields, generated by a stripe domain structure in a ferroelectric grain of a polycrystalline material, is studied within a semiconductor model of ferroelectrics. It is shown that the maximum strength of local depolarization fields is rather determined by the electronic band gap than by the spontaneous polarization magnitude. Furthermore, field screening due to electronic band bending and due to presence of intrinsic defects leads to asymmetric space charge regions near the grain boundary, which produce an effective dipole layer at the surface of the grain. This results in the formation of a potential difference between the grain surface and its interior of the order of 1 V, which can be of either sign depending on defect transition levels and concentrations. Exemplary acceptor doping of BaTiO3 is shown to allow tuning of the said surface potential in the region between 0.1 and 1.3 V.

0 citations0 reviewsNo code linkedOpen access