Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional calculation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Comparing the performance of density functionals in describing the adsorption of atoms and small molecules on Ni(111)

Reaction networks calculated using density functional theory (DFT) play a crucial role in understanding and predicting heterogeneous catalyst performance. However, results obtained using different density functionals can vary significantly. Accordingly, it is important to assess the accuracy of available functionals in capturing adsorption energetics derived from single-crystal adsorption microcalorimetry experiments. Here, we evaluate how different density functionals describe the molecular and dissociative adsorption of CH 3 I, the adsorption of CH 3 , I, and H, and the dissociative adsorption of CH 4 on Ni(111). First, we identify energetically-preferred adsorbed states for each adsorbate at three surface coverages and calculate the corresponding enthalpies of adsorption at 160 K. Then, we compare the results to experimental measurements published in the literature. For each adsorption process we find at least one functional that is quantitatively accurate within experimental error, but no functional is accurate (within experimental error) for all adsorption processes studied. Even when assuming an additional DFT error of ±20 kJ/mol on top of the experimental error, only the PBE-D3 and RPBE-D3 functionals are accurate for all considered adsorption systems. Overall, these results indicate that quantitative agreement between density functional theory calculations and experimental measurements is both system- and functional-dependent. Importantly, single-crystal adsorption microcalorimetry experiments will continue to play an important role for benchmarking density functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Density functional and theoretical study of the temperature and pressure dependency of the plasmon energy of solids

The temperature and pressure dependency of the volume plasmon energy of solids was investigated by density functional theory calculations. The volume change of crystal is the major factor responsible for the variation of valence electron density and plasmon energy in the free electron model. Hence, to introduce the effect of temperature and pressure for the density functional theory calculations of plasmon energy, the temperature and pressure dependency of lattice parameter was used. Also, by combination of the free electron model and the equation of state based on the pseudo-spinodal approach, the temperature and pressure dependency of the plasmon energy was modeled. The suggested model is in good agreement with the results of density functional theory calculations and available experimental data for elements with the free electron behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Zirconia and hafnia polymorphs: Ground-state structural properties from diffusion Monte Carlo

Zirconia (zirconium dioxide) and hafnia (hafnium dioxide) are binary oxides used in a range of applications. Because zirconium and hafnium are chemically equivalent, they have three similar polymorphs, and it is important to understand the properties and energetics of these polymorphs. However, while density functional theory calculations can get the correct energetic ordering, the energy differences between polymorphs depend very much on the specific density functional theory approach, as do other quantities such as lattice constants and bulk modulus. We have used highly accurate quantum Monte Carlo simulations to model the three zirconia and hafnia polymorphs. We compare our results for structural parameters, bulk modulus, and cohesive energy with results obtained from density functional theory calculations. We also discuss comparisons of our results with existing experimental data, in particular for structural parameters where extrapolation to zero temperature can be attempted. We hope our results of structural parameters as well as for cohesive energy and bulk modulus can serve as benchmarks for density-functional theory based calculations and as a guidance for future experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Density functional theory investigation of the electronic structure and thermoelectric properties of layered MoS{sub 2}, MoSe{sub 2} and their mixed-layer compound

First principles density functional theory calculations were carried out for the 2H-MoQ{sub 2} (Q=S and Se) and their hypothetical mixed-layer compound. Due to the different electronegativities of S and Se atoms on MoQ{sub 2}, the band gap size could be adjusted in mixed-layer compound MoS{sub 2}/MoSe{sub 2}. Also, the indirect band gap in pure MoQ{sub 2} compounds is changed to the pseudo direct band gap in mixed-layer MoS{sub 2}/MoSe{sub 2} which is similar to the monolayer compound. The layer mixing enhances the thermoelectric properties because of the increased density of states around the Fermi level and the decreased band gap size. Therefore, we suggest that this layer mixing approach should be regarded as a useful way to modulate their electronic structures and to improve their thermoelectric properties. -- Graphical abstract: On the basis of density functional calculations we predict that the mixed-layer compounds 2H-MoS{sub 2}/2H-MoSe{sub 2}, in which two different layers 2H-MoS{sub 2} and 2H-MoSe{sub 2}, have enhanced thermoelectric properties because of the increased density of states around the Fermi level and the decreased band gap size. Highlights: • We explored a way of improving TE properties of 2H-MoQ{sub 2} on DFT methods. • The mixed-layer compounds MoS{sub 2}/MoSe{sub 2} have enhanced thermoelectric properties. • This is caused by modulated electronic structure of mixed layer compound. • Layer mixing approach should be regarded as a useful way to improve TE properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Anion-mediated electronic effects in reducible oxides: Tuning the valence band of ceria via fluorine doping

Combining experimental spectroscopy and hybrid density functional theory calculations, we show that the incorporation of fluoride ions into a prototypical reducible oxide surface, namely, ceria(111), can induce a variety of nontrivial changes to the local electronic structure, beyond the expected increase in the number of Ce 3+ ions. Our resonant photoemission spectroscopy results reveal new states above, within, and below the valence band, which are unique to the presence of fluoride ions at the surface. With the help of hybrid density functional calculations, we show that the different states arise from fluoride ions in different atomic layers in the near surface region. In particular, we identify a structure in which a fluoride ion substitutes for an oxygen ion at the surface, with a second fluoride ion on top of a surface Ce 4+ ion giving rise to F 2p states which overlap the top of the O 2p band. The nature of this adsorbate F – –Ce 4+ resonant enhancement feature suggests that this bond is at least partially covalent. Here, our results demonstrate the versatility of anion doping as a potential means of tuning the valence band electronic structure of ceria.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Fermi–Löwdin orbital self-interaction correction of adsorption energies on transition metal ions

Density functional theory (DFT)-based descriptions of the adsorption of small molecules on transition metal ions are prone to self-interaction errors. Here, we show that such errors lead to a large over-estimation of adsorption energies of small molecules on Cu + , Zn + , Zn 2+ , and Mn + in local spin density approximation (LSDA) and Perdew, Burke, Ernzerhof (PBE) generalized gradient approximation calculations compared to reference values computed using the coupled-cluster with single, doubles, and perturbative triple excitations method. These errors are significantly reduced by removing self-interaction using the Perdew–Zunger self-interaction correction (PZ-SIC) in the Fermi–Löwdin Orbital (FLO) SIC framework. In the case of FLO-PBE, typical errors are reduced to less than 0.1 eV. Furthermore, analysis of the results using DFT energies evaluated on self-interaction-corrected densities [DFT(@FLO)] indicates that the density-driven contributions to the FLO-DFT adsorption energy corrections are roughly the same size in DFT = LSDA and PBE, but the total corrections due to removing self-interaction are larger in LSDA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Equation of state for 304L stainless steel

Here, a new equation of state for 304L stainless steel has been constructed and entered in the Los Alamos National Laboratory SESAME database with material ID 4273. The equation of state is based on a multiphase approach. The solid and liquid phases were fitted to available experimental data for thermal expansion, heat capacity, and shock Hugoniot, as well as to data from density functional theory calculations of the cold curve and quantum molecular dynamics of the liquid regime. The electronic thermal contribution over the entire equation of state range is based on Tartarus (Green’s functions) average atom calculations. For this low-carbon stainless steel both the density functional and electronic thermal calculations are comprised of 19 wt% Cr, 10 wt% Ni, and the remainder Fe.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A -type antiferromagnetic order in the Zintl-phase insulator EuZn 2 P 2

Zintl phases, containing strongly covalently bonded frameworks with separate ionically bonded ions, have emerged as a critical materials family in which to couple magnetism and strong spin-orbit coupling to drive diverse topological phases of matter. Here we report the single-crystal synthesis, magnetic, thermodynamic, transport, and theoretical properties of the Zintl compound EuZn 2 P 2 that crystallizes in the anti-La 2 O 3 (CaAl 2 Si 2 ) P-3m1 structure, containing triangular layers of Eu 2+ ions. In-plane resistivity measurements reveal insulating behavior with an estimated activation energy of E g = 0.11eV. Specific heat and magnetization measurements indicate antiferromagnetic ordering at T N = 23K. Curie-Weiss analysis of in-plane and out of plane magnetic susceptibility from T = 150 to 300 K yields p eff = 8.61 for μ 0 H⊥c and p eff = 7.74 for μ 0 H//c, close to the expected values for the 4f 7 J = S = 7/2 Eu 2+ ion and indicative of weak anisotropy. Below T N , a significant anisotropy of χ ⊥ /χ // ≈ 2.3 develops, consistent with A-type magnetic order as observed in isostructural analogs and as predicted by the density functional theory calculations reported herein. The positive Weiss temperatures of θ W =19.2K for μ 0 H⊥c and θ W =41.9K for μ 0 H//c show a similar anisotropy and suggest competing ferromagnetic and antiferromagnetic interactions. Comparing Eu magnetic ordering temperatures across trigonal EuM 2 X 2 (M= divalent metal, X= pnictide) shows that EuZn 2 P 2 exhibits the highest ordering temperature, with variations in T N correlating with changes in expected dipolar interaction strengths within and between layers and independent of the magnitude of electrical conductivity. These results provide experimental validation of the crystochemical intuition that the cation Eu 2+ layers and the anionic (M 2 X 2 ) 2– framework can be treated as electronically distinct subunits, enabling further predictive materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Effect of Pressure on Valence and Structural Properties of YbFe 2 Ge 2 Heavy Fermion Compound A Combined Inelastic X-ray Spectroscopy, X-ray Diffraction, and Theoretical Investigation

The crystal structure and the Yb valence of the YbFe 2 Ge 2 heavy fermion compound was measured at room temperature and under high pressures using high-pressure powder X-ray diffraction and X-ray absorption spectroscopy via both partial fluorescence yield and resonant inelastic X-ray emission techniques. Furthermore, the measurements are complemented by first-principles density functional theoretical calculations using the self-interaction corrected local spin density approximation investigating in particular the magnetic structure and the Yb valence. While the ThCr 2 Si 2 -type tetragonal (I4/mmm) structure is stable up to 53 GPa, the X-ray emission results show an increase of the Yb valence from v = 2.72(2) at ambient pressure to v = 2.93(3) at ~9 GPa, where at low temperature a pressure-induced quantum critical state was reported.

Kumar, Ravhi S.↗
DOE OSTI2018

Comparative Study on the Charge-Ordered and Mixed-Valence Phases of LuFe{sub 2}O{sub 4} via First-Principles Calculations

The first-principles density functional theory calculations have been used to investigate the relation between magnetic and structural properties of LuFe{sub 2}O{sub 4} for both the charge-ordered (CO) and the mixed-valence (MV) phases and also to compare the electronic properties for both phases. The computed structural parameters of rhombohedral and monoclinic systems are found to be in good agreement with experimental findings. An analysis of the magnetic properties is provided considering nonmagnetic, ferromagnetic, antiferromagnetic, and ferrimagnetic model orderings. In agreement with experimental results, CO modification of LuFe{sub 2}O{sub 4} is predicted to be ferrimagnetic. The GGA + U electronic structures show a semiconducting behavior of CO phase while a metallic character for MV LuFe{sub 2}O{sub 4} phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Prospects for n -type doping of (Al x Ga 1– x ) 2 O 3 alloys

Here we systematically explore the properties of group-IV (C, Si, Ge, and Sn) and transition metal (Hf, Zr, and Ta) dopants substituting on the cation site in (Al x Ga 1– x ) 2 O 3 (AlGO) alloys using first-principles calculations with a hybrid functional. In Ga 2 O 3 , each of these dopants acts as a shallow donor. In Al 2 O 3 , they are deep defects characterized by the formation of either DX centers or positive- U (+/0) levels. Combining our calculations of dopant charge-state transition levels with information of the AlGO alloy band structure, we estimate the critical Al composition at which each dopant transitions from being a shallow to a deep donor. We identify Si to be the most efficient dopant to achieve n -type conductivity in high Al-content AlGO alloys, acting as a shallow donor over the entire predicted stability range for AlGO solid solution alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of defects in ultra-high gain in fast planar tin gallium oxide UV-C photodetector by MBE

Here we report ultra-high responsivity of epitaxial (Sn x Ga 1–x ) 2 O 3 (TGO) Schottky UV-C photodetectors and experimentally identified the source of gain as deep-level defects, supported by first principles calculations. Epitaxial TGO films were grown by plasma-assisted molecular beam epitaxy on (–201) oriented n-type β-Ga 2 O 3 substrates. Fabricated vertical Schottky devices exhibited peak responsivities as high as 3.5 ×10 4 A/W at –5 V applied bias under 250 nm illumination with sharp cutoff shorter than 280 nm and fast rise/fall time in milliseconds order. Hyperspectral imaging cathodoluminescence (CL) spectra were examined to find the mid-bandgap defects, the source of this high gain. Irrespective of different tin mole fractions, the TGO epilayer exhibited extra CL peaks at the green band (~2.20 eV) not seen in β-Ga 2 O 3 along with enhancement of the blue emission-band (~2.64 eV) and suppression of the UV emission-band. Based on hybrid functional calculations of the optical emission expected for defects involving Sn in β-Ga 2 O 3 , V Ga –Sn complexes are proposed as potential defect origins of the observed green and blue emission-bands. Such complexes behave as acceptors that can efficiently trap photogenerated holes and are predicted to be predominantly responsible for the ultra-high photoconductive gain in the Sn-alloyed Ga 2 O 3 devices by means of thermionic emission and electron tunneling. Regenerating the V Ga –Sn defect complexes by optimizing the growth techniques, we have demonstrated a planar Schottky UV-C photodetector of the highest peak responsivity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Unexpected dependence of the anomalous Hall angle on the Hall conductivity in amorphous transition metal thin films

The anomalous Hall effect (AHE), and magnetic and electronic transport properties were investigated in a series of amorphous transition metal thin films— Fe x S i –x , Fe x Ge 1-x , Co x Ge 1–x , Co x Si 1–x , and Fe 1–y Co y Si. The experimental results are compared with density functional theory calculations of the density of Berry curvature and intrinsic anomalous Hall conductivity. In all samples, the longitudinal conductivity (σ xx ), magnetization (M), and Hall resistivity (ρ xy ) increase with increasing transition metal concentration; due to the structural disorder σ xx is lower in all samples than a typical crystalline metal. In the systems with Fe as the transition metal (including Fe 1–y Co y Si), the magnetization and AHE are large and in some cases greater than the crystalline analog. In all samples, the AHE is dominated by the intrinsic mechanism, arising from a nonzero, locally derived Berry curvature. The anomalous Hall angle (AHA) (= σ xy /σ xx ) is as large as 5% at low temperature. These results are compared with the AHAs reported in a broad range of crystalline and amorphous materials. Previous work has shown that in a typical crystalline ferromagnet the Hall conductivity (σ xy ) and σ xx are correlated and are usually either both large or both small, resulting in an AHA that decreases with increasing σ xy . By contrast, the AHA increases linearly with increasing σ xy in the amorphous systems. This trend is attributed to a generally low σ xx , while σ xy varies and can be large. In the amorphous systems, σ xx and σ xy are not coupled, and there may thus exist the potential to further increase the AHA by increasing σ xy .

36 MATERIALS SCIENCE↗
DOE OSTI2017

Monte Carlo Simulation of Electron Thermalization in Scintillator Materials: Implications for Scintillator Nonproportionality

The lack of reliable quantitative estimates of the length and time scales associated with hot electron thermalization after a gamma-ray induced energy cascade obscures the interplay of various microscopic processes controlling scintillator performance and hampers the search for improved detector materials. We apply a detailed microscopic kinetic Monte Carlo model of the creation and subsequent thermalization of hot electrons produced by gamma irradiation of six important scintillating crystals to determine the spatial extent of the cloud of excitations produced by gamma rays and the time required for the cloud to thermalize with the host lattice. The main ingredients of the model are ensembles of microscopic track structures produced upon gamma excitation (including the energy distribution of the excited carriers), numerical estimates of electron-phonon scattering rates, and a calculated particle dispersion to relate the speed and energy of excited carriers. All these ingredients are based on first-principles density functional theory calculations of the electronic and phonon band structures of the materials. Details of the Monte Carlo model are presented along with results for thermalization time and distance distributions. These results are discussed in light of previous work. It is found that among the studied materials, calculated thermalization distances are positively correlated with measured nonproportionality. In the important class of halide scintillators, the particle dispersion is found to be more influential than the largest phonon energy in determining the thermalization distance.

Phonons, radiation detection, scintillators, thema↗
DOE OSTI2019

Silver clusters tune up electronic properties of graphene nanoflakes: A comprehensive theoretical study

Graphene-silver composites reflect a synergic effect of both components with enhanced electronic, adsorption and catalytic properties. We introduce here the detailed theoretical study on geometric, thermodynamic and electronic properties of silver clusters Ag n (n=2–10) adsorbed on graphene nanoflakes as a model for silver-graphene composites. First, a benchmark study is performed on Ag n …coronene complexes (n=2, 4) by evaluating various functional against reported high level ab-initio method for the selection of best functional for further calculations. Then, density functional theory study at the best chosen level is performed to inquire how clustering of adsorbate attached with carbon surface can impact geometric, electronic and thermodynamic properties of silver-graphene composites. Furthermore, the interaction energies are corrected for dispersion and basis set superposition error in order to get more realistic estimate of interaction energies. In progression, we investigated the dependence of interaction energy upon cluster size and dimensionality. For Ag n …coronene, the binding energies increase monotonically as the number of atoms in the cluster increases. However, for Ag n …dodecabenzocoronene, the interaction energies show a sudden drop at Ag 7 …dodecabenzocoronene. The difference in behavior between Ag n …dodecabenzocoronene and Ag n …coronene for interaction energy is attributed to edge effect present in coronene. Furthermore, coplanar orientations of the silver clusters on polyaromatic hydrocarbons have higher interaction energies. The 2D-3D transition of the silver cluster is observed for Ag 7 …coronene. CDA analysis reveals that backdonation from dodecabenzocoronene to silver clusters is more dominant in Ag 4 –Ag 7 whereas other clusters have dominant charge transfer to dodecabenzocoronene. Charge transfer analysis via constrained density functional theory (CDFT) reveals significant energy lowering from charge donor-acceptor interactions. Further energy decomposition analysis (EDA) using absolutely-localized molecular orbitals (ALMO) suggests that, while charge-transfer is important in these systems, the sum of electrostatics, Pauli repulsion, and London dispersion are far more significant as silver cluster size increases. NCI results also support the presence of van der Waals and electrostatic interactions. The hollow top is most favorite interaction sites compared to bridge top and head top binding sites over coronene lattice. The ionization potential, electron affinity, frontier molecular orbital analysis, chemical hardness, softness, chemical potential, and Fermi levels are deliberated to verify the stability and reactivity of most stable silver-graphene composites. The remarkable outcome of the current findings makes our silver-graphene composite as a potential applicant for high performance electronic and catalytic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Adsorption and activation of CO 2 on Pt/CeOx/TiO 2 (110): Role of the Pt-CeO x interface

The adsorption and dissociation of CO 2 on TiO 2 (110), CeO x /TiO 2 (110) and Pt/CeO x /TiO 2 (110) surfaces has been examined using Ambient Pressure X-ray Photoelectron Spectroscopy (AP-XPS). The substrates under study exhibited different degrees of complexity which were tested for the binding of the adsorbate and the cleavage of C-O bonds. The surfaces were prepared by depositing CeO x (0.1 ML) onto TiO 2 (110) to form a mixed oxide support, onto which Pt nanoparticles (0.2 ML) were deposited. This configuration yields a complex set of interfaces between metal and oxides and we have systematically titrated the active role of each component (Pt 4f, Ce 3d and Ti 2p regions) and the arising surface intermediates (C 1s and O 1s regions). CO 2 barely bonds to stoichiometric TiO 2 (110). It heals oxygen vacancies of this oxide surface (CO 2,gas → CO gas + O a ) and does not form stable carbonates. A stable carbonate was seen upon adsorption of CO 2 on CeO x /TiO 2 (110) and on this type of substrate the adsorbate also removed O vacancies leading to the oxidation of Ti 3+ and Ce 3+ sites. Pt nanoparticles dispersed on CeO x /TiO 2 (110) were highly effective for the binding and dissociation of CO 2 , with the formation of CO 3 , CO, C and CH x species on the Pt/CeO x /TiO 2 (110) system. The results of theoretical calculations based on density-functional theory (DFT) show that Pt/CeOx/TiO2(110) binds CO2 much stronger than surfaces of bulk platinum {(111), (100), (110)} or other late transition metals. On a Pt-CeO x interface, the molecule adsorbs with a bent configuration (~ 130° O-C-O bond angle) and with a substantial elongation (~ 0.1 Å) of the C-O bonds, facilitating its transformation into high value chemicals.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of Eu-Doping in the Electron Transport Behavior in the Zintl Thermoelectric Ca 5-x-y Yb x Eu y Al 2 Sb 6 System

A series of Eu-doped Zintl compounds belonging to theCa 5-x-y Yb x Eu y Al 2 Sb 6 (x = 0, 1.12; 0 ≤ y ≤ 0.63(2)) system have been successfully synthesized by both the arc-melting and the molten Pb-flux methods. All of the five title compounds initially crystallized in the Ca 5 Ga 2 As 6 -type phase (space group Pbam, Z = 2, Pearson code oP26) and maintained their original structure even after the post-heat treatment, unlike the recently reported n-type Zintl analogues in the Ca 5-x-y Yb x RE y Al 2 Sb 6 (RE = Pr, Nd, Sm) systems, which underwent a phase transition from the Ca 5 Ga 2 As 6 -type to the Ca 5 Al 2 Bi 6 -type phase after annealing. This research aimed to understand the origin of the structural preference of the title Ca 5-x-y Yb x Eu y Al 2 Sb 6 system, whether it was affected by the valence electron count or the cationic size. Electrical transport property measurements showed an increase in electrical conductivities and a decrease of Seebeck coefficients for Ca 4.89(1) Eu 0.11 Al 2 Sb 6 , Ca 4.82(1) Eu 0.18 Al 2 Sb 6 , and Ca 4.62(1) Eu 0.38 Al 2 Sb 6 , compared to the parental compound Ca 5 Al 2 Sb 6 . Hole effect measurements proved that these changes should be attributed to the reduced carrier concentration and enhanced carrier mobility. The comprehensive density functional theory calculations including electron density map analysis for the hypothetical model Ca 4.5 Eu 0.5 Al 2 Sb 6 revealed that the polarity between Al and Sb forming the anionic frameworks decreased as the Eu-dopants were introduced, which eventually affected the carrier mobility in the anionic frameworks. Thermal conductivity measurements proved that the Eu-doping successfully lowered the lattice thermal conductivity because of the enhanced atomic disordering. In conclusion, the magnetization measurements for Ca 4.37(2) Eu 0.63 Al 2 Sb 6 showed a typical Curie–Weiss behavior with weak antiferromagnetic nearest-neighbor interactions with θ p = -5.07 K.

36 MATERIALS SCIENCE↗