Wei Wei
Zeng-Xia Zhao1
Wei Li- 1Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun, China
- 2College of Chemistry, Jilin University, Changchun, China
- 3School of Chemistry and Materials Science, Hunan Agricultural University, Changsha, China
Expanded porphyrin systems are capable of binding a variety of substrates due to their increased cavity size and aromatic nature, holding important applications as magnetic resonance imaging contrast agents and as sensitizers for photodynamic therapy. It is there of fundamental interest to know the photoelectrical properties of expanded porphyrins using quantum chemistry calculations. In this work, we theoretically designed and screened a series of expanded porphyrins by incorporating terthiophene (TTH) and dithienothiophene (DTT) moieties. Our calculations showed that all the designed molecules exhibit excellent optoelectronic performance than the reference molecule. It is suggested that the porphyrin molecule with TTH moiety has better stability than the one with DTT moiety. Finally, we demonstrated that molecule 2 features with TTH moiety and the inverted selenophene ring outperform other molecules because it exhibits increased HOMO-LUMO gap, planar geometry, and strengthened aromaticity. We expect that this work can provide theoretical guidelines for the design of novel porphyrin materials.
Introduction
Aromaticity is a fundamental and important concept in the field of chemistry with wide applications (Shen et al., 2018). Many studies have been focusing on the aromatic molecule, benzene (Gil-Guerrero et al., 2018). Aromaticity is directly related to the stability of molecules. In general, aromatic compounds have high thermodynamic stability and are substantially more stable than antiaromatic compounds. Although aromaticity cannot be directly measured experimentally, various parameters based on geometric (Listunov et al., 2018), thermodynamics (Nowroozi1 and Rad, 2017; Rakhi and Suresh, 2016), magnetic (Najmidin et al., 2013; Torrent-Sucarrat et al., 2017), and electronic properties can be used to evaluate aromaticity and antiaromaticity of molecules. In particular, the nucleus-independent chemical shift (NICS) analysis is the most popular method for determining the aromaticity of molecules due to its simplicity and effectivity (Schleyer et al., 1996). In addition to NICS, the highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gap value also has been used to characterize the aromaticity of a molecule (Bédé et al., 2018).
Porphyrin is a good aromatic compound (Cao et al., 2014; Cyrański et al., 2003; Wei et al., 2019; Wei et al., 2020). Porphyrins with 18 π-electrons have been under intense investigation due to their enormous potential for applications in catalysis, magnetic resonance imaging, and photodynamic therapy (Drain et al., 2009; Bryden and Boyle, 2016; Barona-Castaño et al., 2016; Tang et al., 2019). Porphyrin frameworks are composed of four pyrrole subunits that are connected in a coplanar fashion at their α-carbon atoms via a methine carbon bridge. Expanded porphyrins exhibit structures, electronic properties, coordination chemistry, and reactivities that are entirely different from those of porphyrins. Expanded porphyrins exhibit intriguing Hückel and Mobius topologies than tetrapyrrolic porphyrin macrocycles in view of their larger-sized macrocyclic core (Alka et al., 2019), presenting ideal candidates for novel nonlinear optical (NLO) switching devices (Desmedt et al., 2021). Since the first report by Johnson et al., sapphyrins have been well investigated in therapeutic potential due to their optical and material properties (Chatterjee et al., 2017). Compared to porphyrin, sapphyrin contains an additional pyrrole inserted between a meso-carbon and an α-pyrrolic position. There are also modified sapphyrins, that is, oxasapphyrins and dithiabenzisapphyrin (Richter and Lash, 1998; Panda et al., 2005; Jeong et al., 2008), with pyrrolic subunits replaced by other building fragments, exhibiting unusual aromaticity, stability, and absorption properties (Karthik et al., 2014; Szyszko et al., 2017; Rao et al., 2019).
Terthiophene (TTH) has been widely used as the π-conjugated building fragments to substitute pyrrole unit in porphyrin skeleton (Moriaty et al., 1985). Dithienothiophene (DTT), which contains two thiophene rings fused by a bridging sulfur atom, is an electron-rich moiety with the rigid fragment that has been used in many optoelectronic materials (Frey et al., 2002). Both TTH and DTT have three thiophene rings, and the difference is that TTH has the C-C single bond bridge. Therefore, the incorporation of TTH and DTT subunits is a promising way to improve the optoelectronic characteristics of porphyrin molecules. The X-ray analysis of sapphyrins structure reveals that the larger core size and the presence of lighter heteroatoms (i.e., N or O) adjacent to the heterocyclic rleadeads to an inverted structure, while the presence of heavier heteroatoms (i.e., S or Se) leads to the conventional structure (Misra and Chandrashekar, 2008). It would be interesting to investigate how this strategy will change the aromaticity, stability, and optoelectronic properties of sapphyrins.
In this work, we theoretically designed eight macrocycle molecules by incorporating TTH and DTT subunits, coded as 1–8, Scheme 1. The aromaticity, stability, and photophysical properties of the designed molecules are investigated by using density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations. We also investigated the substitutions of the pyrrole N atom by heavier S and Se atoms in expanded porphyrins. This study aims to elucidate the impact of structural modifications on the overall performance of these expanded porphyrins. We expect the current theoretical work can pave the way for the design of novel porphyrin-based optoelectronic materials.
SCHEME 1. Sketch structures of molecules of 1–8.
Computational details
All calculations are performed with DFT/TD-DFT methods, as implemented in the Gaussian 09 software (Frisch et al., 2009). Geometry optimization was carried out using the DFT-B3LYP/6-311G (d,p) level of theory. Vibrational frequencies are calculated for the optimized structure at the same theory level to confirm that there are no imaginary frequencies and that the obtained structure is in local energy minima. Vertical electronic excitations are calculated using TD-DFT/CAM-B3LYP functional together with a 6-311G (d,p) basis set. Based on the optimized ground state structures, the NICS (1) values are calculated by means of the “atoms in molecules” theory, as implemented in the AIM2000 package (Bader et al., 1994). In this work, a few functionals (CAM-B3LYP, M05-2X, M06-2X) are performed on the calculation. The NICS (1) values with CAM-B3LYP are very close to M05-2X and M06-2X functionals. Therefore, CAM-B3LYP functional is the best choice due to its commonly used. Aromatic stabilization energy (ASE) is also used to evaluate the stability and aromaticity of molecules. ASE is an energy-based criterion for analyzing aromaticity. The ASE values can be calculated according to a ring-opening isobond chemical reaction, as we will introduce below. The molecule is aromatic if the reaction energy is positive, otherwise, the molecule is antiaromatic. Readers can refer to ref. 29 for more details.
Results and discussion
The ground state structures of the designed molecules are optimized at the B3LYP/6-311G (d, p) level of theory, as presented in Figure 1. The difference between TTH and DTT units lies in the presence of the C-C single-bond bridge in the former. The designed molecules are constructed by replacing the pyrrole fragment of the porphyrin ring with either a TTH unit (1–4) or a DTT unit (5–8). Molecules 1, 2, 5, and 6 feature the N heteroatom in the pyrrole group, whereas molecules 3, 4, 7 and 8 contain the S heteroatom in the pyrrole group. The key geometrical parameters of all molecules are listed in Table 1.
FIGURE 1. Optimized ground state structures of molecules 1–8.
TABLE 1. Main optimized geometry structure parameters of 1–8.
It can be seen that the N7-N8 bond lengths of molecules 1, 2, 5, and 6 are 7.058, 7.564, 6.106, and 6.172 Å, respectively. The S7-S8 bond lengths of molecules 3, 4, 7, and 8 are 6.520, 6.578, 5.175, and 5.394 Å, respectively. It is obviously that both bond lengths increase in the order of 7 < 8 < 5 < 6 < 3 < 4 < 1 < 2. For bond angles, C4-C5-C6 of molecules 1-8 is 125.5°, 125.6°, 129.2°, 129.7°, 119.8°, 120.8°, 123.2°, and 123.4°, respectively. The C9-C10-C11 bond angles of all molecules change following the same order as the C4-C5-C6 bond angles. The dihedral angles of S1-C4-C5-C6, N7-C9-C10-C11, and S7-C9-C10-C11 are listed in Table 1 as well. The dihedral angles of 1–4 are nearly zero, indicating that the four molecules would have better coplanarity. Figure 1 shows that 1–4 has planner structure, whereas 5 and 6 exhibit bowl structure. For 7 and 8, the DTT unit is perpendicular to the porphyrin core, whereas the three thiophene units of the porphyrin core are nearly coplanar, as evidenced by the values of dihedral angles. The calculated ASE values are listed in Table 2. We choose TTH, DTT, and 5-membered rings as the reference structures that present localized single and double bonds (Scheme 2). Generally, the molecule with a larger (more negative) ASE value would have stronger stability. It is shown from Table 2 that the calculated ASE values decrease in the order of 2 > 1 > 4 > 3 > 6 > 5 > 8 > 7, indicating that stability of molecules decreases in the same order.
TABLE 2. Computed aromatic stabilization energy (ASE, in kcal/mol) of molecules 1–8.
SCHEME 2. The isodesmic reactions of 1–8.
To gain a deep understanding of the stability, the aromaticity of atomic rings of all molecules is calculated. Positive values of NICS indicate antiaromaticity. The more negative (smaller) NICS value suggests better stability. The critical points are analyzed at the geometrical center of a ring or above a molecular plane. To avoid the in-plane components, the NICS values are calculated at 1 Å above the thiophene and pyrrole rings (points a, b, c, d, e, and f), molecular center (point k), and intramolecular (points g, h, i, and j) using the flow around the gauge-independent atomic orbital method. An illustration of the different critical points is detailed in Scheme 1. This approach has been widely used to evaluate the aromaticity and stability of porphyrins, including our and others’ works (Hong and Kim, 2003; Jeon et al., 2004).
The calculated NICS (1) values, as listed in Table 2, show that NICS (1) at 1 Å above the critical points a-k are all negative for molecules 1–8. The NICS (1) values at 1 Å above the critical points a, b, and c for 1–4, and critical points a and c for 5–8 are close. An exception is the NICS (1) value of 5–8 at 1 Å above the critical point of ring b show much smaller absolute values. It shows that the middle thiophene moiety in TTH has larger aromaticity than in DTT. In addition, the NICS (1) values of 1–8 at 1 Å above the critical points of rings d and f are close. The intramolecular 1 Å above critical points g, h, i, and j for 1–8 are similar. Moreover, comparisons of the NICS (1) values at critical point e for molecules with pyrrole N replaced by Se and S suggest that the inverted thiophene ring and selenophene ring are more stable than the conventional structure. As for the molecular center 1 Å above critical point, the NICS (1) values are −20.56, −21.21, −17.43, −18.12, −16.54, −16.91, −15.38, and −15.67 ppm for 1–8, respectively, Table 3. It can be concluded from the aforementioned discussions that the aromaticity of all molecules increase in the order of 7 < 8 < 5 < 6< 3 < 4< 1 < 2. We note that the NICS (1) value of porphyrin is only -14.98 ppm (Wei et al., 2012). This indicates that all designed molecules have strong aromaticity than porphyrin.
TABLE 3. NICS (1) values for 1–8.
Energy levels of HOMO, LUMO, and their gap, ΔH-L, are also used to evaluate the stability of compounds. The Frontier molecular orbital diagrams and the ΔH-L of molecules 1–8 are shown in Figure 2. Both the HOMOs and LUMOs spread over the whole π-conjugated backbones. For 1 and 2, HOMO orbitals exhibit bonding characters, whereas LUMOs exhibit antibonding characters. In contrast to 1 and 2, the HOMO of 3 and 4 show antibonding character, while the LUMO exhibits bonding character. For 5 and 6, the HOMOs display antibonding characters in DTT, while the LUMOs display bonding characters in DTT. For 7 and 8, the HOMOs display antibonding character at the bottom thiophene ring, while the LUMOs display bonding character at the bottom thiophene ring. A larger HOMO-LUMO gap indicates a more stable structure. The HOMO-LUMO gaps of 1 (2.29 eV), 2 (2.30 eV), 3 (2.14 eV), 4 (2.17 eV), 5 (1.94 eV), 6 (1.95 eV), 7(1.60 eV) and 8 (1.64 eV) are listed in Figure 2. For 1–8, molecules 1–4 have larger HOMO-LUMO gaps than 5–8. A positive correlation between NICS (1) (in absolute value) and HOMO-LUMO energy gaps can be observed, that is the larger HOMO-LUMO gaps would have the more negative NICS (1) value. A similar observation has been reported in our previous works [33-34].
FIGURE 2. Molecular orbital distributions and diagrams of the frontier molecular orbitals (from HOMO-3 to LUMO+3) for (A) 1, 2, 3, and 4 and (B) 5, 6, 7, and 8.
The TD-DFT method is used to calculate the vertical electronic excitations, as reported in Table 4. The simulated absorption spectrum is shown below. From Figure 3, it can be seen that, for all molecules, the major absorption bands are located in the visible light region. The dominated transitions for molecules 1 and 2 can be attributed to HOMO→LUMO + 1. For 3–8, the major transitions are assigned to HOMO→LUMO. It is clear that the absorptions are significantly redshifted and intensified with the substitution of the pyrrole N atom by Se atoms as compared to one substitution by S atoms. One can also observe that the absorptions are significantly redshifted and broadened when incorporating DTT moiety. For 2, the lowest-lying excitation is calculated at 707 nm, and the major transition of HOMO→LUMO has a maximum oscillator strength of 1.403.
TABLE 4. Calculated absorption features of 1–8.
FIGURE 3. The simulated absorption spectra of molecules 1–8 were obtained using TD-DFT/CAM-B3LYP/6-311G (d,p) level of theory.
Conclusion
We have investigated the aromaticity, stability, and photophysical properties of 1–8 with TTH and DTT units incorporated in the porphyrins using DFT and TD-DFT calculations. The stability of molecules is directly related to aromaticity, which can be evaluated using the ΔH-L values and the NICS (1) values. The calculation results show that all the designed molecules show larger aromaticity and excellent photophysical properties. It is also revealed that the molecule with TTH unit is more stable than the one with the DTT unit, which can be ascribed to the coplanarity of the TTH unit. For the designed molecules, the one with an inverted thiophene ring or selenophene ring is more stable than the one with a thiophene ring or selenophene ring. Substitution of the pyrrole N atom by Se heteroatom results in a more stable structure than by S heteroatoms. All the absorption bands of the molecules are located in the visible region. The major absorptions are from either HOMO→LUMO + 1 or HOMO→LUMO transitions, which are beneficial for the intramolecular charge transfer process. Compared to the other molecules, compound 2 stands out because it shows unique properties, that is, larger HOMO-LUMO gap, better coplanarity, and stronger aromaticity. We expect these theoretical studies could pave the way for the future development of novel porphyrin molecules.
Data availability statement
All datasets generated for this study are included in the article/supplementary material.
Author contributions
WW conducted the calculation and drafted the manuscript. Z-XB and B-HX assisted the calculation. WW and WL proposed the idea and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This word was supported by the Natural Science Foundation of China (Grant Nos. 21903034 and 21903023), China Postdoctoral Science Foundation (2021M701168), and the Science and Technology Innovation Program of Hunan Province (no. 2021RC3089).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The handling editor FB declared a shared affiliation with the authors WW, ZZ, and BX at the time of review.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Alka, A., Shetti, V. S., and Ravikanth, M. (2019). Coordination chemistry of expanded porphyrins. Coord. Chem. Rev. 401, 213063. doi:10.1016/j.ccr.2019.213063
Bader, R. F. W., Popelier, P. L. A., and Keith, T. A. (1994). Theoretical definition of a functional group and the molecular orbital paradigm. Angew. Chem. Int. Ed. Engl. 33, 620–631. doi:10.1002/anie.199406201
Barona-Castaño, J. C., Carmona-Vargas, C., Brocksom, T. J., and de Oliveira, K. T. (2016). Porphyrins as catalysts in scalable organic reactions. Molecules 21, 310–346. doi:10.3390/molecules21030310
Bédé, A. L., Assoma, A. B., Yapo, K. D., Koné, M. G.-R., Koné, S., Koné, M., et al. (2018). Theoretical study by density functional theory method (DFT) of stability, tautomerism, reactivity and prediction of acidity of quinolein-4-one derivatives. Comput. Chem. 6, 57–70. doi:10.4236/cc.2018.63005
Bryden, F., and Boyle, R. W. (2016). Chapter four - metalloporphyrins for medical imaging applications. Adv. Inorg. Chem. 68, 141–221.
Cao, X. Y., Zhao, Q. Y., Lin, Z. Q., and Ma, H. P. (2014). The chemistry of aromatic osmacycles. Acc. Chem. Res. 47, 341–354. doi:10.1021/ar400087x
Chatterjee, T., Srinivasan, A., Ravikanth, M., and Chandrashekar, T. K. (2017). Smaragdyrins and sapphyrins analogues. Chem. Rev. 117, 3329–3376. doi:10.1021/acs.chemrev.6b00507
Cyrański, M. K., Schleyer, P. v. R., Krygowski, T. M., Jiao, H., and Hohlneicher, G. (2003). Facts and artifacts about aromatic stability estimation. Tetrahedron 59, 1657–1665. doi:10.1016/s0040-4020(03)00137-6
Desmedt, E., Woller, T., Teunissen, J. L., Vleeschouwer, F. D., and Alonso, M. (2021). Fine-tuning of nonlinear optical contrasts of hexaphyrin-based molecular switches using inverse design. Front. Chem. 9, 786036. doi:10.3389/fchem.2021.786036
Drain, C. M., Varotto, A., and Radivojevic, I. (2009). Self-organized porphyrinic materials. Chem. Rev. 109, 1630–1658. doi:10.1021/cr8002483
Frey, J., Bond, A. D., and Holmes, A. B. (2002). Improved synthesis of dithieno[3, 2- b:2′ , 3′ -d]thiophene (DTT) and derivatives for cross coupling. Chem. Commun. 20, 2424–2425. doi:10.1039/b207403f
Frisch, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., et al. (2009). Gaussian 09 revision D.01 Gaussian. Wallingford, CT: Inc.
Gil-Guerrero, S., Ramos-Berdullas, N., and Mandado, M. (2018). Can aromaticity enhance the electron transport in molecular wires? Org. Electron. 61, 177–184. doi:10.1016/j.orgel.2018.05.043
Hong, S. Y., and Kim, S. C. (2003). Towards designing environmentally stable conjugated polymers with very small band-gaps. Bull. Korean Chem. Soc. 24, 1649–1654.
Jeon, K. O., Jun, J. H., Yu, J. S., and Lee, C. K. (2004). Mass spectroscopical properties of benzoyl derivatives of five-membered monoheterocycles and determination of aromaticity indices. Bull. Korean Chem. Soc. 25, 1840–1844.
Jeong, S. D., Sessler, J. L., Lynch, V., and Lee, C. H. (2008). Dithiabenzisapphyrin: a core-modified sapphyrin bearing exocyclic double bonds at the meso-positions. J. Am. Chem. Soc. 130, 390–391. doi:10.1021/ja075800p
Karthik, G., Srinivasan, A., and Chandrashekar, T. K. (2014). Meso-aryl core-modified fused sapphyrins: syntheses and structural diversity. Org. Let. 9, 734–738.
Listunov, D., Duhayon, C., Poater, A., MazHres, S., Saquet, A., Maraval, V., et al. (2018). Steric/π-electronic insulation of the carbo-benzene ring: dramatic effects of tert-butyl versus phenyl crowns on geometric, chromophoric, redox, and magnetic properties. Chem. Eur. J. 24, 10699–10710. doi:10.1002/chem.201800835
Misra, R., and Chandrashekar, T. K. (2008). Structural diversity in expanded porphyrins. Acc. Chem. Res. 41, 265–279. doi:10.1021/ar700091k
Moriaty, R. M., Prakash, O., and Duncan, M. P. (1985). Synthesis of 3, 2′ :5′ , 3′′ - terthiophene and 2, 5-di(3′ -thienyl)furan. Synth. Commun. 15, 789–795. doi:10.1080/00397918508063874
Najmidin, K., Kerim, A., Abdirishit, P., Kalam, H., and Tawar, T. (2013). A comparative study of the aromaticity of pyrrole, furan, thiophene, and their aza-derivatives. J. Mol. Model. 19, 3529–3535. doi:10.1007/s00894-013-1877-x
Nowroozi1, A., and Rad, O. R. (2017). A comparative study of cooperative effects between the intramolecular hydrogen bond and cation π interaction in various complexes of ortho-aminobenzaldehyde with its thio and seleno analogous. Theor. Chem. Acc. 136, 23–31. doi:10.1007/s00214-016-2047-x
Panda, P. K., Kang, Y., and Lee, C. H. (2005). A benzodipyrrole-derived sapphyrin. Angew. Chem. Int. Ed. 44, 4053–4055. doi:10.1002/anie.200500944
Rakhi, R., and Suresh, C. H. (2016). A DFT study on dihydropyrazine annulated linear polyacenes: aromaticity, stability and HOMO–LUMO energy modulation. Phys. Chem. Chem. Phys. 18, 24631–24641. doi:10.1039/c6cp03723b
Rao, Y. T., Zhou, W. J., Xu, L., Zhou, M. B., Yin, B. S., Tanaka, T., et al. (2019). Singly and doubly neo-confused smaragdyrins. J. Am. Chem. Soc. 141, 18836–18844. doi:10.1021/jacs.9b10270
Richter, D. T., and Lash, T. D. (1998). Synthesis of the first expanded carbaporphyrinoid by the “4 + 1” approach. J. Am. Chem. Soc. 120, 9965–9966. doi:10.1021/ja981887q
Schleyer, P. v. R. M., Maerker, C., Dransfeld, A., Jiao, H. J., and Hommes, N. J. R. v. E. (1996). Nucleus-independent chemical shifts: a simple and efficient aromaticity probe. J. Am. Chem. Soc. 118, 6317–6318. doi:10.1021/ja960582d
Shen, C. F., Liu, Z. Z., Liu, H. X., and Zhang, H. Q. (2018). Bond length equalization with molecular aromaticity—a new measurement of aromaticity. Spectrochimica Acta Part A Mol. Biomol. Spectrosc. 201, 392–398. doi:10.1016/j.saa.2018.05.007
Szyszko, B., Białek, M. J., Pacholska-Dudziak, E., and Latos-Grazyṅski, L. ́(2017). Flexible porphyrinoids, Chem. Rev. 117, 2839–2909. doi:10.1021/acs.chemrev.6b00423
Tang, Z., Jiang, Z., Chen, H. J., Su, P. F., and Wu, W. (2019). Energy decomposition analysis based on broken symmetry unrestricted density functional theory. J. Chem. Phys. 151, 244106. doi:10.1063/1.5114611
Torrent-Sucarrat, M., Navarro, S., Cossίo, F. P., Anglada, J. M., and Luis, J. M. (2017). Relevance of the DFT method to study expanded porphyrins with different topologies. J. Comput. Chem. 38, 2819–2828. doi:10.1002/jcc.25074
Wei, W., Bai, F. Q., Xia, B. H., and Zhang, H. X. (2012). Theoretical analysis on molecular magnetic properties of N-confused porphyrins and its derivatives. Bull. Korean Chem. Soc. 33, 2937–2942. doi:10.5012/bkcs.2012.33.9.2937
Wei, W., Li, W., Xia, B. H., Zhang, H. X., and Bai, F. Q. (2019). Theoretical investigations of the aromaticity, stability and photophysical behaviors for expanded porphyrins. Chem. Phys. Lett. 728, 25–31. doi:10.1016/j.cplett.2019.04.067
Keywords: DFT, aromaticity, stability, optoelectronic properties, expanded porphyrins
Citation: Wei W, Zhao Z-X, Xia B-H and Li W (2022) Theoretical analysis of expanded porphyrins: Aromaticity, stability, and optoelectronic properties. Front. Chem. 10:948632. doi: 10.3389/fchem.2022.948632
Received: 20 May 2022; Accepted: 26 July 2022;
Published: 01 September 2022.
Edited by:
Fu-Quan Bai, Jilin University, ChinaReviewed by:
Shamoon Ahmad Siddiqui, Najran University, Saudi ArabiaYuanzuo Li, Northeast Forestry University, China
Copyright © 2022 Wei, Zhao, Xia and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Wei, d2Vpd2Vpd3dAamx1LmVkdS5jbg==; Wei Li, d2VpbGlAaHVuYXUuZWR1LmNu