Structure Low Dimensionality and Lone-Pair Stereochemical Activity: the Key to Low Thermal Conductivity in the Pb-Sn-S System.

Recently, metal sulfides have begun to receive attention as potential cost-effective materials for thermoelectric applications beyond optoelectronic and photovoltaic devices. Herein, based on a comparative analysis of the structural and transport properties of 2D PbSnS2 and 1D PbSnS3, we demonstrate that the intrinsic effects that govern the low lattice thermal conductivity (κL) of these sulfides originate from the combination of the low dimensionality of their crystal structures with the stereochemical activity of the lone-pair electrons of cations. The presence of weak bonds in these materials, responsible for phonon scattering, results in inherently low κL of 1.0 W/m K in 1D PbSnS3 and 0.6 W/m K in 2D PbSnS2 at room temperature. However, the nature of the thermal transport is quite distinct. 1D PbSnS3 exhibits a higher thermal conductivity with a crystalline-like peak at low temperatures, while 2D PbSnS2 demonstrates glassy thermal conductivity in the entire temperature range investigated. First-principles density functional theory calculations reveal that the presence of antibonding states below the Fermi level, especially in PbSnS2, contributes to the very low κL. In addition, the calculated phonon dispersions exhibit very soft acoustic phonon branches that give rise to soft lattices and very low speeds of sounds.

How inversion relates to disordering tendencies in complex oxides.

Complex oxides exhibit great functionality due to their varied chemistry and structures. They are quite flexible in terms of the ordering of cations, which can also impact their functional properties to a large extent. Thus, the propensity for a complex oxide to disorder is a key factor in optimizing and discovering new materials. Here, we show that the propensity to disorder cations in perovskites, pyrochlores, and spinels correlates with the energy to "invert" the structure - to directly swap the cations across the sublattices. This relatively simple metric, involving only two energetic calculations per compound, qualitatively captures disordering trends amongst compounds across these three families of materials and is quantitative in several cases. This provides a fast and robust metric to determine those complex oxides that are easy or hard to disorder, providing new avenues for quick screening of compounds for cation-ordering-dependent functionalities.

Adatom-Driven Oxygen Intermixing during the Deposition of Oxide Thin Films by Molecular Beam Epitaxy.

Thin film deposition from the vapor phase is a complex process involving adatom adsorption, movement, and incorporation into the growing film. Here, we present quantitative experimental data that reveals anion intermixing over long length scales during the deposition of epitaxial Fe2O3 and Cr2O3 films and heterostructures by oxygen-plasma-assisted molecular beam epitaxy. We track this diffusion by incorporating well-defined tracer layers containing 18O and/or 57Fe and measure their redistribution on the nanometer scale with atom probe tomography. Molecular dynamics simulations suggest potential intermixing events, which are then examined via nudged elastic band calculations. We reveal that adatoms on the film surface act to "pull up" subsurface O and Fe. Subsequent ring-like rotation mechanisms involving both adatom and subsurface anions then facilitate their mixing. In addition to film deposition, these intermixing mechanisms may be operant during other surface-mediated processes such as heterogeneous catalysis and corrosion.

The impact of chemical order on defect transport in mixed pyrochlores.

Using temperature accelerated dynamics, an accelerated molecular dynamics method, we examine the relationship between composition and cation ordering and defect transport in the mixed pyrochlore Gd2(Ti1-xZrx)2O7, using the oxygen vacancy as a representative defect structure. We find that the nature of transport is very sensitive to the cation structure, transitioning, as a function of composition, from three-dimensional migration to two-dimensional to pseudo-one-dimensional to becoming essentially immobile before reverting back to three-dimensional as the Zr content is increased. The rates of migration are also affected by the cation structure in the various compositions. This behavior is driven by the connectivity of Ti polyhedra in the material, with more extensive networks of Ti ions leading to a greater ability of the vacancy to traverse the material. Our results indicate that the nature of transport is dictated by the cation structure of the material and that, conversely, the cation structure could be used to control transport and potentially other functionalities in mixed pyrochlores.

Role of Sink Density in Nonequilibrium Chemical Redistribution in Alloys.

Nonequilibrium chemical redistribution in open systems submitted to external forces, such as particle irradiation, leads to changes in the structural properties of the material, potentially driving the system to failure. Such redistribution is controlled by the complex interplay between the production of point defects, atomic transport rates, and the sink character of the microstructure. In this work, we analyze this interplay by means of a kinetic Monte Carlo (KMC) framework with an underlying atomistic model for the Fe-Cr model alloy to study the effect of ideal defect sinks on Cr concentration profiles, with a particular focus on the role of interface density. We observe that the amount of segregation decreases linearly with decreasing interface spacing. Within the framework of the thermodynamics of irreversible processes, a general analytical model is derived and assessed against the KMC simulations to elucidate the structure-property relationship of this system. Interestingly, in the kinetic regime where elimination of point defects at sinks is dominant over bulk recombination, the solute segregation does not directly depend on the dose rate but only on the density of sinks. This model provides new insight into the design of microstructures that mitigate chemical redistribution and improve radiation tolerance.

Efficient Ab initio Modeling of Random Multicomponent Alloys.

We present in this Letter a novel small set of ordered structures (SSOS) method that allows extremely efficient ab initio modeling of random multicomponent alloys. Using inverse II-III spinel oxides and equiatomic quinary bcc (so-called high entropy) alloys as examples, we demonstrate that a SSOS can achieve the same accuracy as a large supercell or a well-converged cluster expansion, but with significantly reduced computational cost. In particular, because of this efficiency, a large number of quinary alloy compositions can be quickly screened, leading to the identification of several new possible high-entropy alloy chemistries. The SSOS method developed here can be broadly useful for the rapid computational design of multicomponent materials, especially those with a large number of alloying elements, a challenging problem for other approaches.

The role of surfaces, chemical interfaces, and disorder on plutonium incorporation in pyrochlores.

Pyrochlores, a class of complex oxides with formula A2B2O7, are one of the candidates for nuclear waste encapsulation, due to the natural occurrence of actinide-bearing pyrochlore minerals and laboratory observations of high radiation tolerance. In this work, we use atomistic simulations to determine the role of surfaces, chemical interfaces, and cation disorder on the plutonium immobilization properties of pyrochlores as a function of pyrochlore chemistry. We find that both Pu(3+) and Pu(4+) segregate to the surface for the four low-index pyrochlore surfaces considered, and that the segregation energy varies with the chemistry of the compound. We also find that pyrochlore/pyrochlore bicrystals A2B2O7/A2'B2'O7 can be used to immobilize Pu(3+) and Pu(4+) either in the same or separate phases of the compound, depending on the chemistry of the material. Finally, we find that Pu(4+) segregates to the disordered phase of an order/disorder bicrystal, driven by the occurrence of local oxygen-rich environments. However, Pu(3+) is weakly sensitive to the oxygen environment, and therefore only slightly favors the disordered phase. This behavior suggests that, at some concentration, Pu incorporation can destabilize the pyrochlore structure. Together, these results provide new insight into the ability of pyrochlore compounds to encapsulate Pu and suggest new considerations in the development of waste forms based on pyrochlores. In particular, the phase structure of a multi-phase pyrochlore composite can be used to independently getter decay products based on their valence and size.

The effects of cation-anion clustering on defect migration in MgAl2O4.

Magnesium aluminate spinel (MgAl2O4), like many other ceramic materials, offers a range of technological applications, from nuclear reactor materials to military body armor. For many of these applications, it is critical to understand both the formation and evolution of lattice defects throughout the lifetime of the material. We use the Speculatively Parallel Temperature Accelerated Dynamics (SpecTAD) method to investigate the effects of di-vacancy and di-interstitial formation on the mobility of the component defects. From long-time trajectories of the state-to-state dynamics, we characterize the migration pathways of defect clusters, and calculate their self-diffusion constants across a range of temperatures. We find that the clustering of Al and O vacancies drastically reduces the mobility of both defects, while the clustering of Mg and O vacancies completely immobilizes them. For interstitials, we find that the clustering of Mg and O defects greatly reduces O interstitial mobility, but has only a weak effect on Mg. These findings illuminate important new details regarding defect kinetics relevant to the application of MgAl2O4 in extreme environments.

Irradiation-induced grain growth and defect evolution in nanocrystalline zirconia with doped grain boundaries.

Grain boundaries are effective sinks for radiation-induced defects, ultimately impacting the radiation tolerance of nanocrystalline materials (dense materials with nanosized grains) against net defect accumulation. However, irradiation-induced grain growth leads to grain boundary area decrease, shortening potential benefits of nanostructures. A possible approach to mitigate this is the introduction of dopants to target a decrease in grain boundary mobility or a reduction in grain boundary energy to eliminate driving forces for grain growth (using similar strategies as to control thermal growth). Here we tested this concept in nanocrystalline zirconia doped with lanthanum. Although the dopant is observed to segregate to the grain boundaries, causing grain boundary energy decrease and promoting dragging forces for thermally activated boundary movement, irradiation induced grain growth could not be avoided under heavy ion irradiation, suggesting a different growth mechanism as compared to thermal growth. Furthermore, it is apparent that reducing the grain boundary energy reduced the effectiveness of the grain boundary as sinks, and the number of defects in the doped material is higher than in undoped (La-free) YSZ.

Carbon-14 decay as a source of non-canonical bases in DNA.

BACKGROUND: Significant experimental effort has been applied to study radioactive beta-decay in biological systems. Atomic-scale knowledge of this transmutation process is lacking due to the absence of computer simulations. Carbon-14 is an important beta-emitter, being ubiquitous in the environment and an intrinsic part of the genetic code. Over a lifetime, around 50 billion (14)C decays occur within human DNA. METHODS: We apply ab initio molecular dynamics to quantify (14)C-induced bond rupture in a variety of organic molecules, including DNA base pairs. RESULTS: We show that double bonds and ring structures confer radiation resistance. These features, present in the canonical bases of the DNA, enhance their resistance to (14)C-induced bond-breaking. In contrast, the sugar group of the DNA and RNA backbone is vulnerable to single-strand breaking. We also show that Carbon-14 decay provides a mechanism for creating mutagenic wobble-type mispairs. CONCLUSIONS: The observation that DNA has a resistance to natural radioactivity has not previously been recognized. We show that (14)C decay can be a source for generating non-canonical bases. GENERAL SIGNIFICANCE: Our findings raise questions such as how the genetic apparatus deals with the appearance of an extra nitrogen in the canonical bases. It is not obvious whether or not the DNA repair mechanism detects this modification nor how DNA replication is affected by a non-canonical nucleobase. Accordingly, (14)C may prove to be a source of genetic alteration that is impossible to avoid due to the universal presence of radiocarbon in the environment.

Competing kinetics and he bubble morphology in W.

The growth process of He bubbles in W is investigated using molecular dynamics and parallel replica dynamics for growth rates spanning 6 orders of magnitude. Fast and slow growth regimes are defined relative to typical diffusion hopping times of W interstitials around the He bubble. Slow growth rates allow the diffusion of interstitials around the bubble, favoring the biased growth of the bubble towards the surface. In contrast, at fast growth rates interstitials do not have time to diffuse around the bubble, leading to a more isotropic growth and increasing the surface damage.

Structure and segregation of dopant-defect complexes at grain boundaries in nanocrystalline doped ceria.

Grain boundaries (GBs) dictate vital properties of nanocrystalline doped ceria. Thus, to understand and predict its properties, knowledge of the interaction between dopant-defect complexes and GBs is crucial. Here, we report atomistic simulations, corroborated with first principles calculations, elucidating the fundamental dopant-defect interactions at model GBs in gadolinium-doped and manganese-doped ceria. Gadolinium and manganese are aliovalent dopants, accommodated in ceria via a dopant-defect complex. While the behavior of isolated dopants and vacancies is expected to depend on the local atomic structure at GBs, the added structural complexity associated with dopant-defect complexes is found to have key implications on GB segregation. Compared to the grain interior, energies of different dopant-defect arrangements vary significantly at the GBs. As opposed to bulk, the stability of oxygen vacancy is found to be sensitive to the dopant arrangement at GBs. Manganese exhibits a stronger propensity for segregation to GBs than gadolinium, revealing that accommodation of dopant-defect clusters depends on the nature of dopants. Segregation strength is found to depend on the GB character, a result qualitatively supported by our experimental observations based on scanning transmission electron microscopy. The present results indicate that segregation energies, availability of favorable sites, and overall stronger binding of dopant-defect complexes would influence ionic conductivity across GBs in nanocrystalline doped ceria. Our comprehensive investigation emphasizes the critical role of dopant-defect interactions at GBs in governing functional properties in fluorite-structured ionic conductors.

Defect interactions with stepped CeO2/SrTiO3 interfaces: implications for radiation damage evolution and fast ion conduction.

Due to reduced dimensions and increased interfacial content, nanocomposite oxides offer improved functionalities in a wide variety of advanced technological applications, including their potential use as radiation tolerant materials. To better understand the role of interface structures in influencing the radiation damage tolerance of oxides, we have conducted atomistic calculations to elucidate the behavior of radiation-induced point defects (vacancies and interstitials) at interface steps in a model CeO2/SrTiO3 system. We find that atomic-scale steps at the interface have substantial influence on the defect behavior, which ultimately dictate the material performance in hostile irradiation environments. Distinctive steps react dissimilarly to cation and anion defects, effectively becoming biased sinks for different types of defects. Steps also attract cation interstitials, leaving behind an excess of immobile vacancies. Further, defects introduce significant structural and chemical distortions primarily at the steps. These two factors are plausible origins for the enhanced amorphization at steps seen in our recent experiments. The present work indicates that comprehensive examination of the interaction of radiation-induced point defects with the atomic-scale topology and defect structure of heterointerfaces is essential to evaluate the radiation tolerance of nanocomposites. Finally, our results have implications for other applications, such as fast ion conduction.

Spatially-varying inversion near grain boundaries in MgAl2O4 spinel.

Complex materials, containing multiple chemical species, often exhibit chemical disorder or inversion. Typically, this disorder is viewed as spatially homogeneous throughout the material. Here, we show, using a simple grain boundary in MgAl2O4 spinel, that this is not the case and that the level of inversion at the grain boundary plane is different than in the bulk. This has ramifications for the energetics of the boundary and how defects interact with it, as exemplified by the relative formation energy of vacancies. Using these results as motivation, we construct a simple model of inversion versus grain size that captures the salient behavior observed in experiments and allows us to extract inversion-relevant properties from those same experiments, suggesting that grain boundaries in the experimentally prepared material are essentially fully inverse. Together, these results highlight the role that microstructure plays on the inversion in the material.

Correction: Spatially-varying inversion near grain boundaries in MgAl2O4 spinel.

[This corrects the article DOI: 10.1039/D0RA00700E.].

First-principles study of fission product (Xe, Cs, Sr) incorporation and segregation in alkaline earth metal oxides, HfO(2), and the MgO-HfO(2) interface.

In order to close the nuclear fuel cycle, advanced concepts for separating out fission products are necessary. One approach is to use a dispersion fuel form in which a fissile core is surrounded by an inert matrix that captures and immobilizes the fission products from the core. If this inert matrix can be easily separated from the fuel, via e.g. solution chemistry, the fission products can be separated from the fissile material. We examine a surrogate dispersion fuel composition, in which hafnia (HfO(2)) is a surrogate for the fissile core and alkaline earth metal oxides are used as the inert matrix. The questions of fission product incorporation in these oxides and possible segregation behavior at interfaces are considered. Density functional theory based calculations for fission product elements (Xe, Sr, and Cs) in these oxides are carried out. We find smaller incorporation energy in hafnia than in MgO for Cs and Sr, and Xe if variation of charge state is allowed. We also find that this trend is reversed or reduced for alkaline earth metal oxides with large cation sizes. Model interfacial calculations show a strong tendency of segregation from bulk MgO to MgO-HfO(2) interfaces.

An Overview of Recent Standard and Accelerated Molecular Dynamics Simulations of Helium Behavior in Tungsten.

One of the most critical challenges for the successful adoption of nuclear fusion power corresponds to plasma-facing materials. Due to its favorable properties in this context (low sputtering yield, high thermal conductivity, high melting point, among others), tungsten is a leading candidate material. Nevertheless, tungsten is affected by the plasma and fusion byproducts. Irradiation by helium nuclei, in particular, strongly modifies the surface structure by a synergy of processes, whose origin is the nucleation and growth of helium bubbles. In this review, we present recent advances in the understanding of helium effects in tungsten from a simulational approach based on accelerated molecular dynamics, which emphasizes the use of realistic parameters, as are expected in experimental and operational fusion power conditions.

Multifidelity Information Fusion with Machine Learning: A Case Study of Dopant Formation Energies in Hafnia.

Cost versus accuracy trade-offs are frequently encountered in materials science and engineering, where a particular property of interest can be measured/computed at different levels of accuracy or fidelity. Naturally, the most accurate measurement is also the most resource and time intensive, while the inexpensive quicker alternatives tend to be noisy. In such situations, a number of machine learning (ML) based multifidelity information fusion (MFIF) strategies can be employed to fuse information accessible from varying sources of fidelity and make predictions at the highest level of accuracy. In this work, we perform a comparative study on traditionally employed single-fidelity and three MFIF strategies, namely, (1) Δ-learning, (2) low-fidelity as a feature, and (3) multifidelity cokriging (CK) to compare their relative prediction accuracies and efficiencies for accelerated property predictions and high throughput chemical space explorations. We perform our analysis using a dopant formation energy data set for hafnia, which is a well-known high-k material and is being extensively studied for its promising ferroelectric, piezoelectric, and pyroelectric properties. We use a dopant formation energy data set of 42 dopants in hafnia-each studied in six different hafnia phases-computed at two levels of fidelities to find merits and limitations of these ML strategies. The findings of this work indicate that the MFIF based learning schemes outperform the traditional SF machine learning methods, such as Gaussian process regression and CK provides an accurate, inexpensive and flexible alternative to other MFIF strategies. While the results presented here are for the case study of hafnia, they are expected to be general. Therefore, materials discovery problems that involve huge chemical space explorations can be studied efficiently (or even made feasible in some situations) through a combination of a large number of low-fidelity and a few high-fidelity measurements/computations, in conjunction with the CK approach.

Dissociated vacancies and screw dislocations in MgO and UO2: atomistic modeling and linear elasticity analysis.

Understanding the effect of dislocations on the mass transport in ionic ceramics is important for understanding the behavior of these materials in a variety of contexts. In particular, the dissociated nature of vacancies at screw dislocations, or more generally, at a wide range of low-angle twist grain-boundaries, has ramifications for the mechanism of defect migration and thus mass transport at these microstructural features. In this paper, a systematic study of the dissociated vacancies at screw dislocations in MgO is carried out. The important role of stress migration in the atomistic modeling study is identified. Another aspect of the current work is a rigorous treatment of the linear elasticity model. As a result, good agreement between the atomistic modeling results and the linear elasticity model is obtained. Furthermore, we demonstrate that the proposed vacancy dissociation mechanism can also be extended to more complicated ionic ceramics such as UO2, highlighting the generality of the mechanism.