
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39232091
70500
10.1038/s41598-024-70500-9
Article
Chemical short-range order increases the phonon heat conductivity in a refractory high-entropy alloy
Mora-Barzaga Geraudys 1
Urbassek Herbert M. urbassek@rhrk.uni-kl.de

2
Deluigi Orlando R. 1
Pasinetti P. Marcelo 3
Bringa Eduardo M. 14
1 https://ror.org/01es6dz53 grid.441701.7 0000 0001 2163 0608 CONICET and Faculty of Engineering, University of Mendoza, Mendoza, 5500 Argentina
2 grid.519840.1 Physics Department, University of Kaiserslautern-Landau, Erwin-Schrödinger-Straße, 67663 Kaiserslautern, Germany
3 grid.412115.2 0000 0001 2309 1978 INFAP-CONICET and National University of San Luis, San Luis, 5700 Argentina
4 https://ror.org/00pn44t17 grid.412199.6 0000 0004 0487 8785 Center for Applied Nanotechnology, Faculty of Sciences, Universidad Mayor, Santiago, 8580745 Chile
4 9 2024
4 9 2024
2024
14 2062814 6 2024
16 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
We study the effects of the chemical short-range order (SRO) on the thermal conductivity of the refractory high-entropy alloy HfNbTaTiZr using atomistic simulation. Samples with different degrees of chemical SRO are prepared by a Monte Carlo scheme. With increasing SRO, a tendency of forming HfTi and TiZr clusters is found. The phonon density of states is determined from the velocity auto-correlation function and chemical SRO modifies the high-frequency part of the phonon density of states. Lattice heat conductivity is calculated by non-equilibrium molecular dynamics simulations. The heat conductivity of the random alloy is lower than that of the segregated binary alloys. Phonon scattering by SRO precipitates might be expected to reduce scattering times and, therefore, decrease thermal conductivity. We find that, in contrast, due to the increase of the conductivity alongside SRO cluster percolation pathways, SRO increases the lattice heat conductivity by around 12 %. This is expected to be a general result, extending to other HEAs.

Keywords

Vibrational density of states
High-entropy alloy
Heat conductivity
Chemical short-range order
Percolation
Subject terms

Condensed-matter physics
Information theory and computation
Atomistic models
http://dx.doi.org/10.13039/501100005961 Universidad Nacional de Cuyo UNCuyo 06/M008-T1 UNCuyo 06/M008-T1 UNCuyo 06/M008-T1 Mora-Barzaga Geraudys Urbassek Herbert M. Pasinetti P. Marcelo http://dx.doi.org/10.13039/501100021778 Agencia Nacional de Promoción de la Investigación, el Desarrollo Tecnológico y la Innovación PICTO-UUMM-2019-00048 PICTO-UUMM-2019-00048 PICTO-UUMM-2019-00048 Mora-Barzaga Geraudys Urbassek Herbert M. Pasinetti P. Marcelo http://dx.doi.org/10.13039/501100002923 Consejo Nacional de Investigaciones Científicas y Técnicas PIP 2021-2023 11220200102578CO PIP 2021-2023 11220200102578CO PIP 2021-2023 11220200102578CO Mora-Barzaga Geraudys Urbassek Herbert M. Pasinetti P. Marcelo Rheinland-Pfälzische Technische Universität Kaiserslautern-Landau (6375)Open Access funding enabled and organized by Projekt DEAL.

issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

High-entropy alloys (HEAs) are a class of metals with a high configurational entropy due to the large number of elements in their composition. Given the combination of different elements with different atomic radii, the HEA lattice is intrinsically distorted1, which in turn influences the structure and thermal properties of the alloy. HEAs are known for their superior mechanical and thermal properties such as a good mechanical performance at different temperature, radiation resistance, and fracture toughness2–5.

For a number of applications, the thermal conductivity of HEAs is relevant. Thermal energy is transported by phonons, electrons, and magnons. Generally, magnon conductivity is much lower than the others and can hence be neglected. The thermal conductivity of monoatomic metals and binary/ternary alloys is mainly due to electrons. However, a recent study6 finds that in transition metals and transition-intermetallic compounds the phononic contributions are non-negligible.

Due to chemical disorder, the electronic contribution to heat conductivity can decrease by orders of magnitude in HEAs7, reaching values similar to those of phonon contribution8. In the kinetic theory of gases, the thermal conductivity depends on the product of the group velocity and the mean free path of the carrier. The mean free path of the carriers decreases as a result of disorder. For instance, irradiation typically leads to point defects which cause a significant decrease of thermal conductivity with radiation dose9. In a HEA, chemical disorder leads to a significant distortion of the crystal lattice, even in the absence of point defects or other defects. Therefore, it has been proposed that the thermal conductivity of HEAs could be adjusted by varying the composition and elemental concentration of the alloy10. Computational simulations are an efficient tool in the search for optimized HEAs with desired properties. An important area of application is the design of thermoelectric materials, where the reduction of the phonon lattice heat conductivity is even a design principle and has been exploited for HEAs11.

The phonon part of the heat conduction is of intrinsic interest for studies of keV-ion induced collision cascades in solids, since here the projectile ion imparts its energy to the atomic motion (i.e., phonons) of the target material, and electrons only provide a small contribution to the heat conduction in later phases, after electron-phonon coupling conveyed energy to the electron system4,12,13. Reductions in thermal conductivity would lead to the cascade staying hot longer, favoring defect recombination and lowering the final defect content, resulting in enhanced radiation resistance.

Recent studies showed that the atomic composition in HEAs need not be random on an atomic scale, but chemical effects can generate a short-range order (SRO)14–16. This SRO is often associated with nm-sized ordered precipitates and influences many properties of the material17, including mechanical18–20 and thermal properties17,21–23. The presence of several alloying elements per se implies a relatively low phonon thermal conductivity24,25, but the presence of precipitates would generally be expected to reduce thermal conductivity even further25.

Although chemical complexity has been emphasized as an important factor in reducing thermal conductivity,26, studies of the phonon thermal conductivity of concentrated alloys are relatively rare, in contrast to elemental materials and binary alloys6,27,28. A notable exception is the model study by Caro et al.10 which showed how changes in atomic mass and force constants decrease the thermal conductivity of fcc binary alloys – and also fcc multi-component alloys – considerably below that of the pure elements. Cheng et al.29 investigated the phonon conductivity of fcc CrFeCoNiCu alloys as a function of temperature and found that it deviates increasingly strongly from that of elemental Cu with decreasing temperature. Jin et al.30 discussed the contributions of electrons and phonons to the heat conductivity in HEAs for the special example of fcc-structured Ni-based alloys and demonstrate that, in particular at low temperatures, phonons dominate the heat conductivity. Samples of (FeNiCrCo)1-xAlx have been studied as a function of Al content and temperature, but considering random samples, without the appropriate SRO that arises from the effective repulsion between Al-Al pairs31. Samples of (FeNiCrCo)0.7Al0.3 were simulated without SRO, but compared satisfactorily with experiments32. DFT simulation of the elastic constants coupled to analytical models to obtain the thermal conductivity was used to study 12 different bcc HEAs, finding that the relative phonon contribution can vary significantly33. Another ab-initio study found various behaviors in the contribution of electrons to the thermal conductivity of bcc HEAs, while the phonon conductivity was estimated with an analytical model requiring the Debye temperature34. Finally, experiments on both thermal conductivity and electrical resistivity were used to estimate the phonon contributions to thermal transport for several bcc HEAs below 200 K35. All the simulation studies above were for ‘perfect’ bulk crystals, without defects. Finally, a recent study36 investigated the effect of mass disorder on the thermal phonon conductivity, and found that correlations in the spatial mass distribution, i.e., mass SRO, decreased the heat conductivity.

Alternative approaches help the study of defective samples. Recently, the topology of a disordered network was related to its thermal conductivity37, including the generation of a dataset for prediction and tuning thermal transport. On the other hand, percolation theory has been applied to a variety of phenomena in physics, chemistry, biology, and materials science, where connectivity and clustering play an important role38,39. This theory has also provided insights into the behavior of more complicated models exhibiting phase transitions and critical phenomena. A recent work in this line is Ref.40 where the authors study the percolative phase transitions of a NbMoTaW HEA thin film during the process of filling the lattice on interstitial sites with small non-metal atoms, O or N. The phenomenon is explained by transforming a multi-color percolation into a two-color classical percolation problem where clusters with high and low electrical resistivity eventually form a percolation path contributing to a sudden reduction of the resistivity. They find a critical concentration of the conductive clusters similar to the percolative thresholds of the fcc and amorphous lattices.

Available analytical approaches to predict the phonon thermal conductivity of HEAs such as in Ref.33,41 are not able to include the chemical SRO. SRO in refractory HEAs is complex, and still not fully understood23. The low thermal conductivity of HfNbTaTiZr is desirable for thermal insulation technologies42. The low-temperature thermal conductivity, from nearly 0 K to 200 K was measured recently alongside the electrical resistivity, leading to estimated values of the phonon thermal conductivity similar to the ones due to the electronic thermal conductivity35. For high temperatures, there are measurements for a sample with 1 wt.% oxygen, giving a relatively low conductivity in the range 10–17 W/(Km) at 300–1200 K42.

In the present study, we focus on the effect of chemical SRO on the vibrational density of states and the phonon heat conductivity in a prototypical refractory HEA, HfNbTaTiZr, which crystallizes in the bcc structure3,43,44. Using an atomistic model of this alloy, we are able to study, besides the random alloy, various degrees of chemical SRO in the alloy. As a result, we show that chemical SRO can modify the heat conductivity in a sizable way.

Methods

In this section, the main methods used for performing and analyzing the simulations in this study are presented. These are the preparation of the random HfNbTaTiZr sample, the creation of SRO, methods for the sample analysis, models to assess the percolation of clusters in the sample, the calculation of the phonon density of states, and finally the calculation of the heat conductivity.

The molecular dynamics (MD) simulations are performed with the LAMMPS code45. We use the embedded-atom-method (EAM) type potential by Xu et al.46. This potential was developed explicitly in order to describe the short-range order in HfNbTaTi-based quinary refractory multi-principal element alloys such as we employ in the present study.

Sample preparation

As a first step, we create a random HfNbTaTiZr alloy. The single-crystalline body-centered cubic (bcc) sample has a cubic shape with {100} faces and an edge length of around 13.5 nm (40 lattice constants); it contains 128 000 atoms. We simulate this specimen using periodic boundary conditions.

In the random crystal, the atoms are distributed randomly over the lattice points with equal fractions of 20%. After relaxation to a final temperature of around 10 mK and low stress components around 10-4 GPa, it has a lattice constant of 3.389 Å, in good agreement with experimental data43,47,48.

SRO samples

The unrelaxed random sample was evolved using a Metropolis Monte Carlo (MC), with Kawasaki-type evolution, such that atoms are swapped49,50. The temperature for the Metropolis rejection factor was 100 K. 10 different swaps were attempted per MC step, one for each possible chemical pair AB, with A different from B. A similar scheme was recently employed to evolve a similar sample later studied under nanoindentation51. The simulation was run for 4.7e6 steps. The evolution of the potential energy during the Monte Carlo run is shown in the Supplementary Material (SM), Fig. S17, and demonstrates that after 4.7e6 steps, the sample has nearly reached equilibrium.

SRO is usually quantified using the Warren-Cowley (WC) parameters52,53,1 αij=1-pijcj,

to quantify the chemical SRO. Here, pij is the probability of finding an atom of type j in the first-neighbor shell of an atom of type i and cj is the (average) concentration of atoms of type j. Thus, for a random alloy it is αij=0; αij>0 characterizes atomic repulsion and while αij<0 denotes attraction.

At this low temperature, the chemical SRO of the alloy increases with MC steps until reaching a “steady” state, according to the evolution of both the potential energy and the Warren-Cowley parameters, at few million MC steps. Therefore, in what follows, the level of SRO is directly associated to the MC step number.

The WC parameters provide likely segregation pairs, but do not provide information about the structure of the aggregation. For bcc crystals, a common structure arising due to SRO is the B2 structure23, where an equiatomic binary alloy accommodates one chemical type in the center of the cube and the other type in the vertices.

Recent simulations of MC evolution of the HfNbTaTiZr HEA show the rise of such B2 structures for the most negative WC parameters51. In a B2 structure made with atoms A and B, atoms of type A only have neighbors of type B and vice versa. Therefore, to obtain clusters for a given B2 structure we select only atoms of type A which are linked to type B at the nearest neighbor distance, using the Create bonds tool in OVITO. Interconnected atoms form clusters and their cluster size can also be obtained with OVITO.

Sample analysis

Atomistic snapshots are generated with OVITO54. OVITO was also employed to obtain coordination, clusters, and solid volume fraction of selected clusters.

The volume fraction of SRO clusters can be evaluated in different manners. It can be obtained assuming that all atoms occupy the same volume in the regular bcc lattice used for the SRO MC simulations. Then the volume fraction is given by Na/Nt, where Na is the total number of atoms in the cluster, and Nt is the total number of atoms in the whole sample. An alternative estimation could be obtained from the ConstructSurfaceMesh tool54 in OVITO. In this method, a virtual probe sphere of radius r is rolled over the cluster surface to define a smooth geometric surface, and the volume inside that surface is calculated based on Delaunay tessellation of that volume. A spherical probe size of 0.29 nm is used, unless otherwise noted. For this case, to distinguish from the previous estimations using atom numbers, volume fractions are termed as cn, where n=AB refers to clusters the binary B2 alloy with atom types A and B.

SRO clusters are irregular. We could consider them as formed by interconnected filaments, and the freeware software FoamExplorer55 was employed to obtain an estimate of the size of those filaments. The wording “filament” is borrowed from bicontinuous foams, which form a single connected porous cluster, formed by interconnected filaments, and FoamExplorer has been used to characterize such structures. FoamExplorer works in the following way: surface atoms are identified according to their coordination, and then a linear “chord” is obtained, going from one surface atom to another, and the different chord lengths generate a histogram which represents the filament size distribution. A nanowire will give a very narrow histogram with the diameter of the nanowire, but irregular filaments will give some wider distributions. The following parameters were employed for FoamExplorer: for external (surface) atoms, a cutoff of 0.3 nm and a maximum coordination of 7. For internal atoms, a cutoff of 0.3 nm for first neighbors, cutoff of 0.49 nm for third neighbors, and maximum coordination at third neighbors of 26 were used.

Percolation theory and models

Usually, the percolation model in a lattice is classified into two categories, namely, the site model and the bond model. The percolation transition is a geometrical phase transition where the critical concentration of sites (or bonds) separates a phase of finite clusters from a phase where a macroscopic, spanning, or infinite cluster (in the thermodynamic limit) is present38. In our case we will assume a 3D site percolation model, considering that all atomic species occupy the same volume and that they are arranged on a regular bcc lattice.

There are five atomic species involved (multicolor percolation), however, as we are interested in studying the conductivity given by the HfTi association in a B2 structure, this association will be considered as the conductive (white) species, while any other species as non-conductive (black), reducing the problem to a two-color percolation40.

The substrate is initially generated completely at random. Then, a Kawasaki MC simulation brings the system closer to a condition of thermodynamic equilibrium, increasing the SRO of the system. This equilibrium is given at a certain temperature and keeping constant the number of atoms of each species (canonical ensemble). Therefore, this is not a strictly classical random percolation involving a random sequential adsorption process, which is essentially an out-of-equilibrium process, where each site has the same probability p of being occupied. Instead, this is a sort of thermal percolation, as in Giménez et al.56, where the probability for each site of being occupied depends on the temperature and the potential used in the Kawasaki exchange process. In our case, the percolation threshold, pc, will be associated with the number of MC steps necessary for the formation of a percolating cluster. It will be shown that the relevant percolating cluster in this case is HfTi with B2 structure.

At the percolation transition, the cluster size distribution of the system, n(s), is expected to have a power law behavior38,2 n(s)∝s-τ,

where τ is known as the Fisher exponent, taking the value τ≈2.19 for site percolation in the bcc lattice.

Phonon density of states

Previous research has underscored the significance of the vibrational density of states (VDOS) in evaluating the thermal characteristics of diverse nanostructured materials57–61. Molecular dynamics simulation allows to calculate easily the velocity autocorrelation function3 Z(t)=1N∑i=1Nvi(t)·vi(0)vi(0)·vi(0),

where vi(t) denotes the velocity vector of atom i at time t, and the summation is performed over all atoms within the system. The VDOS is then given by the Fourier transformation,4 G(ω)=∫0∞Z(t)cos(ωt)dt.

To obtain the dependency on the wavelength λ, the approximate relation λ=vs/ω can be used, where the sound velocity vs= 3.67 km/s can be obtained as vs=B/ρ, using the bulk modulus B = 134 GPa and the mass density ρ = 9.93 gcm-3 of the HfNbTaTiZr HEA as obtained by Alhafez et al.51.

Heat conductivity

There are several approaches to calculate the thermal conductivity62. For instance, both equilibrium and non-equilibrium approaches were shown to give similar results for amorphous C63. In this study, we took the approach of non-equilibrium molecular dynamics, that has frequently been used to calculate the thermal conductivity κ10,64. It obtains the thermal conductivity from Fourier’s law:5 κ=JΔzΔT,

where J is the heat flux and ΔT/Δz the temperature gradient across the system.

In detail, first a Langevin thermostat is applied for 0.1 ns to the entire crystal, so that it reaches the target temperature T=300 K. Then, two regions at opposite boundaries of the HEA, 0.7 nm thick, are fixed during the rest of the simulation. Adjacent to these fixed regions, the atoms in a strip of width 0.7 nm are subjected to Langevin thermostats to keep their temperatures at 0.8T and 1.2T, respectively. This induces a thermal gradient, ΔT/Δz in the system, which is applied for 0.6 ns. During the following 0.4 ns the thermodynamic properties of the system are calculated every 0.5 ps. This is in particular the heat Q which is calculated as the average between the energy added and subtracted to the hot and cold thermostat, respectively; the heat flux J is determined from Q by dividing by the lateral cross section of the crystal and the time during which the heat is delivered.

Results

In this section, the results of this study are presented. We discuss the chemical SRO established in the HfNbTaTiZr sample, the atomic displacements due to lattice distortion, the growth of clusters in the sample and their percolation as the SRO increases, the phonon density of states, and the thermal conductivity in the sample. Finally, models for the thermal conductivity and its temperature dependence are discussed.

Chemical SRO in HfNbTaTiZr

The evolution of the WC parameters with the number of MC steps is plotted in Fig. 1. Only the atomic pairs showing the largest deviations from 0 are shown. A strong tendency for attraction (α<0) is observed for several atom pairs, in particular for HfNb, HfTi, and TiZr; this attraction leads to the build-up of clusters ordered in the B2 structure. We note that this finding is in agreement with previous studies46,51. Also, Ta-Ta shows a strongly negative Warren-Cowley parameter, indicating the formation of Ta clusters. Strong positive WC parameter (α>0) show up for instance for Hf-Hf pairs; this demonstrates that Hf prefers not to be nearest neighbor of its own species, but rather second-nearest neighbor such as in the B2 structure. Other repelling atom pairs include Zr-Zr, Nb-Ti, Hf-Zr, and Ta-Ti.

In addition to the EAM potential65 used in the present study, there is an alternative modified embedded atom model (MEAM) potential for HfNbTaTiZr, developed by Huang et al.15,16. In the SM, Fig. S16, the WC parameters obtained for these two potentials are compared to DFT results65. The EAM potential65 employed here shows good agreement with DFT at “low” temperatures, identifying the same two pairs with the lowest WC, HfTi and TiZr, while MEAM provides positive WC values for those pairs66. We note that experiments at high temperatures, above 800 K, provide a complex picture, including TaNb-rich precipitates and HfZr-rich precipitates with a slightly different bcc phase47,67–69, and a HfZr-rich hexagonal phase47,68,69. Further studies, beyond the scope of this work, would be required to asses the applicability at high temperatures of the EAM potential used here..

Atomic displacements due to lattice distortion

In general, large atomic displacements are related to anharmonicity and lowering of the thermal conductivity70. In addition to chemical disorder, we note that lattice distortion has been assigned a large role in decreasing the thermal conductivity in bcc HEAs33, and it can be included in simulation supercells within the Coherent Potential Approximation (CPA) method, before calculating transport properties35. Within the Debye approximation, thermal displacements are 0.007 nm at 300 K. For HfNbTaTiZr, an ab-initio lattice distortion leading to atomic displacements with a mean value of 0.018 nm, much larger than expected thermal values, was reported in Ref.35.

In the SM, we display the pair correlation function, g(r), for a perfect bcc crystal and both the random and SRO alloys, with SRO after 4.7e6 steps, after relaxation, at 300 K as Fig. S1. The width at half-maximum of the nearest-neighbor peak is 0.025 nm for random and 0.022 nm for SRO, within the range of values for other bcc HEAs33. As expected, SRO reduces the width of the g(r) peaks. Partial g(r) are also shown, indicating a clear narrowing of the peaks for SRO, and the formation of HfTi pairs.

We also obtain the histogram of atomic displacements, as detailed in the SM, which extends to ∼0.05 nm, as in other bcc HEAs35. This is shown in Fig. S2a and gives mean displacement values of 0.0173 nm and 0.0153 nm, for random and SRO samples. The value for the random sample is slightly lower than the ab-initio value. SRO reduces the lattice distortion thanks to the formation of the B2 ordered structures and, therefore, this would help increasing the thermal conductivity. For NbTaTiVZr, atomic displacements were discriminated by chemical species, with Ta having the lowest value and Ti the largest one. Similar results are found here, as shown in Figs. S2b–c. Ti would be a “rattler”, i.e. the species with the largest displacements, reducing heat conductivity70.

A recent experimental study on Si nanoribbons found that strain gradients increase phonon scattering and significantly decrease thermal conductivity71. They had strains ∼0.05 and strain gradients ∼0.1/nm. In our HEA there are huge strain fluctuations along the sample due to size mismatch, which are akin to strain gradients. The strain and strain gradient can be roughly estimated as the root-mean-square displacement over the nearest-neighbor distance, and the nearest-neighbor distance squared, respectively, to give strains ∼0.06 and strain gradients ∼0.2/nm, consistent with the low thermal conductivity obtained in the simulations.

Cluster growth and percolation

Considering that the heat flow might prefer to travel across SRO structures, it is important to analyse the evolution of the SRO clusters in detail. The more relevant results are presented here while a more detailed analysis of the cluster structure, including coordination histograms, surface-to-volume evolution and fractal dimension is presented in the SM.

Fig. 2 presents the volume fractions c of the 4 relevant phases in the samples as a function of the MC steps. These have been determined using the ConstructSurfaceMesh tool in OVITO software54,72 with a probe sphere diameter of 4.0 Å. Only the volume fraction of the largest cluster was considered because, from percolation theory, the largest cluster will dominate transport across the sample. Note that these volume fractions are at most of the order of a few percent, for the last frame with the largest SRO. Their values can be found in Table S1. All these clusters display a jump in the associated volume fraction at the same time, slightly before 1e5 MC steps, which can be associated with a percolation threshold.

Precipitates in alloys can be complex, with coherent or incoherent interfaces and complex geometries, including lamellar or fibrous structures73. We illustrate the appearance and evolution of short-range-ordered clusters for the example of HfTi in Fig. 3, which shows the largest HfTi B2 cluster at a given MC step; the corresponding snapshots for the other atom pairs discussed above are contained in the SM, Fig. S3. Similar B2 clusters, with highly irregular shape, were recently reported in simulations of the NbMoTa alloy, due to diffusive evolution74.

Cluster topology appears quite complex. Given the “tortuosity” of the SRO cluster shown in Fig. 3, some atoms with low coordination Z could be eliminated from the estimated cluster volume fraction, and this is also shown in Fig. 4. It can be seen that the volume fraction cHfTi estimated with ConstructSurfaceMesh is only slightly lower than the volume fraction determined by counting the number of atoms in the cluster and eliminating all atoms with Z≤6. This is equivalent to taking bulk atoms on atoms from flat surfaces, neglecting lower-coordinated surface atoms at kinks, steps, etc. The coordination histograms at different MC steps are shown in Fig. S4, with the mean coordination shown in Fig. S5.

Fig. S6a is similar to Fig. 4, but includes not only the largest cluster but all clusters containing more than 8 atoms; these are equivalent to one atom and all its bcc neighbors in a bulk sample. Fig. S6b quantifies the evolution of the distribution of cluster sizes with the number of MC steps for the example of HfTi B2 clusters. The power law size distribution is steep, with a Fisher exponent τ=2.5 for low MC steps, then decreases and reaches the value expected for the bcc lattice site percolation, τ=2.1938, after the percolation threshold, as it will be discussed later. Figure S7 shows clusters ordered by size, with a particularly strong increase in the largest size observed between 10 000 and 50 000 MC steps. This is also related to the percolation transition. Note, also, that the number of clusters shrinks from initially > 2000 to less than 200 in the sample with the largest SRO.

Fig. 5 plots the distribution of the diameter D of Hf-Ti filaments present in the largest HfTi cluster shown in Fig. 3, after several MC steps, as obtained by using the software FoamExplorer55. The figure shows that SRO features are of the order of 1 nm, in agreement with experimental studies for HEA22,75. Aggregation of these filaments forms the entire SRO cluster, which eventually traverses the sample. The number of filaments with that size displays a steep increase just before 50 000 MC steps, consistently with the result above.

Fig. 6 shows the size of the largest HfTi cluster (which is also shown in the upper curve of Fig. 4) versus MC steps. This size can be interpreted as the order parameter of the percolation transition38. There is a clear jump and saturation above 5e4 MC steps, which can be associated with the percolation phase transition. Fig. 6 also shows that the number of filaments with diameters smaller than 1.5 nm increases with the number of MC steps, following the same trend as the cluster size. As SRO increases, there is an increase in the number of small, interconnected filaments, which is also connected to the large jump in atoms with bulk coordination, as shown in Fig. S3 in the SM.

Figs. S8–S11 display a deeper analysis of the cluster topology, including surface to volume evolution and fractal dimension estimates. Given that the fractal dimension is close to 2 for the largest SRO cluster, it can be concluded that these are not typical precipitates as found in alloys like steel.

Vibrational density of states analysis

We now discuss the VDOS for the different alloys appearing due to SRO and the clusters described previously. The VDOS of the random alloy and its change with the number of MC steps is displayed in Fig. 7a. At low frequencies f (large wavelengths λ), the alloy features the universal ∝f2 behavior; it is virtually unchanged by the effect of the SRO. For wavelengths above ∼2.5 nm, the random alloy features a rather structureless plateau. Similar VDOS shapes, without sharp peaks, have been observed for amorphous C63. At smaller wavelengths λ≲1 nm, this plateau becomes affected by the SRO in that modes at around 6 THz are depleted and modes at 7 THz are populated. Beyond 8 THz, no more phonon modes exist.

A similar behavior was found when comparing the experimental VDOS for Ni versus fcc FeCoCrMnNi76. Our results are also consistent with ab-initio results for a 54–128-atom bcc alloy where SRO was modeled by special quasi random structures77. This study showed narrow peaks for the binary alloy MoTa and much broader structures for the quinary alloy MoNbTaVW. They found that the broadening in the VDOS is related to phonon broadening, caused by the reduced phonon lifetime which appears due to mass disorder and force-constant disorder in the HEA.

Fig. 7b compares the VDOS of the random alloy and of the short-range-ordered alloy after 1 million MC steps with several ordered binary alloys for comparison. HfTi, HfNb and TiZr are chosen since these compounds are found as clusters in the short-range-ordered alloys, in agreement with the WC parameters in Fig. 1. On the contrary, NbTa is an alloy which is not favored in the short-range ordering process for this interatomic potential. The B2 phase was considered for all of these binary alloys.

Clearly, the long-range ordered binaries have a strongly different VDOS than the quinary HEA. HfTi and also TiZr feature a prominent gap between the acoustic and optical phonons. In the optical part of the spectrum, the van-Hove singularities show up as strong peaks at the edges of the optical band. The comparison with the random HEA shows that the gap between acoustic and optical phonons is completely wiped out. Due to the strong lattice distortions, the structures within the optical band are destroyed and only a featureless VDOS is left over. With the onset of SRO, the highest peak of the optical band – present both in HfTi and TiZr at around 7 THz – starts to show up. The VDOS depletion in the SRO HEA around 6 THz is reminiscent of the gap between optical and acoustical phonons and is roughly situated at the corresponding gap of both HfTi and TiZr.

Finally, we display in Fig. 7c the partial (element-specific) VDOS, i.e., the contribution of each of the constituting elements to the total VDOS. It is obtained from Eq. (3) by restricting the velocity autocorrelation to the atoms of the element investigated. The sum of all partial VDOS equals the total VDOS. We observe that the high-frequency part of the VDOS (6–8 THz), which is influenced by short-range order, is populated mostly by Ti atoms, and to a lesser degree by Zr and Hf; Nb and Ta do not contribute to this frequency region. This is in agreement with our finding that the short-range order manifests mostly in the creation of HfTi and TiZr clusters, and that the features of these clusters are mirrored in the VDOS. Recently, it was found that also in the fcc HEA AlCrCoNiFe, SRO develops sharper peaks in the VDOS as compared to the random alloy21, in agreement with our finding for a bcc HEA.

In general, nanoscale features, like the ones for SRO, will affect thermal conductivity due to high-frequency phonons interacting with those features25. Figure 5 indicates that phonons with wavelength < 1.25 nm will have no problem crossing the HfTi B2 substructure after percolation occurs. This includes most relevant phonons according to the VDOS (see top axis in Fig. 7b and c).

Lattice thermal conductivity

The influence of the chemical SRO on the lattice thermal conductivity κ (at a temperature of 300 K) is displayed in Fig. 8. An increase of κ with SRO is observed and amounts to 12% for the case of maximum SRO, after 4.7 million MC steps. It thus appears that the effect of chemical SRO on heat conductivity should be measurable in dedicated experiments of heat conductivity. The increase in heat conductivity appears to be small as long as only small ordered clusters are present in the sample (≤104 MC steps), but increases more strongly when large clusters – and in particular percolating clusters ≥105 MC steps) – have been created.

There are a multitude of models for determining the effective heat conductivity in inhomogeneous media consisting of two components78–80, and comparison between our simulated values and models for a binary system, composed only of the random alloy matrix and a HfTi B2 phase, appear in Fig. S12, showing reasonable agreement with the Maxwell model ME-278.

A more realistic model for our simulated system would include many distinct components. This is readily achieved in the ‘parallel’ model of heat conduction in which the heat flux passes through either of the media i,6 κp=∑iciκi,

and the ‘serial’ model in which the heat flux passes sequentially through all media i,7 1κs=∑iciκi.

Here, ci and κi denote the volume fraction and the heat conductivity of medium i, respectively. Wang et al.81 discuss the use of the harmonic mean of these extremes,8 1κh=βκs+1-βκp,

where the parameter β weighs the contributions of κs and κp to the total conductivity. We denote this model as the ‘mixed model’ in the following.

Table 1 lists the thermal conductivities of the various individual and binary components encountered in the short-range-ordered alloy. We note that for hcp Ti, we can compare our MD value of κ=5.10 W/(Km) with the value of 5.32 W/(Km) determined by Tong et al.6 using first-principles calculations; the agreement is satisfactory.

In order to estimate very roughly alloy effects for the random alloy, one can follow this simple approach of calculating an effective thermal conductivity using a serial or parallel model and the conductivities of all the individual elements. We obtain 5.85 or 6.24 W/(Km) for the random alloy. This is a decrease of 35% or 31% compared to pure Ta, which has the largest thermal conductivity at around 9 W/(Km). However, the decrease due to alloying effects (chemical disorder) in the random alloy has a simulated value around 0.765 W/(Km). This decrease of 91% is much larger than the estimate above, and also much larger than the maximum decrease of 70% in simulations of 4-element random alloys10. This can be compared to experimental results for sintered HfNbTaTiZr with 1 wt.% O content, giving 10-17 W/(Km) between 300 and 1200 K, while the parallel average gives about 35-40 W/(Km)42.

As SRO increases, thermal conductivity has to include different ordered B2 phases in a random matrix, and the simulated value for maximum SRO is an increase of 12% giving 0.86 W/(Km). For the model, besides the random HEA as matrix, (at least) 4 further phases are considered: HfNb, HfTi, TiZr, and Ta clusters.

Fig. 9 shows the results of the parallel and serial modes of heat conduction, Eqs. (6) and (7), respectively, calculated with these values for the B2 phases. They bracket the simulated heat conductivities as the short-range-ordered HEA evolves with MC steps. The mixed model, Eq. (8), allows to describe the data accurately, using a coefficient of β=0.65. The data are plotted vs. the volume fraction of the random phase, cr, which is determined from the volume fractions, ci of the 4 ordered phases considered (Fig. 2) as cr=1-∑ici. Good agreement of the mixed model, Eq. (8), with the simulation data is observed. The values for these plots, are available in the Supplementary Material (Table S1).

Recently, a somewhat similar model was proposed for homogeneous alloys with arbitrary stoichiometry, where the conductivity is evaluated as the sum over a mixture of solid solutions with integer composition, including the configurational entropy as a dividing prefactor82.

Thermal conductivity from elastic properties

There are simple models of the phonon thermal conductivity, using the elastic constants to calculate the Cahill lower limit to the conductivity83, the Slack conductivity84, and the mixed-conductivity26. The elastic constants can be evaluated with DFT for a perfect crystal, but this presents limitations for a SRO sample33. The interatomic potential for HfNbTaTiZr used in this work provides excellent agreement with ab-initio elastic constants46, and values for both random and SRO samples were provided in Ref.51. All relevant parameters can be calculated from those, the lattice parameter and atomic masses.

Values for the random and SRO alloys are included in Table S2, showing almost no difference between the samples. The Grüneisen parameter is fairly high, 2.2, due to the large difference between elastic and shear moduli, as reported for many bcc HEAs33. Regarding the thermal conductivity, our MD results are slightly above the Cahill limit, indicating that HEAs are relatively effective as insulators. The mixed model and the Slack model values are about 7–8 times larger than the MD values. Other refractory HEAs give similar lower limits, but significantly smaller values for the other models33. This could be traced to the large average mass and Debye temperature for HfNbTaTiZr seen in Table S2. The VDOS for our material is far from the one for a Debye model and quantitative agreement with MD is lacking.

These analytical models do not take into account SRO except for possible elasticity changes, failing to predict an increase in conductivity due to SRO. Recent studies33,34 used the Slack model to estimate the phonon thermal conductivity for many HEAs, but the large discrepancy found here with MD results points to possible pitfalls in this approach.

Temperature dependence of the thermal conductivity at high temperatures

For a model alloy with simple cubic structure, it was found that spatial correlations in the atomic mass distribution might suppress the decrease of phonon thermal conductivity with temperature for a random alloy36. Experiments and MD simulations of random fcc HEAs show a 50% decrease of thermal conductivity with temperature going from 300 K to 1100 K32.

For the case of the random alloy used here we have calculated κ at 300 K, 1000 K, and 1300 K. The decrease at high temperatures is less than 7% compared to room temperature, unlike the calculations for the fcc HEA32. The sample with SRO from 4.7e6 MC steps behaves similarly, with a negligible thermal conductivity decrease with temperature. Therefore, there is a nearly constant conductivity, even for the very low level of correlated disorder in the random alloy.

The experimental thermal conductivities of individual components in HfNbTaTiZr change only moderately at high temperatures, sometimes with a slight increment42,85. This is different from fcc metals and alloys in the same range, typically showing some thermal conductivity decrease. Experiments for a sintered HfNbTaTiZr sample show about the same thermal conductivity for 300 K and 1000 K42, around 10 W/(Km).

We note that, for the time scale of the simulations, near 1 ns, the sample did not show any phase change at such a high temperature. However, there are experiments for this alloy at high temperature showing a fraction of hcp phase86, with calculated transition temperature values of 800 K42. Because of this, our simulation results at high temperatures might be considered as results only for a model single-phase bcc HEA.

Regarding collision cascades, they typically reach temperatures well above melting within a nano-sized volume4. The cooling time of this volume depends on the thermal conductivity at high temperatures. Therefore, it would be of interest to design HEAs not only with low conductivity, but also with a decrease in conductivity with temperature. Ab-initio simulations for several bcc HEAs do show a lowering of the electronic thermal conductivity with temperature34. This decrease would augment defect recombination and help radiation resistance. From this work, SRO causes an increase in the conductivity, together with possible contribution to avoid a decrease at high temperatures. From the above, this would be detrimental to radiation resistance. However, there is a competing effect, given that SRO phases might help trapping defects and facilitating their recombination. The final net outcome would have to be studied in detail for the relevant HEA.

Discussion

Caro et al.10 simulated fcc random metallic alloys with Lennard-Jones pair potentials. They found that the thermal conductivity in random binary alloys decreases by 30–70% with respect to the conductivity of individual components, with a minimum roughly around equiatomic composition. For random 4-component alloys the decrease depends on stoichiometry, but it can be as large as 70%.

Lennard-Jones potentials were used in early studies of HEAs10 because the construction of reliable interatomic potentials for HEAs is challenging. Even if they work well for a random alloy, they might not offer correct results for SRO evolution. As an example, the AlFeNiCoCr potential by Farkas and Caro87 has been used for some SRO studies88,89, including conductivity simulations21, despite the fact that the paper describing the potential warns about the error in describing preferred formation of FeAl instead of FeNi pairs, as found in experiments and DFT simulations. For the material described in this study, there are experiments in the temperature range 750–1300 K, which find bcc precipitates rich in NbTa67–69,90,91, bcc precipitates rich in HfZr91, and hcp precipitates rich in HfZr68,69,91. Regarding atomistic simulations of the same alloy, using a MEAM potential, Huang et al.15 observe negative WC parameters for NbTa, TiZr, HfZr and HfTi. After MC/MD at 800 K, they observed HfTiZr-rich clusters, which gave rise to hcp precipitates; at 300 K they did not observe any hcp precipitates. Wu et al.15,92 used a charge-transfer ionic potential (CTIP) optimized with machine learning to anneal a sample above 2000 K and then cool down to room temperature. This leads to significant chemical segregation, and the most negative WC parameter is HfTi, as in our low-temperature results, without annealing. They also observe negative WC for NbTa and ZrTi pairs. Given the lack of experimental and ab-initio simulation results for SRO in this alloy at room temperature, additional studies would be required to ascertain the appropriate description of SRO evolution.

A recent analysis of conductivity in SRO samples presented the conductivity increase in terms of the Warren-Cowley parameter of the pair with lowest WC value, for an fcc HEA21. They used a cubic box with 21×21×21 fcc cells (37044 atoms), which would limit the SRO cluster sizes. A similar analysis for our sample, focusing only on the pair with lowest WC parameter, which in our case is HfTi, is shown in the SM, also showing the correlation between conductivity and WC.

We note that for fcc HEAs, there are several favored pairs, as it happens for our bcc HEA (seeFig. 2) and for most HEAs. Therefore, we consider that our analysis in terms of the fraction of various B2 phases, as shown in Fig. 9, provides an improved description. We show αHfTi versus MC step in the SM (Fig. S13), which displays a rapid increase until about 1e5 MC steps, in agreement with the observed percolation threshold, but the percolation transition is not as clear for the WC parameter as for the largest cluster size.

For the most preferred pair, HfTi, the WC parameter αHfTi increases linearly with the fraction of B2 HfTi phase, cHfTi, after an initial stage with negligible B2 phase, as shown in Fig. S14 of the SM. The fraction of B2 phase might be measured experimentally from diffraction studies, while the WC parameters appear as more challenging to obtain from experiments. The lattice heat conductivity κ has a fairly linear dependence with αHfTi, as shown in Fig. S15.

Returning to bcc HEAs, there is a recent study on 12 different alloys, using ab-initio simulations to obtain their elastic constants, and models to obtain their thermal properties from those33. This results in thermal conductivities in the range of 0.3–1.8 W/(Km), with the exception of two quaternary HEAs with higher conductivity. For instance, HfNbTiZrV has κ=0.6 W/(Km). Using the Cahill model83 for the minimum thermal conductivity gives values of around 0.5 W/(Km) for all those alloys. It can be seen that the conductivity values found here are well within the expected conductivity range.

Machine learning was used to predict thermal conductivity versus temperature in alloys, based on chemical composition and a large training dataset with experimental conductivity values93. Such an approach would have to be enhanced for SRO alloys, given that chemical composition alone is not sufficient to determine κ.

Percolation theory and models applied to thermal conductivity

The fact that the growth of HfTi B2 clusters with the number of MC steps is related to an increase in thermal conductivity is explained here by the presence of a percolative transition. This transition is not driven by a non-equilibrium deposition process, as is usual in percolation of RSA-type processes, but by a Kawasaki-type exchange process which is responsible for the clustering of some atomic species.

Recently, the magnetization of finite Ising systems was evaluated within the numerical mean spin method assuming that there is a fraction of non-magnetic atoms occupying a regular lattice, with magnetic atoms interacting with first neighbors94. The critical temperature for the ferromagnetic/paramagnetic transition was associated to a percolation transition of the sites occupied by magnetic atoms, i.e., above the percolation transition the average magnetization drops to zero. In our study, the increase in conductivity (magnetization) could also be related to the percolation transition for conducting/non-conducting sites.

Summary and conclusions

By performing MD simulations on bcc single-crystalline equiatomic HfNbTaTiZr alloys with various degrees of chemical SRO, we obtained the following findings on the influence of chemical SRO on the lattice thermal conductivity. Increasing chemical SRO in HfNbTaTiZr consists in the build-up of nm-sized HfTi and TiZr clusters, in agreement with the evolution of Warren-Cowley parameters.

Evolution of the size of the largest cluster indicates a percolation transition, which is also reflected in other metrics, including number of filaments composing the cluster, atomic coordination and fractal dimension. The cluster size distribution follows a power law, with a Fisher exponent τ=2.19, as expected from site percolation in a bcc lattice38.

The cluster structure is complex, with fractal dimension around 2.4 at the transition and less than 2 at the end of our simulations, indicating a highly irregular cluster structure.

The percolative transition is observed in samples generated with increasing SRO through a MC Kawasaki exchange process. Therefore, this is not a strictly “classical” random percolation but a “thermal” percolation, as the one studied by Giménez et al.56. The agreement in the values of both Fisher exponent and fractal dimension at the transition, within the numerical errors, suggest, as in56, that the problem still belongs to the classical random percolation universality class.

Lattice distortion has been associated to a reduction of the thermal conductivity33,35. We observe that SRO leads to a 12.5% reduction of the atomic displacements due to the formation of ordered phases.

Thermal conductivity is related to the vibrational density of states (VDOS). In contrast to long-range chemically ordered crystals (such as B2-phase binary alloys), showing pronounced localized peaks, the VDOS of the random HfNbTaTiZr alloy is characterized by an almost structureless plateau. Acoustic and optic bands are merged, without the appearance of a gap. Van-Hove singularities, which feature band edges, are absent.

With increasing SRO, the VDOS of HfNbTaTiZr shows structure in the high-frequency part of the spectrum. A detailed analysis correlates this structure with the optical bands in HfTi and TiZr.

The MD simulated phonon thermal conductivity, κMD, can be compared to models which use the elastic constants for the samples. We find that κMD is slightly above the Cahill lower limit83, but significantly lower than the estimates provided by the mixed model26 and the Slack model84.

Calculation of the phonon thermal conductivity from MD simulations shows a continuous increase with chemical SRO. The increase amounts to roughly 12 % above that of the random alloy and might be experimentally measurable. This increase is not related to changes in elastic constants.

The increase can be quantitatively explained by a simple composite model, where the material is assumed to be a combination of several phases, each with different heat conductivity. The phases correspond to the various SRO binary alloys, which have higher thermal conductivity than the random alloy matrix, leading to the net conductivity increase.

Conductivity does not decrease significantly with temperature for the random or SRO samples, unlike the strong decrease observed for a random fcc HEA32.

We note that some materials, including semiconductors like Si, include important contributions to the heat conduction by long-wavelength phonons of tens of nm95,96. There are relevant phonon mean-free paths of up to 1 μm, as demonstrated by experiments and simulations97. This leads to strong reduction of the thermal conductivity with porosity25, but might also lead to thermal conductivity reduction with SRO in semiconductor alloys, due to the impossibility of phonon transmission along SRO structures with a width smaller than their wavelength.

Phonon thermal conductivity will be reduced by alloy effects in most materials24, and precipitates within the alloy will further reduce the thermal conductivity25. This is because the boundary between matrix and precipitate can contribute to phonon scattering, but also because in many alloys of technological interest, precipitates have lower heat conduction that the matrix, as is the case for Si precipitates in Al alloys73. Despite this quite general behavior, we observe an increase in thermal conductivity in HEA with SRO, because ordered SRO precipitates have higher thermal conductivity than the random HEA matrix. We expect this behavior to be typical for any HEA with SRO, including fcc HEAs. It might be also found in other multicomponent materials with similar behavior of the thermal conductivity. In Mg alloys, addition of alloying elements usually decreases heat conductivity, but aging the alloy can lead to large thermal conductivity increments partly due to the growth of second-phase precipitates98. A recent study in W finds that phonon conductivity dominates over electronic conductivity for pressures above 20 GPa99. Some W alloys might display similar behavior and the findings in this paper might be relevant to those cases.

Besides phonons, electrons also contribute to the thermal conductivity of HEAs. Using ab-initio theory for the example of the Cantor alloy100,101, it has been found that – while the electron density is comparable to that of elemental metals – the electron mean-free path decreases strongly in HEAs, due to the strong lattice and chemical disorder. In the Wiedemann-Franz law the Lorenz number relates the ratio of electrical conductivity and the electronic contribution of heat conductivity to temperature102. Recently, the Wiedemann-Franz law was applied to experimental measurements of the electrical resistivity in HfNbTaTiZr, in order to separate phonon and electron contributions to the thermal conductivity at low temperatures35. Although there are many limitations to this approximate law, given a typical short mean free path for electrons in solids, we expect that electrical conductivity might also increase with SRO for HEA.

Since the alloys and phases arising with SRO have a large influence in the thermodynamic properties of HEAs, atomistic simulations, including both Monte Carlo and molecular dynamics simulations will require in the future improved potential development, including machine learning potentials, which needs to include the thermodynamics of several possible binary phases and their interfaces, in order to model SRO evolution and its effect on different properties of technological relevance, alongside machine learning approaches for improved prediction and design of thermal conductivity in complex materials37,103.Table 1 Thermal conductivities κ of the relevant species and B2 alloys found in the short-range-ordered HEA, as determined by MD in the present study.

Crystal	Structure	κ [W/(Km)]	
HEA(r)	BCC	0.765	
HEA(SRO)	BCC	0.86	
HfNb	B2	3.27	
HfTi	B2	2.66	
TiZr	B2	4.31	
Hf	BCC	4.71	
Nb	BCC	6.42	
Ta	BCC	9.59	
Zr	BCC	5.39	
Ti	BCC	8.24	
Ti	HCP	5.10	
HEA(r) indicates the random alloy, while HEA(SRO) indicates the alloy after 4.7e6 MC steps

Figure 1 Evolution of the Warren–Cowley parameters, Eq. (1), with the number of MC steps. Only the atomic pairs showing the largest deviations from 0 are shown. After 1e3 steps, the WC parameters are nearly zero, indicating that SRO is almost negligible.

Figure 2 Volume fractions cn of the different B2 phases found in the short-range-ordered HEA after a certain number of MC steps, obtained with SurfaceMesh from OVITO54, for the largest clusters.

Figure 3 Series of snapshots illustrating the presence of B2-structured HfTi clusters in samples with the SRO established after an increasing number of MC steps: (a) 1 000, (b) 10 000, (c) 50 000, (d) 100 000, (e) 1 000 000. Only the largest cluster is shown in each snapshot; due to the periodic boundary conditions it may appear disconnected, but it is not.

Figure 4 Volume fraction of the largest B2 HfTi cluster, calculated as Na/Nt, filtering the number of atoms (Na) according to their coordination Z. Z=8 corresponds to bulk B2 atoms. Only Hf and Ti atoms with at least one HfTi bond are included. cHfTi, indicating the fraction obtained from SurfaceMesh in OVITO54, is also shown.

Figure 5 Evolution of the distribution of filament sizes, i.e., diameters D, for the largest HfTi cluster, for different MC steps, obtained from FoamExplorer55.

Figure 6 Evolution of the fraction of atoms contained in the largest HfTi cluster, together with the number of HfTi filaments with diameters smaller than 1.5 nm (arbitrary units), with the number of MC steps. Here the percolation threshold is associated with a critical number of MC steps of around 105, from which a spanning cluster emerges, and it is marked by a vertical dashed line. At this threshold, just below 105 steps, the number of filaments shows a sudden increase.

Figure 7 (a) Evolution of the VDOS for various MC steps indicated in the legend. (b) VDOS of the short-range-ordered alloy (after 1 million MC steps), compared to that of several binary B2 alloys. (c) Contribution of the constituting elements to the VDOS of the fully SRO alloy (after 1 million MC steps). In all panels, the VDOS is normalized to area 1 and shown versus frequency f and wavelength λ.

Figure 8 Evolution of the thermal conductivity with the number of MC steps in the SRO alloy. Lines before (dotted) and after (dashed) the percolation threshold, have been added to guide the eye.

Figure 9 Thermal conductivity of the short-range-ordered HEA. The number of MC steps is coded in the volume fraction of the random phase, cr. Simulation data are compared to the parallel, serial and mixed models of heat conduction, Eqs. (6), (7) and (8), respectively.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-70500-9.

Acknowledgements

EMB would like to thank Francisco Muñoz for pointing out related research which motivated this study. This work used the TOKO Cluster from FCEN-UNCuyo and the Serafín Cluster from CCAD-UNC, which are part of the SNCAD-MinCyT, Argentina.

Author contributions

All authors read and approved the final manuscript. G.M.B. performed and analyzed the calculations. O.R.D. performed and analyzed the calculations to obtain the virtual samples. M.P. analyzed percolation. H.M.U. wrote the initial manuscript. E.M.B. conceived the study and analyzed some of the results. All authors discussed the results and contributed to the final manuscript.

Funding

Open Access funding enabled and organized by Projekt DEAL. EMB, GMB and ORD thank support from SIIP-UNCuyo 06/M008-T1 grant, the Agencia Nacional de Promoción Científica y Tecnológica (ANPCyT) PICTO-UUMM-2019-00048 grant, and CONICET PIP 2021-2023 11220200102578CO grant.

Data availability

The data that supports the findings of this study are available from the corresponding author upon reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Zhang Y Zhou YJ Lin JP Chen GL Peter K Liaw Solid-solution phase formation rules for multi-component alloys Adv. Eng. Mater. 2008 10 534 538 10.1002/adem.200700240
Zhang, Y., Zhou, Y. J., Lin, J. P., Chen, G. L. & Peter, K. Liaw Solid-solution phase formation rules for multi-component alloys. Adv. Eng. Mater. 10, 534–538. 10.1002/adem.200700240 (2008).10.1002/adem.200700240
2. Li W Xie D Li D Zhang Y Gao Y Peter KL Mechanical behavior of high-entropy alloys Prog. Mater Sci. 2021 118 100777 10.1016/j.pmatsci.2021.100777
Li, W. et al. Mechanical behavior of high-entropy alloys. Prog. Mater Sci. 118, 100777 (2021).10.1016/j.pmatsci.2021.100777
3. Liu F Liaw PK Zhang Y Recent progress with bcc-structured high-entropy alloys Metals 2022 12 501 10.3390/MET12030501
Liu, F., Liaw, P. K. & Zhang, Y. Recent progress with bcc-structured high-entropy alloys. Metals 12, 501. 10.3390/MET12030501 (2022).10.3390/MET12030501
4. Deluigi OR Pasianot RC Valencia FJ Caro A Farkas D Bringa EM Simulations of primary damage in a high entropy alloy: Probing enhanced radiation resistance Acta Mater. 2021 213 116951 10.1016/j.actamat.2021.116951
Deluigi, O. R. et al. Simulations of primary damage in a high entropy alloy: Probing enhanced radiation resistance. Acta Mater. 213, 116951. 10.1016/j.actamat.2021.116951 (2021).10.1016/j.actamat.2021.116951
5. Gludovatz B Hohenwarter A Catoor D Chang EH George EP Ritchie RO A fracture-resistant high-entropy alloy for cryogenic applications Science 2014 345 1153 1158 10.1126/science.1254581 25190791
Gludovatz, B. et al. A fracture-resistant high-entropy alloy for cryogenic applications. Science 345, 1153–1158. 10.1126/science.1254581 (2014).25190791 10.1126/science.1254581
6. Tong Z Li S Ruan X Bao H Comprehensive first-principles analysis of phonon thermal conductivity and electron-phonon coupling in different metals Phys. Rev. B 2019 100 144306 10.1103/PhysRevB.100.144306
Tong, Z., Li, S., Ruan, X. & Bao, H. Comprehensive first-principles analysis of phonon thermal conductivity and electron-phonon coupling in different metals. Phys. Rev. B 100, 144306. 10.1103/PhysRevB.100.144306 (2019).10.1103/PhysRevB.100.144306
7. Jin K Mu S An K Porter WD Samolyuk GD Stocks GM Bei H Thermophysical properties of Ni-containing single-phase concentrated solid solution alloys Mater. Des. 2017 117 185 192 10.1016/j.matdes.2016.12.079
Jin, K. et al. Thermophysical properties of Ni-containing single-phase concentrated solid solution alloys. Mater. Des. 117, 185–192. 10.1016/j.matdes.2016.12.079 (2017).10.1016/j.matdes.2016.12.079
8. Tsai M-H Physical properties of high entropy alloys Entropy 2013 15 5338 5345 10.3390/e15125338
Tsai, M.-H. Physical properties of high entropy alloys. Entropy 15, 5338–5345. 10.3390/e15125338 (2013).10.3390/e15125338
9. Zhao Y Liu D Chen J Zhu L Belianinov A Ovchinnikova OS Unocic RR Burch MJ Kim S Hao H Pickard DS Li B Thong JTL Engineering the thermal conductivity along an individual silicon nanowire by selective helium ion irradiation Nat. Commun. 2017 8 15919 10.1038/ncomms15919 28653663
Zhao, Y. et al. Engineering the thermal conductivity along an individual silicon nanowire by selective helium ion irradiation. Nat. Commun. 8, 15919. 10.1038/ncomms15919 (2017).28653663 10.1038/ncomms15919
10. Caro M Béland LK Samolyuk GD Stoller RE Caro A Lattice thermal conductivity of multi-component alloys J. Alloy. Compd. 2015 648 408 413 10.1016/j.jallcom.2015.06.035
Caro, M., Béland, L. K., Samolyuk, G. D., Stoller, R. E. & Caro, A. Lattice thermal conductivity of multi-component alloys. J. Alloy. Compd. 648, 408–413. 10.1016/j.jallcom.2015.06.035 (2015).10.1016/j.jallcom.2015.06.035
11. Fan Z Wang H Wu Y Liu XJ Lu ZP Thermoelectric high-entropy alloys with low lattice thermal conductivity RSC Adv. 2016 6 52164 52170 10.1039/C5RA28088E
Fan, Z., Wang, H., Wu, Y., Liu, X. J. & Lu, Z. P. Thermoelectric high-entropy alloys with low lattice thermal conductivity. RSC Adv. 6, 52164–52170. 10.1039/C5RA28088E (2016).10.1039/C5RA28088E
12. Nordlund K Wei L Zhong Y Averback RS Role of electron-phonon coupling in collision cascade development in Ni, Pd, and Pt Phys. Rev. B 1998 57 R13965 10.1103/PhysRevB.57.R13965
Nordlund, K., Wei, L., Zhong, Y. & Averback, R. S. Role of electron-phonon coupling in collision cascade development in Ni, Pd, and Pt. Phys. Rev. B 57, R13965. 10.1103/PhysRevB.57.R13965 (1998).10.1103/PhysRevB.57.R13965
13. Zhang Y Stocks GM Jin K Chenyang L Bei H Sales BC Wang L Beland LK Stoller RE Samolyuk GD Caro M Caro A Weber WJ Influence of chemical disorder on energy dissipation and defect evolution in concentrated solid solution alloys Nat. Commun. 2015 6 8736 10.1038/ncomms9736 26507943
Zhang, Y. et al. Influence of chemical disorder on energy dissipation and defect evolution in concentrated solid solution alloys. Nat. Commun. 6, 8736 (2015).26507943 10.1038/ncomms9736
14. Chen S Aitken ZH Pattamatta S Zhaoxuan W Zhi Gen Yu Banerjee R Srolovitz DJ Liaw PK Zhang Y-W Chemical-affinity disparity and exclusivity drive atomic segregation, short-range ordering, and cluster formation in high-entropy alloys Acta Mater. 2021 206 116638 10.1016/j.actamat.2021.116638
Chen, S. et al. Chemical-affinity disparity and exclusivity drive atomic segregation, short-range ordering, and cluster formation in high-entropy alloys. Acta Mater. 206, 116638. 10.1016/j.actamat.2021.116638 (2021).10.1016/j.actamat.2021.116638
15. Huang X Liu L Duan X Liao W Huang J Sun H Chunyan Yu Atomistic simulation of chemical short-range order in HfNbTaZr high entropy alloy based on a newly-developed interatomic potential Mater. Des. 2021 202 109560 10.1016/j.matdes.2021.109560
Huang, X. et al. Atomistic simulation of chemical short-range order in HfNbTaZr high entropy alloy based on a newly-developed interatomic potential. Mater. Des. 202, 109560. 10.1016/j.matdes.2021.109560 (2021).10.1016/j.matdes.2021.109560
16. Huang X Liu L Liao W Huang J Sun H Chunyan Yu Characterization of nucleation behavior in temperature-induced BCC-to-HCP phase transformation for high entropy alloy Acta Metall. Sinica. 2021 34 1546 1556 10.1007/s40195-021-01282-6
Huang, X. et al. Characterization of nucleation behavior in temperature-induced BCC-to-HCP phase transformation for high entropy alloy. Acta Metall. Sinica. 34, 1546–1556. 10.1007/s40195-021-01282-6 (2021).10.1007/s40195-021-01282-6
17. Ferrari A Körmann F Asta M Neugebauer J Simulating short-range order in compositionally complex materials Nat. Comput. Sci. 2023 3 221 229 10.1038/s43588-023-00407-4 38177884
Ferrari, A., Körmann, F., Asta, M. & Neugebauer, J. Simulating short-range order in compositionally complex materials. Nat. Comput. Sci. 3, 221–229. 10.1038/s43588-023-00407-4 (2023).38177884 10.1038/s43588-023-00407-4
18. Lokman Ali Md Shinzato S Wang V Shen Z Du J Ogata S An atomistic modeling study of the relationship between critical resolved shear stress and atomic structure distortion in fcc high entropy alloys - relationship in random solid solution and chemical-short-range-order alloys - Mater. Trans. 2020 61 605 609 10.2320/matertrans.MT-MK2019007
Lokman Ali, Md. et al. An atomistic modeling study of the relationship between critical resolved shear stress and atomic structure distortion in fcc high entropy alloys - relationship in random solid solution and chemical-short-range-order alloys -. Mater. Trans. 61, 605–609. 10.2320/matertrans.MT-MK2019007 (2020).10.2320/matertrans.MT-MK2019007
19. Antillon E Woodward C Rao SI Akdim B Chemical short range order strengthening in bcc complex concentrated alloys Acta Mater. 2021 215 117012 10.1016/j.actamat.2021.117012
Antillon, E., Woodward, C., Rao, S. I. & Akdim, B. Chemical short range order strengthening in bcc complex concentrated alloys. Acta Mater. 215, 117012. 10.1016/j.actamat.2021.117012 (2021).10.1016/j.actamat.2021.117012
20. Liu X Curtin WA Atomistic simulations reveal strength reductions due to short-range order in alloys Acta Mater. 2024 263 119471 10.1016/j.actamat.2023.119471
Liu, X. & Curtin, W. A. Atomistic simulations reveal strength reductions due to short-range order in alloys. Acta Mater. 263, 119471. 10.1016/j.actamat.2023.119471 (2024).10.1016/j.actamat.2023.119471
21. Hasan, M. A. A., Shin, S. & Liaw, K. L. Short-range order effects on the thermodynamic behavior of AlCoCrFeNi high-entropy alloys. Comput. Mater. Sci. 239, 112980. 10.1016/j.commatsci.2024.112980 (2024).
22. Wu Xiaolei Chemical short-range orders in high-/medium-entropy alloys J. Mater. Sci. Technol. 2023 147 189 196 10.1016/J.JMST.2022.10.070
Wu, X. Chemical short-range orders in high-/medium-entropy alloys. J. Mater. Sci. Technol. 147, 189–196. 10.1016/J.JMST.2022.10.070 (2023).10.1016/J.JMST.2022.10.070
23. Chen S Aitken ZH Pattamatta S Zhaoxuan W Zhi Gen Yu Srolovitz DJ Liaw PK Zhang YW Short-range ordering alters the dislocation nucleation and propagation in refractory high-entropy alloys Mater. Today 2023 65 14 25 10.1016/J.MATTOD.2023.03.009
Chen, S. et al. Short-range ordering alters the dislocation nucleation and propagation in refractory high-entropy alloys. Mater. Today 65, 14–25. 10.1016/J.MATTOD.2023.03.009 (2023).10.1016/J.MATTOD.2023.03.009
24. Seyf, H R., Yates, L., Bougher, T L., Graham, S., Cola, B. A., Detchprohm, T., Ji, M-H., Kim, J., Dupuis, R., LV, W. et al., Rethinking phonons: The issue of disorder npj Comput. Mater. 3, 49 (2017) 10.1038/s41524-017-0052-9
25. Qian X Zhou J Chen G Phonon-engineered extreme thermal conductivity materials Nat. Mater. 2021 20 1188 1202 10.1038/s41563-021-00918-3 33686278
Qian, X., Zhou, J. & Chen, G. Phonon-engineered extreme thermal conductivity materials. Nat. Mater. 20, 1188–1202. 10.1038/s41563-021-00918-3 (2021).33686278 10.1038/s41563-021-00918-3
26. Toberer ES Zevalkink A Jeffrey Snyder G Phonon engineering through crystal chemistry J. Mater. Chem. 2011 21 15843 15852 10.1039/C1JM11754H
Toberer, E. S., Zevalkink, A. & Jeffrey Snyder, G. Phonon engineering through crystal chemistry. J. Mater. Chem. 21, 15843–15852. 10.1039/C1JM11754H (2011).10.1039/C1JM11754H
27. Giri A Braun JL Tomko JA Hopkins PE Reducing the thermal conductivity of chemically ordered binary alloys below the alloy limit via the alteration of phonon dispersion relations Appl. Phys. Lett. 2017 110 233112 10.1063/1.4985204
Giri, A., Braun, J. L., Tomko, J. A. & Hopkins, P. E. Reducing the thermal conductivity of chemically ordered binary alloys below the alloy limit via the alteration of phonon dispersion relations. Appl. Phys. Lett. 110, 233112. 10.1063/1.4985204 (2017).10.1063/1.4985204
28. McGaughey AJH Jain A Kim H-Y Fu B Phonon properties and thermal conductivity from first principles, lattice dynamics, and the Boltzmann transport equation J. Appl. Phys. 2019 125 011101 10.1063/1.5064602
McGaughey, A. J. H., Jain, A., Kim, H.-Y. & Fu, B. Phonon properties and thermal conductivity from first principles, lattice dynamics, and the Boltzmann transport equation. J. Appl. Phys. 125, 011101. 10.1063/1.5064602 (2019).10.1063/1.5064602
29. Cheng C Ma S Wang S The role of phonon anharmonicity on the structural stability and phonon heat transport of CrFeCoNiCux high-entropy alloys at finite temperatures J. Alloy. Compd. 2023 935 168003 10.1016/j.jallcom.2022.168003
Cheng, C., Ma, S. & Wang, S. The role of phonon anharmonicity on the structural stability and phonon heat transport of CrFeCoNiCu high-entropy alloys at finite temperatures. J. Alloy. Compd. 935, 168003. 10.1016/j.jallcom.2022.168003 (2023).10.1016/j.jallcom.2022.168003
30. Jin K Sales BC Stocks GM Samolyuk GD Daene M Weber WJ Zhang Y Bei H Tailoring the physical properties of Ni-based single-phase equiatomic alloys by modifying the chemical complexity Sci. Rep. 2016 6 20159 10.1038/srep20159 26832223
Jin, K. et al. Tailoring the physical properties of Ni-based single-phase equiatomic alloys by modifying the chemical complexity. Sci. Rep. 6, 20159. 10.1038/srep20159 (2016).26832223 10.1038/srep20159
31. Al Hasan MA Wang J Shin S Gilbert DA Liaw PK Tang N Namila WL Liyanage C Santodonato L DeBeer-Schmitt L Butch NP Butch, effects of aluminum content on thermoelectric performance of AlxCoCrFeNi high-entropy alloys J. Alloys Compd. 2021 883 160811 10.1016/j.jallcom.2021.160811
Al Hasan, M. A. et al. Effects of aluminum content on thermoelectric performance of AlCoCrFeNi high-entropy alloys. J. Alloys Compd. 883, 160811. 10.1016/j.jallcom.2021.160811 (2021).10.1016/j.jallcom.2021.160811
32. Sun Z Shi C Gao L Lin S Li W Thermal physical properties of high entropy alloy Al0.3CoCrFeNi at elevated temperatures J. Alloys Compd. 2022 901 163554 10.1016/j.jallcom.2021.163554
Sun, Z., Shi, C., Gao, L., Lin, S. & Li, W. Thermal physical properties of high entropy alloy Al0.3CoCrFeNi at elevated temperatures. J. Alloys Compd. 901, 163554 (2022).10.1016/j.jallcom.2021.163554
33. San S Hasan S Adhikari P Ching W-Y Designing quaternary and quinary refractory-based high-entropy alloys: Statistical analysis of their lattice distortion, mechanical, and thermal properties Metals 2023 13 1953 10.3390/met13121953
San, S., Hasan, S., Adhikari, P. & Ching, W.-Y. Designing quaternary and quinary refractory-based high-entropy alloys: Statistical analysis of their lattice distortion, mechanical, and thermal properties. Metals 13, 1953. 10.3390/met13121953 (2023).10.3390/met13121953
34. Singh, P. et al. Alloying effects on the transport properties of refractory high-entropy alloys. Acta Materialia 276, 120032. 10.1016/j.actamat.2024.120032 (2024).
35. Dong Y Sai M Guo X Han J Duan J Jia N Xue Y Jin K Transport properties of refractory high-entropy alloys with single-phase body-centered cubic structure Scripta Mater. 2023 231 115464 10.1016/j.scriptamat.2023.115464
Dong, Y. et al. Transport properties of refractory high-entropy alloys with single-phase body-centered cubic structure. Scripta Mater. 231, 115464. 10.1016/j.scriptamat.2023.115464 (2023).10.1016/j.scriptamat.2023.115464
36. Thébaud S Lindsay L Berlijn T Breaking Rayleigh’s law with spatially correlated disorder to control phonon transport Phys. Rev. Lett. 2023 131 026301 10.1103/PhysRevLett.131.026301 37505967
Thébaud, S., Lindsay, L. & Berlijn, T. Breaking Rayleigh’s law with spatially correlated disorder to control phonon transport. Phys. Rev. Lett. 131, 026301. 10.1103/PhysRevLett.131.026301 (2023).37505967 10.1103/PhysRevLett.131.026301
37. Zhu C Luo T Li B Shen X Zhu G Machine learning aided understanding and manipulating thermal transport in amorphous networks J. Appl. Phys. 2024 135 195103 10.1063/5.0200779
Zhu, C., Luo, T., Li, B., Shen, X. & Zhu, G. Machine learning aided understanding and manipulating thermal transport in amorphous networks. J. Appl. Phys. 135, 195103. 10.1063/5.0200779 (2024).10.1063/5.0200779
38. Stauffer D Aharony A Introduction to Percolation Theory 2018 New York Taylor & Francis
Stauffer, D. & Aharony, A. Introduction to Percolation Theory (Taylor & Francis, New York, 2018).
39. Sahimi M Applications of Percolation Theory 1994 New York CRC Press
Sahimi, M. Applications of Percolation Theory (CRC Press, New York, 1994).
40. Wang X Li X Li Z Zheng Y Dong C Liaw PK Utilizing cluster percolation theory to analyze electrical resistivity jumps in high-entropy alloy thin films containing small atoms Scripta Mater. 2024 238 115727 10.1016/j.scriptamat.2023.115727
Wang, X. et al. Utilizing cluster percolation theory to analyze electrical resistivity jumps in high-entropy alloy thin films containing small atoms. Scripta Mater. 238, 115727. 10.1016/j.scriptamat.2023.115727 (2024).10.1016/j.scriptamat.2023.115727
41. Abere MJ Ziade E Ping L Saltonstall CB Xiaojun G Wright WJ Argibay N Kustas AB A predictive analytical model of thermal conductivity for aluminum/transition metal high-entropy alloys Scripta Mater. 2022 208 114330 10.1016/j.scriptamat.2021.114330
Abere, M. J. et al. A predictive analytical model of thermal conductivity for aluminum/transition metal high-entropy alloys. Scripta Mater. 208, 114330. 10.1016/j.scriptamat.2021.114330 (2022).10.1016/j.scriptamat.2021.114330
42. Zhang D Hsu Y-C Dunand DC Ink-extrusion 3D printing and silicide coating of HfNbTaTiZr refractory high-entropy alloy for extreme temperature applications Adv. Sci. 2024 11 2309693 10.1002/advs.202309693
Zhang, D., Hsu, Y.-C. & Dunand, D. C. Ink-extrusion 3D printing and silicide coating of HfNbTaTiZr refractory high-entropy alloy for extreme temperature applications. Adv. Sci. 11, 2309693. 10.1002/advs.202309693 (2024).10.1002/advs.202309693
43. Senkov ON Scott JM Senkova SV Miracle DB Woodward CF Microstructure and room temperature properties of a high-entropy TaNbHfZrTi alloy J. Alloy. Compd. 2011 509 6043 6048 10.1016/j.jallcom.2011.02.171
Senkov, O. N., Scott, J. M., Senkova, S. V., Miracle, D. B. & Woodward, C. F. Microstructure and room temperature properties of a high-entropy TaNbHfZrTi alloy. J. Alloy. Compd. 509, 6043–6048. 10.1016/j.jallcom.2011.02.171 (2011).10.1016/j.jallcom.2011.02.171
44. Senkov ON Miracle DB Chaput KJ Couzinie J-P Development and exploration of refractory high entropy alloys—a review J. Mater. Res. 2018 33 3092 3128 10.1557/jmr.2018.153
Senkov, O. N., Miracle, D. B., Chaput, K. J. & Couzinie, J.-P. Development and exploration of refractory high entropy alloys—a review. J. Mater. Res. 33, 3092–3128. 10.1557/jmr.2018.153 (2018).10.1557/jmr.2018.153
45. Thompson, A. P., Aktulga, H.M., Berger, R., Bolintineanu, D. S., Brown, W.M., Crozier, P. S., Pieter J. in ’t Veld, Axel Kohlmeyer, Stan G. Moore, Trung Dac Nguyen, Ray Shan, Mark J. Stevens, Julien Tranchida, Christian Trott, and Steven J. Plimpton, LAMMPS—a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales, 10.1016/j.cpc.2021.108171Computer Physics Communications 271, 108171 (2022)
46. Xu S Jian W-R Beyerlein IJ Ideal simple shear strengths of two HfNbTaTi-based quinary refractory multi-principal element alloys APL Mater. 2022 10 111107 10.1063/5.0116898
Shuozhi, X., Jian, W.-R. & Beyerlein, I. J. Ideal simple shear strengths of two HfNbTaTi-based quinary refractory multi-principal element alloys. APL Mater. 10, 111107. 10.1063/5.0116898 (2022).10.1063/5.0116898
47. Malek J Zyka J Lukac F Monika V Vlasak T Cizek J Melikhova O Adela M Kim H-S The effect of processing route on properties of HfNbTaTiZr high entropy alloy Materials 2019 12 4022 10.3390/ma12234022 31817034
Malek, J. et al. The effect of processing route on properties of HfNbTaTiZr high entropy alloy. Materials 12, 4022. 10.3390/ma12234022 (2019).31817034 10.3390/ma12234022
48. Dirras G Gubicza J Heczel A Lilensten L Couzinie J-P Perriere L Guillot I Hocini A Microstructural investigation of plastically deformed Ti20Zr20Hf20Nb20Ta20 high entropy alloy by X-ray diffraction and transmission electron microscopy Mater. Charact. 2015 108 1 7 10.1016/j.matchar.2015.08.007
Dirras, G. et al. Microstructural investigation of plastically deformed TiZrHfNbTa high entropy alloy by X-ray diffraction and transmission electron microscopy. Mater. Charact. 108, 1–7. 10.1016/j.matchar.2015.08.007 (2015).10.1016/j.matchar.2015.08.007
49. Kawasaki K Diffusion constants near the critical point for time-dependent Ising models I Phys. Rev. 1966 145 224 10.1103/PhysRev.145.224
Kawasaki, K. Diffusion constants near the critical point for time-dependent Ising models I. Phys. Rev. 145, 224. 10.1103/PhysRev.145.224 (1966).10.1103/PhysRev.145.224
50. Fey LTW Beyerlein IJ Random generation of lattice structures with short-range order Integr. Mater. Manuf. Innov. 2022 11 382 390 10.1007/S40192-022-00269-0/TABLES/1
Fey, L. T. W. & Beyerlein, I. J. Random generation of lattice structures with short-range order. Integr. Mater. Manuf. Innov. 11, 382–390. 10.1007/S40192-022-00269-0/TABLES/1 (2022).10.1007/S40192-022-00269-0/TABLES/1
51. Alhafez IA Deluigi OR Tramontina D Merkert N Urbassek HM Bringa EM Nanoindentation into a bcc high-entropy HfNbTaTiZr alloy—an atomistic study of the effect of short-range order Sci. Rep. 2024 14 9112 10.1038/s41598-024-59761-6 38643297
Alhafez, I. A. et al. Nanoindentation into a bcc high-entropy HfNbTaTiZr alloy—an atomistic study of the effect of short-range order. Sci. Rep. 14, 9112 (2024).38643297 10.1038/s41598-024-59761-6
52. de Fontaine D The number of independent pair-correlation functions in multicomponent systems J. Appl. Crystallogr. 1971 4 15 19 10.1107/S0021889871006174
de Fontaine, D. The number of independent pair-correlation functions in multicomponent systems. J. Appl. Crystallogr. 4, 15–19. 10.1107/S0021889871006174 (1971).10.1107/S0021889871006174
53. Zhang R Chen Y Fang Y Yu Q Characterization of chemical local ordering and heterogeneity in high-entropy alloys MRS Bull. 2022 47 186 193 10.1557/s43577-022-00277-7
Zhang, R., Chen, Y., Fang, Y. & Yu, Q. Characterization of chemical local ordering and heterogeneity in high-entropy alloys. MRS Bull. 47, 186–193. 10.1557/s43577-022-00277-7 (2022).10.1557/s43577-022-00277-7
54. Stukowski A Computational analysis methods in atomistic modeling of crystals JOM 2014 66 399 407 10.1007/s11837-013-0827-5
Stukowski, A. Computational analysis methods in atomistic modeling of crystals. JOM 66, 399–407. 10.1007/s11837-013-0827-5 (2014).10.1007/s11837-013-0827-5
55. Aparicio E Millán EN Ruestes CJ Bringa EM Foamexplorer: Automated measurement of ligaments and voids for atomistic systems Comput. Mater. Sci. 2020 185 109942 10.1016/j.commatsci.2020.109942
Aparicio, E., Millán, E. N., Ruestes, C. J. & Bringa, E. M. Foamexplorer: Automated measurement of ligaments and voids for atomistic systems. Comput. Mater. Sci. 185, 109942. 10.1016/j.commatsci.2020.109942 (2020).10.1016/j.commatsci.2020.109942
56. Gimenez MC Nieto F Ramirez-Pastor AJ Thermal percolation for interacting monomers adsorbed on square lattices J. Phys. A Math. Gen. 2005 38 3253 10.1088/0305-4470/38/15/002
Gimenez, M. C., Nieto, F. & Ramirez-Pastor, A. J. Thermal percolation for interacting monomers adsorbed on square lattices. J. Phys. A Math. Gen. 38, 3253 (2005).10.1088/0305-4470/38/15/002
57. Stevens RJ Zhigilei LV Norris PM Effects of temperature and disorder on thermal boundary conductance at solid-solid interfaces: Nonequilibrium molecular dynamics simulations Int. J. Heat Mass Transf. 2007 50 3977 3989 10.1016/j.ijheatmasstransfer.2007.01.040
Stevens, R. J., Zhigilei, L. V. & Norris, P. M. Effects of temperature and disorder on thermal boundary conductance at solid-solid interfaces: Nonequilibrium molecular dynamics simulations. Int. J. Heat Mass Transf. 50, 3977–3989. 10.1016/j.ijheatmasstransfer.2007.01.040 (2007).10.1016/j.ijheatmasstransfer.2007.01.040
58. Shafai G Ortigoza MA Rahman TS Vibrations of Au13 and FeAu12 nanoparticles and the limits of the Debye temperature concept J. Phys. Condens. Matter 2012 24 104026 10.1088/0953-8984/24/10/104026 22354075
Shafai, G., Ortigoza, M. A. & Rahman, T. S. Vibrations of Au and FeAu nanoparticles and the limits of the Debye temperature concept. J. Phys. Condens. Matter 24, 104026. 10.1088/0953-8984/24/10/104026 (2012).22354075 10.1088/0953-8984/24/10/104026
59. Seyf HR Gordiz K DeAngelis F Henry A Using Green–Kubo modal analysis (GKMA) and interface conductance modal analysis (ICMA) to study phonon transport with molecular dynamics J. Appl. Phys. 2019 10.1063/1.5081722
Seyf, H. R., Gordiz, K., DeAngelis, F. & Henry, A. Using Green–Kubo modal analysis (GKMA) and interface conductance modal analysis (ICMA) to study phonon transport with molecular dynamics. J. Appl. Phys. 125, 081101 10.1063/1.5081722 (2019).10.1063/1.5081722
60. Li M Deng T Zheng B Zhang Y Liao Y Zhou H Effect of defects on the mechanical and thermal properties of graphene Nanomaterials 2019 9 347 10.3390/nano9030347 30832437
Li, M. et al. Effect of defects on the mechanical and thermal properties of graphene. Nanomaterials 9, 347. 10.3390/nano9030347 (2019).30832437 10.3390/nano9030347
61. Yang X Tiwari J Feng T Reduced anharmonic phonon scattering cross-section slows the decrease of thermal conductivity with temperature Mater. Today Phys. 2022 24 100689 10.1016/j.mtphys.2022.100689
Yang, X., Tiwari, J. & Feng, T. Reduced anharmonic phonon scattering cross-section slows the decrease of thermal conductivity with temperature. Mater. Today Phys. 24, 100689. 10.1016/j.mtphys.2022.100689 (2022).10.1016/j.mtphys.2022.100689
62. Cahill, D. G., Braun, P. V., Chen, G., Clarke, D. R., Fan, S., Goodson, K. E., Keblinski, P., King, W. P., Mahan, G. D., Majumdar, A., Maris, H. J., Phillpot, S. R., Pop, E., and  Shi, L. Nanoscale thermal transport. II. 2003-2012, Appl. Phys. Rev. 1, 011305 (2014) 10.1063/1.4832615
63. Giri A Dionne CJ Hopkins PE Atomic coordination dictates vibrational characteristics and thermal conductivity in amorphous carbon npj Comput. Mater. 2022 8 55 10.1038/s41524-022-00741-7
Giri, A., Dionne, C. J. & Hopkins, P. E. Atomic coordination dictates vibrational characteristics and thermal conductivity in amorphous carbon. npj Comput. Mater. 8, 55. 10.1038/s41524-022-00741-7 (2022).10.1038/s41524-022-00741-7
64. Esfahani MN Jabbari M Xu Y Soutis C Effect of nanoscale defects on the thermal conductivity of graphene Mater. Today. Commun. 2021 26 101856 10.1016/j.mtcomm.2020.101856
Esfahani, M. N., Jabbari, M., Xu, Y. & Soutis, C. Effect of nanoscale defects on the thermal conductivity of graphene. Mater. Today. Commun. 26, 101856. 10.1016/j.mtcomm.2020.101856 (2021).10.1016/j.mtcomm.2020.101856
65. Xun K Zhang B Wang Q Zhang Z Ding J Ma E Local chemical inhomogeneities in TiZrNb-based refractory high-entropy alloys J. Mater. Sci. Technol. 2023 135 221 230 10.1016/j.jmst.2022.06.047
Xun, K. et al. Local chemical inhomogeneities in TiZrNb-based refractory high-entropy alloys. J. Mater. Sci. Technol. 135, 221–230. 10.1016/j.jmst.2022.06.047 (2023).10.1016/j.jmst.2022.06.047
66. Luis César R. Aliaga, Alexandre Melhorance Barboza, Loena Marins de Couto, and Ivan Napoleão Bastos, Molecular dynamics on Hf-Nb-Ta-Ti-Zr high entropy alloy, in 10.5772/intechopen.1004372High Entropy Alloys – Composition and Microstructure Design, edited by Yu Yin, Han Huang, Mingxing Zhang, and Libo Zhou (IntechOpen, Rijeka, 2024) Chap. 3
67. Senkov ON Pilchak AL Semiatin SL Effect of cold deformation and annealing on the microstructure and tensile properties of a HfNbTaTiZr refractory high entropy alloy Metall. Mater. Trans. A. 2018 49 2876 2892 10.1007/s11661-018-4646-8
Senkov, O. N., Pilchak, A. L. & Semiatin, S. L. Effect of cold deformation and annealing on the microstructure and tensile properties of a HfNbTaTiZr refractory high entropy alloy. Metall. Mater. Trans. A. 49, 2876–2892. 10.1007/s11661-018-4646-8 (2018).10.1007/s11661-018-4646-8
68. Schuh B Völker B Todt J Schell N Perrière L Li J Couzinié JP Hohenwarter A Thermodynamic instability of a nanocrystalline, single-phase TiZrNbHfTa alloy and its impact on the mechanical properties Acta Mater. 2018 142 201 212 10.1016/j.actamat.2017.09.035
Schuh, B. et al. Thermodynamic instability of a nanocrystalline, single-phase TiZrNbHfTa alloy and its impact on the mechanical properties. Acta Mater. 142, 201–212. 10.1016/j.actamat.2017.09.035 (2018).10.1016/j.actamat.2017.09.035
69. Chen SY Tong Y Tseng K-K Yeh J-W Poplawsky JD Wen JG Gao MC Kim G Chen W Ren Y Feng R Li WD Liaw PK Phase transformations of HfNbTaTiZr high-entropy alloy at intermediate temperatures Scripta Mater. 2019 158 50 56 10.1016/j.scriptamat.2018.08.032
Chen, S. Y. et al. Phase transformations of HfNbTaTiZr high-entropy alloy at intermediate temperatures. Scripta Mater. 158, 50–56. 10.1016/j.scriptamat.2018.08.032 (2019).10.1016/j.scriptamat.2018.08.032
70. Zeng Z Shen X Cheng R Perez O Ouyang N Fan Z Lemoine P Raveau B Guilmeau E Chen Y Pushing thermal conductivity to its lower limit in crystals with simple structures Nat. Commun. 2024 15 3007 10.1038/s41467-024-46799-3 38589376
Zeng, Z. et al. Pushing thermal conductivity to its lower limit in crystals with simple structures. Nat. Commun. 15, 3007. 10.1038/s41467-024-46799-3 (2024).38589376 10.1038/s41467-024-46799-3
71. Yang L Yue S Tao Y Qiao S Li H Dai Z Song B Chen Y Du J Li D Gao P Suppressed thermal transport in silicon nanoribbons by inhomogeneous strain Nature 2024 10.1038/s41586-024-07390-4 39223371
Yang, L. et al. Suppressed thermal transport in silicon nanoribbons by inhomogeneous strain. Nature 629, 1021 10.1038/s41586-024-07390-4 (2024).39223371 10.1038/s41586-024-07390-4
72. Stukowski A Visualization and analysis of atomistic simulation data with OVITO—the open visualization tool Model. Simul. Mater. Sci. Eng. 2010 18 015012 10.1088/0965-0393/18/1/015012
Stukowski, A. Visualization and analysis of atomistic simulation data with OVITO—the open visualization tool. Model. Simul. Mater. Sci. Eng. 18, 015012. 10.1088/0965-0393/18/1/015012 (2010).10.1088/0965-0393/18/1/015012
73. Zhang A Li Y Thermal conductivity of aluminum alloys-a review Materials 2023 16 2972 10.3390/ma16082972 37109807
Zhang, A. & Li, Y. Thermal conductivity of aluminum alloys-a review. Materials 16, 2972. 10.3390/ma16082972 (2023).37109807 10.3390/ma16082972
74. Xing B Rupert TJ Pan X Cao P Neural network kinetics for exploring diffusion multiplicity and chemical ordering in compositionally complex materials Nat. Commun. 2024 15 3879 10.1038/s41467-024-47927-9 38724515
Xing, B., Rupert, T. J., Pan, X. & Cao, P. Neural network kinetics for exploring diffusion multiplicity and chemical ordering in compositionally complex materials. Nat. Commun. 15, 3879 (2024).38724515 10.1038/s41467-024-47927-9
75. Zhang M Walsh F Asta M Ritchie RO Minor AM Diffuse electron diffraction intensities in concentrated solid solutions do not necessarily come from short-range order Microsc. Microanal. 2023 29 437 438 10.1093/MICMIC/OZAD067.207 37613268
Zhang, M., Walsh, F., Asta, M., Ritchie, R. O. & Minor, A. M. Diffuse electron diffraction intensities in concentrated solid solutions do not necessarily come from short-range order. Microsc. Microanal. 29, 437–438. 10.1093/MICMIC/OZAD067.207 (2023).37613268 10.1093/MICMIC/OZAD067.207
76. Turner SR Pailhès S Bourdarot F Ollivier J Sidis Y Castellan J-P Zanotti J-M Berrod Q Porcher F Bosak A Phonon behavior in a random solid solution: A lattice dynamics study on the high-entropy alloy FeCoCrMnNi Nat. Commun. 2022 13 7509 10.1038/s41467-022-35125-4 36473859
Turner, S. R. et al. Phonon behavior in a random solid solution: A lattice dynamics study on the high-entropy alloy FeCoCrMnNi. Nat. Commun. 13, 7509 (2022).36473859 10.1038/s41467-022-35125-4
77. Körmann F Ikeda Y Grabowski B Sluiter MHF Phonon broadening in high entropy alloys Npj Comput. Mater. 2017 3 36 10.1038/s41524-017-0037-8
Körmann, F., Ikeda, Y., Grabowski, B. & Sluiter, M. H. F. Phonon broadening in high entropy alloys. Npj Comput. Mater. 3, 36 (2017).10.1038/s41524-017-0037-8
78. Hashin Z Shtrikman S A variational approach to the theory of the effective magnetic permeability of multiphase materials J. Appl. Phys. 1962 33 3125 3131 10.1063/1.1728579
Hashin, Z. & Shtrikman, S. A variational approach to the theory of the effective magnetic permeability of multiphase materials. J. Appl. Phys. 33, 3125–3131. 10.1063/1.1728579 (1962).10.1063/1.1728579
79. Kirkpatrick S Percolation and conduction Rev. Mod. Phys. 1973 45 574 588 10.1103/RevModPhys.45.574
Kirkpatrick, S. Percolation and conduction. Rev. Mod. Phys. 45, 574–588. 10.1103/RevModPhys.45.574 (1973).10.1103/RevModPhys.45.574
80. Carson JK Lovatt SJ Tanner DJ Cleland AC Thermal conductivity bounds for isotropic, porous materials Int. J. Heat Mass Transf. 2005 48 2150 2158 10.1016/j.ijheatmasstransfer.2004.12.032
Carson, J. K., Lovatt, S. J., Tanner, D. J. & Cleland, A. C. Thermal conductivity bounds for isotropic, porous materials. Int. J. Heat Mass Transf. 48, 2150–2158. 10.1016/j.ijheatmasstransfer.2004.12.032 (2005).10.1016/j.ijheatmasstransfer.2004.12.032
81. Wang J Carson JK North MF Cleland DJ A new structural model of effective thermal conductivity for heterogeneous materials with co-continuous phases Int. J. Heat Mass Transf. 2008 51 2389 2397 10.1016/j.ijheatmasstransfer.2007.08.028
Wang, J., Carson, J. K., North, M. F. & Cleland, D. J. A new structural model of effective thermal conductivity for heterogeneous materials with co-continuous phases. Int. J. Heat Mass Transf. 51, 2389–2397. 10.1016/j.ijheatmasstransfer.2007.08.028 (2008).10.1016/j.ijheatmasstransfer.2007.08.028
82. Li M Cao G Luo Y Sheng C Predicting the lattice thermal conductivity of alloyed compounds from the perspective of configurational entropy npj Comput. Mater. 2022 8 75 10.1038/s41524-022-00771-1
Li, M., Cao, G., Luo, Y. & Sheng, C. Predicting the lattice thermal conductivity of alloyed compounds from the perspective of configurational entropy. npj Comput. Mater. 8, 75. 10.1038/s41524-022-00771-1 (2022).10.1038/s41524-022-00771-1
83. Cahill DG Watson SK Pohl RO Lower limit to the thermal conductivity of disordered crystals Phys. Rev. B 1992 46 6131 6140 10.1103/PhysRevB.46.6131
Cahill, D. G., Watson, S. K. & Pohl, R. O. Lower limit to the thermal conductivity of disordered crystals. Phys. Rev. B 46, 6131–6140. 10.1103/PhysRevB.46.6131 (1992).10.1103/PhysRevB.46.6131
84. Morelli, D. T., Slack, G. A. High lattice thermal conductivity solids, in High Thermal Conductivity Materials (Springer, 2006) pp. 37–68
85. Ho, C. Y., Powell, R. W. & Liley, P. E. Thermal Conductivity of the Elements: A Comprehensive Review (American Institute of Physics, New York, 1974).
86. Mills LH Emigh MG Frey CH Philips NR Murray SP Shin J Gianola DS Pollock TM Temperature-dependent tensile behavior of the HfNbTaTiZr multi-principal element alloy Acta Mater. 2023 245 118618 10.1016/j.actamat.2022.118618
Mills, L. H. et al. Temperature-dependent tensile behavior of the HfNbTaTiZr multi-principal element alloy. Acta Mater. 245, 118618 (2023).10.1016/j.actamat.2022.118618
87. Farkas D Caro A Model interatomic potentials for Fe–Ni–Cr–Co–Al high-entropy alloys J. Mater. Res. 2020 35 3031 3040 10.1557/JMR.2020.294
Farkas, D. & Caro, A. Model interatomic potentials for Fe–Ni–Cr–Co–Al high-entropy alloys. J. Mater. Res. 35, 3031–3040. 10.1557/JMR.2020.294 (2020).10.1557/JMR.2020.294
88. Babicheva R Jarlöv A Zheng H Dmitriev S Korznikova E Nai MLS Ramamurty U Zhou K Effect of short-range ordering and grain boundary segregation on shear deformation of CoCrFeNi high-entropy alloys with Al addition Comput. Mater. Sci. 2022 215 111762 10.1016/J.COMMATSCI.2022.111762
Babicheva, R. et al. Effect of short-range ordering and grain boundary segregation on shear deformation of CoCrFeNi high-entropy alloys with Al addition. Comput. Mater. Sci. 215, 111762. 10.1016/J.COMMATSCI.2022.111762 (2022).10.1016/J.COMMATSCI.2022.111762
89. Li G Yuan L Li J Zhang M Li D Effect of al related chemical short-range order on the irradiation resistance of Al0.5CoCrFeNi multi-principal element alloy: A molecular dynamics simulation Mater. Today Commun. 2024 38 108302 10.1016/J.MTCOMM.2024.108302
Li, G., Yuan, L., Li, J., Zhang, M. & Li, D. Effect of al related chemical short-range order on the irradiation resistance of AlCoCrFeNi multi-principal element alloy: A molecular dynamics simulation. Mater. Today Commun. 38, 108302. 10.1016/J.MTCOMM.2024.108302 (2024).10.1016/J.MTCOMM.2024.108302
90. Moschetti M Alan X Hohenwarter A Wei T Davis J Short K Thorogood GJ Kong C Couzinié J-P Bhattacharyya D Kruzic JJ Gludovatz B The influence of phase formation on irradiation tolerance in a nanocrystalline TiZrNbHfTa refractory high-entropy alloy Adv. Eng. Mater. 2024 26 2300863 10.1002/adem.202300863
Moschetti, M. et al. The influence of phase formation on irradiation tolerance in a nanocrystalline TiZrNbHfTa refractory high-entropy alloy. Adv. Eng. Mater. 26, 2300863. 10.1002/adem.202300863 (2024).10.1002/adem.202300863
91. Rong C Yang J Zhao X Huang K Liu Y Wang X Zhu D Chen R Microstructure recrystallization and mechanical properties of a cold-rolled TiNbZrTaHf refractory high-entropy alloy Acta Metall. Sinica 2024 37 633 647 10.1007/S40195-023-01649-X/METRICS
Rong, C. et al. Microstructure recrystallization and mechanical properties of a cold-rolled TiNbZrTaHf refractory high-entropy alloy. Acta Metall. Sinica 37, 633–647. 10.1007/S40195-023-01649-X/METRICS (2024).10.1007/S40195-023-01649-X/METRICS
92. Wu Y Yu W Shen S Developing a variable charge potential for Hf/Nb/Ta/Ti/Zr/O system via machine learning global optimization Mater. Des. 2023 230 111999 10.1016/J.MATDES.2023.111999
Wu, Y., Yu, W. & Shen, S. Developing a variable charge potential for Hf/Nb/Ta/Ti/Zr/O system via machine learning global optimization. Mater. Des. 230, 111999. 10.1016/J.MATDES.2023.111999 (2023).10.1016/J.MATDES.2023.111999
93. Bhandari U Chen Y Ding H Zeng C Emanet S Gradl PR Guo S Machine-learning-based thermal conductivity prediction for additively manufactured alloys J. Manuf. Mater. Process. 2023 7 160 10.3390/jmmp7050160
Bhandari, U. et al. Machine-learning-based thermal conductivity prediction for additively manufactured alloys. J. Manuf. Mater. Process. 7, 160. 10.3390/jmmp7050160 (2023).10.3390/jmmp7050160
94. Afremov LL Kirienko YV Petrov AA Chepak AK Size effect’s influence on the magnetic phase transitions in the nanosized magnets J. Supercond. Novel Magn. 2023 36 587 600 10.1007/s10948-022-06482-4
Afremov, L. L., Kirienko, Y. V., Petrov, A. A. & Chepak, A. K. Size effect’s influence on the magnetic phase transitions in the nanosized magnets. J. Supercond. Novel Magn. 36, 587–600 (2023).10.1007/s10948-022-06482-4
95. Minnich AJ Advances in the measurement and computation of thermal phonon transport properties J. Phys. Condens. Matter 2015 27 053202 10.1088/0953-8984/27/5/053202 25603881
Minnich, A. J. Advances in the measurement and computation of thermal phonon transport properties. J. Phys. Condens. Matter 27, 053202. 10.1088/0953-8984/27/5/053202 (2015).25603881 10.1088/0953-8984/27/5/053202
96. Dong H Shi Y Ying P Ke X Liang T Wang Y Zeng Z Xin W Zhou W Xiong S Chen S Fan Z Molecular dynamics simulations of heat transport using machine-learned potentials: A mini-review and tutorial on GPUMD with neuroevolution potentials J. Appl. Phys. 2024 135 161101 10.1063/5.0200833
Dong, H. et al. Molecular dynamics simulations of heat transport using machine-learned potentials: A mini-review and tutorial on GPUMD with neuroevolution potentials. J. Appl. Phys. 135, 161101. 10.1063/5.0200833 (2024).10.1063/5.0200833
97. Hu Y Zeng L Minnich AJ Dresselhaus MS Chen G Spectral mapping of thermal conductivity through nanoscale ballistic transport Nat. Nanotechnol. 2015 10 701 706 10.1038/nnano.2015.109 26030656
Hu, Y., Zeng, L., Minnich, A. J., Dresselhaus, M. S. & Chen, G. Spectral mapping of thermal conductivity through nanoscale ballistic transport. Nat. Nanotechnol. 10, 701–706. 10.1038/nnano.2015.109 (2015).26030656 10.1038/nnano.2015.109
98. Li S Yang X Hou J Du W A review on thermal conductivity of magnesium and its alloys J. Mag. Alloys 2020 8 78 90 10.1016/j.jma.2019.08.002
Li, S., Yang, X., Hou, J. & Du, W. A review on thermal conductivity of magnesium and its alloys. J. Mag. Alloys 8, 78–90. 10.1016/j.jma.2019.08.002 (2020).10.1016/j.jma.2019.08.002
99. Bhatt N Karna P Thakur S Giri A Transition from electron-dominated to phonon-driven thermal transport in tungsten under extreme pressures Phys. Rev. Mater. 2023 7 115001 10.1103/PhysRevMaterials.7.115001
Bhatt, N., Karna, P., Thakur, S. & Giri, A. Transition from electron-dominated to phonon-driven thermal transport in tungsten under extreme pressures. Phys. Rev. Mater. 7, 115001. 10.1103/PhysRevMaterials.7.115001 (2023).10.1103/PhysRevMaterials.7.115001
100. Mu S Uncovering electron scattering mechanisms in NiFeCoCrMn derived concentrated solid solution and high entropy alloys npj Comput. Mater. 2019 5 1 10.1038/s41524-018-0138-z
Mu, S. et al. Uncovering electron scattering mechanisms in NiFeCoCrMn derived concentrated solid solution and high entropy alloys. npj Comput. Mater. 5, 1. 10.1038/s41524-018-0138-z (2019).10.1038/s41524-018-0138-z
101. Raghuraman V Wang Y Widom M An investigation of high entropy alloy conductivity using first-principles calculations Appl. Phys. Lett. 2021 119 121903 10.1063/5.0065239
Raghuraman, V., Wang, Y. & Widom, M. An investigation of high entropy alloy conductivity using first-principles calculations. Appl. Phys. Lett. 119, 121903. 10.1063/5.0065239 (2021).10.1063/5.0065239
102. Kittel C McEuen P Introduction to Solid State Physics 2018 9 New Jersey Wiley
Kittel, C. & McEuen, P. Introduction to Solid State Physics 9th edn. (Wiley, New Jersey, 2018).
103. Luo Y Li M Yuan H Liu H Fang Y Predicting lattice thermal conductivity via machine learning: A mini review npj Comput. Mater. 2023 10.1038/s41524-023-00964-2
Luo, Y., Li, M., Yuan, H., Liu, H. & Fang, Y. Predicting lattice thermal conductivity via machine learning: A mini review. npj Comput. Mater. 9, 4 10.1038/s41524-023-00964-2 (2023).10.1038/s41524-023-00964-2
