BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Paul Erhart

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

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.comp-phcond-mat.mes-hallphysics.chem-phcond-mat.stat-mechcond-mat.str-el

Catalog footprint

What is connected

41works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

41 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

Polaronic contributions to oxidation and hole conductivity in acceptor-doped BaZrO$_3$

Acceptor-doped perovskite oxides like BaZrO$_3$ are showing great potential as materials for renewable energy technologies where hydrogen acts an energy carrier, such as solid oxide fuel cells and hydrogen separation membranes. While ionic transport in these materials has been investigated intensively, the electronic counterpart has received much less attention and further exploration in this field is required. Here, we use density functional theory (DFT) to study hole polarons and their impact on hole conductivity in Y-doped BaZrO$_3$. Three different approaches have been used to remedy the self-interaction error of local and semi-local exchange-correlation functionals: DFT$+U$, pSIC-DFT and hybrid functionals. Self-trapped holes are found to be energetically favorable by about $-$0.1 eV and the presence of yttrium results in further stabilization. Polaron migration is predicted to occur through intraoctahedral transfer and polaron rotational processes, which are associated with adiabatic barriers of about 0.1 eV. However, the rather small energies associated with polaron formation and migration suggest that the hole becomes delocalized and band-like at elevated temperatures. These results together with an endothermic oxidation reaction [A. Lindman et al., Phys. Rev. B 91, 245114 (2015)] yield a picture that is consistent with experimental data for the hole conductivity. The results we present here provide new insight into hole transport in acceptor-doped BaZrO$_3$ and similar materials, which will be of value in the future development of sustainable technologies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The role of interstitial binding in radiation induced segregation in W-Re alloys

Due to their high strength and advantageous high-temperature properties, tungsten-based alloys are being considered as plasma-facing candidate materials in fusion devices. Under neutron irradiation, rhenium, which is produced by nuclear transmutation, has been found to precipitate in elongated precipitates forming thermodynamic intermetallic phases at concentrations well below the solubility limit. Recent measurements have shown that Re precipitation can lead to substantial hardening, which may have a detrimental effect on the fracture toughness of W alloys. This puzzle of sub-solubility precipitation points to the role played by irradiation induced defects, specifically mixed solute-W interstitials. Here, using first-principles calculations based on density functional theory, we study the energetics of mixed interstitial defects in W-Re, W-V, and W-Ti alloys, as well as the heat of mixing for each substitutional solute. We find that mixed interstitials in all systems are strongly attracted to each other with binding energies of -2.4 to -3.2 eV and form interstitial pairs that are aligned along parallel first-neighbor <111> strings. Low barriers for defect translation and rotation enable defect agglomeration and alignment even at moderate temperatures. We propose that these elongated agglomerates of mixed-interstitials may act as precursors for the formation of needle-shaped intermetallic precipitates. This interstitial-based mechanism is not limited to radiation induced segregation and precipitation in W-Re alloys but is also applicable to other body-centered cubic alloys.

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
preprint2015arXiv

Implications of the band gap problem on oxidation and hydration in acceptor-doped barium zirconate

Charge carrier concentrations in acceptor-doped proton-conducting perovskites are to a large extent determined by the hydration and oxidation of oxygen vacancies, which introduce protons and holes, respectively. First-principles modeling of these reactions involves calculation of formation energies of charged defects, which requires an accurate description of the band gap and the position of the band edges. Since density-functional theory (DFT) with local and semi-local exchange-correlation functionals (LDA and GGA) systematically fails to predict these quantities this can have serious implications on the modeling of defect reactions. In this study we investigate how the description of band gap and band edge positions affects the hydration and oxidation in acceptor-doped BaZrO$_3$. First-principles calculations are performed in combination with thermodynamic modeling in order to obtain equilibrium charge carrier concentrations at different temperatures and partial pressures. Three different methods have been considered: DFT with both semi-local (PBE) and hybrid (PBE0) exchange-correlation functionals, and many-body perturbation theory within the $G_0W_0$-approximation. All three methods yield similar results for the hydration reaction, which are consistent with experimental findings. For the oxidation reaction, on the other hand, there is a qualitative difference. PBE predicts the reaction to be exothermic while the two others predict an endothermic behavior. Results from thermodynamic modeling are compared with available experimental data, such as enthalpies, concentrations and conductivities, and only the results obtained with PBE0 and $G_0W_0$, with an endothermic oxidation behavior, give a satisfactory agreement with experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficacy of the DFT+U formalism for modeling hole polarons in perovskite oxides

We investigate the formation of self-trapped holes (STH) in three prototypical perovskites (SrTiO3, BaTiO3, PbTiO3) using a combination of density functional theory (DFT) calculations with local potentials and hybrid functionals. First we construct a local correction potential for polaronic configurations in SrTiO3 that is applied via the DFT+U method and matches the forces from hybrid calculations. We then use the DFT+U potential to search the configuration space and locate the lowest energy STH configuration. It is demonstrated that both the DFT+U potential and the hybrid functional yield a piece-wise linear dependence of the total energy on the occupation of the STH level suggesting that self-interaction effects have been properly removed. The DFT+U model is found to be transferable to BaTiO3 and PbTiO3, and formation energies from DFT+U and hybrid calculations are in close agreement for all three materials. STH formation is found to be energetically favorable in SrTiO3 and BaTiO3 but not in PbTiO3, which can be rationalized by considering the alignment of the valence band edges on an absolute energy scale. In the case of PbTiO3 the strong coupling between Pb 6s and O 2p states lifts the valence band minimum (VBM) compared to SrTiO3 and BaTiO3. This reduces the separation between VBM and STH level and renders the STH configuration metastable with respect to delocalization (band hole state). We expect that the present approach can be adapted to study STH formation also oxides with different crystal structures and chemical composition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Origin of resolution enhancement by co-doping of scintillators: Insight from electronic structure calculations

It was recently shown that the energy resolution of Ce-doped LaBr$_3$ scintillator radiation detectors can be crucially improved by co-doping with Sr, Ca, or Ba. Here we outline a mechanism for this enhancement on the basis of electronic structure calculations. We show that (i) Br vacancies are the primary electron traps during the initial stage of thermalization of hot carriers, prior to hole capture by Ce dopants; (ii) isolated Br vacancies are associated with deep levels; (iii) Sr doping increases the Br vacancy concentration by several orders of magnitude; (iv) $\text{Sr}_\text{La}$ binds to $V_\text{Br}$ resulting in a stable neutral complex; and (v) association with Sr causes the deep vacancy level to move toward the conduction band edge. The latter is essential for reducing the effective carrier density available for Auger quenching during thermalization of hot carriers. Subsequent de-trapping of electrons from $\text{Sr}_\text{La}-V_\text{La}$ complexes then can activate Ce dopants that have previously captured a hole leading to luminescence. This mechanism implies an overall reduction of Auger quenching of free carriers, which is expected to improve the linearity of the photon light yield with respect to the energy of incident electron or photon.

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
preprint2013arXiv

Formation and switching of defect dipoles in acceptor doped lead titanate: A kinetic model based on first-principles calculations

The formation and field-induced switching of defect dipoles in acceptor doped lead titanate is described by a kinetic model representing an extension of the well established Arlt-Neumann model [Ferroelectrics {\bf 76}, 303 (1987)]. Energy barriers for defect association and reorientation of oxygen vacancy-dopant (Cu and Fe) complexes are obtained from first-principles calculations and serve as input data for the kinetic coefficients in the rate equation model. The numerical solution of the model describes the time evolution of the oxygen vacancy distribution at different temperatures and dopant concentrations in the presence or absence of an alternating external field. We predict the characteristic time scale for the alignment of all defect dipoles with the spontanenous polarization of the surrounding matrix. In this state the defect dipoles act as obstacles for domain wall motion and contribute to the experimentally observed aging. Under cycling conditions the fully aligned configuration is perturbed and a dynamic equilibrium is established with defect dipoles in parallel and anti-parallel orientation relative to the spontaneous polarization. This process can be related to the deaging behavior of piezoelectric ceramics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Low-temperature criticality of martensitic transformations of Cu nanoprecipitates in α-Fe

Nanoprecipitates form during nucleation of multiphase equilibria in phase segregating multicomponent systems. In spite of their ubiquity, their size-dependent physical chemistry, in particular at the boundary between phases with incompatible topologies, is still rather arcane. Here we use extensive atomistic simulations to map out the size-temperature phase diagram of Cu nanoprecipitates in α-Fe. The growing precipitates undergo martensitic transformations from the body-centered cubic (BCC) phase to multiply-twinned 9R structures. At high temperatures, the transitions exhibit strong first-order character and prominent hysteresis. Upon cooling the discontinuities become less pronounced and the transitions occur at ever smaller cluster sizes. Below 300 K the hysteresis vanishes while the transition remains discontinuous with a finite but diminishing latent heat. This unusual size-temperature phase diagram results from the entropy generated by the soft modes of the BCC-Cu phase, which are stabilized through confinement by the α-Fe lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermodynamic and mechanical properties of copper precipitates in alpha-iron from atomistic simulations

The Fe-Cu system has attracted much attention over the last several decades due to its technological importance as a model alloy for Cu steels. In spite of these efforts several aspects of its phase diagram remain unexplained. Here we use atomistic simulations to characterize the polymorphic phase diagram of Cu precipitates in body-centered cubic (BCC) Fe and establish a consistent link between their thermodynamic and mechanical properties in terms of thermal stability, shape, and strength. The size at which Cu precipitates transform from BCC to a close-packed 9R structure is found to be strongly temperature dependent, ranging from approximately 4 nm in diameter (~2,700 atoms) at 200 K to about 8 nm (~22,800 atoms) at 700 K. These numbers are in very good agreement with the interpretation of experimental data given Monzen et al. [Phil. Mag. A 80, 711 (2000)]. The strong temperature dependence originates from the entropic stabilization of BCC Cu, which is mechanically unstable as a bulk phase. While at high temperatures the transition exhibits first-order characteristics, the hysteresis, and thus the nucleation barrier, vanish at temperatures below approximately 300\,K. This behavior is explained in terms of the mutual cancellation of the energy differences between core and shell (wetting layer) regions of BCC and 9R nanoprecipitates, respectively. The proposed mechanism is not specific for the Fe--Cu system but could generally be observed in immiscible systems, whenever the minority component is unstable in the lattice structure of the host matrix. Finally, we also study the interaction of precipitates with screw dislocations as a function of both structure and orientation. The results provide a coherent picture of precipitate strength that unifies previous calculations and experimental observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A first-principles study of helium storage in oxides and at oxide--iron interfaces

Density-functional theory calculations based on conventional as well as hybrid exchange-correlation functionals have been carried out to study the properties of helium in various oxides (Al2O3, TiO2, Y2O3, YAP, YAG, YAM, MgO, CaO, BaO, SrO) as well as at oxide-iron interfaces. Helium interstitials in bulk oxides are shown to be energetically more favorable than substitutional helium, yet helium binds to existing vacancies. The solubility of He in oxides is systematically higher than in iron and scales with the free volume at the interstitial site nearly independently of the chemical composition of the oxide. In most oxides He migration is significantly slower and He--He binding is much weaker than in iron. To quantify the solubility of helium at oxide-iron interfaces two prototypical systems are considered (Fe|MgO, Fe|FeO|MgO). In both cases the He solubility is markedly enhanced in the interface compared to either of the bulk phases. The results of the calculations allow to construct a schematic energy landscape for He interstitials in iron. The implications of these results are discussed in the context of helium sequestration in oxide dispersion strengthened steels, including the effects of interfaces and lattice strain.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A scalable parallel Monte Carlo algorithm for atomistic simulations of precipitation in alloys

We present an extension of the semi-grandcanonical (SGC) ensemble that we refer to as the variance-constrained semi-grandcanonical (VC-SGC) ensemble. It allows for transmutation Monte Carlo simulations of multicomponent systems in multiphase regions of the phase diagram and lends itself to scalable simulations on massively parallel platforms. By combining transmutation moves with molecular dynamics steps structural relaxations and thermal vibrations in realistic alloys can be taken into account. In this way, we construct a robust and efficient simulation technique that is ideally suited for large-scale simulations of precipitation in multicomponent systems in the presence of structural disorder. To illustrate the algorithm introduced in this work, we study the precipitation of Cu in nanocrystalline Fe.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Calculations of excess free energies of precipitates via direct thermodynamic integration across phase boundaries

We describe a technique for constraining macroscopic fluctuations in thermodynamic variables well-suited for Monte Carlo (MC) simulations of multiphase equilibria. In particular for multicomponent systems this amounts to a statistical ensemble that implements constraints on both the average composition as well as its fluctuations. The variance-constrained semi-grandcanonical (VC-SGC) ensemble allows for MC simulations, in which single-phase systems can be reversibly switched into multiphase equilibria allowing the calculation of excess free energies of precipitates of complex shapes by thermodynamic integration. The basic features as well as the scaling and convergence properties of this technique are demonstrated by application to an Ising model. Finally, the VC-SGC MC simulation technique is used to calculate alpha/alpha' interface free energies in Fe-Cr alloys as a function of orientation and temperature taking into account configurational, vibrational, and structural degrees of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Contributions of point defects, chemical disorder, and thermal vibrations to electronic properties of Cd(1-x)Zn(x)Te alloys

We present a first principles study based on density functional theory of thermodynamic and electronic properties of the most important intrinsic defects in the semiconductor alloy Cd(1-x)Zn(x)Te with x<0.13. The alloy is represented by a set of supercells with disorder on the Cd/Zn sublattice. Defect formation energies as well as electronic and optical transition levels are analyzed as a function of composition. We show that defect formation energies increase with Zn content with the exception of the neutral Te vacancy. This behavior is qualitatively similar to but quantitatively rather different from the effect of volumetric strain on defect properties in pure CdTe. Finally, the relative carrier scattering strengths of point defects, alloy disorder, and phonons are obtained. It is demonstrated that for realistic defect concentrations carrier mobilities are limited by phonon scattering for temperature above approximately 150 K.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Electronic structure of LaBr3 from quasi-particle self-consistent GW calculations

Rare-earth based scintillators in general and lanthanum bromide (LaBr_3) in particular represent a challenging class of materials due to pronounced spin-orbit coupling and subtle interactions between d and f states that cannot be reproduced by standard density functional theory (DFT). Here a detailed investigation of the electronic band structure of LaBr_3 using the quasi-particle self-consistent GW (QPscGW) method is presented. This parameter-free approach is shown to yield an excellent description of the electronic structure of LaBr_3. Specifically it is able to reproduce the band gap, the correct level ordering and spacing of the 4f and 5d states, as well as the spin-orbit splitting of La-derived states. The QPscGW results are subsequently used to benchmark several computationally less demanding techniques including DFT+U, hybrid exchange-correlation functionals, and the G_0W_0 method. Spin-orbit coupling is included self-consistently at each QPscGW iteration and maximally localized Wannier functions are used to interpolate quasi-particle energies. The QPscGW results provide an excellent starting point for investigating the electronic structure of excited states, charge self-trapping, and activator ions in LaBr_3 and related materials.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modeling the electrical conductivity in BaTiO3 on the basis of first-principles calculations

The dependence of the electrical conductivity on the oxygen partial pressure is calculated for the prototypical perovskite $\Ba\TiØ_3$ based on data obtained from first-principles calculations within density functional theory. The equilibrium point defect concentrations are obtained via a self-consistent determination of the electron chemical potential. This allows to derive charge carrier concentrations for a given temperature and chemical environment and eventually the electrial conductivity. The calculations are in excellent agreement with experimental data if an accidental acceptor dopant level of $10^{17}\,\cm^{-3}$ is assumed. It is shown that doubly charged oxygen vacancies are accountable for the high-temperature $n$-type conduction under oxygen-poor conditions. The high-temperature $p$-type conduction observed at large oxygen pressures is due to barium vacancies and titanium-oxygen di-vacancies under Ti and Ba-rich conditions, respectively. Finally, the connection between the present approach and the mass-action law approach to point defect thermodynamics is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermodynamics of mono and di-vacancies in barium titanate

The thermodynamic and kinetic properties of mono and di-vacancy defects in cubic (para-electric) barium titanate are studied by means of density-functional theory calculations. It is determined which vacancy types prevail for given thermodynamic boundary conditions. The calculations confirm the established picture that vacancies occur in their nominal charge states almost over the entire band gap. For the dominating range of the band gap the di-vacancy binding energies are constant and negative. The system, therefore, strives to achieve a state in which under metal-rich (oxygen-rich) conditions all metal (oxygen) vacancies are bound in di-vacancy clusters. The migration barriers are calculated for mono-vacancies in different charge states. Since oxygen vacancies are found to readily migrate at typical growth temperatures, di-vacancies can be formed at ease. The key results of the present study with respect to the thermodynamic behavior of mono and di-vacancies influence the initial defect distribution in the ferroelectric phases and therefore the conditions for aging.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Are there stable long-range ordered Fe(1-x)Cr(x) compounds?

The heat of formation of Fe-Cr alloys undergoes an anomalous change of sign at small Cr concentrations. This observation raises the question whether there are intermetallic phases present in this composition range. Here we report the discovery of several long-range ordered structures that represent ground state phases at zero Kelvin. In particular we have identified a structure at 3.7% Cr with an embedding energy which is 49 meV/Cr atom below the solid solution. This implies there is an effective long-range attractive interaction between Cr atoms. We propose that the structures found in this study complete the low temperature-low Cr region of the phase diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Association of oxygen vacancies with impurity metal ions in lead titanate

Thermodynamic, structural and electronic properties of isolated copper and iron atoms as well as their complexes with oxygen vacancies in tetragonal lead titanate are investigated by means of first principles calculations. Both dopants exhibit a strong chemical driving force for the formation of M_Ti-V_O (M = Cu, Fe) defect associates. The most stable configurations corresponds to a local dipole aligned along the tetragonal axis parallel to the spontaneous polarization. Local spin moments are obtained and the calculated spin densities are discussed. The calculations provide a simple and consistent explanation for the experimental findings. The results are discussed in the context of models for degradation of ferroelectric materials.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Calculation of the substitutional fraction of ion-implanted He in an Fe target

Ion-implantation is a useful technique to study irradiation damage in nuclear materials. To study He effects in nuclear fusion conditions, He is co-implanted with damage ions to reproduce the correct He/dpa ratios in the desired or available depth range. However, the short-term fate of these He ions, i.e over the time scales of their own collisional phase, has not been yet unequivocally established. Here we present an atomistic study of the short-term evolution of He implantation in an Fe substrate to approximate the conditions encountered in dual ion-implantation studies in ferritic materials. Specifically, we calculate the fraction of He atoms that end up in substitutional sites shortly after implantation, i.e. before they contribute to long-term miscrostructural evolution. We find that fractions of at most 3% should be expected for most implantation studies. In addition, we carry out an exhaustive calculation of interstitial He migration energy barriers in the vicinity of matrix vacancies and find that they vary from approximately 20 to 60 meV depending on the separation and orientation of the He-vacancy pair.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Extrinsic point defects in aluminum antimonide

We investigate thermodynamic and electronic properties of group IV (C, Si, Ge, Sn) and group VI (O, S, Se, Te) impurities as well as P and H in aluminum antimonide (AlSb) using first-principles calculations. To this end, we compute the formation energies of a broad range of possible defect configurations including defect complexes with the most important intrinsic defects. We also obtain relative scattering cross strengths for these defects to determine their impact on charge carrier mobility. Furthermore, we employ a self-consistent charge equilibration scheme to determine the net charge carrier concentrations for different temperatures and impurity concentrations. Thereby, we are able to study the effect of impurities incorporated during growth and identify optimal processing conditions for achieving compensated material. The key findings are summarized as follows. Among the group IV elements, C, Si, and Ge substitute for Sb and act as shallow acceptors, while Sn can substitute for either Sb or Al and displays amphoteric character. Among the group VI elements, S, Se, and Te substitute for Sb and act as deep donors. In contrast, O is most likely to be incorporated as an interstitial and predominantly acts as an acceptor. As a group V element, P substitutes for Sb and is electrically inactive. C and O are the most detrimental impurities to carrier transport, while Sn, Se, and Te have a modest to low impact. Therefore, Te can be used to compensate C and O impurities, which are unintentionally incorporated during the growth process, with minimal effect on the carrier mobilities.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Intrinsic point defects in aluminum antimonide

Calculations within density functional theory on the basis of the local density approximation are carried out to study the properties of intrinsic point defects in aluminum antimonide. Special care is taken to address finite-size effects, band gap error, and symmetry reduction in the defect structures. The correction of the band gap is based on a set of GW calculations. The most important defects are identified to be the aluminum interstitial $Al_{i,Al}^{1+}$, the antimony antisites $Sb_{Al}^0$ and $Sb_{Al}^{1+}$, and the aluminum vacancy $V_{Al}^{3-}$. The intrinsic defect and charge carrier concentrations in the impurity-free material are calculated by self-consistently solving the charge neutrality equation. The impurity-free material is found to be n-type conducting at finite temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Pressure-induced phase transition in the electronic structure of palladium nitride

We present a combined theoretical and experimental study of the electronic structure and equation of state (EOS) of crystalline PdN2. The compound forms above 58 GPa in the pyrite structure and is metastable down to 11 GPa. We show that the EOS cannot be accurately described within either the local density or generalized gradient approximations. The Heyd-Scuseria-Ernzerhof exchange-correlation functional (HSE06), however, provides very good agreement with experimental data. We explain the strong pressure dependence of the Raman intensities in terms of a similar dependence of the calculated band gap, which closes just below 11 GPa. At this pressure, the HSE06 functional predicts a first-order isostructural transition accompanied by a pronounced elastic instability of the longitudinal-acoustic branches that provides the mechanism for the experimentally observed decomposition. Using an extensive Wannier function analysis, we show that the structural transformation is driven by a phase transition of the electronic structure, which is manifested by a discontinuous change in the hybridization between Pd-d and N-p electrons as well as a conversion from single to triple bonded nitrogen dimers. We argue for the possible existence of a critical point for the isostructural transition, at which massive fluctuations in both the electronic as well as the structural degrees of freedom are expected.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Short-range order and precipitation in Fe-rich Fe-Cr alloys: Atomistic off-lattice Monte Carlo simulations

Short-range order (SRO) in Fe-rich Fe-Cr alloys is investigated by means of atomistic off-lattice Monte Carlo simulations in the semi-grand canonical ensemble using classical interatomic potentials. The SRO parameter defined by Cowley [Phys. Rev. B 77, 669 (1950)] is used to quantify the degree of ordering. In agreement with experiments a strong ordering tendency in the Cr distribution at low Cr concentrations (~< 5%) is observed, as manifested in negative values of the SRO parameters. For intermediate Cr concentrations (5% ~< c_Cr ~< 15%) the SRO parameter for the alpha-phase goes through a minimum, but at the solubility limit the alpha-phase still displays a rather strong SRO. In thermodynamic equilibrium for concentrations within the two-phase region the SRO parameter measured over the entire sample therefore comprises the contributions from both the alpha and alpha-prime phases. If both of these contributions are taken into account, it is possible to quantitatively reproduce the experimental results and interpret their physical implications. It is thereby shown that the inversion of the SRO observed experimentally is due to the formation of stable (supercritical) alpha-prime precipitates. It is not related to the loss of SRO in the alpha-phase or to the presence of unstable (subcritical) Cr precipitates in the alpha-phase.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Urbach tail in silica glass from first principles

We present density-functional theory calculations of the optical absorption spectra of silica glass for temperatures up to 2400 K. The calculated spectra exhibit exponential tails near the fundamental absorption edge that follow the Urbach rule, in good agreement with experiments. We also discuss the accuracy of our results by comparing to hybrid exchange correlation functionals. By deriving a simple relationship between the exponential tails of the absorption coefficient and the electronic density-of-states, we establish a direct link between the photoemission and the absorption spectra near the absorption edge. This relationship is subsequently employed to determine the lower bound to the Urbach frequency regime. Most interestingly, in this frequency interval, the optical absorption is Poisson distributed with very large statistical fluctuations. Finally, We determine the upper bound to the Urbach frequency regime by identifying the frequency at which transition to Poisson distribution takes place.

0 citations0 reviewsNo code linkedOpen access