Vanessa Meschkea,
Andrew Novick
a,
Jen Rogersb,
Claire Porter
a,
Remco Changb,
Thomas Proffen
c and
Eric S. Toberer
*a
aPhysics Department, Colorado School of Mines, Golden, CO 80401, USA. E-mail: etoberer@mines.edu
bDepartment of Computer Science, Tufts University, Medford, MA 02115, USA
cNeutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
First published on 2nd August 2024
Characterizing local bonding environments in complex materials is essential for understanding and optimizing their properties. Equally as important is the ability to predict local motifs as a function of synthesis conditions, enhancing chemists’ ability to design properties into materials. In this study, we present an approach to leverage statistical mechanics to generate temperature- and energy-informed ensemble averaged pair distribution functions (PDFs). This method, which we have named Thermodynamic Ensemble Averages of PDFs for Ordering and Transformations (TEAPOT), utilizes density functional theory (DFT) to relax supercells while incorporating energetic penalties for local order, enabling accurate and computationally efficient analysis of local structure. We apply this method to the neutron PDF measurements of the pseudobinary MnTe–GeTe (MGT) alloy, demonstrating its capability to resolve complex local distortions and chemical ordering. Our results reveal detailed insights into phase transformations and local distortions driven by Mn substitution. For compositions that globally present as rock salt, our analysis reveals that Ge coordination geometry is heavily impacted by synthesis temperature. We propose that high temperature synthesis conditions promote a lowered Ge polyhedra distortion, promoting high charge carrier mobility due to the alignment of local and global structure. Incorporating statistical mechanics and computation into experimental analysis thus guides synthesis of tailored local structure.
Experimentally, probing local structure has been challenging. Spectroscopic techniques such as X-ray absorption fine structure (XAFS) and X-ray absorption near-edge structure (XANES) can provide valuable qualitative information about the first and second nearest neighbor chemistries, but quantitative descriptions in complex alloys are generally lacking.9 Alternatively, electron diffraction can probe the overall structure of a hyperlocal section of a sample, but fail to pick up on the nuances of local coordination. To measure both local chemistry and bond distributions, though, one technique particularly shines, especially for complex materials: the pair distribution function (PDF).10,11 The PDF provides valuable insight into both the local bonding environment and mesoscale structure for bulk materials that isn’t accessible with other techniques.12 PDF is thus, in principle, an excellent technique for probing the impacts of local order and distortion in complex materials. However, analyzing the results of PDF measurements in materials with significant local order and distortions remains challenging. In this paper, we develop a thermodynamics-aware first-principles approach to analyzing disorder and distortion and apply it to the high-perfoming thermoelectric alloy MnxGe1−xTe (MGT).13–17 MGT is not only interesting due to its high thermoelectric performance; it also offers a temperature and composition driven phase change that significantly alters the local bonding environment.18–21
When analyzing a PDF using common small box approaches, a useful first fitting attempt is to approximate perfect chemical disorder by creating atoms that are fractionally occupied according to the overall stoichiometry.22 This average structure representation is computationally inexpensive and, for highly disordered but undistorted materials, fits reasonably well. Due to its simplicity, the average structure is unable to capture any degree of chemical ordering; however, it can capture a small amount of distortion when used in conjunction with the box-car refinement strategy. By refining the average structure's atomic positions and lattice constants over small r-intervals, a box-car refinement can provide insight into the change in average bond distributions at different length scales. However, this approach is more useful in extracting information about the length scales of order and distortion rather than truly capturing local distortions.3,23,24 As such, the averaged structure approach will ultimately fail to capture local order and distortions in more complex systems.
To more accurately capture local distortions, the special quasirandom structures (SQS) method is an attractive alternative to the average structure for materials with minimal local order.25,26 In SQS, a single supercell that mimics perfect chemical disorder is generated, and local distortions are accounted for by relaxing the supercell with molecular dynamics (MD) or density functional theory (DFT).27–29 However, for high-dimensional alloy spaces, these relaxations become prohibitively expensive, making this method ineffective for complex chemistries. Additionally, the SQS will still fail to model local order due to its overestimation of entropic contributions.
To incorporate local order as well as distortions, the Reverse Monte Carlo (RMC) technique is very popular in PDF analysis.30,31 RMC minimizes a goodness of fit parameter between calculated and experimental PDFs by swapping or hopping atoms from an initial structure, with local distortions again being realized from piping results to MD or DFT.32 RMC alone has already been used with fair results in disordered materials such as MGT at the end of its cubic phase transition,18 while RMC + MD method has been particularly effective for modeling dynamic disorder in superionic conductors.33 Similarly to SQS, RMC again becomes exceedingly expensive as more elements are incorporated, since complex spaces require the use of large supercells. Additionally, while RMC allows for geometric constraints, it has no concept of the energetic penalty of the various structural or chemical motifs motifs, only fit quality. As such, any order captured with RMC risks being nonphysical and simply one of many possible solutions that can describe the PDF pattern well.
To account for both local order and distortion while considering energetic penalties, building a cluster expansion (CE) to use with a Monte Carlo simulation (MC) and first principles relaxation is an elegant workflow. With this method, the energetic penalties of local order are calculated when fitting the model Hamiltonian for the CE. Entropy contributions are then accounted for by annealing a starting structure using MC at a given temperature. Finally, local distortions are incorporated by using MD or DFT to relax the annealed structure from the CE + MC simulation.
Szymanski, et. al. demonstrate an excellent example of the CE + MC workflow by analyzing the local order and distortions in two disordered oxyfluoride alloys.29 In their work, Szymanski, et. al. explicitly show the degree of complexity that can be added to a model by separately investigating the effects of ordering and distortion using the various methods outlined in this introduction. Ultimately, the energy-informed CE + MC approach they use captures the experimentally measured PDFs better than any energy-unaware approximations to structure for their complex alloy system. However, approaches that rely on CE + MC have constraints. First, developing a CE can be challenging for systems that exhibit strong local distortions, including materials undergoing phase transitions. Second, MC-based annealing is necessary at each temperature of interest. Finally, the chemical complexity of the alloy is constrained by the need to ultimately relax the supercell with first principles calculations.
In the work discussed above, the goal of PDF analysis is to use computation to resolve local structure, but not necessarily provide guidance on how to experimentally proceed in designing tailored local structure. Indeed, the overarching philosophy is to find a singular representation of a structure that is instead best viewed statistically. Here, we present an alternative approach to PDF analysis that takes inspiration from the success of modeling alloys’ structures and properties through the independent cell approximation.34–44 By leveraging statistical mechanics, we are able to create thermodynamic ensemble averages of PDFs for ordering and transformations (TEAPOT). By developing a statistical ensemble of supercells relaxed with DFT, we minimize the number of PDF fitting parameters while remaining energy-aware. Further, the ensemble approach of TEAPOT allows for intensive variables such as temperature to be swept gracefully without requiring additional calculations of local structure. As such, having determined the thermodynamic state of a given sample, we are able to hypothesize how alternative synthetic strategies could be employed to impact local structure.
We present a use case for TEAPOT through a comprehensive study of neutron PDF measurements of MGT alloys. To begin, we develop a benchmark by invoking perfect chemical disorder in the PDF refinements. Following this initial fitting, we delve into the results of TEAPOT, highlighting the impact of Mn on the local coordination through the phase transformation. Given the large amount of data generated by this technique, we also demonstrate complex dimension reduction techniques to increase understanding of the local bonding environments. Through this study, we pave the way for a more thorough comprehension of complex alloy structures and showcase the potential of our approach to accelerate PDF analysis.
The total energies of the supercells were calculated using density functional theory48,49 (DFT) as implemented in the Vienna Ab Initio Simulation Package (VASP).50 The exchange-correction energy functional was approximated with the Perdew–Burke–Ernzerhof (PBE) functional51 within the projector-augmented wave formalism.52 A planewave cutoff of 340 eV and a Γ-centered k-point mesh were used such that the energy converged to <3 meV per atom. Each Mn atom was initialized with a starting spin of ±5μB, and spin configurations were chosen such that each supercell had a net spin of zero. Spin polarized relaxations were conducted such that the volume, cell shape, atom positions, and magnetic moments were all optimized.
Following the calculation of the supercell energies, thermodynamic averages of the PDFs were fit to the experimental data. First, structures were sorted into ending in rhombohedral or rock salt structure by assessing the calculated PDFs. Fewer than 5 specimens per composition were unable to be sorted into either rhombohedral or rock salt and were left out of the ensemble. The ensemble averaged PDF () is calculated as:
![]() | (1) |
![]() | (2) |
In calculating the probability of any PDF, Ei is the relaxed energy of supercell i, kb is the Boltzmann constant, and T is the temperature in K. The most significant portion of the probability is the normalization factor, Z, otherwise known as the partition function. The partition function for a given phase α is calculated as:
![]() | (3) |
As shown in Fig. S10 (ESI†), the individual PDFs have significant variation. To assess the uncertainty in the high temperature PDF for a given ensemble size, the central limit theorem was employed. The resulting uncertainty across the PDF is shown in Fig. S7 (ESI†) and is sufficiently small. While not employed herein, the uncertainty for other temperatures can be determined following the methods outlined by Novick, et al.40
The resulting structures from the DFT relaxation were converted to the CIFs using Pymatgen, and PDFs were calculated from these CIFs using PDFFit2.22,53 PDF refinements were carried out from 2–15 Å for each composition in the alloy space at ensemble calculation temperatures ranging from 100 to 10000 K. Before temperature was considered, each supercell's lattice was isotropically dilated or compressed to match the experimental pattern's peak locations. This dilation constant was then used to scale the supercell's lattice across all refinements at all temperatures. Following individual lattice refinements, temperature dependent ensemble refinements commenced. For these refinements, a singular isotropic U was refined for all supercells at each temperature considered. Additionally, the relative fractions of each phase (rhombohedral and cubic) were refined for the ensemble at each temperature considered. SciPy was used to perform a non-linear least squares refinement using the curve fit method with the Levenberg–Marquardt algorithm.54
In addition to calculating the ensemble PDFs, data were collected on the bond lengths and chemical coordination of individual atoms within the supercells. The coordination shells of all cations were collected using the Python library PyLaDa Light.55 All coordination shells were forced to have the same number of atoms as the corresponding coordination shells of a perfectly cubic structure (6, 12, 8, and 6 atoms for the first, second, third and fourth coordination shells, respectively).
Previous studies of the MGT phase diagram report a continuous rhombohedral to cubic phase transition as the Mn content is increased.57 This phase transition can be modeled as a growth of the ratios of the a and c lattice parameters, which vary from 4.165Å:
10.663 Å = 0.391 (rhombohedral) to 4.164
Å:
10.204 Å = 0.408 (rock salt). This transition is shown in Fig. 2(b). Previous studies have reported the alloy to be fully cubic for x ∼ 0.18 in MnxGe1−xTe from powder XRD, though x as high as 0.25 has been reported.19,20,57,58 Our experimental data shows agreement with the cubic phase transition being complete no earlier than x = 25%, and individual fits for each composition are shown in the ESI.† In Fig. 2(d), inspecting the peaks at ∼3 and 4.8 Å−1 show the merging of the doublet present for rhombohedral GeTe (bottom, black) to a single peak that is a hallmark of the rocksalt structure for cubic MGT. Additionally, by inspecting the a to c lattice parameter ratios, shown as filled circles in Fig. 2(c), a trend line can be approximated from the pure GeTe to x = 0.25 Mn content sample, showing the lattice parameters changing continuously throughout the phase transition until the rock salt a:c ratio (0.408) is reached.
In addition to the lattice transition, the atomic site positions move from the distorted, non-symmetric rhombohedral positions (0.236 and 0.762 for (Mn,Ge) and Te, respectively) to high symmetry 0.25 and 0.75 positions, shown as the solid squares in Fig. 2(f). The shifting atomic sites fully completes the transition from rhombohedral to rock salt and merges all peaks to the expected space group 225 pattern by the end of the alloy series studied, with a similar slowing rate of change as x ∼ 25% is reached. In contrast to the lattice parameters, the atomic site positions show continued evolution after x = 25%. In summary, we see that the rhombohedral to cubic phase transition occurs slowly and continuously until 25% Mn, at which point atomic positions and lattice parameters largely converge to the expected cubic diffraction pattern.
Despite the longer range structure being well fit across the sample space, the modeled nearest neighbors peak appears much less asymmetric than experiment suggests. The extracted a:
c lattice parameters ratio and fractional atomic positions from the PDF refinements show interesting results. At a low-r box car range, increasing the Mn content results in increased disagreement for both parameters as compared to the diffraction pattern. However, good agreement between diffraction and PDF analysis is found for all boxcar ranges beyond 2–15 Å and is shown in ESI,† Fig. S9.
Beyond the local coordination range, the averaged structures show varying fit quality from the local to the long range order. A summary of the fit quality over all selected box-car r ranges for all compositions is shown in the ESI.† Overall, the 25% Mn sample shows the best fit quality at all r ranges inspected, while the 35% Mn sample shows the worst. This is likely due to the degree of disorder present in the system at each of the compositions. The highest quality fits for all compositions lie in the 20–30 Å range, indicating the alloy shows local orderings that are not present at the intermediate range.
As discussed in the Methods section, 50 supercells each were generated in the cubic and rhombohedral structures for each stoichiometry studied and relaxed using DFT. While this approach will not be able to discern information regarding the widths of the energetic wells of the possible structures in an alloy space, it is significantly less computationally intensive than fully sampling configuration space.59 No constraints were placed on the relaxations to prevent bias towards either end of the structural transition. Given the fraught nature of space group classification, relaxed structures were then sorted based on the shape of the first nearest neighbors peak of the calculated PDF. Supercells with a doublet as the first nearest neighbors peak were sorted into the rhombohedral category, while single peaks were considered cubic. Fewer than five structures at each stoichiometry were unable to be sorted into either category and were discarded. The sorted structures energies’ were then used to create two structure-based TDOS for each composition, shown in Fig. 3, that determine the probability of a structure.
For the MGT system, the TDOS for the rhombohedral (green) and cubic (teal) structures are shown for each composition in Fig. 3. All distributions were referenced to the lowest energy for that composition and structure to provide a consistent reference point for comparison in the plots. Within these distributions, the 15% and 25% Mn samples have particularly heavy samplings of the lowest energy states, particularly on the cubic side of the distribution. The remaining samples have relatively Gaussian samplings across the energy range. The spread of the energies within any TDOS is relatively tight, spanning less than 10 meV per atom from the most to least favorable configuration. While this energy spread is larger than the error of a DFT calculation, it is relatively small and may lead to issues with distinguishing between energetic configurations in this space.
Incorporating the ensemble temperature via eqn (1)–(3) turns the above energy distributions into probability distributions. With Boltzmann averaging, the prevalence of a single representation is dictated by the supercell energies and temperature, as discussed in the Methods. An example of the relationship between an individual representation's ensemble probability and temperature is shown in Fig. 4(a). At low temperatures, only the few lowest energy structures from the ensemble contribute to the average structure significantly. For the rhombohedral and cubic structures, this manifests as the singular green and teal curves with probabilities near unity at 100 K. At high temperatures, however, all structures contribute equally to the structure with 1/n weight, where n is the number of samples in the ensemble. Thus, the high temperature limit approximates perfect disorder similar to the average structure or SQS technique, indicating high degrees of disorder in the system. In Fig. 4(a), this limit is approached near 7000 K and 10000 K for the rhombohedral and cubic ensembles, respectively. The temperature dependent distributions for the remaining compositions are shown in the ESI,† in Fig. S8.
The impact of temperature is realized as changes in the TEAPOT PDF, an example of which is shown in Fig. 4(b) for the 35% Mn ensemble. The averaged PDFs for Mn0.35Ge0.65Te were calculated from 100 K (blue) to 10000 K (red) with all calculation parameters were kept constant. For this ensemble, the temperature mostly changes peak heights and broadness, but not peak location. This strong variation is imperative for the efficacy of the TEAPOT method. Without the temperature dependence, any nuances of local coordination are unable to be distinguished from the high temperature limit, providing little information for guiding synthetic conditions to change alloy structure.
As such, considering the minimally achievable variation in TEAPOT PDFs is an important step in determining if the method can be successful in predicting structure changes. To quantify the potential impact of temperature on the ensemble structure, we consider two hypothetical situations. First, we construct a base dataset using the ensemble at 700 K. This pattern will be referred to as To. The degree of structural diversity compared to this base PDF can then be quantified by calculating the mean squared difference (MSD) between the PDF pattern at To and a PDF pattern at an arbitrary temperature T:
![]() | (4) |
![]() | ||
Fig. 6 Experimental PDF patterns (color) and TEAPOT fits (black) with difference curves below (grey) at best temperature fit from 2–15 Å show quality fits with Rw ranging from 0.09–0.12 across the entire alloy series. The fit qualities are shown for all temperatures in Fig. 7. |
Considering this fitting space, the refined phase fraction transitioned from predominantly rhombohedral to rock salt with increasing Mn content. This evolution is consistent with prior understanding of the composition-driven phase transition in the MGT chemical space. The optimal atomic displacement value was consistent across the composition space (U ∼ 0.01 Å2); we’ll return to this in the context of the bond lengths from the DFT-relaxed supercells. The optimal temperature for each composition (Fig. 7) varied from the low to the high temperature limit. The highest Mn content samples showed the strongest temperature dependence, with the high temperature limit being a clear best fit. In contrast, the low Mn samples showed minima at finite temperature, though these minima were less well defined than the high Mn samples. In the following paragraphs, we will delve into the dependency of each composition across the entire fitting space.
For each of the best fits, the upper end of the 2–15 Å box-car fit exceptionally well, though the quality of the first nearest neighbors peak is again lacking. However, previous studies have shown that the left side of this doublet is associated with inelastic scattering associated with correlated motion between the sheets.60 We are delighted to see that with increasing Mn content, this correlated motion rapidly fades away, likely due to Mn's propensity for 6-fold coordination.
The two most Mn-rich samples show the Rw curve flattens out for temperatures above 1000 K. This result is unsurprising when considering the variation in the TEAPOT PDFs considered in Fig. 5(b), where differences in the PDF are nearly imperceptible from 1000 K to 10000 K. As such, none of these samples have a singular best fit temperature, though each of the PDFs is nearly equivalent to the high temperature limit. However, the change in fit and fit quality from the low to high energy PDFs is very evident and is consistent with our expectation from Fig. 5. In the high Mn compositions, the combination of a low enthalpy of mixing (Fig. 3) and high configurational entropy support a randomly disordered cation sublattice.
Next considering the 25% Mn sample, this composition has a high degeneracy of low energy states (Fig. 3), and as such significant structural averaging has already occurred even at low temperature. The low temperature ensemble structure proves an excellent fit. However, a further improvement in fit occurs with temperature, with a shallow minimum centered at 428 K. Inspecting the variation in the fits with temperature, shown in Fig. S12–S16 (ESI†), the quality of this fit is largely dictated by changes in fourth nearest neighbor peak (∼6 A) intensities. The impact of such longer-range ordering is typically outside the scope of SQS and related modelling techniques. Little variation is observed in peak locations, so we attribute most of the ensemble variation to differences in the cation bonding environments.
For the 20% and 15% Mn content refinements, the fit quality is again quite good at all temperatures, though little variation is observed in the PDF below 10 Å (Fig S12 and S13, ESI†). This is consistent with the dilute nature of these compositions. Given that ordering at this length scale is unlikely to be a primary driver in the system energetics, the variation in Rw with temperature is unlikely to be predictive of local structure for these two samples. These samples highlight the importance of inspecting the temperature-dependence of the PDF rather than simply focusing on the optimal Rw, especially in more dilute systems.
A common trend throughout the composition space is the critical role that peak intensity, rather than peak location, played in determining the optimal temperature. Since all TEAPOT PDFs were modeled with the same number and ratio of cations in all supercells, changes in peak intensity are caused by changes in the number of Mn–Mn, Ge–Ge, and Mn–Ge interactions in the second coordination shells. Given the signs of the scattering lengths for Ge, Mn, and Te (8.185 fm, −3.73 fm, and 5.80 fm, respectively),61 the cation pairwise interactions can either result in a net positive (Ge–Ge, Mn–Mn) or negative (Mn–Ge, Ge–Mn) intensity. For example, the second and fourth nearest neighbor peaks at ∼4.2 and 6 Å show the opposite trend in peak intensity with increasing temperature (Fig. S12–S16, ESI†). This trend is consistent with the evolution of the Ge–Mn fractional coordination counts with temperature (Fig. S17 and S18, ESI†). The tightening of the bond distributions can also contribute to the increased peak intensity; we primarily see this in the first peak at ∼3 Å.
The energy-derived probabilities associated with each supercell provide direct insight into the statistics of bonding (Fig. 8). For comparison, we provide the bond lengths obtained using the average structure PDF fits from 2–15 Å as dashed lines. Beginning with low Mn content, Fig. 8 shows that the expected broad range (2.8–3.3 Å) of nearest neighbor bonds is present. However, the strong presence of ∼3 Å bonds in Ge-centered coordination shells indicates the presence of local octahedral motifs, even at this low Mn content. Compared the Ge-centered coordination shells, the distribution for Mn is significantly narrower. On the other compositional extreme, the Ge-centered bond length distribution continue to show a broad range of bond lengths, with peaks near 2.8 and 3.0 Å remaining through the structure transition. In contrast, the Mn distribution grows narrower as the Ge content decays. As such, we see the preference for Mn to retain its six-fold coordination driving the change in the overall structure from rhombohedral to cubic. While binary MnTe does not adopt the rock salt structure, the NiAs structure that it adopts is built from face sharing MnTe6 octahedra. The role of temperature is also considered in the bond length distributions of Fig. 8. Across all compositions, the largest impact of temperature is to smooth the bond length distributions and tighten the extrema inward. Little impact is seen on the average bond lengths, consistent with the temperature dependent PDFs shown in the ESI.†
Beyond the temperature dependent bond length distributions, the overall distortion of the coordination polyhedra can also be assessed. The degree of distortion of each polyhedron can be calculated using our bond distortion index, which is based off the quadratic elongation.62 We calculate the distortion as:
![]() | (5) |
Histograms of the temperature dependent distortions for the first nearest neighbors coordination shell are shown for each composition at each studied temperature and composition in Fig. 9. Even at 15%, the Mn polyhedra are, on average, close to a regular octahedra and the Ge polyhedra do not resemble the parent distortion value (solid line) or the average structure distortion value (dashed line). With increased Mn content, the Ge polyhedra become less distorted.
![]() | ||
Fig. 9 Distortion metric as defined in eqn (5) evolve across all compositions (panels) and temperatures (colors) for Ge-centered polyhedra (left) and Mn-centered (right). For comparison, the distortion values are also shown for pure rhombohedral GeTe (solid line, left panel) and the results from average structure refinements (dashed line). As Mn content increases, Ge-centered polyhedra distortions in the system unsurprisingly drop due to Mn's rigidity in its coordination preference, seen in each panel. |
Considering the temperature dependencies of the distortions shown in Fig. 9, Mn-centered polyhedra retain a regular octahedral coordination for all composition/temperature pairs. For Ge polyhedra in the Ge-rich samples, no clear trends emerge. However, at Ge-poor compositions, the Ge polyhedra show an inverse relationship between distortion and temperature. Chemically, this result is rather intuitive. Ge has previously been known to resist bonding in a perfect octahedron due to the presence of a lone pair that distorts its local bonding environment.63–65 Indeed, the distortion index distribution at low temperature is narrow and approximately the value of rhombohedral GeTe. At high temperatures, though, the entire configurational density of states (Fig. 3) is activated, thereby bringing in configurations that place the Ge into regular octahedral environments. We note that this is distinct from the impact of vibrations, as these calculations are relaxations, and thus the distortion index is strictly driven by next-nearest neighbor and beyond chemistry. In an effort to elucidate the driving force behind the temperature dependent trends, the UMAP projection method was used to find correlations in distortion metrics, partial PDFs, and coordination chemistry.66 However, selecting the appropriate featurization of the data proved rather challenging in the high dimensional alloy space, and no clear energy-dependent trends emerged from the projections.
In light of the TEAPOT results (Fig. 7), the Mn-rich samples considered here are best fit at the high temperature limit. As such, the fully random tiling of Ge and Mn lead to a preponderance of regular Ge octahedra. Using this result to intentionally design properties into this alloy space is slightly less intuitive. In one case, we might seek to maximize the electronic mobility to improve the thermoelectric efficiency. Here, annealing at high temperatures would more effectively align the local structure with the global rock salt structure; we hypothesize this will minimize charge carrier scattering and maximize the mobility. Alternatively, should a reduced lattice thermal conductivity be desired, we expect a lower annealing temperature would promote local site asymmetry and drive larger lattice anharmonicity. Thus, we expect that there is, unusually, anti-correlation between mobility and lattice thermal conductivity that can be tuned via synthetic temperature. In our Mn-rich samples, we expect they lie on one extreme based on the TEAPOT refinements shown in Fig. 7.
For the MGT system, we benchmark our thermodynamic average PDFs against the average structure method, showcasing our ability to accurately model local structures and to extract temperature-dependent knowledge about bonding environments. The knowledge extracted from these trends allows us to make informed extrapolations about synthetically controlled structures, highlighting the potential of this method for designing materials with tailored structures, and thus, properties. In conclusion, we show that thermodynamic averages offer an advancement in analyzing the PDFs of disordered, distorted materials, enabling accurate and efficient characterization of local structures.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4tc02896a |
This journal is © The Royal Society of Chemistry 2024 |