Search

Understanding molten salt chemistry is essential in ongoing research of the molten salt nuclear reactor (MSR). In this context, detailed understanding of the mechanisms underlying selective oxidation of metal species, such as Cr, is required to guide the design of effective corrosion mitigation strategies in molten salts. An important starting point for such mechanistic understanding is knowledge of the solvation structure and its role in controlling metal speciation. In this work, we use ab initio molecular dynamics simulations to study the short-range (on the scale of the nearest-neighbor bond lengths) and medium-range (over length scales of several neighbor spacings) structure in three different fluoride melts with and without Cr addition; namely, 2KF-NaF, 2LiF-BeF2, and 3LiF-AlF3. We find that Cr0,Cr2+,Cr3+ can each be coordinated by different numbers of F-, with the variance in coordination number decreasing as oxidation state increases, and that these coordination geometries are largely independent of solvent. The manner by which Cr changes the medium-range structure, however, is found to be solvent-dependent. While 2KF-NaF melts show short and medium range order that is highly dynamic, 2LiF-BeF2 and 3LiF-AlF3 are characterized by molecular associates that are relatively long-lived that organize into oligomer structures on larger length scales. Rather than being solvated by F- ions alone, we find that Cr can incorporate into and be solvated within this oligomer structure. Fluoroacidity, alone, may therefore prove too simple a metric for assessing the corrosivity of molten fluorides. As our work suggests, the ability of Cr to solvate must be understood in the context of the short- and medium-range structure of the solvent. more...

In traditional body-centered cubic (bcc) metals, the core properties of screw dislocations play a critical role in plastic deformation at low temperatures. Recently, much attention has been focused on refractory high-entropy alloys (RHEAs), which also possess bcc crystal structures. However, unlike face-centered cubic high-entropy alloys (HEAs), there have been far fewer investigations into bcc HEAs, specifically on the possible effects of chemical short-range order (SRO) in these multiple principal element alloys on dislocation mobility. Here, using density functional theory, we investigate the distribution of dislocation core properties in MoNbTaW RHEAs alloys, and how they are influenced by SRO. The average values of the core energies in the RHEA are found to be larger than those in the corresponding pure constituent bcc metals, and are relatively insensitive to the degree of SRO. However, the presence of SRO is shown to have a large effect on narrowing the distribution of dislocation core energies and decreasing the spatial heterogeneity of dislocation core energies in the RHEA. It is argued that the consequences of the mechanical behavior of HEAs is a change in the energy landscape of the dislocations, which would likely heterogeneously inhibit their motion. more...

Lead halide perovskites are strong candidates for high-performance low-cost photovoltaics, light emission, and detection applications. A hot-phonon bottleneck effect significantly extends the cooling time of hot charge carriers, which thermalize through carrier-optic phonon scattering, followed by optic phonon decay to acoustic phonons and finally thermal conduction. To understand these processes, we adjust the lattice dynamics independently of electronics by changing isotopes. We show that doubling the mass of hydrogen in methylammonium lead iodide by replacing protons with deuterons causes a large 20 to 50% softening of the longitudinal acoustic phonons near zone boundaries, reduces thermal conductivity by ~50%, and slows carrier relaxation kinetics. Phonon softening is attributed to anticrossing with the slowed libration modes of the deuterated molecules and the reduced thermal conductivity to lowered phonon velocities. Our results reveal how tuning the organic molecule dynamics enables control of phonons important to thermal conductivity and the hot-phonon bottleneck. more...

The formation of complex hierarchical nanostructures has attracted a lot of attention from both the fundamental science and potential applications point of view. Spherulite structures with radial fibrillar branches have been found in various solids; however, their growth mechanisms remain poorly understood. Here, we report real time imaging of the formation of two-dimensional (2D) iron oxide spherulite nanostructures in a liquid cell using transmission electron microscopy (TEM). By tracking the growth trajectories, we show the characteristics of the reaction front and growth kinetics. Our observations reveal that the tip of a growing branch splits as the width exceeds certain sizes (5.5–8.5 nm). The radius of a spherulite nanostructure increases linearly with time at the early stage, transitioning to nonlinear growth at the later stage. Furthermore, a thin layer of solid is accumulated at the tip and nanoparticles from secondary nucleation also appear at the growing front which later develop into fibrillar branches. The spherulite nanostructure is polycrystalline with the co-existence of ferrihydrite and Fe3O4 through-out the growth. A growth model is further established, which provides rational explanations on the linear growth at the early stage and the nonlinearity at the later stage of growth. [Figure not available: see fulltext.]. more...

While rhenium has proven to be an ideal material in fast-cycling high-temperature applications such as rocket nozzles, its prohibitive cost limits its continued use and motivates a search for viable cost-effective substitutes. We show that a simple design principle that trades off average valence electron count and cost considerations proves helpful in identifying a promising pool of candidate substitute alloys: The Mo-Ru-Ta-W quaternary system. We demonstrate how this picture can be combined with a computational thermodynamic model of phase stability, based on high-throughput ab initio calculations, to further narrow down the search and deliver alloys that maintain rhenium's desirable hcp crystal structure. This thermodynamic model is validated with comparisons to known binary phase diagram sections and corroborated by experimental synthesis and structural characterization demonstrating multiprinciple-element hcp solid-solution samples selected from a promising composition range. more...

Double perovskite halides are a class of materials with diverse chemistries that are amenable to solution-based synthesis routes, and display a range of properties for a variety of potential applications. Starting from a consideration of the octahedral and tolerance factors of ∼2000 candidate double perovskite compounds, we compute structural, electronic, and transport properties of ∼1000 using first-principles calculations based on density-functional-theory methods. The computational results have been assembled in a database that is accessible through the Materials Project online. As one potential application, double perovskites are candidates in the search for lead-free halide photovoltaic absorbers. We present the application of our database to aid the discovery of new double perovskite halide photovoltaic materials, by combining the results with optical absorption and phonon stability calculations. From three distinct classes of chemistries, 11 compounds were identified as promising solar absorbers and the complex chemical trends for band gap within each of these are analyzed, to provide guidelines for the use of substitutional alloying as a means of further tuning the electronic structure. Other possible applications of the database are also discussed. more...

Lead toxicity has sparked interest into alternative halide nanomaterials with properties similar to CsPbX 3 perovskites. A promising alternative suggested from bulk studies is the family of double perovskites of the form Cs 2 AgMX 6 . Here, we report the synthesis of colloidal Cs 2 AgInCl 6 and Cs 2 AgSbCl 6 nanocrystals via injection of acyl halides into a metal acetate solution under atmospheric conditions and relatively mild temperatures. We demonstrate the synthesis of single-crystalline cubic nanocrystals of ca. 10 nm side length and their morphological similarities to other double perovskite nanostructures in terms of their [200] facet termination and decoration with Ag (0) smaller nanocrystallites. To compare the stabilities of the synthesized materials, we develop a titration assay based on the degradation of nanocrystals with amines as a proxy for degradation by humidity, which provides a quantifiable stability metric. This measurement shows that Cs 2 AgSbCl 6 releases more than twice the decomposition energy compared to Cs 2 AgInCl 6 or CsPbCl 3 and degrades in the presence of approximately one molar equivalent of amine, whereas the other two materials require more than a 100-fold excess. Using facile chemical titration to quantitatively determine chemical stability provides an additional tool to aid in the basic understanding of what makes some of these materials more environmentally stable than others. more...

Metallic glasses (MGs) possess remarkably high strength but often display only minimal tensile ductility due to the formation of catastrophic shear bands. Purposely enhancing the inherent heterogeneity to promote distributed flow offers new possibilities in improving the ductility of monolithic MGs. Here, we report the effect of the spatial heterogeneity of elasticity, resulting from the inherently inhomogeneous amorphous structures, on the deformation behavior of MGs, specifically focusing on the ductility using multiscale modeling methods. A highly heterogeneous, Gaussian-type shear modulus distribution at the nanoscale is revealed by atomistic simulations in Cu64Zr36 MGs, in which the soft population of the distribution exhibits a marked propensity to undergo the inelastic shear transformation. By employing a mesoscale shear transformation zone dynamics model, we find that the organization of such nanometer-scale shear transformation events into shear-band patterns is dependent on the spatial heterogeneity of the local shear moduli. A critical spatial correlation length of elastic heterogeneity is identified for the simulated MGs to achieve the best tensile ductility, which is associated with a transition of shear-band formation mechanisms, from stress-dictated nucleation and growth to structure-dictated strain percolation, as well as a saturation of elastically soft sites participating in the plastic flow. This discovery is important for the fundamental understanding of the role of spatial heterogeneity in influencing the deformation behavior of MGs. We believe that this can facilitate the design and development of new ductile monolithic MGs by a process of tuning the inherent heterogeneity to achieve enhanced ductility in these high-strength metallic alloys. more...

Recent experimental and computational studies have demonstrated pressure and epitaxial stabilization of polar PbVO3 phases with perovskite-derivative crystal structures. In this study, we demonstrate, by density functional theory (DFT) calculations, the stability of similar perovskite-derivative structures in the KVO3 and NaVO3 systems when subjected to compressive biaxial strain. The electronic structure and polar properties of these compounds are computed as a function of biaxial strain, and the results are compared to those obtained for experimentally observed PbVO3 structures. It is demonstrated that the substitution of Pb with monovalent K or Na cations increases the strength of the vanadyl bond due to the removal of the spatially extended Pb 6p states. Both KVO3 and NaVO3 exhibit epitaxially stabilized perovskite-derivative phases having large polarizations and only small total energy increases relative to their unstrained bulk structures. The calculated epitaxial phase diagram for KVO3 predicts a strain-energy driving force for a phase separation from ∼% to 1.5% misfit strain into a polar Cm phase, having square-pyramidal coordination of the B-site, and a paraelectric Pbcm phase, having tetrahedral coordination of the B-site. The results show that strain-stabilized polar vanadate compounds may occur for other compositions in addition to PbVO3 and that changes in the A-site species can be used to tune bonding, structure, and functional properties in these systems. more...

As materials data sets grow in size and scope, the role of data mining and statistical learning methods to analyze these materials data sets and build predictive models is becoming more important. This manuscript introduces matminer, an open-source, Python-based software platform to facilitate data-driven methods of analyzing and predicting materials properties. Matminer provides modules for retrieving large data sets from external databases such as the Materials Project, Citrination, Materials Data Facility, and Materials Platform for Data Science. It also provides implementations for an extensive library of feature extraction routines developed by the materials community, with 47 featurization classes that can generate thousands of individual descriptors and combine them into mathematical functions. Finally, matminer provides a visualization module for producing interactive, shareable plots. These functions are designed in a way that integrates closely with machine learning and data analysis packages already developed and in use by the Python data science community. We explain the structure and logic of matminer, provide a description of its various modules, and showcase several examples of how matminer can be used to collect data, reproduce data mining studies reported in the literature, and test new methodologies. more...

At elevated temperatures, heavy liquid metals and their alloys are known to create a highly corrosive environment that causes irreversible degradation of most iron-based materials. It has been found that an appropriate concentration of oxygen in the liquid alloy can significantly reduce this issue by creating a passivating oxide scale that controls diffusion, especially if Al is present in Fe-based materials (by Al-oxide formation). However, the increase of the temperature and of oxygen content in liquid phase leads to the increase of oxygen diffusion into bulk, and to promotion of the internal Al oxidation. This can cause a strain in bulk near the oxide layer, due either to mismatch between the thermal expansion coefficients of the oxides and bulk material, or to misfit of the crystal lattices (bulk vs. oxides). This work investigates the strain induced into proximal bulk of a Fe-Cr-Al alloy by oxide layers formation in liquid lead-bismuth eutectic utilizing synchrotron X-ray Laue microdiffraction. It is found that internal oxidation is the most likely cause for the strain in the metal rather than thermal expansion mismatch as a two-layer problem. more...

Point defects have a strong impact on the performance of semiconductor and insulator materials used in technological applications, spanning microelectronics to energy conversion and storage. The nature of the dominant defect types, how they vary with processing conditions, and their impact on materials properties are central aspects that determine the performance of a material in a certain application. This information is, however, difficult to access directly from experimental measurements. Consequently, computational methods, based on electronic density functional theory (DFT), have found widespread use in the calculation of point-defect properties. Here we have developed the Python Charged Defect Toolkit (PyCDT) to expedite the setup and post-processing of defect calculations with widely used DFT software. PyCDT has a user-friendly command-line interface and provides a direct interface with the Materials Project database. This allows for setting up many charged defect calculations for any material of interest, as well as post-processing and applying state-of-the-art electrostatic correction terms. Our paper serves as a documentation for PyCDT, and demonstrates its use in an application to the well-studied GaAs compound semiconductor. We anticipate that the PyCDT code will be useful as a framework for undertaking readily reproducible calculations of charged point-defect properties, and that it will provide a foundation for automated, high-throughput calculations. Program summary: Program title: PyCDT Program Files doi: http://dx.doi.org/10.17632/7vzk5gxzh3.1 Licensing Provisions: MIT License. Programming language: Python External routines/libraries: NumPy [1], matplotlib [2], and Pymatgen [3], Nature of problem: Computing the formation energies and stable point defects with finite size supercell error corrections for charged defects in semiconductors and insulators. Solution method: Automated setup, and parsing of defect calculations, combined with local use of finite size supercell corrections. All combined into a code with a standard user-friendly command line interface that leverages a core set of tools with a wide range of applicability. Additional comments: This article describes version 1.0.0. Program obtainable from https://bitbucket.org/mbkumar/pycdt more...

The structural properties, energetics, and polarizations of perovskite-based thin-film oxide systems are computed as a function of biaxial strain state and epitaxial orientation, employing an automated computational workflow based on density functional theory. A total of 14 compositions are considered, of the form ABO3, with A = Ba, K, Na, Pb, and Sr and B = Hf, Sn, Ti, Zr, Nb, Ta, and V site cations chosen to yield tolerance factors with values ranging between 0.95 and 1.1. Three biaxial strain states corresponding to epitaxial growth of (100)-, (110)-, and (111)-oriented films are considered, with misfit strains ranging between -4% to 4%. Results are presented for the series of perovskite-derived phases, and their corresponding symmetries, which are energetically favorable as a function of misfit strain, along with their corresponding equilibrium atomic positions, lattice parameters, and electric polarizations. The results demonstrate robust trends of in-plane polarization enhancement under tensile strain for all epitaxial orientations, and out-of-plane polarization enhancement with compression for the (100)- and (110)-oriented films. Strains corresponding to the (111)-growth orientation lead to a wider variety of out-of-plane polarization behavior, with BaTiO3 showing anomalous diminishing polarization with compression. Epitaxial orientation is shown to have a strong effect on the nature of strain-induced phase transitions, with (100)-oriented systems tending to have smooth, second-order transitions and (110)- and (111)-oriented systems more commonly exhibiting first-order transitions. The significance of this effect for device applications is discussed, and a number of systems are identified as potentially interesting for ferroelectric thin-film applications based on energetic stability and polarization behavior. Analysis of polarization behavior across different orientations reveals distinct groups into which compositions can be organized, some of which have polarization dependencies on misfit strain that have not been reported previously. more...

Semiconductor compounds are widely used for photocatalytic hydrogen production applications, where photogenerated electron-hole pairs are exploited to induce catalysis. Recently, powders of a metallic oxide (Sr1-xNbO3, 0.033+δ and find that their bandgaps are ∼4.1 eV. Surprisingly, the carrier density of the conducting phase exceeds 1022 cm-3 and the carrier mobility is only 2.47 cm2 V-1 s-1. Contrary to earlier reports, the visible light absorption at 1.8 eV (∼688 nm) is due to the plasmon resonance, arising from the large carrier density. We propose that the hot electron and hole carriers excited via Landau damping (during the plasmon decay) are responsible for the photocatalytic property of this material under visible light irradiation. more...

We present a combination of machine learning and high throughput calculations to predict the points defects behavior in binary intermetallic (A-B) compounds, using as an example systems with the cubic B2 crystal structure (with equiatomic AB stoichiometry). To the best of our knowledge, this work is the first application of machine learning-models for point defect properties. High throughput first principles density functional calculations have been employed to compute intrinsic point defect energies in 100 B2 intermetallic compounds. The systems are classified into two groups: (i) those for which the intrinsic defects are antisites for both A and B rich compositions, and (ii) those for which vacancies are the dominant defect for either or both composition ranges. The data was analyzed by machine learning-techniques using decision tree, and full and reduced multiple additive regression tree (MART) models. Among these three schemes, a reduced MART (r-MART) model using six descriptors (formation energy, minimum and difference of electron densities at the Wigner-Seitz cell boundary, atomic radius difference, maximal atomic number and maximal electronegativity) presents the highest fit (98 %) and predictive (75 %) accuracy. This model is used to predict the defect behavior of other B2 compounds, and it is found that 45 % of the compounds considered feature vacancies as dominant defects for either A or B rich compositions (or both). The ability to predict dominant defect types is important for the modeling of thermodynamic and kinetic properties of intermetallic compounds, and the present results illustrate how this information can be derived using modern tools combining high throughput calculations and data analytics. more...

Multifunctional, complex oxides capable of exhibiting highly-coupled electrical, mechanical, thermal, and magnetic susceptibilities have been pursued to address a range of salient technological challenges. Today, efforts are focused on addressing the pressing needs of a range of applications and identifying, understanding, and controlling materials with the potential for enhanced or novel responses. In this prospective, we highlight important developments in theoretical and computational techniques, materials synthesis, and characterization techniques. We explore how these new approaches could revolutionize our ability to discover, probe, and engineer these materials and provide a context for new arenas where these materials might make an impact. more...

In this work, we investigated several titanates with lepidocrocite-type structures (general formula AxTi1-yMyO4, with A = Na and M = Li or Mg), having potential utility as anode materials for sodium-ion batteries. First-principles calculations were used to determine key battery metrics, including potential profiles, structural changes during sodiation, and sodium diffusion energy barriers for several compositions, and were compared to experimental results. Site limitations were found to be critical determinants of the gravimetric capacities, which are also affected both by the stacking arrangement of the corrugated layers and the identity of M (Li or Mg). To explain the experimentally observed lattice parameter changes observed as a function of the state of charge, it was necessary to assume the participation of water/solvent during the sodium intercalation process. Sodium diffusion barriers were also found to vary as a function of state of charge and diffusion direction, with a spread of 0.06-1.3 eV at low sodium contents, narrowing to 0.3-0.5 eV at higher sodium contents. Based on these results, strategies for selecting and improving the performance of these electrode materials are suggested. more...

Correction for ‘Computational and experimental investigation of TmAgTe2 and XYZ2 compounds, a new group of thermoelectric materials identified by first-principles high-throughput screening’ by Hong Zhu et al., J. Mater. Chem. C, 2015, 3, 10554–10565. more...

We present an overview and preliminary analysis of computed thermoelectric properties for more than 48:000 inorganic compounds from the Materials Project (MP). We compare our calculations with available experimental data to evaluate the accuracy of different approximations in predicting thermoelectric properties. We observe fair agreement between experiment and computation for the maximum Seebeck coefficient determined with MP band structures and the BoltzTraP code under a constant relaxation time approximation (R2 = 0.79). We additionally find that scissoring the band gap to the experimental value improves the agreement. We find that power factors calculated with a constant and universal relaxation time approximation show much poorer agreement with experiment (R2 = 0.33). We test two minimum thermal conductivity models (Clarke and Cahill-Pohl), finding that both these models reproduce measured values fairly accurately (R2 = 0.82) using parameters obtained from computation. Additionally, we analyze this data set to gain broad insights into the effects of chemistry, crystal structure, and electronic structure on thermoelectric properties. For example, our computations indicate that oxide band structures tend to produce lower power factors than those of sulfides, selenides, and tellurides, even under the same doping and relaxation time constraints. We also list families of compounds identified to possess high valley degeneracies. Finally, we present a clustering analysis of our results. We expect that these studies should help guide and assess future high-throughput computational screening studies of thermoelectric materials. more...

Intrinsically doped samples of YCuTe2 were prepared by solid state reaction of the elements. Based on the differential scanning calorimetry and the high temperature X-ray diffraction analyses, YCuTe2 exhibits a first order phase transition at ∼440 K from a low-temperature-phase crystallizing in the space group P3m1 to a high-temperature-phase in P3. Above the phase transition temperature, partially ordered Cu atoms become completely disordered in the crystal structure. Small increases to the Cu content are observed to favour the formation of the high temperature phase. We find no indication of superionic Cu ions as for binary copper chalcogenides (e.g., Cu2Se or Cu2Te). All investigated samples exhibit very low thermal conductivities (as low as ∼0.5 W m-1 K-1 at 800 K) due to highly disordered Cu atoms. Electronic structure calculations are employed to better understand the high thermoelectric efficiency for YCuTe2. The maximum thermoelectric figure of merit, zT, is measured to be ∼0.75 at 780 K for Y0.96Cu1.08Te2, which is promising for mid-temperature thermoelectric applications. more...

A new group of thermoelectric materials, trigonal and tetragonal XYZ2 (X, Y: rare earth or transition metals, Z: group VI elements), the prototype of which is TmAgTe2, is identified by means of high-throughput computational screening and experiment. Based on density functional theory calculations, this group of materials is predicted to attain high zT (i.e. ∼1.8 for p-type trigonal TmAgTe2 at 600 K). Among approximately 500 chemical variants of XYZ2 explored, many candidates with good stability and favorable electronic band structures (with high band degeneracy leading to high power factor) are presented. Trigonal TmAgTe2 has been synthesized and exhibits an extremely low measured thermal conductivity of 0.2-0.3 W m-1 K-1 for T > 600 K. The zT value achieved thus far for p-type trigonal TmAgTe2 is approximately 0.35, and is limited by a low hole concentration (∼1017 cm-3 at room temperature). Defect calculations indicate that TmAg antisite defects are very likely to form and act as hole killers. Further defect engineering to reduce such XY antisites is deemed important to optimize the zT value of the p-type TmAgTe2. more...

Point defects play an important role in determining the structural stability and mechanical behavior of intermetallic compounds. To help quantitatively understand the point defect properties in these compounds, we developed PyDII, a Python program that performs thermodynamic calculations of equilibrium intrinsic point defect concentrations and extrinsic solute site preferences in intermetallics. The algorithm implemented in PyDII is built upon a dilute-solution thermodynamic formalism with a set of defect excitation energies calculated from first-principles density-functional theory methods. The analysis module in PyDII enables automated calculations of equilibrium intrinsic antisite and vacancy concentrations as a function of composition and temperature (over ranges where the dilute solution formalism is accurate) and the point defect concentration changes arising from addition of an extrinsic substitutional solute species. To demonstrate the applications of PyDII, we provide examples for intrinsic point defect concentrations in NiAl. more...

Our program on capacitive energy storage is a comprehensive one that combines experimental and computational components to achieve a fundamental understanding of charge storage processes in redox-based materials, specifically transition metal oxides. Some of the highlights of this program are the identification of intercalation pseudocapacitance in Nb2O 5, which enables high energy density to be achieved at high rates, and the development of a new route for synthesizing mesoporous films in which preformed nanocrystal building blocks are used in combination with polymer templating. The resulting material architectures have large surface areas and enable electrolyte access to the redox active pore walls, while the interconnected mesoporous film provides good electronic conductivity. Select first-principles density-functional theory studies of prototypical pseudocapacitor materials are reviewed, providing insight into the key physical and chemical features involved in charge transfer and ion diffusion. Rigorous multiscale physical models and numerical tools have been developed and used to reproduce electrochemical properties of carbon-based electrochemical capacitors with the ultimate objective of facilitating the optimization of electrode design. For the organic photovoltaic (OPV) program, our focus has been ongoing beyond the trial-and-error Edisonian approaches that have been responsible for the increase in power conversion efficiency of blend-cast (BC) bulk heterojunction blends of polymers and fullerenes. Our first approach has been to use molecular self-assembly to create the ideal nanometer-scale architecture using thermodynamics rather than relying on the kinetics of spontaneous phase segregation. We have created fullerenes that self-assemble into one-dimensional stacks and have shown that use of these self-assembled fullerenes lead to dramatically enhanced OPV performance relative to fullerenes that do not assemble. We also have created self-assembling conjugated polymers that form gels based on electrically continuous cross-linked micelles in solution, opening the possibility for water-processable "green" production of OPVs based on these materials. Our second approach has been to avoid kinetic control over phase separation by using a sequential processing (SqP) technique to deposit the polymer and fullerene materials in separate deposition steps. The polymer layer is deposited first, using solvents and deposition conditions that optimize the polymer crystallinity for swelling and hole mobility. The fullerene layer is then deposited in a second step from a solvent that swells the polymer but does not dissolve it, allowing the fullerene to penetrate into the polymer underlayer to the desired degree. Careful comparison of composition- and thickness-matched BC and SqP devices shows that SqP not only produces more efficient devices but also leads to devices that behave more consistently. © 2014 American Chemical Society. more...

The energetics of mixing and defect ordering in solid solutions of fluorite-structured ThO2 with oxides of trivalent cations (Sc, In, Y, Nd, and La) are investigated by electronic density-functional theory (DFT). Through DFT calculations of structures enumerated by lattice-algebra techniques, we identify the lowest energy patterns for defect clustering for four separate dopant concentrations. The most stable structures are characterized by a repulsive interaction between nearest-neighbor vacancies on the oxygen sublattice. The enthalpies of formation with respect to constituent oxides are positive for all dopants considered, and show a tendency to decrease in magnitude as the size and electronegativity of the trivalent dopant decrease. Due to the small positive formation enthalpies and low oxygen-vacancy binding energy with La dopants, La2 O3 -ThO2 solid solutions are predicted to have relatively high ionic conductivities relative to those for the other aliovalent dopants considered. Our results are compared with those for the more widely studied ZrO2 and CeO2 fluorite-structured solid solutions with trivalent cations. © 2010 The American Physical Society. more...