Literature DB >> 33718829

High-entropy materials for energy-related applications.

Maosen Fu1, Xiao Ma1, Kangning Zhao2, Xiao Li2, Dong Su2.   

Abstract

High-entropy materials (HEMs), including high-entropy alloys (HEAs), high-entropy oxides (HEOs), and other high-entropy compounds, have gained significant interests over the past years. These materials have unique structures with the coexistence of antisite disordering and crystal periodicity, which were originally investigated as structural materials. Recently, they have emerged for energy-related applications, such as catalysis, energy storage, etc. In this work, we review the research progress of energy-related applications of HEMs. After an introduction on the background, theory, and syntheses of HEMs, we survey their applications including electrocatalysis, batteries, and others, aiming to retrieve the correlations between their structures and performances. In the end, we discussed the challenges and future directions for developing HEMs.
© 2021 The Author(s).

Entities:  

Keywords:  Energy Materials; Materials Science; Physics

Year:  2021        PMID: 33718829      PMCID: PMC7921604          DOI: 10.1016/j.isci.2021.102177

Source DB:  PubMed          Journal:  iScience        ISSN: 2589-0042


Introduction

High-entropy alloys (HEAs), also known as multi-principal element alloys, emerged into the academic literature in 2004 and have grown in applications and understanding since. They are multicomponent solid solution materials that contain five or more elements in near-equiatomic proportions to stabilize their structures by maximizing the configurational entropy (Cantor et al., 2004; Yeh et al., 2004). These unconventional structures provide opportunities for achieving unprecedented combinations of phase stability and mechanical performance (Ding et al., 2019; Gludovatz et al., 2014; Lei et al., 2018), especially overcoming the strength-ductility trade-off (Li et al., 2016). Although there are still arguments on if these multicomponent alloys are really stabilized by entropy or not (Löffler et al., 2019), we here call it HEAs for its popularity in the literature (George et al., 2019). It is believed that the excellent performances of HEAs are promoted by four “core effects”: the high-entropy effect, the lattice distortion effect, the sluggish diffusion effect, and the “cocktail” effect, all of which are pivotal to structure-property studies in this field (Miracle and Senkov, 2017a). HEAs are defined from a signature concept, the high-entropy effect, which highlights that when five or more elements are combined they may favor forming of a solid solution over an intermetallic compound (Yeh et al., 2004). For our view, the most unique structural character of HEAs is the antisite disordering of atomic types within an ordered crystal. This character could be extended to compounds and other applications. Rost et al. have synthesized rock-salt (NiCoCuZnMg)Ox solid solution oxide (Rost et al., 2015) and these oxides are named as high-entropy oxides (HEOs) because they are also stabilized by high configurational entropy (Batchelor et al., 2019; Yan and Zhang, 2020). Up to now, numerous functional HEOs, high-entropy metallic glass (HEMGs) nanoparticles (NPs) (Glasscott et al., 2019), high-entropy borides (HEBs) (Gild et al., 2016), high-entropy sulfides (HESs) (Zhang et al., 2018b), high-entropy carbides (HECs) (Sarker et al., 2018), and high-entropy nitrides (HENs) (Jin et al., 2018) have also been developed. High-entropy materials (HEMs) including above compounds, are defined as any solid solution materials that consist of quasi-equimolar multicomponent (Dai, 2020). Recent works show that HEMs can be utilized for various energy-related applications, as catalysts for ammonia oxidation and decomposition (Xie et al., 2019; Yao et al., 2018); as electrocatalysts for oxidation (methanol [Wang et al., 2014)], carbon monoxide [Chen et al., 2018a; Dai, 2020]), evolution (oxygen [Qiu et al., 2019a], hydrogen [Zhang et al., 2018a]), reduction (carbon dioxide [Nellaiappan et al., 2020], oxygen [Löffler et al., 2018]), and degradation (azo dye [Lv et al., 2016]); as electrode for batteries (Dragoe and Bérardan, 2019; Oses et al., 2020); as hydrogen storage (Shen et al., 2019); as supercapacitors (Kong et al., 2019); etc. The structural tunability of HEMs provides them a great possibility for energy-related applications. For example, the large number of atomic sites on the surface of HEMs provides huge opportunity to tune the surface electronic structures, adsorption energies, and then the catalytic performances. HEMs for energy-related applications are usually in the form of nanoporous alloys or NPs to maximize their surface areas to prompt the reaction kinetics. Nanosized HEMs can be synthesized with wet-chemistry (solvothermal, ultrasonicated-assisted wet-chemistry, sol-gel auto-combustion, etc.), pulsed laser deposition (PLD), electrodeposition, carbothermal shock (CTS) methods, and so on, where processing temperature, heating, and cooling rates are key factors to form stable nanosized HEMs. In addition to thermodynamic stability, the uniform configuration of elements with different sizes through a single-phase solid solution also leads to severe lattice distortion, which then increases the activation barrier of diffusion (the sluggish diffusion effect) and improves kinetic stability (Usharani et al., 2020). The applications of HEMs for catalysis or energy storage have been investigated for less than one decade, and lots of explorations are needed for understanding well their structure-performance correlations. In this work, we summarize the research progress of energy-related applications of HEMs. After a survey on the syntheses of HEMs, we introduce the structure and theory of HEMs and then the applications of nanosized HEMs, correlating with their structures, as shown in the roadmap of Figure 1. To keep the article brief, we skip the background knowledge of each application and focus on the structure-performance correlations of HEMs. In the last section, we present our perspectives on the challenges and future research opportunities.
Figure 1

Roadmap of HEMs for recent energy-related applications

Before 2015, HEMs were mainly used as structural materials (Reproduced with permission from Gludovatz et al., 2014. Copyright © 2014, American Association for the Advancement of Science). Rost et al. discovered the entropy-stabilized oxides and named them as HEOs (Reproduced with permission from Rost et al., 2015. Copyright © 2015, Springer Nature). Then many other HEMs have been synthesized and utilized in solid electrolyte (Reproduced with permission from Bérardan et al., 2016a, 2016b. Copyright © 2016, The Royal Society of Chemistry), exchange coupling (Reproduced with permission from Meisenheimer et al., 2017. Copyright © 2017, Springer Nature), water decomposition (Reproduced with permission from Zhai et al., 2018. Copyright © 2018, The Royal Society of Chemistry), Li-ion battery (Reproduced with permission Wang et al., 2019. Copyright © 2019, Elsevier), COOR/CO2RR catalysts (Reproduced with permission from Chen et al., 2019a. Copyright © 2019, American Chemical Society), and Na-ion battery (Reproduced with permission from Zhao et al., 2020a. Copyright © 2020, Wiley-VCH). Simultaneously, predictive models have been developed to guide the development of HEMs (Reproduced with permission from Löffler et al., 2020. Copyright © 2020, Wiley-VCH).

Roadmap of HEMs for recent energy-related applications Before 2015, HEMs were mainly used as structural materials (Reproduced with permission from Gludovatz et al., 2014. Copyright © 2014, American Association for the Advancement of Science). Rost et al. discovered the entropy-stabilized oxides and named them as HEOs (Reproduced with permission from Rost et al., 2015. Copyright © 2015, Springer Nature). Then many other HEMs have been synthesized and utilized in solid electrolyte (Reproduced with permission from Bérardan et al., 2016a, 2016b. Copyright © 2016, The Royal Society of Chemistry), exchange coupling (Reproduced with permission from Meisenheimer et al., 2017. Copyright © 2017, Springer Nature), water decomposition (Reproduced with permission from Zhai et al., 2018. Copyright © 2018, The Royal Society of Chemistry), Li-ion battery (Reproduced with permission Wang et al., 2019. Copyright © 2019, Elsevier), COOR/CO2RR catalysts (Reproduced with permission from Chen et al., 2019a. Copyright © 2019, American Chemical Society), and Na-ion battery (Reproduced with permission from Zhao et al., 2020a. Copyright © 2020, Wiley-VCH). Simultaneously, predictive models have been developed to guide the development of HEMs (Reproduced with permission from Löffler et al., 2020. Copyright © 2020, Wiley-VCH).

Syntheses and structures

There are many physical or chemical synthetic routes being developed in the past 20 years to fabricate HEMs (Amiri and Shahbazian-Yassar, 2021; Tomboc et al., 2020; Xin et al., 2020). In the following discussion, we will briefly introduce separately the synthetic approaches of HEAs, HEOs, and other HEMs, correlating to their structural features.

Syntheses and structures of HEAs

Syntheses of HEAs

For energy-related applications, HEAs usually adopt nanostructures as nanoporous or NPs to maximize the surface area. Dealloying of bulk alloy is a common method to fabricate nanoporous HEA (such as AlMoCuPdAu, AlMnCoIrMo, AlNiCuPtPdAu, etc.) with enhanced surface area and uniform pore structures, ensuring high catalytic activity (Qiu et al., 2019b). Wet-chemistry is the general method to fabricate NPs through reduction of metallic precursors in solution. The general solution-based common chemical synthesis might not be adequate for obtaining homogeneous nano-HEA with multi-components, and often severe decomposition or shock syntheses or reduction reactions are required (Bondesgaard et al., 2019). Liu et al. take advantage of ultrasonication to develop a facile ultrasonication-assisted wet-chemistry for preparing HEA NPs. The local temperature of the ultrasonication cavities is supposed to exceed extremely high temperature up to ∼ 5,000°C, directly accelerating the metallic ions reduction to form entropy-maximized state, facilitating the formation of HEAs (Liu et al., 2019b). In another case, Yao et al. developed the CTS method to prepare multicomponent HEA NPs. Precursor of metal salt mixtures was thermally shocked at ∼ 2,000 K within 55 ms (as shown in Figure 2A), alloying up to eight dissimilar elements into single-phase solid-solution NPs (PtPdCoNiFeCuAuSn, etc.) loaded onto carbon support. This method can synthesize a wide range of multicomponent NPs with desired compositions and proper sizes by controlling the CTS parameters as well as supports, which may be a general route toward high performance and cost-effective catalysts (Yao et al., 2018). However, productivity of HEA NPs via CTS is very low. Future work to enhance the productivity for industrial applications is needed. Furthermore, Gao et al. developed a fast-moving bed pyrolysis (FMBP) method to prepare HEAs with diverse supports. Mixed metal precursors were rapidly heated to ∼ 923 K (higher than their pyrolysis temperature) within 5 s. With this strategy, quinary FeCoPdIrPt, senary AuPdPtCuNiSn, octonary AuIrPdPtCoCuNiSn, and denary AuIrPdPtRhCoCuMnNiSn HEA NPs are successfully synthesized, exhibiting high mass activity and stability (Gao et al., 2020).
Figure 2

Syntheses of nanosized HEMs

(A) Carbothermal synthesis of HEA-NPs on carbon supports. Reproduced with permission from Yao et al. (2018). Copyright © 2018, American Association for the Advancement of Science.

(B) Nanodroplet-mediated electrodeposition synthesis controlling nanoparticle stoichiometry and microstructure. Reproduced with permission from Glasscott et al. (2019). Copyright © 2019, Springer Nature.

Syntheses of nanosized HEMs (A) Carbothermal synthesis of HEA-NPs on carbon supports. Reproduced with permission from Yao et al. (2018). Copyright © 2018, American Association for the Advancement of Science. (B) Nanodroplet-mediated electrodeposition synthesis controlling nanoparticle stoichiometry and microstructure. Reproduced with permission from Glasscott et al. (2019). Copyright © 2019, Springer Nature. Table 1 has summarized the methods for preparing HEAs. For example, solvothermal (Bondesgaard et al., 2019), co-sputtering (Löffler et al., 2018), and deposition (Glasscott et al., 2019) have also been successfully developed to prepare HEA NPs. Electrodeposition can be used to produce multi-component HEMGs as PtCoFeLaNi, by confining multiple metal salt precursors to water nanodroplets emulsified in dichloroethane within 100 ms electroshock as shown in Figure 2B (Glasscott et al., 2019). The contents of nanodroplet are entirely reduced during collision, enabling disordered co-deposition of different metal precursors. CrMnFeCoNi NPs with an averaged diameter of 1.7 nm were synthesized by combinatorial co-sputtering into the [Bmim]-[Tf2N] (Löffler et al., 2018). Bondesgaard et al. developed a general solvothermal synthesis method for PtPdIrRhRu HEA nanocatalysts (Bondesgaard et al., 2019).
Table 1

Synthesis, structure, and application of HEMs

ApplicationsCompositionStructureFormSynthesisReferences
Ammonia oxidationPtPdRhRuCe, PtCoNiFeCuAu, PtPdCoNiCuAu, PtPdCoNiFeCuAuSnfccNanoparticleCarbothermal shock(Yao et al., 2018)
Ammonia decompositionCoxMoyFe10Ni10Cu10 (x + y = 70)fccNanoparticleCarbothermal shock(Xie et al., 2019)
HERCoFeLaNiPtfccNanoparticleNanodroplet-mediated electrodeposition(Glasscott et al., 2019)
PtAuPdRhRufccNanoparticleUltrasonication(Liu et al., 2019b)
IrPdPtRhRufccNanoparticleFacile one-pot polyol synthesis(Wu et al., 2020a)
AlMoCuPdAufccNanoporous alloyDealloying(Qiu et al., 2019b)
FeCoPdIrPtfccNanoparticleFast-moving bed pyrolysis(Gao et al., 2020)
Cr15Fe20Co35Ni20Mo10fccBulkArc-melting(Zhang et al., 2018a)
(FeMnCoNi)xORock-salt + spinelMicro-powdersSol-gel and solid-state sintered(Zhai et al., 2018)
(CrMnFeCoNi)xPHexagonalNanosheetEutectic solvent method(Zhao et al., 2020b)
OERCrMnFeCoNifccNanoparticleLaser(Waag et al., 2019)
AlMoCuPdAufccNanoporous alloyDealloying(Qiu et al., 2019b)
CoFeLaNiPtfccNanoparticleNanodroplet-mediated electrodeposition(Glasscott et al., 2019)
AlMoCoIrMofccNanoporous alloyDealloying(Jin et al., 2019)
(FeMnCoNi)Ox (x ~ 1.2)Rock-salt + spinelMicro-powdersSol-gel and solid-state sintered(Zhai et al., 2018)
FeMnCoNi (+(FeMnCoNi)Ox)fccMicro-powdersMechanical alloying(Dai et al., 2019)
(CrMnFeCoNi)xPHexagonalNanosheetEutectic solvent method(Zhao et al., 2020b)
K(MgMnFeCoNi)F3, K(MgMnCoNiZn)F3PerovskiteNanoparticleBoiled solution by combining a hydrothermal method with mechanochemistry(Wang et al., 2020b)
ORRCrMnFeCoNiNb, CrMnFeCoNiMofccNanoparticleCombinatorial co-sputtering(Löffler et al., 2018)
AlNiCuPtPdAufccNanoporousDealloying(Qiu et al., 2019b)
Al-Cu-Ni-Pt-MnfccNanoporous alloySelective dealloying(Li et al., 2020b)
COORAlMoCuPdAufccNanoporous alloyDealloying(Qiu et al., 2019b)
PtNiMgCuZnCoOxRock-saltMicro-powderSolid-state sintered(Chen et al., 2018a)
Ru/BaSrBi(ZrHfTiFe)O3PerovskiteNanoparticleSonochemical-based synthesis(Okejiri et al., 2020)
(CuNiFeCoMg)Ox−Al2O3Rock-saltMesoporousMechanochemical solid state(Zhang et al., 2019b)
CO2RRAuAgPtPdCufccNanocrystallineCast cum cryo-milling(Nellaiappan et al., 2020)
Pt/Ru-(NiMgCuZnCo)OfccMicro-powderMechanochemical(Chen et al., 2019a)
Methanol oxidation reactionIr0.19Os0.22Re0.21Rh0.20Ru0.19hcpBulkThermal decomposition(Yusenko et al., 2017)
AlMoCuPdAufccNanoporous alloyDealloying(Qiu et al., 2019b)
Ethanol oxidation reactionRuRhPdOsIrPtfccNanoparticleWet-chemistry(Wu et al., 2020b)
Water oxidation(Co,Cu,Fe,Mn,Ni)3O4SpinelNanoparticleFacile methodology(Wang et al., 2019)
Oxidation of aromatic alcohols(Mg0.2Co0.2Ni0.2Cu0.2Zn0.2)ORock-saltHoley lamellarAnchoring-merging process(Feng et al., 2020a)
Li-ion battery(Co0.2Cu0.2Mg0.2Ni0.2Zn0.2)ORock-saltNanoparticleFlame spray pyrolysis, Nebulized spray pyrolysis, Reverse co-precipitation(Sarkar et al., 2018)
(Co0.2Cu0.2Mg0.2Ni0.2Zn0.2)ORock-saltNanoparticleSolid-state sintered & ball milling(Chen et al., 2019b)
(MgCoNiZn)1−xLixO (x = 0.05, 0.15, 0.25, and 0.35)Rock-saltBulkSolid-state sintered(Lökçü et al., 2020)
(Mg0.2Ti0.2Zn0.2Cu0.2Fe0.2)3O4SpinelPowderSolid-state sintered(Chen et al., 2020b)
(Mg,Ti,Zn,Cu,Fe)3O4SpinelNanoparticleSurfactant-assisted hydrothermal(Nguyen et al., 2020)
[(Bi,NA)1/5(La,Li)1/5(Ce,K)1/5Ca1/5Sr1/5]TiO3PerovskiteBulkSolid-state sintered(Yan et al., 2020)
Na-ion batteryNaNi0.12Cu0.12Mg0.12Fe0.15Co0.15Mn0.1Ti0.1Sn0.1Sb0.04O2O3Micro-powderSolid-state sintered(Zhao et al., 2020a)
Solid oxide fuel cellLa1-xSrx(CoCrFeMnNi)O3-δ (x = 0, 0.1, 0.2, 0.3)PerovskiteBulkSol-gel & solid-state sintered(Dąbrowa et al., 2020)
Solid-state electrolytes(MgCoNiCuZn)1-x-yGayAxO (A = Li+, Na+, K+)Rock-saltBulkSolid-state sintered(Bérardan et al., 2016b)
Semiconductor(Cr,Fe,Mg,Mn,Ni)3O4SpinelBulkSolid-state sintered(Stygar et al., 2020)
Hydrogen storageTiZrHfMoNbbccMicro-powderArc-melting & grinding(Shen et al., 2019)
SupercapacitorsFeNiCoMnMgfccNanoparticleCarbothermal shock(Xu et al., 2020)
AlCoCrFeNifccNanoporous alloySelective dissolution(Kong et al., 2019)
(TiNbTaZrHf)CfccNanoparticleFacile electrochemical(Sure et al., 2020)
(CrMoVZrNb)NfccMicro-powderMechanochemical-assisted soft urea strategy(Jin et al., 2018)
Colossal dielectric(Mg,Co,Ni,Cu,Zn)1-xLixORock-saltBulkSolid-state sintered(Bérardan et al., 2016a)
Exchange coupling(Mg0.25(1-x)CoxNi0.25(1-x)Cu0.25(1-x)Zn0.25(1-x))ORock-saltFilmSolid-state sintered & pulsed laser deposition(Meisenheimer et al., 2017)
Wear(Al0.17Ti0.41V0.14Cr0.04Nb0.24)N, (Al0.31Ti0.34V0.12Cr0.06Nb0.17)NfccFilmCathodic arc deposition(Yalamanchili et al., 2017)
High-temperature oxidation resistance(Hf0.2Zr0.2Ta0.2Mo0.2Ti0.2)B2Hexagonal layeredBulkSpark plasma sintering(Gild et al., 2016)
Synthesis, structure, and application of HEMs It was proposed that the key factor to synthesize HEA NPs is to control the cooling rate (Yao et al., 2018). A proper cooling rate, ranging from several seconds to several minutes depending on the system or synthesis method, facilitated the formation of single-phase solid solutions with certain crystal structures, enabling HEA structures. However, we believe that in most cases, reduction of all components at the same time in solutions is a key issue instead of the cooling rate, which is evidenced by the work by Gao et al. and Wu et al., showing that a slow cooling rate may also prompt the formation of HEAs (Gao et al., 2020; Wu et al., 2020a).

Structure of HEAs

HEAs can have cubic crystal structure with fcc lattice (Figure 3A) and bcc lattice (Figure 3B), hexagonal crystal structure (Figure 3C), and other structures (Ostovari Moghaddam and Trofimov, 2021; Yan and Zhang, 2020; Ye et al., 2016). Fcc lattice is by far the most extensively investigated structure of HEAs (George et al., 2019; Gorsse et al., 2017; Li et al., 2019; Miracle and Senkov, 2017b). For energy-related applications, nano-HEAs such as CoMoFeNiCu and CoFeLaNiPt with fcc lattice have a lattice parameter of a ≈ 3.6 Å. HEAs containing heavy hcp metals (mainly Os, Re, Ru, and Zr) could form hcp structure. Yusenko et al. synthesized Ir0.19Os0.22Re0.21Rh0.20Ru0.19 HEA with space group P63/mmc and lattice parameters а = 2.728 Å and c = 4.338 Å (Yusenko et al., 2017). Reports on the syntheses and structures of HEA NPs with different crystal structures are summarized in Table 2. The fundamental correlations between the compositions and phases, e.g., a phase diagram for HEAs, are not clear at present, which is highly needed in future.
Figure 3

Typical structures of HEAs, HEOs, and other HEMs

(A) fcc lattice structure, (B) bcc lattice structure, (C) hexagonal crystal structure, (D) rock-salt structure, (E) spinel structure, and (F) perovskite structure.

Table 2

Synthesis and structure of HEA NPs

StructureCompositionSynthetic methodFormCatalytic applicationReferences
fccPtRuCuOsIrMechanical alloying (chemical dealloying)NanoporousMethanol oxidation & oxygen reduction(Chen et al., 2015)
AlNiCuPtPdAuAlNiCuPtPdAuCoFeAlNiCuMoCoFeArc-melting & dealloyingCO oxidation(Qiu et al., 2019b)
AlNiCoIrMe (Me: Mo, Cr, Cu, Nb, V)Arc-melting & dealloyingOxygen evolution reaction(Jin et al., 2019)
AlNiCoFeMe (Me: Mo, Cr, V, Nb)Arc-melting & dealloyingOxygen evolution reaction(Qiu et al., 2019a)
AlCuNiPtMe (Me: Mn, Co, Mo, Pd, Au)Arc-melting & dealloyingOxygen reduction reaction(Li et al., 2020b)
MnFeCoNiBall millingPowderOxygen evolution reaction(Dai et al., 2019)
CuAgAuPtPdArc-melting & cryogrindingFormic acid & methanol electro-oxidation(Katiyar et al., 2020)
AuAgPtPdCuArc-melting & cryogrindingCO2 & CO reduction reactions(Nellaiappan et al., 2020)
CoCrFeMnNiMechanical alloying (ball milling)Degrading azo dye solutions(Hu et al., 2019)
CoMoFeNiCuRuRhCoNiIrCarbothermal shockNanoparticleNH3 decomposition(Xie et al., 2019; Yao et al., 2020b)
PtPdRhRuCeCarbothermal shockNH3 oxidation(Yao et al., 2018)
RuIrCeNiWCuCrCoCarbothermal shockOxygen electrocatalysis(Lacey et al., 2019)
PtPdRhNiPtPdFeCoNiCarbothermal shockElectrocatalysts(Yao et al., 2020a)
NiCoCuFePtAerosol synthesis(Yang et al., 2020)
IrOsReRhRuThermal decompositionMethanol oxidation(Yusenko et al., 2017)
CrMnFeCoNiCombinatorial co-sputteringOxygen reduction reaction(Löffler et al., 2018)
PtAuPdRhPtAuPdRhRuUltrasonication-assisted wet chemistryHydrogen evolution reaction(Liu et al., 2019b)
MnFeCoNiCuMOF-template methodOxygen evolution reaction(Huang et al., 2020)
IrPdPtRhRuOne-pot polyol processHydrogen evolution reaction(Wu et al., 2020a)
CoCrCuNiAlSol-gel autocombustion synthesisSuperparamagnetic behavior(Niu et al., 2017)
CoCrFeMnNiLaser ablationOxygen evolution reaction(Waag et al., 2019)
PtIrPdRhRuAuto-catalytic formation(Broge et al., 2020)
CuPdSnPtAuNiCuPdSnAuPtMnCoNiCuRhPdSnIrPtAu, etc.Fast moving bed pyrolysisHydrogen evolution reaction(Gao et al., 2020)
PtNiFeCoCuOne-pot oil phase synthesis methodHydrogen evolution reaction & methanol oxidation reaction(Li et al., 2020a)
CuPdPtAuElectrochemical synthesisOxygen reduction reaction(Park and Ahn, 2020)
IrPdPtRhRuOne-pot oil phase synthesis methodEthanol oxidation reaction(Wu et al., 2020b)
RuFeCoNiCuOne-pot oil phase synthesis methodElectrocatalytic nitrogen reduction reaction(Zhang et al., 2020)
AuAgCuPdAuNiAgCrCoMo, etc.Oscillatory spark methodEthanol oxidation reaction, ethanol oxidation reaction, & formic acid oxidation reaction(Feng et al., 2020b)
PtFeCoNiCuAgRadio-frequency sputter depositionMethanol oxidation reaction(Tsai et al., 2009)
bccAlCoCrTiZnAlCoCrFeNiMechanical alloying (ball milling)PowderDegrading azo dye solutions(Hu et al., 2019)
CoCrCuNiAlSol-gel autocombustion synthesisNanoparticleSuperparamagnetic behavior(Niu et al., 2017)
hcpIrOsReRhRuThermal decompositionNanoparticleMethanol oxidation(Yusenko et al., 2017)
AmorphousCoFeLaNiPtCoCrMnNiVCoFeLaMnNiNanodroplet-mediated electrodepositionNanoparticleOxygen reduction reaction & hydrogen evolution reaction(Glasscott et al., 2019)
Typical structures of HEAs, HEOs, and other HEMs (A) fcc lattice structure, (B) bcc lattice structure, (C) hexagonal crystal structure, (D) rock-salt structure, (E) spinel structure, and (F) perovskite structure. Synthesis and structure of HEA NPs

Syntheses and structures of HEOs

Syntheses of HEOs

High-entropy compounds can be synthesized from binary metal compound precursors. For HEO in bulk form, traditional methods including solid-state sintering (Rost et al., 2015), spark plasma sintering (Biesuz et al., 2019; Gild et al., 2019), and reactive flash sintering (Liu et al., 2020; Wang et al., 2020a) have been successfully developed. The sintering temperature and soaking time are key factors for the syntheses of bulk HEOs. With proper solvent, heating/cooling rates, and reaction temperatures, shocking methods as CTS, FMBP, and pulsed laser deposition (PLD) might be adequate to fabricate HEO NPs from bulk HEOs similar with HEA. Films of (Mg0.2(1-x)CoxNi0.2(1-x)Cu0.2(1-x)Zn0.2(1-x))O (x = 0.2, 0.27, 0.33) were fabricated by solid-state sintering combined with PLD method (Meisenheimer et al., 2017). This method was also proved to be practicable in the fabrication of (MgNiCoCuZnSc)xO (x ∼ 0.167) and (LaPrNdSmEu)0.2NiO3 (Kotsonis et al., 2018; Patel et al., 2020). Wet-chemistry methods (Biesuz et al., 2018) such as solution combustion (Xiang et al., 2020), flame spray pyrolysis, nebulized spray pyrolysis (NSP), reverse co-precipitation (Sarkar et al., 2017), etc., have been successfully applied to the fabrication of (CoCuMgNiZn)Ox HEO NPs. Among these methods, the NSP method might have broad prospect because it can directly synthesize nanocrystalline single phase. Wang et al. have used Co0.2Cu0.2Mg0.2Ni0.2Zn0.2O as the precursor and LiF or NaCl as the reactant to prepare Li-containing entropy-stabilized oxyfluoride (HEOF) of Lix(Co0.2Cu0.2Mg0.2Ni0.2Zn0.2)OFx (Wang et al., 2019b).

Structure of HEOs

HEOs with fcc lattice constitute a large portion of energy-related applications such as catalysis and battery electrode materials. Among fcc structures, rock-salt HEO of (MgCoNiCuZn)Ox was synthesized by Rost et al. (Figure 3D) (Rost et al., 2015). Extended X-ray absorption fine structure and scanning transmission electron microscopy (STEM) analyses reveal that the first-near-neighbor cation-to-anion distances of this compound are identical (Rost et al., 2017). With a substitution of Li+ up to 16.6%, this HEO still keeps its rock-salt structure, shedding potential for energy-storage applications (Bérardan et al., 2016b). Recently, Jimenez-Segura et al. reported that this kind of rock-salt HEO exhibits a long-range magnetic order coexisted with the component disorder, and Zhang et al. found that the magnetic order in this HEO might correlate with Cu2+ Jahn-Teller distortions through neutron diffraction (Jimenez-Segura et al., 2019; Zhang et al., 2019a). Berardan et al. studied the distortion and atomic disordering of this HEO using X-band electron paramagnetic resonance spectroscopy. They proposed that post-annealing treatment can adjust the local environment of Cu2+, from octahedral in air-quenched Cu-sub-stoichiometric samples to rhombic in Cu-enriched HEOs, and then adjust the Jahn-Teller distortion and the dielectric property (Berardan et al., 2017). The results agreed with the density functional theory (DFT) calculations by Rák et al. Splitting of the Cu d-states near the Fermi level after relaxation suggests a Jahn-Teller distortion involving O-atom, which offers new possibilities of tuning the dielectric, magnetic, superconductivity, etc., properties of HEOs (Rák et al., 2018). HEOs composed of multicomponent rare earth elements as (Ce,La,Pr,Sm,Y)2O3+x synthesized by NSP method have CaF2-type structures. Interestingly, Djenadic et al. suggested that they cannot maintain single phase without Ce (such as (Ga, La, Nd, Pr, Sm, Y)O) (Djenadic et al., 2017). In addition, Chen et al. found that (HfZrCeTiSn)O2 with CaF2 structure exhibits a single phase at high temperature and multiphase at low temperature, which is believed as an evidence of entropy stabilization (Chen et al., 2018b). It seems that not only the entropy effect but also the compositions affect the structural stabilization of HEOs in addition to the synthesis procedure and the crystal size. Spinel HEOs (as shown in Figure 3E) such as (Cr,Fe,Mg,Mn,Ni)3O4 and (Co,Mg,Mn,Ni,Zn) Al,Co,Cr,Fe,Mn)2O4 with a space group of Fd3¯m can be synthesized by solid-state sintering (Dąbrowa et al., 2018; Fracchia et al., 2020b; Stygar et al., 2020). Based on NiFe2O4, (Co0.2Cr0.2Fe0.2Mn0.2Ni0.2)Fe1.9(Dy0.02Er0.02Gd0.02Ho0.02Tb0.02)O4 were synthesized by sol-gel method enhanced by entropy stabilization (Parida et al., 2020). (Cr0.2Fe0.2Mn0.2Ni0.2Zn0.2)3O4 nanocrystalline spinel powders have been synthesized by solution combustion with good magnetic properties (Mao et al., 2020). HEMs with perovskite structure (ABX3) have fruitful physical properties with a broad range of energy-related applications, which contains 12-fold coordinated A-site cations, 6-fold coordinated B-site cations, and X anion octahedra (Figure 3F). To form a single high-entropy perovskite phase, the “tolerance factor” of t has to be close to unity (t ≈ 1.0) (Patel et al., 2020).

Syntheses and structures of other HEMs

Methods similar to preparation of HEOs can be used to synthesize other high-entropy compounds. Sure et al. synthesized nanoscale HEC of (TiNbTaZrHf)Cx with fcc lattice via a facile electrochemical process (Sure et al., 2020). Liu et al. developed a facile borothermal reduction method to produce HEB of (Hf0.2Zr0.2Ta0.2Nb0.2Ti0.2)B2 in hexagonal crystal structure with average particle size of 310 nm (Liu et al., 2019a). Similar method as molten salt-mediated magnesiothermic reduction can prepare (Zr0.25Ta0.25Nb0.25Ti0.25)B2 HEB NPs with size of 28 – 56 nm (Ye et al., 2020). Table 1 summarizes the synthetic methods for the construction of efficient nano-HEMs with well-controlled composition, size, and uniformity, including both top-down and bottom-up methods. For the top-down methods, proper cooling/heating rate and appropriate duration time are required to be well controlled to prepare HEMs with adequate size. For the bottom-up methods based on wet-chemistry, low cost and high productivity are key factors for industrial production. Here we listed the reported structures of HEMs. More explorations will be emerging as according to our view there is no fundamental rule excluding HEMs from other structures.

Surface structure and activity

HEMs provide a large amount of surface sites with different atomic environments associated with various adsorption energies. For example, heterogeneous catalysts may speed up chemical reactions at active sites on the surface. HEM-based catalysts provide plentiful possible combinations of surface structure (then electronic structure) and the opportunity to achieve optimized catalytic performances. Batchelor et al. used theory to understand the activity of HEA catalysts and optimize the use of IrPdPtRhRu HEAs for catalyzing the oxygen reduction reaction (ORR). Figure 4A shows a schematic of ∗OH adsorbed on (111) plane of IrPdPtRhRu using nearest neighbor model. Through the DFT calculations of ∗OH adsorption energies on surface sites (Figure 4B), they predicted that near-continuous adsorption energies can be spanned out by tuning surface configurations (Batchelor et al., 2019; Pedersen et al., 2020). Then it is possible to engineer the surface structure to maximize the catalytic activity with peak value of adsorption energies (Figure 4C). The IrPdPtRhRu HEA and binary IrPt alloy identified in this study show significant activity enhancement over Pt(111) (Wang and Xin, 2019). Besides the surface structure, the variety of local electronic structure will intrinsically affect the other basic properties of HEMs, such as ionic/electronic conductivity, thermal conductivity, charge distribution magnetic properties, and dielectric properties. These properties may affect their catalytic properties through mass transfer, charge transfer, etc.
Figure 4

Predictive model of HEA for electrocatalysis

(A) Schematic of oxygen reduction reaction on the surface of IrPdPtRhRu HEA.

(B) Comparison of ∗OH adsorption energies machine learned on 871 samples (blue dots) and DFT calculation tested on 76 samples (red crosses).

(C) Distribution of adsorption energies for predicted Ir10.2Pd32.0Pt9.30Rh19.6Ru28.9 (upper) and Ir17.5Pt82.5 (lower), which show excellent activity for ORR. Reproduced with permission from Batchelor et al. (2019). Copyright © 2018, Elsevier.

Predictive model of HEA for electrocatalysis (A) Schematic of oxygen reduction reaction on the surface of IrPdPtRhRu HEA. (B) Comparison of ∗OH adsorption energies machine learned on 871 samples (blue dots) and DFT calculation tested on 76 samples (red crosses). (C) Distribution of adsorption energies for predicted Ir10.2Pd32.0Pt9.30Rh19.6Ru28.9 (upper) and Ir17.5Pt82.5 (lower), which show excellent activity for ORR. Reproduced with permission from Batchelor et al. (2019). Copyright © 2018, Elsevier. The structure of HEMs is believed to be stabilized by the high level of atomic disorder throughout the simple crystal structure and maximized entropy. Ratio of formation enthalpy to entropy and the differences in atomic radii are the core factors that affect their stability (Sarkar et al., 2020a). It is not clear if the surface structures of HEMs can also be preserved at working conditions.

Structural characterizations of HEMs

The structure of HEMs is featured with a single solid solution where the atoms disorderly distributed in lattice. The structural characterization demands the crystallography characterization and the chemical configuration including chemical composition and distribution. The former one can be realized with analytic methods such as X-ray diffraction (XRD), selected area electron diffraction (SAED), electron back-scattered diffraction, and high-resolution transmission electron microscopy (HRTEM) (Bondesgaard et al., 2019). The last of these is normally done with energy-dispersive spectrometry (EDS), wavelength-dispersive spectrometry, X-ray absorption spectroscopy (XAS), electron energy loss spectrometry, atomic probe tomography (APT), and inductively coupled plasma-optical emission spectrometry (Xie et al., 2019). In addition, scanning electron microscopy (SEM), atomic force microscopy, TEM, and high-angle annular dark-field (HAADF) STEM are useful for microstructure analysis (Liu et al., 2019b). HRTEM and high-resolution STEM have been used to clarify the atomic configuration uniformity (Gao et al., 2020). Hard X-ray can detect the electron structure and bond length of HEM (Sarkar et al., 2020b). The combination of the aforementioned methods may have advantages to reveal the local structure: for example, Ding et al. proposed that elements adopted greater aggregation with a wavelength of incipient concentration waves in CrFeCoNiPd bulk alloy, with the assistant of aberration-corrected TEM and EDS mapping (Ding et al., 2019).

Applications

Catalysis

Ammonia oxidation

Yao et al. synthesized PtPdRhRuCe, PtCoNiFeCuAu, PtPdCoNiCuAu, and PtPdCoNiFeCuAuSn nanocrystals with fcc lattice via the CTS method. They found that PtPdRhRuCe NPs exhibits excellent activity for ammonia oxidation with ∼ 100% conversion and > 99% nitrogen oxide selectivity toward NOx (NO + NO2) (Figure 5A) (Yao et al., 2018). The enhanced catalytic selectivity is suggested due to the highly atomically homogeneous nature of the solid-solution NPs, compared with the lower activity of the phase separated heterostructures. In addition, PtPdRhRuCe NPs show no activity degradation owing to their high-entropy nature preventing phase separation or elemental segregation (as shown in Figure 5B).
Figure 5

HEA NPs for ammonia oxidation and decomposition reactions

(A) Detailed steps of ammonia oxidation process over HEA NP catalysts highlighting how the apparent catalytic activity differences between the PtPdRhRuCe HEA NPs and phase-separated PtPdRhRuCe are due to different structure.

(B) Time-dependent catalytic performance with consideration of selectivity of PtPdRhRuCe HEA NPs at 700 °C. Reproduced with permission from Yao et al. (2018). Copyright © 2018, American Association for the Advancement of Science.

(C) Schematic of the rate-limiting factors in NH3 decomposition over HEA NP catalysts.

(D) NH3 decomposition efficiency in different systems. Reproduced with permission from Xie et al. (2019). Copyright © 2019, Springer Nature.

HEA NPs for ammonia oxidation and decomposition reactions (A) Detailed steps of ammonia oxidation process over HEA NP catalysts highlighting how the apparent catalytic activity differences between the PtPdRhRuCe HEA NPs and phase-separated PtPdRhRuCe are due to different structure. (B) Time-dependent catalytic performance with consideration of selectivity of PtPdRhRuCe HEA NPs at 700 °C. Reproduced with permission from Yao et al. (2018). Copyright © 2018, American Association for the Advancement of Science. (C) Schematic of the rate-limiting factors in NH3 decomposition over HEA NP catalysts. (D) NH3 decomposition efficiency in different systems. Reproduced with permission from Xie et al. (2019). Copyright © 2019, Springer Nature.

Ammonia decomposition

Co-Mo nanoalloys show high activity for the reaction of ammonia decomposition. Xie et al. found that quinary CoMoFeNiCu NPs (CoxMoyFe10Ni10Cu10 (x + y = 70)) on carbon nanofiber exhibit an even higher NH3 conversion rate than either binary Co-Mo or Ru catalysts (Xie et al., 2019). The HEA-CoxMoy catalysts are stable and highly efficient for NH3 decomposition, achieving a more than 20 times improvement versus Ru catalyst. It is well known that in a Co-rich system, the rate is limited by NH3 activation, whereas in a Mo-rich system, the rate is limited by the ∗N desorption. Importantly, the catalytic activity of these HEA NPs is found regulated by the Co/Mo ratio (Figure 5C). Thus, their catalytic activities and catalytic kinetics exhibited a volcano-type profile described by the nitrogen adsorption energy (Figure 5D) (Xie et al., 2019). Guided by the DFT calculations, another quinary system, RuRhCoNiIr exhibits 5 times higher efficiency for NH3 decomposition than the separated phases at 550 K. Similarly, its outstanding catalytic performance and high stability were attributed to the synergistic effect of the ultrafine size, the uniform dispersion, the multi-elemental composition, and the fcc lattice structure (Yao et al., 2020b).

Oxidation of aromatic alcohols

Feng et al. synthesized a holey lamellar HEO of Co0.2Ni0.2Cu0.2Mg0.2Zn0.2O with mesoporous structure by anchoring and merging. This compound has ultra-high catalytic activity for the solvent-free aerobic oxidation of benzyl alcohol, with 98% conversion achieved in only 2 h. It has been considered that the highly catalytic activity could be enhanced by abundant oxygen vacancies and holey lamellar framework in the HEO (Feng et al., 2020a). Understanding this example, it is clear that HEO catalysts show promise in a variety of catalytic reactions including many examples of thermocatalysis.

Electrocatalysis

Electrochemical devices for energy-related applications were normally utilized at relatively low temperature of ∼ 80 °C and liquid environment. HEM catalysts have been developed for the following hydrogen evolution reaction (HER), oxygen evolution reaction (OER), ORR, CO/CO2 reduction reaction (CORR/CO2RR), CO oxidation reaction (COOR), methanol/ethanol oxidation reaction (MOR/EOR), etc. We introduced the representative works in following text.

Hydrogen evolution reaction

Normally HEAs show better corrosion resistance than normal alloys due to their locally disordered chemistry. Therefore, HEAs were considered as highly active and stable electrocatalysts for HER in both acidic and basic electrolytes. It is reported that the Ni20Fe20Mo10Co35Cr15 catalyst has an overpotential of 107 mV in acidic solutions and 172 mV in the basic solution, respectively, showing much better activity than the dual-phase catalyst. This activity is even comparable to that of Pt sheet. The high HER activity of HEA is believed to be enhanced by their partially filled d-orbitals fitting for the gain and loss of electrons, promoting the adsorption of hydrogen (Zhang et al., 2018a). As shown in Figures 6A and 6B, the CoFeLaNiPt HEMGs can be electro synthesized with equimolar components by confining multiple metal salt precursors to water nanodroplets (Glasscott et al., 2019). This system shows a high HER overpotential of 555 ± 2 mV and a high OER overpotential of 377 ± 4 mV, respectively, benefitted from synergistic activity from Pt sites and transition metal sites. Wu et al. also successfully synthesized a highly active RuRhPdIrPt HER catalyst. Hard X-ray photoelectron spectroscopy (XPS) analysis reveals that this kind of HEA has a broad and featureless valence band spectrum. The authors believe that the featureless XPS spectra imply random atomic configurations with a variety of local electronic structures and unique local density states in HEA, and different atomic configurations have distinct local electronic structures and spectra characteristics in general phases (Wu et al., 2020a).
Figure 6

HEM electrocatalysts for HER and OER

(A) Elemental mappings and electron diffraction pattern show structural uniformity of CoFeNiLaPt HEMP NP.

(B) Electrocatalytic evaluation of a CoFeLaNiPt electrocatalyst. Reproduced with permission from Glasscott et al. (2019). Copyright © 2019, Springer Nature.

(C) Schematic of two-step thermochemical water splitting using a (FeMnCoNi)xO HEO. Reproduced with permission from Zhai et al. (2018). Copyright © 2018, The Royal Society of Chemistry.

(D) HAADF image and structure of nanoporous AlNiCoFeMo HEA covered with (AlNiCoFeMo)xO HEO surface.

(E) Linear scanning voltammetry curves showing the oxygen evolution activity on the different nanoporous HEAs and nanoporous AlNiCoFe. Reproduced with permission from Qiu et al. (2019a). Copyright © 2019, American Chemical Society.

HEM electrocatalysts for HER and OER (A) Elemental mappings and electron diffraction pattern show structural uniformity of CoFeNiLaPt HEMP NP. (B) Electrocatalytic evaluation of a CoFeLaNiPt electrocatalyst. Reproduced with permission from Glasscott et al. (2019). Copyright © 2019, Springer Nature. (C) Schematic of two-step thermochemical water splitting using a (FeMnCoNi)xO HEO. Reproduced with permission from Zhai et al. (2018). Copyright © 2018, The Royal Society of Chemistry. (D) HAADF image and structure of nanoporous AlNiCoFeMo HEA covered with (AlNiCoFeMo)xO HEO surface. (E) Linear scanning voltammetry curves showing the oxygen evolution activity on the different nanoporous HEAs and nanoporous AlNiCoFe. Reproduced with permission from Qiu et al. (2019a). Copyright © 2019, American Chemical Society. In addition to HEA catalysts, Zhai et al. discovered HEO mixture of (FeMgCoNi)Ox (x ∼ 1.2) composed of rock-salt structure and spinel structure can thermochemically split water and produce H2 and O2 in a two-step cycle. During the reaction, the HEO mixture experiences a reversible phase swing between rock-salt and spinel phases (Figure 6C). The rock-salt oxide became thermally reduced at a higher temperature and released O2, and then the reduced oxide cooled to a lower temperature and was oxidized by water to produce H2 during thermochemical cycling. Interestingly, XAS analysis shows Fe component, the only redox-active species in the HEO. However, Fe ions in well-studied spinel ferrites are purely Fe3+, which is believed to have a larger redox capacity than a spinel/rock-salt two-phase oxide system where both Fe2+ and Fe3+ are present (Zhai et al., 2018). Furthermore, Zhao et al. prepared (CrMnFeCoNi)xP HEMP nanosheets with hexagonal crystal structure, which presents a current density of 10 mA/cm2 at an overpotential of 136 mV for HER. The enhanced electrocatalytic activity might be benefited from the synergistic activity of the components (Zhao et al., 2020b).

Oxygen evolution reaction

Developed from NiFe or NiCoFe nano-alloys, nanoporous AlNiCoFeX (X = Mo, Nb, Cr) HEAs covered with a naturally oxidized surface exhibit low impedance resistance and excellent electrochemical durability as OER catalysts (Figure 6D). These catalysts were synthesized with chemical etching, a top-down strategy that is typically considered to be highly reproducible and applicable to the syntheses of other HEA catalysts. Cr, Nb, V, Zr, Mn, Cu, Mo, etc., have been added and formed nanoporous AlNiCoFeX HEAs, enhancing the OER activity. Among them, enhancement from Cr, Nb, and Mo is strong and that from V, Zr, and Cu is weak, as shown by the linear scanning voltammetry curves in Figure 6E (Qiu et al., 2019a). Jin et al. synthesized a quinary nanoporous AlNiCoIrMo HEA with 2 nm well-aligned nanopores via dealloying. This system exhibits good electrocatalytic activity (overpotential of 233 mV at 10 mA/cm2) and cycling stability than those Ir-based catalysts in the literature (highest is the result published by Jiang et al. on Ir nanosheets with overpotential of 240 mV at 10 mA/cm2) (Jiang et al., 2018) that are enhanced by the high-entropy effect and sluggish diffusion effect of HEA (Jin et al., 2019). It was proposed by Wang et al. that the high-entropy perovskite fluorides of K(MgMnFeCoNi)F3, etc., can act as efficient OER catalysts in alkaline media. The substitution of K with Na can further improve the OER performance. For example, K0.8Na0.2(MgMnFeCoNi)F3 exhibits a low overpotential of 314 mV at a current density of 10 mA/cm2, which is much lower than that of either IrO2 or Ba0.5Sr0.5Co0.8Fe0.2O3−d (Wang et al., 2020b). The authors suggest that the superior performance might be benefited from its highly dispersed active sites, low charge transfer resistance, and fast mass transfer enhanced by high entropy.

Oxidation reduction reaction

Pt-based HEA nanostructures have been explored for electrocatalysts for ORR. Li et al. have synthesized nanoporous AlCuNiPtX (X = Pd, V, Co, Mo, etc.) by a selective dealloying method, which provides tunable ORR catalytic performance. Among these quinary alloys, the nanoporous AlCuNiPtMn exhibited high ORR activity and good electrochemical cycling durability, with a half-wave potential of 0.945 V in acidic media which is a ∼ 16-fold mass activity of Pt/C. The authors believe that the Mn component enables modulation of the electronic properties and enhances the ORR activity (Li et al., 2020b). Löffler et al. investigated the variation of adsorption energy and ORR activity when adding or replacing components into CrMnFeCoNi NPs (as shown in Figure 7A). Upon adding Nb or Mo, the first two adsorption peaks split and the intrinsic ORR activity improved (Figure 7B). However, adding Cu or replacing Ni by V will reduce the ORR activity. The enhancement of ORR activity might be attributed to the continuous coverage of adsorption energies and the presence of highly active sites at the surface (Figure 7C) (Löffler et al., 2020). The authors also suggested that the high activity not only correlated with synergistic combination of right choice of elements but also was affected by their atomic ratios (Löffler et al., 2018). In another case, Qiu et al. put immiscible elements together into a single-phase nanoporous structure to form senary AlNiCuPtPdAu, octonary AlNiCuPtPdAuCoFe, and senary AlNiCuMoCoFe by a dealloying strategy. The nanoporous AlNiCuPtPdAu exhibited greatly enhanced high-temperature stability (up to 600 °C) and good ORR activity of 10 times mass activity than that of Pt/C (Qiu et al., 2019b). So far, the strategy of Pt-based HEAs for ORR mainly lies in their tunable lattice parameters and electronic structures, which help to optimize the absorption energy. However, more work will be needed to find out the best recipe of elements and more importantly, the rules of optimizing the activity and stability.
Figure 7

Nanosized HEM electrocatalysts for ORR

(A) Changes of adsorption energy distribution patterns when adding or replacing elements from CrMnFeCoNi.

(B) Schematic illustration of senary HEA nanoparticle with near-continuous adsorption energy distribution.

(C) Comparison of intrinsic activity of different catalysts toward ORR highlighting importance of HEM structure on activity. Reproduced with permission from Löffler et al. (2020). Copyright © 2020, Wiley-VCH.

Nanosized HEM electrocatalysts for ORR (A) Changes of adsorption energy distribution patterns when adding or replacing elements from CrMnFeCoNi. (B) Schematic illustration of senary HEA nanoparticle with near-continuous adsorption energy distribution. (C) Comparison of intrinsic activity of different catalysts toward ORR highlighting importance of HEM structure on activity. Reproduced with permission from Löffler et al. (2020). Copyright © 2020, Wiley-VCH.

CO/CO2 reduction reaction

Recently, Pedersen et al. proposed an approach combining DFT calculations and supervised machine learning to discover selective and active CORR/CO2RR catalysts (Figure 8A). This work calculated the CO and H adsorption energies on the (111) surfaces of CoCuGaZnNi and AgAuCuPdPt HEAs. Based on the calculation, they proposed that HEAs with weak hydrogen adsorption could suppress the molecular hydrogen formed, and those with strong CO adsorption on the (111) surfaces could favor the reduction of CO (Pedersen et al., 2020). Nellaiappan et al. found that AuAgPtPdCu exhibited a faradaic efficiency of about 100% toward gaseous products at a low applied potential (−0.3 V versus reversible hydrogen electrode), adopting high CO2RR activity. After comparing with pristine Cu(111) surface, the authors believe that the high activity might be caused by the reversal in adsorption trends for ∗OCH3 and ∗O on Cu(111) and HEA surfaces (Figure 8B). Although many elements are present in this HEA catalysts, the electrocatalytic activity might be predominantly enhanced by the redox-active Cu metal (Cu2+/Cu0), and other metals provide synergistic effect (Nellaiappan et al., 2020). Chen et al. used a low-temperature mechanochemical method-synthesized (NiMgCuZnCo)O catalysts decorated with noble metal Pt/Ru single atoms or nanoclusters. This mixed HEO shows a good stability at high temperature and a high catalytic activity in the hydrogenation of atmospheric CO2 to CO (44.9% at 500 °C) (Figure 8C). The HEOs are supposed to be thermodynamically unstable at low temperature; however, the high entropy here lowers Gibbs free energy and helps the homogeneous distribution of cation species and improves the stability (Chen et al., 2019a).
Figure 8

HEM electrocatalysts for CO2RR/CORR

(A) Plots of the CO2RR/CORR selectivity and CORR activity space achievable over AgAuCuPdPt NP catalysts. Reproduced with permission from Pedersen et al. (2020). Copyright © 2020, American Chemical Society.

(B) Free-energy diagram of CO2RR on the AuAgPtPdCu HEA surface. Reproduced with permission from Nellaiappan et al. (2020). Copyright © 2020, American Chemical Society.

(C) Schematic illustration of Pt-(NiMgCuZnCo)O and CO2 hydrogenation activity. Reproduced with permission from Chen et al. (2019a). Copyright © 2019, American Chemical Society.

HEM electrocatalysts for CO2RR/CORR (A) Plots of the CO2RR/CORR selectivity and CORR activity space achievable over AgAuCuPdPt NP catalysts. Reproduced with permission from Pedersen et al. (2020). Copyright © 2020, American Chemical Society. (B) Free-energy diagram of CO2RR on the AuAgPtPdCu HEA surface. Reproduced with permission from Nellaiappan et al. (2020). Copyright © 2020, American Chemical Society. (C) Schematic illustration of Pt-(NiMgCuZnCo)O and CO2 hydrogenation activity. Reproduced with permission from Chen et al. (2019a). Copyright © 2019, American Chemical Society.

CO oxidation reaction

Chen et al. discovered that NiMgCuZnCoOx can be used as COOR catalyst, which exhibits a complete CO to CO2 conversion at 305 °C. To obtain considerable activity for the CO conversion at lower temperature, Pt was introduced into the catalyst to form rock-salt PtNiMgCuZnCoOx. After being reduced in H2, this catalyst shows a 100% conversion of CO even at 155 °C. The high activity is attributed to the highly dispersed active sites and entropy stabilization (Chen et al., 2018a). On the other hand, the authors prepared a heterostructure of (NiMgCuZnCo)O-CuCeOx via entropy-driven mechanochemistry. CuCeOx NPs were decorated on the surface of HEOs, providing catalytic sites for COOR (Chen et al., 2020a). It has also been reported that Cu2+ in Mg0.2Co0.2Ni0.2Cu0.2Zn0.2O can be rapidly reduced to Cu+ by CO at a temperature higher than 130 °C. Operando soft X-ray absorption spectroscopy suggests that Cu is the active metal, whereas Co and Ni do not play any role in the COOR. Moreover, Cu2+ oxidation state can be easily but slowly recovered by annealing under O2, which allows the catalytic cycle to continue (Fracchia et al., 2020a). Okejiri et al. used a sonochemical-based method to synthesize monodispersed perovskite HEO NPs with an average crystallite size ∼ 5.9 nm. It is shown that Ru/BaSrBi(ZrHfTiFe)O3 with excellent dispersion of Ru adopts good catalytic activity for COOR. For example, 100% conversion can be achieved at 118°C (Okejiri et al., 2020).

Methanol/ethanol oxidation reaction

HEAs composed of Pt-group metals may form not only fcc lattice but also hcp lattice, varying with the component and concentration. Yusenko et al. synthesized hcp structured Ir0.19Os0.22Re0.21Rh0.20Ru0.19 and fcc structured Ir0.26Os0.05Pt0.31Rh0.23Ru0.15 HEAs via the thermal decomposition of single-source precursor. The former shows a high thermal stability at 1500 K as well as exceptional MOR activity, which was suggested to be due to both electronic and cocktail effects (Yusenko et al., 2017).

Energy storage

Recently, HEOs have been explored for the applications of electrode materials for Li/Na ion batteries, solid electrolytes, and electrode materials for Li-sulfur battery. These applications demand the materials for different properties, such as high capacity, either high electronic conductivity, or high ionic conductivity, etc. The cation tunability of HEOs may fulfill these requests.

Electrode materials in Li ion battery

As shown in Figure 9A, HEOs have been used as both cathode and anode materials for lithium-ion batteries (LIBs) (Oses et al., 2020). Then rational design of HEOs by combining both active and inactive cations is the common strategy for exploring high-performance electrode materials. Most of the investigated HEO electrodes have a rock-salt structure. Lun et al. design high-entropy LIB cathodes with rock-salt structure guided by the DFT calculations. 30% Li excess has been added to ensure good Li transport without severely compromising the capacity. By comparing a group of rock-salt cathodes containing two, four, and six transition metal species (TM6: Li1.3Mn2+0.1Co2+0.1Mn3+0.1Cr3+0.1Ti0.1Nb0.2O1.7F0.3), they found both energy density and rate capability increased with the number of metal species. TM6 exhibits an excellent rate capability, achieving ∼307 mAh/g at a rate of 20 mA/g and ∼170 mAh/g at a rate of 2000 mA/g, respectively, as shown in Figure 9B (Lun et al., 2020). Wang et al. synthesized Li-containing HEOF of Li(Co0.2Cu0.2Mg0.2Ni0.2Zn0.2)OF via a facile mechanochemistry method. The HEOF exhibits a potential of 3.4 V versus Li+/Li, enabling its application as cathode active material due to entropy stabilization (Wang et al., 2019b).
Figure 9

HEM electrodes for lithium storage

(A) Schematic illustration of a rechargeable battery with HEOs as the electrode. Reproduced with permission Oses et al. (2020) and Wang et al. (2019a). Copyright © 2019, Elsevier. Copyright © 2020, Springer Nature.

(B) The first cycle charge-discharge curves of rock-salt TM6 (Li1.3Mn2+0.1Co2+0.1Mn3+0.1Cr3+0.1Ti0.1Nb0.2O1.7F0.3) at different currents. Reproduced with permission from Lun et al. (2020). Copyright © 2020, Springer Nature.

(C) Operando XRD studies in the first cycle and (C′) SAED pattern from the lithiated sample of (Co0.2Cu0.2Mg0.2Ni0.2Zn0.2)O HEO anode.

(D) Schematics of the proposed de-/lithiation mechanism during the conversion reaction of HEO. M stands only for the transition metal cations. Reproduced with permission from Sarkar et al. (2018). Copyright © 2018, Springer Nature.

HEM electrodes for lithium storage (A) Schematic illustration of a rechargeable battery with HEOs as the electrode. Reproduced with permission Oses et al. (2020) and Wang et al. (2019a). Copyright © 2019, Elsevier. Copyright © 2020, Springer Nature. (B) The first cycle charge-discharge curves of rock-salt TM6 (Li1.3Mn2+0.1Co2+0.1Mn3+0.1Cr3+0.1Ti0.1Nb0.2O1.7F0.3) at different currents. Reproduced with permission from Lun et al. (2020). Copyright © 2020, Springer Nature. (C) Operando XRD studies in the first cycle and (C′) SAED pattern from the lithiated sample of (Co0.2Cu0.2Mg0.2Ni0.2Zn0.2)O HEO anode. (D) Schematics of the proposed de-/lithiation mechanism during the conversion reaction of HEO. M stands only for the transition metal cations. Reproduced with permission from Sarkar et al. (2018). Copyright © 2018, Springer Nature. HEO with transition metal cations can serve as anode for LIBs. Sarkar, Chen, Qiu, and Wang et al. studied the rock-salt (Co0.2Cu0.2Mg0.2Ni0.2Zn0.2)O for anode of LIBs, respectively (Sarkar et al., 2018; Qiu et al., 2019; Wang et al., 2019a; Chen et al., 2019b). Operando XRD and ex situ TEM analysis suggest the lithiation reaction of this HEO might be MOx + 2xLi → M + xLi2O, with a capacity of 700–1,000 mAh/g (Figure 9C). The authors suggest that Co2+ and Cu2+ cations provide capacity with a conversion reaction, whereas the other cations act as matrix to maintain the rock-salt structure. In addition, it is suggested that Mg2+ ions are vital in making the HEO structure intact and the rock-salt structure can be preserved after fully lithiation (Figures 9D, [Sarkar et al., 2018]). It has also been proposed by Qiu et al. that a small amount of inactive MgO would be formed after the discharge of (Co0.2Cu0.2Mg0.2Ni0.2Zn0.2)O HEO, which may help buffer the volume change of the HEO anode and keep the separation of active nanosized crystals (Qiu et al., 2019). Wang et al. has built a full cell with an anode of (Co0.2Cu0.2Mg0.2Ni0.2Zn0.2)O and a cathode of LiNi1/3Co1/3Mn1/3O2. The full batteries show pretty good performances, delivering an initial capacity of 446 mAh/g and a capacity of 256 mAh/g after 100 cycles (Wang et al., 2019a). Lökçü et al. synthesized rock-salt (MgCoNiZn)1−xLixO (x = 0.05, 0.15, 0.25, 0.35) anode via a conventional solid-state sintering. XPS analysis revealed that the oxidation of Co and the concentration of oxygen vacancies enhanced by the Li content. The (MgCoNiZn)0.65Li0.35O anode exhibited an initial capacity of 1,930 mAh/g and a discharge capacity of 610 mAh/g after 130 cycles (Lökçü et al., 2020). HEOs with perovskite and spinel structures have also been developed. For example, (Bi,Na)1/5(La,Li)1/5(Ce,K)1/5Ca1/5Sr1/5]TiO3 with perovskite structure shows an initial discharge capacity of 125.9 mAh/g and good stability (Yan et al., 2020). A spinel anode of (Mg0.2Ti0.2Zn0.2Cu0.2Fe0.2)3O4 shows a large reversible capacity of 504 mAh/g at a current density of 100 mA/g after 300 cycles (Chen et al., 2020b).

Electrode materials in Na ion battery

Besides the application in LIBs, it is not surprising to see HEOs being used as electrodes of Na-ion batteries. Zhao et al. have synthesized O3-type NaNi0.12Cu0.12Mg0.12Fe0.15Co0.15Mn0.1Ti0.1Sn0.1Sb0.04O2 via solid-state reaction as cathode (Figure 10A). It is believed that different cations have different roles. For example, cations of Ni2+, Cu2+, Fe3+, and Co3+ are expected to provide capacity; Mg2+ and Ti4+ may stabilize the host's structure during cycling; Mn2+ can help in the formation of the structure; Sn4+ and Sb5+ may raise the average voltage (Sathiya et al., 2018; Yuan et al., 2014). The HEO exhibits a long cycling stability (83% of capacity retention after 500 cycles) as well as a good rate capability, enhanced by the highly reversible O3-P3 phase transition (as shown in Figures 10B and 10C). The authors propose that the multiple transition-metal components may accommodate the changes of local bonds during (de)intercalation, resulting in varied interactions between TMO2 and NaO2 slabs, as shown in Figure 10D. During the cycle processes, the size and oxidation states of transition metals that participate in the charge compensation can change. The high entropy composition promotes homogeneous distribution of these redox elements and accommodates the local changes. Consequently, it may help to maintain its layered O3-type structure, supporting the long-term cycling stability as well as better rate performance (Zhao et al., 2020a).
Figure 10

Sodium storage of HEMs

(A) XRD pattern shows the as-prepared NaNi0.12Cu0.12Mg0.12Fe0.15Co0.15Mn0.1Ti0.1Sn0.1Sb0.04O2 HEO with O3-type crystal structure.

(B) Charge-discharge curves cycled at rates from 0.1 to 5.0 C in the voltage range 2.0–3.9 V.

(C) Capacity retention and coulombic efficiency of the O3-type HEO at 3.0 C up to 500 cycles.

(D) Schematic of possible Na+ (de)intercalation process of the O3-type HEO cathode (left) with multiple transition metals (right). Reproduced with permission from Zhao et al., 2020a. Copyright © 2020, Wiley-VCH.

Sodium storage of HEMs (A) XRD pattern shows the as-prepared NaNi0.12Cu0.12Mg0.12Fe0.15Co0.15Mn0.1Ti0.1Sn0.1Sb0.04O2 HEO with O3-type crystal structure. (B) Charge-discharge curves cycled at rates from 0.1 to 5.0 C in the voltage range 2.0–3.9 V. (C) Capacity retention and coulombic efficiency of the O3-type HEO at 3.0 C up to 500 cycles. (D) Schematic of possible Na+ (de)intercalation process of the O3-type HEO cathode (left) with multiple transition metals (right). Reproduced with permission from Zhao et al., 2020a. Copyright © 2020, Wiley-VCH.

All-solid-state electrolytes

Bérardan et al. synthesized a solid-state HEO electrolyte of (MgCoNiCu)1-x-yGayAxO (A = Li, Na, K) via solid-state reaction. Impedance spectroscopy measurements reveal that this rock-salt structure has a superionic Li+ mobility of ∼ 10−3 S/cm at room temperature as well as a fast ionic Na+ mobility of 5 × 10−6 S/cm. The authors believe that high ionic conductivity occurring via the oxygen vacancies is created by a charge compensation (Bérardan et al., 2016b). So far, there are not many reports on the application of HEOs on solid electrolyte. To make this kind of electrolyte useful, the chemical stabilities and electrochemical compatibility between electrolyte and electrodes are important for further investigations.

Chemical anchor/sulfur host material/cathode

(NiMgCuZnCo)Ox HEO has also been utilized as the cathode catalyst to anchor the lithium polysulfides (LIPSs) in Li-S batteries. The homogeneously dispersive multiple metal active species in HEO and the synergistic contribution of Li-O and S-Ni bonds favor the restriction of LIPSs and facilitate the redox reaction in the cathode. The authors suggest that this HEO can act as a promoter for polysulfide immobilizing in Li-S batteries with a competitive reversible capacity and a low capacity decay (Zheng et al., 2019).

Hydrogen storage

Shen et al. fabricate an equiatomic HEA of TiZrHfMoNb with the arc-melting method for hydrogen storage. They found this HEA changes its structure from bcc to fcc during hydrogen absorption. The metastable fcc lattice can transform back to bcc during hydrogen desorption at around 302°C. The reversible transformation of this alloy may lead to an improved hydrogen recycling rate and an enhanced energy efficiency for application of solar thermal energy storage (Shen et al., 2019).

Supercapacitor

Kong et al. synthesize nanoporous natured HEAs by dissolving the selective HEA alloys into sulfuric acid solution. Utilized as a binder-free electrode for supercapacitor, nanoporous HEA of AlCoCrFeNi exhibits a high volumetric capacitance of 700 F/cm3 and an excellent cycling stability of over 3,000 cycles. The good performance may be due to the large surface area from its interconnected large pore channels (Kong et al., 2019). In addition, other HEMs with high surface area were synthesized for the application of supercapacitor. For example, Jin et al. synthesize (VCrNbMoZr)Nx HEN via a mechanochemical-assisted soft urea method. They found that the HEN exhibits a specific capacitance of 78 F/g at a scan rate of 100 mV/s in 1 m KOH solution (Jin et al., 2018). Xu et al. synthesize uniform quinary FeNiCoMnMg (∼ 30 nm) and FeNiCoMnCu (∼ 50 nm) NPs on superaligned electrospun carbon nanofibers via CTS. The latter one exhibits a high capacitance of 203 F/g and a specific energy density of 21.7 Wh/kg (Xu et al., 2020).

Other energy-related applications

Dielectric materials for energy applications

The dielectric compound may store the electric energy via its polarization in the presence of external electric field and therefore be used for various applications such as capacitors, energy harvesting and storage devices, and high-power electronic transducers. Bérardan et al. suggest that (MgCoNiCuZn)O can be substituted by aliovalent ions (+1 or +3) with a charge compensation mechanism, to obtain HEOs as (MgCoNiCu)0.8(LiGa)0.2O, (MgNiCoCuZn)0.95Li0.05O, and (MgNiCoCuZn)0.84Li0.16O. These HEOs exhibit promising colossal dielectric properties in a wide frequency range (2.3 MHz –100 Hz) (Bérardan et al., 2016a). Moreover, high-entropy perovskite (Na0.2Bi0.2Ba0.2Sr0.2Ca0.2)TiO3 is reported to produce a discharge energy density of 1.02 J/cm3 under an electric field of 145 kV/cm, showing promising applications as solid-state refrigeration and energy storage in circuit (Pu et al., 2019).

Magnetic devices

The magnetic properties of functional oxides are dependent on their structure, both atomically and electronically. Meisenheimer et al. take advantage of the inherent chemical and structural disorder of HEO to boost the exchange bias in Permalloy/(Mg0.25(1-x)CoxNi0.25(1-x)Cu0.25(1-x)Zn0.25(1-x))O heterostructures. The heterostructures are integration of ferromagnetic and antiferromagnetic compounds, which can achieve a 10-fold increase in the exchange interaction depending on the concentration of Co. The authors propose that both local disorder of cations and the concentration of magnetic ions contribute to an improved exchange coupling (Meisenheimer et al., 2017). HEMs containing magnetic ions may serve as a tunable platform for magnetic properties. With proper integration with other functional compounds, HEM-based devices can be explored for applications of magnetic memories, magnetometers, and devices for magneto-optics.

Perspectives

The past few years have witnessed the emergence of HEMs for energy-related applications, especially catalysis and energy storage. HEMs provide a rich platform for tuning the antisite disordering of atoms within a crystalline structure, for achieving an optimized performance. However, the understanding between the performance and surface structure at atomic level is still not clear. The essential task for developing HEMs is to achieve this structure-property relationship and precisely control their structure and local electronic structure. We briefly list several challenges below. With the help of advanced calculations, syntheses, and characterization tools, these challenges also provide huge opportunities to explore HEMs. Strategies of designing proper HEMs. Theoretical investigations based on state-of-the-art DFT, molecular dynamics, and others may provide useful guidelines for designing proper HEMs. However, the atomic type disordering may give a huge difficulty in building the atomic model and discovering special constituent combinations for special energy-related applications. Lots of work needs to be done for this goal, and machine learning and other big data approaches may help to deal with the enormous dataset. Precise structural control of HEMs. Experimentally, synthetic method with controllable manners, simpler operation, and larger scale are still much needed. Strategies of well size, shape, and composition control of HEMs, which are analogous to the case of well-defined metal (including alloy and intermetallic) and semiconductor nanomaterials, remain to be extended. Advanced characterization. Various advanced techniques, such as X-ray-based technique, SEM, and TEM have been used for structural characterization, which may provide atomically accurate structural and chemical information of HEMs. Developing a technique to address the atomic structure of HEMs quantitatively in 3D would be a challenging task. 3D electron tomography and APT can be promising for this goal. Thermal stability of HEMs. Although the term “high-entropy” is prevalent in the literature and the idea of using configurational entropy to stabilize a single-phase solid solution is well accepted, the actual contribution of entropy (including configurational, vibrational, electronic, and magnetic entropy) and enthalpy to the total Gibbs free energy are difficult to determine. Attentions are also needed to be paid to the surface strain. Surface stability at work environment. The surface structure and stability of HEMs during reaction are essential for catalytic applications. This issue is important for catalysis research, but the surface structure of HEMs can be more complicated than the current used catalysts. Exploration of new types of HEMs. The emergence of HEMs and their energy-related application are still in their infancy. The idea of configuring more than five elements into one lattice may provide infinite combinations. It is also important to investigate how the antisite disordering interplays with other physical parameters, such as strain, 2-dimentional geometry, magnetism, etc. We believe that HEMs will work as important materials platform for future energy and other functional devices.

Resource availability

Lead contact

Further information and requests for resources should be directly to and will be fulfilled by the Lead Contact, Dong Su (dongsu@iphy.ac.cn).

Materials availability

This study did not generate any new unique reagents.

Data and code availability

The data that support the findings of this study are available from the corresponding authors on reasonable request.

Methods

All methods can be found in the accompanying Transparent Methods supplemental file.
  41 in total

1.  High-Entropy Perovskite Fluorides: A New Platform for Oxygen Evolution Catalysis.

Authors:  Tao Wang; Hao Chen; Zhenzhen Yang; Jiyuan Liang; Sheng Dai
Journal:  J Am Chem Soc       Date:  2020-02-27       Impact factor: 15.419

2.  Order emerging from disorder.

Authors:  Nita Dragoe; David Bérardan
Journal:  Science       Date:  2019-11-01       Impact factor: 47.728

3.  Tuning element distribution, structure and properties by composition in high-entropy alloys.

Authors:  Qingqing Ding; Yin Zhang; Xiao Chen; Xiaoqian Fu; Dengke Chen; Sijing Chen; Lin Gu; Fei Wei; Hongbin Bei; Yanfei Gao; Minru Wen; Jixue Li; Ze Zhang; Ting Zhu; Robert O Ritchie; Qian Yu
Journal:  Nature       Date:  2019-10-09       Impact factor: 49.962

4.  Eutectic Synthesis of High-Entropy Metal Phosphides for Electrocatalytic Water Splitting.

Authors:  Xinhui Zhao; Zhimin Xue; Wenjun Chen; Yaqing Wang; Tiancheng Mu
Journal:  ChemSusChem       Date:  2020-03-03       Impact factor: 8.928

5.  Across the Board: Sheng Dai on Catalyst Design by Entropic Factors.

Authors:  Sheng Dai
Journal:  ChemSusChem       Date:  2020-02-27       Impact factor: 8.928

6.  Stable Multimetallic Nanoparticles for Oxygen Electrocatalysis.

Authors:  Steven D Lacey; Qi Dong; Zhennan Huang; Jingru Luo; Hua Xie; Zhiwei Lin; Dylan J Kirsch; Vivek Vattipalli; Christopher Povinelli; Wei Fan; Reza Shahbazian-Yassar; Dunwei Wang; Liangbing Hu
Journal:  Nano Lett       Date:  2019-07-17       Impact factor: 11.189

7.  Data-Driven Design of Ecofriendly Thermoelectric High-Entropy Sulfides.

Authors:  Rui-Zhi Zhang; Francesco Gucci; Hongyu Zhu; Kan Chen; Michael J Reece
Journal:  Inorg Chem       Date:  2018-09-26       Impact factor: 5.165

8.  Highly efficient decomposition of ammonia using high-entropy alloy catalysts.

Authors:  Pengfei Xie; Yonggang Yao; Zhennan Huang; Zhenyu Liu; Junlei Zhang; Tangyuan Li; Guofeng Wang; Reza Shahbazian-Yassar; Liangbing Hu; Chao Wang
Journal:  Nat Commun       Date:  2019-09-05       Impact factor: 14.919

9.  Entropy-stabilized oxides.

Authors:  Christina M Rost; Edward Sachet; Trent Borman; Ali Moballegh; Elizabeth C Dickey; Dong Hou; Jacob L Jones; Stefano Curtarolo; Jon-Paul Maria
Journal:  Nat Commun       Date:  2015-09-29       Impact factor: 14.919

View more
  1 in total

Review 1.  Connecting battery technologies for electric vehicles from battery materials to management.

Authors:  Gang Zhao; Xiaolin Wang; Michael Negnevitsky
Journal:  iScience       Date:  2022-01-07
  1 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.