收藏切换
High-throughput screening of SrxA1−xFeyB1−yO3 perovskites for low-temperature chemical looping air separation using graph neural networks
收藏切换
PDF
Jie ZHAO1, Changqing DONG1, 2, *, Junjie XUE1, Xiaoying HU1, Junjiao ZHANG3
Journal of Fuel Chemistry and Technology | 2026, 54(7) : 20250302-1423 - 20250302-1429
Less
收藏切换
Journal of Fuel Chemistry and Technology | 2026, 54(7): 20250302-1423-20250302-1429
• RESEARCH PAPER •
High-throughput screening of SrxA1−xFeyB1−yO3 perovskites for low-temperature chemical looping air separation using graph neural networks
Full
Jie ZHAO1, Changqing DONG1, 2, *, Junjie XUE1, Xiaoying HU1, Junjiao ZHANG3
Affiliations
  • 1National Engineering Laboratory for Biomass Power Generation Equipment (NCEPU), School of New Energy, North China Electric Power University, Beijing 102206, China
  • 2New Energy Generation National Engineering Research Center (NCEPU), School of New Energy, North China Electric Power University, Beijing 102206, China
  • 3School of Energy Power and Mechanical Engineering, North China Electric Power University, Beijing 102206, China
Published: 2026-07-10 doi: 10.1016/S1872-5813(26)60680-9
Outline
收藏切换

Low-temperature Chemical Looping Air Separation (CLAS) is a promising technology for producing oxygen-enriched gas streams. It utilizes the redox properties of solid oxygen carriers to selectively capture and release oxygen from the air. Oxygen vacancy formation energy (Eovf) can be used to evaluate the ease of oxygen release by the oxides. In this study, the applicable range of Eovf for CLAS oxygen carriers was determined to be <2.3 eV using the thermodynamic calculations. A graph neural network (GNN) was trained and used to predict the Eovf. By predicting Eovf with the GNN model, a screening from the 3649 compositions of SrxA1−xFeyB1−yO3 perovskites was conducted to discover the CLAS oxygen carriers. The prediction results indicate that the doping of Ba and Ca ions at the A site and Co ions at the B-site can effectively reduce the Eovf of the perovskite structure. The feasibility of screening the CLAS oxygen carriers using the Eovf method was validated through the rediscovery some known materials. These results support the effectiveness of GNNs in accelerating the discovery of CLAS oxygen carriers.

chemical looping air separation  /  oxygen vacancy formation energy  /  graph neural network  /  material discovery
赵杰, 董长青, 薛俊杰, 胡笑颖, 张俊姣. 使用图神经网络高通量筛选用于低温化学链空气分离的SrxA1−xFeyB1−yO3钙钛矿. 燃料化学学报(中英文), 2026 , 54 (7) : 20250302-1423 -20250302-1429 . DOI: 10.1016/S1872-5813(26)60680-9
Jie ZHAO, Changqing DONG, Junjie XUE, Xiaoying HU, Junjiao ZHANG. High-throughput screening of SrxA1−xFeyB1−yO3 perovskites for low-temperature chemical looping air separation using graph neural networks[J]. Journal of Fuel Chemistry and Technology, 2026 , 54 (7) : 20250302-1423 -20250302-1429 . DOI: 10.1016/S1872-5813(26)60680-9
Chemical looping air separation (CLAS) is an attractive technology for industrial applications, owing to its high efficiency, low energy consumption, and the potential to integrate with thermoelectric power generation cycles to reduce CO2 emissions[1−3]. CLAS operates using oxide-based oxygen carriers that can reversibly absorb and release O2 in response to changes in the oxygen partial pressure or temperature. As shown in Figure 1, when the oxygen partial pressure is low, the lattice oxygen in the oxide is released in the form of gaseous oxygen, accompanied by the reduction of the oxide. Subsequently, the reduced oxide can be regenerated in air or oxygen-rich stream. Repeating this process in a cyclic manner enables a continuous air separation and oxygen production.
Perovskite materials possess superior oxygen transport capabilities, high chemical stability, and tunable redox properties, which render them ideal oxygen-carrier candidates for chemical looping (CL) processes[1,4−6]. SrFeO3−δ with a perovskite structure has been studied for low-temperature CLAS reactions (400−650 °C) due to its excellent oxygen release and absorption capabilities[2,7]. The perovskite structure of SrFeO3−δ provides the possibility of doping substitution atoms at the A-site or B-site to further regulate and enhance the material´s redox properties[7−8]. For instance, the doping of A-site with Ba[9] or Ca[8,10] and B-site with Co[11−12] of SrFeO3 has demonstrated an enhanced CLAS performance. However, the search for optimal doping elements and concentrations still relies largely on conventional heuristics and trial-and-error methods.
Oxygen vacancy formation energy (Eovf) can be used to evaluate the redox performance of CLAS materials[13−15]. Oxygen carrier materials with lower Eovf are more likely to release lattice oxygen and have a better performance in low-temperature CLAS[16−17]. Density functional theory (DFT) method is commonly used to calculate the Eovf of oxides in crystal structures[18−19]. The Eovf can be calculated through Equation (1):
$ \begin{matrix}{E}_{{\mathrm{ovf}}}={E}_{{\mathrm{def}}}-{E}_{{\mathrm{perf}}}+{\mu }_{{\mathrm{O}}} \end{matrix} $
where $ {E}_{{\mathrm{def}}} $ is the total energy of the supercell with an O vacancy, $ {E}_{{\mathrm{perf}}} $ is the total energy of the pristine supercell, $ {\mu }_{{\mathrm{O}}} $ is the O chemical potential. However, calculating the Eovf using DFT is expensive and time-consuming.
Machine learning methods have been researched for rapid prediction of Eovf. Linear models based on DFT descriptors such as band gap and formation energy can predict Eovf with a mean absolute error (MAE) of 0.39 eV, without the need to construct a supercell with an oxygen vacancy[20−21]. The unified atomistic line graph neural network force-field (ALIGNN-FF), trained with formation energy data of the defect-free crystals, can be directly applied to predict the Eovf with an MAE of 0.87 eV[22−23]. The defect graph neural network (dGNN) model, trained with Eovf data calculated by DFT, predicts the Eovf of oxides (containing elements present in the training set) with an MAE of 0.52 eV[14].
In this study, an ALIGNN model was trained using a curated dataset of Eovf and used for the prediction of Eovf. The range of Eovf suitable for CLAS oxygen carriers was estimated using thermodynamic calculations. High-throughput screening for CLAS was performed on 3649 SrxA1−xFeyB1−yO3 perovskites through stability judgment and Eovf prediction. The feasibility of the Eovf method for screening CLAS oxygen carriers was confirmed through materials rediscovery.
The charge-neutral Eovf dataset[14,24] utilized in this article was based on DFT-computed data using the Generalized Gradient Approximation (GGA) with Hubbard U correction. The dataset comprises 48 cationic elements, 102 space groups, 901 compounds, 2329 distinct oxygen vacancy structures, 9342 oxygen sites (see Section S1 in Supplemental Information). The distribution of elements in the dataset is shown in Figure 2(a).
The Atomistic Line Graph Neural Network (ALIGNN) demonstrated superior performance over other graph neural networks in multiple atomistic prediction tasks by incorporating bond angle information[23,25]. As shown in Figure 2(b), the ALIGNN was trained to predict the Eovf at all oxygen sites within the crystal structures. Details on the ALIGNN can be found in the related paper[23,25]. In brief, each node in the atomistic graph is assigned 9 input node features based on its atomic species: electronegativity, group number, covalent radius, valence electrons, first ionization energy, electron affinity, block and atomic volume. The inter-atomic bond distances are used as edge features with a radial basis function up to a cut-off of 8 Å. A periodic 12-nearest-neighbor graph construction was used. This atomistic graph is then used for constructing the corresponding line graph using interatomic bond-distances as nodes and bond-angles as edge features. The dataset was partitioned into training, validation, and test sets at a ratio of 80∶10∶10, and an ALIGNN model was trained to predict the Eovf.
The Eovf predicted by the ALIGNN corresponds to the zero-temperature limit. Therefore, thermodynamic calculations are required to determine the range of Eovf values that are suitable for various chemical looping applications. In light of the comparable Eovf across distinct sites in the ideal perovskite structure[14,26], it is reasonable to assume that the properties of oxygen vacancies have a degree of similarity among these sites. At a specified temperature T, the Gibbs free energy for the formation of oxygen vacancies in the crystalline structure is expressed as follows[27]:
$ \Delta G=n*{h}_{{\mathrm{ovf}}}-T{\Delta S}_{{\mathrm{cf}}} $
where n represents the number of vacancies, hovf denotes the formation enthalpy per defect (including contributions from vibrational enthalpy and entropy for each defect), and ΔScf is the configurational entropy. Assuming a total number of N atoms in the system (where N ≫ n), the structural entropy is derived via Stirling’s approximation as follows[27]:
$ {\Delta S}_{{\mathrm{cf}}}={k}_{{\mathrm{B}}}\times n \times (\ln N-\ln n)$
In the context of solids within the zero-temperature limit, the volume change is negligible, and thus the hovf can be reasonably approximated by the Eovf[28]. The chemical potential Δμ is introduced to determine the Gibbs free energy of oxygen vacancy formation under the relevant thermodynamic constraints[14]:
$ \Delta G=n \times \left({E}_{{\mathrm{ovf}}}+{\Delta }\mu \right)-T\times {k}_{{\mathrm{B}}}\times n \times (\ln N-\ln n)$
when the system reaches equilibrium, the condition $ \Delta G=0 $ leads to the following expression for the dimensionless oxygen vacancy concentration:
$ {c}_{v}=\frac{n}{N}={{\mathrm{e}}}^{-\left({E}_{{\mathrm{ovf}}}+{\Delta }\mu \right)/\left({k}_{{\mathrm{B}}}T\right)} $
Given a specific oxygen vacancy concentration, the corresponding Δμ at the thermodynamic constraints can be determined by inverting the above expression[14]. By comparing the Δμ curves corresponding to different Eovf values on the Ellingham diagram, the applicable ranges of Eovf for different chemical looping applications can be identified[13−14].
The SrxA1−xFeyB1−yO3 perovskite structures were constructed by substituting A- and/or B-site cations in SrFeO3 using the pymatgen package[29]. Table 1 summarizes the dopant types and concentrations investigated in this study. By exhausting all possible combinations of A- and/or B-site dopants in the selected cation network, 3520 perovskite models were constructed. Through a preliminary screening process, some of these combinations were filtered out due to their inability to form stable perovskite phases. Firstly, the total valence of the compositions was examined to determine whether it can be zero (see Section S2 in Supplemental Information). Besides charge neutrality, a modified tolerance factor $ \tau =\dfrac{{r}_{{\mathrm{O}}}}{{r}_{{\mathrm{B}}}}-{n}_{{\mathrm{A}}}\left({n}_{{\mathrm{A}}}-\dfrac{{r}_{{\mathrm{A}}}/{r}_{{\mathrm{B}}}}{\ln\left({r}_{{\mathrm{A}}}/{r}_{{\mathrm{B}}}\right)}\right) $[30] was applied to estimate the stabilities of the perovskite structures, where rO, rA, rB represent the radii of oxygen anion, and the A- and B-site cations (see Table S1 in Supplemental Information). nA and nB are the oxidation states of the A- and B-site cations, respectively. The modified tolerance factor (τ) carries more chemical information and exhibits better predictive ability than the classical Goldschmidt tolerance factor[4,30].
After the initial construction of the doped crystal structures, the fast structure optimization was performed using the ALIGNN-FF[23] model, for which the model parameters were directly adopted from the literature[23] and shared publicly. The ALIGNN-FF was developed using both energies and forces for applications, such as, in fast structure optimization. Then, the optimized crystal structures were input into the ALIGNN model which is then trained for Eovf to predict the Eovf for each oxygen site. The average of all the Eovf values for each crystal structure was calculated and used as an indicator to evaluate the redox performance. The process is illustrated in Figure 3.
As shown in Figure 4, the Gibbs free energy of several key reactions is used to distinguish different chemical looping reactions[13,31]. The data of reactions was from calculations with the FactSage software. For a target oxygen vacancy concentration of [VO] = 1%, Figure 4 shows the oxygen vacancy formation chemical potential Δμ at various temperatures for different values of Eovf. It has been shown that the trends of Eovf across different perovskite oxides remain the same for different extents of oxygen vacancy[32]. Therefore, it is reasonable to estimate the performance of the oxygen carrier by using Eovf of the chemical potential corresponding to [VO] = 1%. From Figure 4, it can be inferred that for CLAS, the suitable Eovf values should be approximately less than 2.3 eV when the temperature is higher than 500 ℃. Previous studies have indicated that materials with Eovf less than 2.5 eV are suitable for thermochemical O2 storage and pumping/separation[13,15], which is in agreement with the estimated results. Additionally, this estimation method can also be applied to assess the suitability of materials in fields such as chemical looping full oxidation (CLFO), chemical looping with oxygen uncoupling (CLOU), and chemical looping syngas production.
The ALIGNN model trained for predicting Eovf was tested. The oxides in the test set did not exist in the training set or the validation set. As shown in Figure 5, the ALIGNN model achieved a mean absolute error (MAE) of 0.26 eV in predicting the Eovf values at different oxygen sites of the oxides in the test set, and the coefficient of determination (R2) reached 0.94. The closer the MAE is to 0 and the R2 is to 1, the higher the accuracy of the model. This result outperforms the linear model using the DFT descriptors (0.39 eV)[20−21] and the dGNN model based on a simpler convolutional architecture relative to the ALIGNN (0.52 eV)[14].
The screening results from the 3649 compositions of SrxA1−xFeyB1−yO3 perovskites are presented in Figure 6. Most of the compositions containing ions of Na, K, Mg and Cu were filtered out due to the inability to achieve electrical neutrality. It is noted that the nonstoichiometric cases of perovskites (SrxA1−xFeyB1−yO3−δ) were not considered in the preliminary screening process. The different grids in the Figure 6 may represent the same chemical composition.
The slight differences in the average Eovf for the same material, such as SrxBa1−xFe1Ti0O3 and SrxBa1−xFe1V0O3, may be attributed to random errors introduced during the crystal structure optimization process. For the Sr1−xCaxFeO3, it is known that the Eovf first decreases and then increases with the increase of the Ca2+ concentration from 0 to 0.5[10]. This is consistent with the predictions. In certain compositions, such as Sr0.125Ba0.875CoO3, the doping of the Ba ion at the A-sites also leads to a decrease in Eovf. The doping of the Co ion at the B-sites shows a pronounced trend of decreasing Eovf with an increasing concentration. Among all the compositions with a predicted average Eovf below 2.0 eV, the Co ion is consistently present (see Section S4 in Supplemental Information).
Through the screening criterion Eovf < 2.3 eV, a series of known low-temperature CLAS oxygen carriers were rediscovered, as shown in Table 2. This indicates that the Eovf can be used to accelerate the material discovery for the low-temperature CLAS process. Still, this method has certain limitations. The Eovf can be used to assess the ease of oxygen release from the oxides, but it cannot be used to evaluate the material´s re-oxidation extent. For example, although the predicted Eovf for BaCoO3 is as low as 1.44 eV, it is difficult to be fully oxidized after oxygen loss, making it unsuitable for low-temperature CLAS[12].
Based on the graph neural network for predicting the Eovf, this work proposes a fast method for the screening of low-temperature CLAS oxygen carriers, and conducts a high-throughput screening from the 3649 compositions of SrxA1−xFeyB1−yO3 perovskites. Through the thermodynamic estimation, it was found that the Eovf value for the CLAS material should be less than 2.3 eV. An the ALIGNN model was trained for the prediction of the Eovf and the test set shows that the model has a MAE of 0.26 eV. In the high-throughput screening, the doping of the Co ions at the B-site was found to be able to effectively reduce the Eovf of the perovskite structure. And the doping of specific proportions of Ba and Ca ions at the A site can also contribute to lowering the Eovf. Under the screening criterion of Eovf < 2.3 eV, the model rediscovered several previously reported oxygen carriers for the low-temperature CLAS. This work demonstrates that using GNN to predict Eovf can accelerate the discovery process of CLAS materials. In a broader context, the presented methods may also accelerate the discovery of materials for other related chemical looping applications, such as chemical looping full oxidation and chemical looping syngas production.
1
KRZYSTOWCZYK E, HARIBAL V, DOU J, et al. Chemical looping air separation using a perovskite-based oxygen sorbent: System design and process analysis[J]. ACS Sustainable Chem Eng, 2021, 9(36): 12185−12195.
2
CAI R, DOU J, KRZYSTOWCZYK E, et al. Chemical looping air separation with Sr0.8Ca0.2Fe0.9Co0.1O3−δ perovskite sorbent: Packed bed modeling, verification, and optimization[J]. Chem Eng J, 2022, 429: 132370.
3
FAN Q, LI H, SAQLINE S, et al. An investigation of the structural and electronic origins of enhanced chemical looping air separation performance of B-site substituted SrFe1−xCoxO3−δ perovskites[J]. Phys Chem Chem Phys, 2024, 26(30): 20511−20521.
4
WANG X, GAO Y, KRZYSTOWCZYK E, et al. High-throughput oxygen chemical potential engineering of perovskite oxides for chemical looping applications[J]. Energy Environ Sci, 2022, 15(4): 1512−1528.
5
LIU C, PARK J, DE SANTIAGO H A, et al. Perovskite oxide materials for solar thermochemical hydrogen production from water splitting through chemical looping[J]. ACS Catal, 2024, 14(19): 14974−15013.
6
WANG Y, LI S, ZOU D, et al. Multifunctional SrFeO3-@Al composite macrocapsules for synergistic high-temperature thermal storage and oxygen separation[J]. Chem Eng J, 2025, 520: 166053.
7
DOU J, KRZYSTOWCZYK E, WANG X, et al. Sr1−xCaxFe1−yCoyO3−δ as facile and tunable oxygen sorbents for chemical looping air separation[J]. J Phys-Energy, 2020, 2(2): 025007.
8
CAPSTICK S, BULFIN B, NAIK J M, et al. Oxygen separation via chemical looping of the perovskite oxide Sr0.8Ca0.2FeO3 in packed bed reactors for the production of nitrogen from air[J]. Chem Eng J, 2023, 452: 139289.
9
BEKTAS H, CAI R, BRODY L, et al. Structural and thermodynamic assessment of Ba and Ba/Mg substituted SrFeO3−δ for “Low-Temperature” chemical looping air separation[J]. Energy Fuels, 2024, 38(12): 11107−11118.
10
JIA T, POPCZUN E J, LEKSE J W, et al. Effective Ca2+-doping in Sr1−xCaxFeO3−δ oxygen carriers for chemical looping air separation: A theoretical and experimental investigation[J]. Appl Energy, 2021, 281: 116040.
11
DOU J, KRZYSTOWCZYK E, MISHRA A, et al. Perovskite promoted mixed cobalt–iron oxides for enhanced chemical looping air separation[J]. ACS Sustainable Chem Eng, 2018, 6(11): 15528−15540.
12
EZBIRI M, ALLEN K M, GàLVEZ M E, et al. Design principles of perovskites for thermochemical oxygen separation[J]. ChemSusChem, 2015, 8(11): 1966−1971.
13
ZENG L, CHENG Z, FAN J A, et al. Metal oxide redox chemistry for chemical looping processes[J]. Nat Rev Chem, 2018, 2(11): 349−364.
14
WITMAN M D, GOYAL A, OGITSU T, et al. Defect graph neural networks for materials discovery in high-temperature clean-energy applications[J]. Nat Comput Sci, 2023, 3(8): 675−686.
15
WEXLER R B, GAUTAM G S, STECHEL E B, et al. Factors governing oxygen vacancy formation in oxide perovskites[J]. J Am Chem Soc, 2021, 143(33): 13212−13227.
16
MISHRA A, LI T, LI F, et al. Oxygen vacancy creation energy in Mn-containing perovskites: An effective indicator for chemical looping with oxygen uncoupling[J]. Chem Mater, 2018, 31(3): 689−698.
17
DOU J, KRZYSTOWCZYK E, WANG X, et al. A-and B-site codoped SrFeO3 oxygen sorbents for enhanced chemical looping air separation[J]. ChemSusChem, 2019, 13(2): 385−393.
18
WANG Z, LIN R, HUO Y, et al. Formation, detection, and function of oxygen vacancy in metal oxides for solar energy conversion[J]. Adv Funct Mater, 2021, 32(7): 2109503.
19
ZHOU M, WU J, LU Y, et al. Non-metallic phosphorus-induced lattice engineering in LaFeO3 perovskite: Oxygen vacancy modulation for enhanced methane conversion in chemical looping partial oxidation[J]. Chem Eng J, 2025, 520: 165963.
20
DEML A M, HOLDER A M, O’HAYRE R P, et al. Intrinsic material properties dictating oxygen vacancy formation energetics in metal oxides[J]. J Phys Chem Lett, 2015, 6(10): 1948−1953.
21
DEML A M, STEVANOVIĆ V, MUHICH C L, et al. Oxide enthalpy of formation and band gap energy as accurate descriptors of oxygen vacancy formation energetics[J]. Energ Environ Sci, 2014, 7(6): 1996.
22
CHOUDHARY K, SUMPTER B G. Can a deep-learning model make fast predictions of vacancy formation in diverse materials?[J]. Aip Adv, 2023, 13(9): 095109.
23
CHOUDHARY K, DECOST B, MAJOR L, et al. Unified graph neural network force-field for the periodic table: Solid state applications[J]. Digit Discov, 2023, 2(2): 346−355.
24
KUMAGAI Y, TSUNODA N, TAKAHASHI A, et al. Insights into oxygen vacancies from high-throughput first-principles calculations[J]. Phys Rev Mater, 2021, 5(12): 123803.
25
CHOUDHARY K, DECOST B. Atomistic line graph neural network for improved materials property predictions[J]. Npj Comput Mater, 2021, 7(1): 185.
26
EMERY A A, WOLVERTON C. High-throughput DFT calculations of formation energy, stability and oxygen vacancy formation energy of ABO3 perovskites[J]. Sci Data, 2017, 4(1): 1−10.
27
LANY S. Communication: The electronic entropy of charged defect formation and its impact on thermochemical redox cycles[J]. J Chem Phys, 2018, 148(7): 071101.
28
STEVANOVIĆ V, LANY S, ZHANG X, et al. Correcting density functional theory for accurate predictions of compound enthalpies of formation: Fitted elemental-phase reference energies[J]. Phys Rev B, 2012, 85(11): 641.
29
ONG S P, RICHARDS W D, JAIN A, et al. Python materials genomics (pymatgen): A robust, open-source python library for materials analysis[J]. Comp Mater Sci, 2013, 68: 314−319.
30
BARTEL C J, SUTTON C, GOLDSMITH B R, et al. New tolerance factor to predict the stability of perovskite oxides and halides[J]. Sci Adv, 2019, 5(2): eaav0693.
31
LUO S, ZENG L, FAN L-S. Chemical looping technology: oxygen carrier characteristics[J]. Annu Rev Chem Biomol, 2015, 6(1): 53−75.
32
MAITI D, DAZA Y A, YUNG M M, et al. Oxygen vacancy formation characteristics in the bulk and across different surface terminations of La(1−x)SrxFe(1−y)CoyO(3−δ) perovskite oxides for CO2 conversion[J]. J Mater Chem A, 2016, 4(14): 5137−5148.
33
MIURA N, IKEDA H, TSUCHIDA A. Sr1–xCaxFeO3−δ as a new oxygen sorbent for the high-temperature pressure-swing adsorption process[J]. Ind Eng Chem Res, 2016, 55(11): 3091−3096.
Year 2026 volume 54 Issue 7
PDF
105
46
Cite this Article
BibTeX
Article Info
doi: 10.1016/S1872-5813(26)60680-9
  • Receive Date:2025-10-23
  • Online Date:2026-09-23
  • Published:2026-07-10
Article Data
Affiliations
History
  • Received:2025-10-23
  • Revised:2025-11-21
  • Accepted:2025-12-14
Affiliations
    1National Engineering Laboratory for Biomass Power Generation Equipment (NCEPU), School of New Energy, North China Electric Power University, Beijing 102206, China
    2New Energy Generation National Engineering Research Center (NCEPU), School of New Energy, North China Electric Power University, Beijing 102206, China
    3School of Energy Power and Mechanical Engineering, North China Electric Power University, Beijing 102206, China

Corresponding:

References
Share
https://castjournals.cast.org.cn/joweb/rlhxxb/EN/10.1016/S1872-5813(26)60680-9
Share to
QR

Scan QR to access full text

Cite this article
BibTeX
Citations
表12种不同金属材料的力学参数
科
Family
属数
Number of
genus
种数
Number of
species
占总种数比例
Percentage of
total species (%)
属
Genus
种数
Number of
species
占总种数比例
Percentage of total
species (%)
鹅膏菌科Amanitaceae 2 11 5.26 鹅膏菌属 Amanita 10 4.78
小菇科 Mycenaceae 2 12 5.74 丝盖伞属 Inocybe 5 2.39
多孔菌科 Polyporaceae 8 14 6.70 蜡蘑属 Laccaria 5 2.39
红菇科 Russulaceae 3 23 11.00 小皮伞属 Marasmius 6 2.87
小菇属 Mycena 11 5.26
光柄菇属 Pluteus 5 2.39
红菇属 Russula 17 8.13
栓菌属 Trametes 5 2.39
关闭全屏
  • BibTeX
  • EndNote
  • RefWorks
  • TxT