Skip to main content

ORIGINAL RESEARCH article

Front. Public Health, 11 March 2021
Sec. Radiation and Health
This article is part of the Research Topic Medical Application and Radiobiology Research of Particle Radiation View all 16 articles

Design of Beam Shaping Assemblies for Accelerator-Based BNCT With Multi-Terminals

  • 1Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, China
  • 2College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
  • 3Key Laboratory of Heavy Ion Beam Radiation Biomedicine, Chinese Academy of Sciences, Lanzhou, China
  • 4Gansu Provincial Key Laboratory of Heavy Ion Beam Radiation Medical Application, Lanzhou, China

To moderate fast neutrons produced by accelerator to appropriate therapeutic energies for boron neutron capture therapy (BNCT), beam shaping assembly (BSA) is required definitely. In this work, based on a model of 2.5 MeV/30 mA proton accelerator, the Monte Carlo simulation software MCNPX was employed to design multi-terminal BSAs. All parameters for both the thermal and epithermal neutron beams at the exit ports of the designed BSAs meet the treatment recommendation values proposed by the International Atomic Energy Agency (IAEA). The clinical parameters of the thermal and epithermal neutron beams were also calculated for clinical indication consideration.

Introduction

Boron neutron capture therapy (BNCT) is a novel modality of radiation cancer therapy. Boron compounds are administered to tumor cells, then the tumor is irradiated with neutrons, inducing the 10B (n, α)7Li nuclear reaction:

 10B+nth11B4He(1.78 MeV)+7Li(1.01 MeV) (6.3%) 10B+nth11B4He(1.47 MeV)+7Li(0.84 MeV)+γ(0.48 MeV) (93.7%)    (1)

Both α particles and 7Li nuclei deposit their energies along their very short paths, which are comparable to the size of cells. As a result, tumor cells are destroyed accurately without harming healthy tissues (1).

BNCT was firstly proposed by G. Locher in 1936 (2) and firstly practiced by W. Sweet in 1951 for the clinical trial of glioma (3). Over the past two decades, many research groups around the world have continued the work of W. Sweet and the others, particularly the pioneering clinical work of Hatanaka (4). Subsequently, clinical trials of BNCT were conducted in the United States, Sweden, Finland, the Czech Republic, Argentina, the European Union (centered on Finland) and Japan.

Nuclear reactors were firstly used to produce neutrons for BNCT. However, although they could provide high-intensity neutron beam, they have numerous shortcomings: most of them are located far from hospitals, and are also very expensive. Besides, nuclear reactors have too huge size to be suitable for being used in hospital. Accelerator-based BNCT (AB-BNCT) facilities therefore are being developed to replace nuclear reactors. In AB-BNCT, fast neutrons are obtained by bombarding lithium or beryllium target with protons. However, the fast neutrons produced by this method cannot be used for BNCT treatment directly and need to be moderated by beam shaping assembly (BSA). The functions of BSA are: (1) slow fast neutrons down to thermal neutrons (<0.5 eV) or epithermal neutrons (0.5 eV−10 keV), (2) reduce the composition of fast neutron, thermal neutron and γ ray as much as possible, and (3) collimate neutron beam. The thermal neutron is suitable for treating superficial lesions while the epithermal neutron is for treating deep ones. BSA is mainly composed of the following components: moderator, reflector, gamma filter, collimator, etc., and thermal neutron filters are also required if neutrons need to be moderated to the energy range of epithermal neutrons.

Currently, the proposed BSA designs around the world mainly focus on generating epithermal neutrons which are essential for the treatment of deep-seated tumors, such as the BSAs in Tsukuba University (5), Nagoya University (6, 7), and Kyoto University (8). However, thermal neutron beam cannot be ignored anyway. It is applicable to the treatment of superficial tumors, such as melanoma, as well as cell and animal pre-clinical experiments. Therefore, multiple BSAs were designed for the generation of thermal and epithermal neutrons, respectively, and both of them fulfill the IAEA recommended values (9) which are listed in Table 1. In this work, the Monte Carlo simulation program MCNPX was used to design multiple BSAs based on an AB-BNCT model and the clinical parameters of the thermal and epithermal neutrons generated from the multiple BSAs were calculated, aiming at providing reference for the construction of AB-BNCT facility.

TABLE 1
www.frontiersin.org

Table 1. Neutron beam parameters and IAEA recommended values.

Methods and Materials

Initial BSA Model

An initial BSA model was proposed and shown in Figure 1. The whole BSA has a cylinder structure with lithium target and beam channel located in the central axis of the cylinder. The lithium target is 10 cm in diameter and 100 μm in thickness, and there is a copper holder of 2.3 cm in thickness below the target, which plays a role in heat dissipation and structure support. Thirty milliampere protons would generate lots of heat in target indeed and Li target has a low melt pointing. So, it is necessary to consider a cooling system for the Li target. In this work, copper was used to roughly represent the cooling system for target. The beam pipe is made of 316L stainless steel with a thickness of 1 cm for 2.5 MeV proton transport, and the thickness of the stainless steel above the target is 4 cm. In this way, the radiation damage of recoiling neutrons and protons can be reduced. Boron containing polyethylene (10 wt% natural B) is used in the outer side of the BSA as an absorption shield for neutrons. On the bottom of the collimator, a gamma shield is designed in the inner side of the collimator to further reduce the gamma component in the beam. The diameter of the BSA beam port is set to 14 cm. Finally, we used the MCNPX software to calculate the physical and clinical parameters of thermal and epithermal neutrons.

FIGURE 1
www.frontiersin.org

Figure 1. Structure of initial design of BSA.

In addition, in the optimization process of BSA for thermal neutrons, the energy of thermal neutron is close to the kinetic energy of nuclear thermal motion in materials of reflector, moderator etc. So, the library of the thermal scattering law data S (α, β) may be required. This library plays an important role in describing the transport of thermal neutrons (10). This work used the ENDF/B-VII.0 cross section library for the simulation of thermal neutron scattering.

Neutron Source

The neutron beam generated from bombarding lithium target by 2.5 MeV protons was simulated using the MCNPX software, whose spectrum and angular distribution are shown in Figure 2. The neutron beam was made as a dumb data file, which acted as a neutron source used in the BSA optimization. In this way, the computation time was greatly reduced.

FIGURE 2
www.frontiersin.org

Figure 2. Spectrum (A) and angular distribution (B) of the neutron beam generated from bombarding Li target by 2.5 MeV/30 mA protons.

Moderator and Reflector

The most important part of BSA is moderator. Its role is to moderate the energy of neutrons produced by protons into the energy range of thermal or epithermal neutrons without producing excessive gamma rays. Thus, the moderator should have a high scattering cross section at desired energies (thermal or epithermal energy), low one for undesired energies (thermal or fast energy) and absorption cross section, avoiding loss of neutron intensity and producing large quantities of gamma-rays (11). In the optimization process, different moderator materials (Fluental, TiF3, CaF2, Al, AlF3 and MgF2 for epithermal neutron beam BSA, and D2O, normal polyethylene and graphite for thermal neutron BSA, respectively) were considered (1214).

Another important part is reflector which is used to reflect scattered neutrons back into the beam. Reflector should have a low absorption cross section, a high elastic scattering cross section for thermal or epithermal neutrons and also a large mass number in which less loss of energy with elastic collision. We considered Teflon, Pb, 316L stainless steel, BeO, and Al2O3 for thermal neutron beam BSA, and Teflon, Al2O3, Pb for epithermal neutron beam BSA, respectively.

Neutron and Gamma Filters

To minimize the damage to healthy tissue around the tumor, beam filters are necessary for reducing contaminations of fast neutrons, thermal neutrons and gamma rays. For the thermal neutron beam BSA, Pb, and Bi were compared as gamma filter. As for the epithermal neutron beam BSA, Ni, and 7LiF were used as fast neutron filter and thermal neutron filter, respectively.

Collimator

The collimator can limit divergence of the neutron beam and, reduce undesired irradiation and focus neutrons to patient position. And we calculated J/Φ to measure the beam divergence variation. A high ratio means that the neutron beam is close to the beam port and change slightly with distance from the port. A target value for this ratio should be > 0.7 (9).

Clinical Parameters

Under clinical conditions, it is vital to investigate the dosimetry performance in the patients. So, in-phantom parameters were calculated. These parameters are advantage depth (AD), advantage ratio (AR), AD dose rate (ADDR), and treatment time (TT), where AD is the depth in phantom at which the total therapeutic dose in tumor equals the maximum dose of the normal tissue. AD indicates the depth of effective beam penetration. The AR is the ratio of the total therapeutic dose in tumor to the total normal tissue dose over a given depth (usually from the surface to AD). It is a measure of the therapeutic gain. ADDR is defined as the maximum dose rate for the normal tissue (15, 16).

Four components contributing to the absorbed dose in BNCT as follows:

(1) The gamma dose Dγ: the dose due to gamma rays in the neutron beam as well as gamma rays induced in the tissue from nuclear reactions, like 1H (n,γ) 2H reaction;

(2) The hydrogen dose DH: the dose due to recoil protons from 1H (n,n′) 1H reaction;

(3) The nitrogen dose DN: the dose due to energetic proton and the recoiling 14C nucleus from 14N (n,p)14C reaction;

(4) The boron dose DB: the dose due to α particles and the recoiling 7Li nuclei from 10B (n,α) 7Li reaction.

The total RBE-weighted dose, DT is expressed in the unit of RBE-Gy, as a sum of physical dose components multiplied by appropriate weighting-factors (RBE or CBE) for each dose component. It can be calculated using Equation (2) as below:

DT=CB×DB+ωN×DN+ωH×DH+ωγ×Dγ    (2)

where ωγ, ωH, ωN, and CB are the weighting factors for gamma rays, hydrogen, nitrogen and boron, respectively. The values of ωH and ωN were taken as 3.0, ωγ was considered as 1, while ωB was 1.35 for boron in the normal tissue and 3.8 for boron in the tumor. A simple phantom was considered to be a cylinder with a simplified composition of soft tissue as shown in Figure 3. The elemental compositions for the material of the cylinder phantom were also listed in Figure 3 (17). 10B was added to the phantom directly at the tumor concentration, as BT, of 30 ppm. The normal tissue concentration, BN, was chosen as 9 ppm, so that the ratio of BT to BN, or T/N was 3.33 (18).

FIGURE 3
www.frontiersin.org

Figure 3. Geometry and composition of the phantom used in dose calculation.

Results and Discussion

Optimization Design of the Thermal Neutron Beam BSA

Moderator for the Thermal Neutron Beam BSA

The thickness of the material varied from 10 to 50 cm with a step of 2 cm in the calculations, and the radius of the moderator was set to be 22 cm. The results are shown in Figure 4. D2O gave the highest thermal neutron flux (Φth), and when the thickness exceeded 30 cm, it gave low epithermal and fast neutron component (Depi−fastth), low γ ray component (Dγth), and high proportion of thermal neutron (Φthtotal). Therefore, D2O was chosen as the moderator for the thermal neutron beam BSA with a thickness of 40 cm.

FIGURE 4
www.frontiersin.org

Figure 4. Thermal neutron parameters of moderator, (A–D) are thermal neutron flux, thermal neutron ratio, epithermal, and fast neutron component, Gamma component, respectively.

Reflector for the Thermal Neutron Beam BSA

As the next step of the optimization, BeO, Al2O3, 316L stainless steel, Teflon and Pb were compared for the reflector. The radius of the reflector was changed with a range of 25–75 cm and a step length of 5 cm. The results are shown in Figure 5. First of all, among the five kinds of materials, BeO presented the highest thermal neutron flux (Φth), the lowest epithermal and fast neutron component (Depi−fastth) and the highest proportion of thermal neutron (Φthtotal). The γ ray component is also lower than 316L Stainless steel, Al2O3 and Teflon, so BeO was selected as the reflector material. In addition, it shows that when the radius is > 45 cm, the reflector radius has no significant influence on the neutron beam parameters. Considering the purpose of saving material and reducing the size of BSA, the reflector radius is chosen as 50 cm.

FIGURE 5
www.frontiersin.org

Figure 5. Thermal neutron parameters of reflector, (A–D) are thermal neutron flux, thermal neutron ratio, epithermal, and fast neutron component, Gamma component, respectively.

Gamma Filter for Thermal Neutron Beam BSA

After determining the materials and sizes of the moderator and reflector, it is necessary to optimize the design of gamma filter to reduce γ rays' component because that γ rays which produced during moderation process cause unnecessary dose to normal tissue. The commonly used gamma filter materials are Pb and Bi. We made a comparison between these two materials. The thickness of gamma filter varied from 5–10 cm and the step length was 1 cm. Results as shown in Figure 6, Bi is a better choice for BNCT because it provides high photon elimination and low loss of thermal neutron flux. Therefore, we choose 9 cm Bi as gamma filter.

FIGURE 6
www.frontiersin.org

Figure 6. Thermal neutron parameters of gamma filter, (A–D) are thermal neutron flux, thermal neutron ratio, epithermal, and fast neutron component, Gamma component, respectively.

Collimator

We choose Bi as material of collimator. And it is designed as a cone. It is 16 cm long (in height), and radius of bottoms are 22 and 7 cm. The J/Φ is 0.654 which is basically meets the target value. The details of collimator are shown in Figure 7A.

FIGURE 7
www.frontiersin.org

Figure 7. Structures of proposed epithermal and thermal neutron beam BSAs, (A) is thermal neutron BSA, and (B) is epithermal neutron BSA.

In conclusion, the final design of thermal neutron beam BSA can be obtained by using D2O as moderator, BeO as reflector, and Bi as gamma filter. The structure of thermal neutron beam BSA is shown in Figure 7A, and the neutron beam parameters were listed in Table 2 with other designs in the world. The flux of proposed thermal neutron beam BSA is higher than INFN and not much less than IHNI. However, IHNI is based on reactor, so it's difficult to build in hospital. So, the proposed thermal neutron beam BSA has its advantages.

TABLE 2
www.frontiersin.org

Table 2. Thermal neutron beam parameters of proposed design and other designs in the world.

Optimization Design of Epithermal Neutron Beam BSA

The process of the optimization design is basically the same as that of the thermal neutron beam BSA, so we show the final optimization results directly. We choose 45 cm MgF2 as moderator, 50 cm Pb as reflector, 0.4 mm 6LiF as the thermal neutron filter. The collimator is as same as thermal neutron BSA except its material is Pb. The structure of epithermal neutron beam BSA is shown in Figure 7B, and the parameters are listed in Table 3 with other designs in the world.

TABLE 3
www.frontiersin.org

Table 3. Epithermal neutron beam parameters of proposed design and other designs in the world.

It is obvious that all the parameters of proposed BSAs fulfill the IAEA recommended values. And the flux of proposed epithermal neutron beam BSA is higher than other designs.

Calculation of Clinical Parameters

Firstly, we calculated neutron flux distribution in the phantom of thermal neutron beam and epithermal neutron beam generated by proposed BSAs, respectively. As shown in Figures 8A,B, the maximum depth of flux in the phantom of thermal neuron beam is 5 cm, and the maximum depths of thermal and epithermal flux of epithermal neutron beam are 12 and 8 cm, respectively. The components of epithermal and fast neutron are too low so that they are not shown in the Figure 8A.

FIGURE 8
www.frontiersin.org

Figure 8. Neutron flux distribution in the phantom, (A) is thermal neutron beam, and (B) is epithermal neutron beam.

As for the dose in tumor and normal tissue, we also calculated them in the phantom. As shown in Figure 9, the AD is 7.93 cm for epithermal neutron beam and 3.52 cm for thermal neutron beam. The ADDR for thermal neutron beam is 0.476RBE-Gy/min, and for epithermal neutron beam it is 0.208RBE-Gy/min. According to the results, we found that the large amount of dose was delivered to skin and superficial normal tissue. If we defined the time that dose delivered to normal tissue exceed the maximum tolerated dose (12.5RBE-Gy) as the treatment time (TT) which is the maximum value. Then, TTs are 60.1 min for epithermal neutron beam and 26.3 min for thermal neutron beam. During this time, the maximum DTs of tumor are 51.87 RBE-Gy for epithermal neutron beam and 65.75 RBE-Gy for epithermal neutron beam. And maximum DRs are 0.863RBE-Gy/min for epithermal neutron beam and 2.5 RBE-Gy/min for thermal neutron beam, respectively.

FIGURE 9
www.frontiersin.org

Figure 9. Dose distribution in the phantom, (A) is thermal neutron beam, and (B) is epithermal neutron beam.

Conclusions

In this paper, MCNPX is used to design the BSAs of an AB-BNCT using 7Li(p,n)7Be reaction. The optimized BSAs for thermal and epithermal neutron beams can ensure that the beam parameters at the exit all meet the recommended values of IAEA. And the clinical parameters are also calculated so that it can give a reference for clinical condition. So the final optimal design of thermal and epithermal BSA can be an important reference for the BSA engineering scheme of multi-terminal AB-BNCT device. In the follow-up work, the accelerator-based multi-terminal BSA conversion device will be further designed to facilitate the flexible switching and replacement of BSA, ensure the safe operation of the device, and give better play to the advantages of multi-terminal devices, which will promote the development of AB-BNCT in China.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

GL adapted the neutron source in the Li target (with WJ), performed the simulations on BSA, calculated clinical parameters, and wrote the manuscript. WJ and LZ helped execute MCNP and adjust the input files. WC participated in the study design and reviewed the results. QL supervised this work and contributed to the study design. All authors contributed to the article and approved the submitted version.

Funding

This study was funded by National Key R&D Program of China (2018YFC0115700), National Natural Science Foundation of China (11875299), Major science and technology project of Gansu Province (1602FKDA005), The key research program of the Chinese Academy of Sciences (KFZD-SW-222), West Light Foundation of The Chinese Academy of Sciences (29Y86205), and the project of Youth Innovation Promotion Association, CAS (2017454).

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.

References

1. Barth RF, Soloway AH. Boron neutron capture therapy of primary and metastatic brain tumors. Mol Chem Neuropathol. (1994) 21:139–54. doi: 10.1007/BF02815348

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Locher GL. Biological effects and therapeutic possibilities of neutrons. Am J Roentgenol Radium Ther. (1936) 36:1–13.

Google Scholar

3. Sweet WH. The uses of nuclear disintegration in the diagnosis and treatment of brain tumor. N Engl J Med. (1951) 245:875–8. doi: 10.1056/NEJM195112062452301

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Hatanaka H, Nakagawa Y. Clinical results of long-surviving brain tumor patients who underwent boron neutron capture therapy. Int J Radiat Oncol. (1994) 28:1061–6. doi: 10.1016/0360-3016(94)90479-0

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Kumada H, Naito F, Hasegawa K, Kobayashi H, Kurihara T, Takada K, et al. Development of LINAC-based neutron source for boron neutron capture therapy in University of Tsukuba. Plasma Fusion Res. (2018) 13:2406006. doi: 10.1585/pfr.13.2406006

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Uritani A, Menjo Y, Watanabe K, Yamazaki A, Kiyanagi Y, Tsuchida K. Design of beam shaping assembly for an accelerator-driven BNCT system in Nagoya University. J Phys Soc Jpn. (2018) 011002:1–7. doi: 10.7566/JPSCP.22.011002

CrossRef Full Text | Google Scholar

7. Sato K, Uritani A, Watanabe K, Yoshihashi S, Yamazaki A, Kiyanagi Y, et al. Improved design of the exit of a beam shaping assembly for an accelerator-driven BNCT system in Nagoya University. J Phys Soc Jpn. (2018) 011003:1–7. doi: 10.7566/JPSCP.22.011003

CrossRef Full Text | Google Scholar

8. Tanaka H, Sakurai Y, Suzuki M, Masunaga S, Mitsumoto T, Fujita K, et al. Experimental verification of beam characteristics for cyclotron-based epithermal neutron source (C-BENS). Appl Radiat Isot. (2011) 69:1642–5. doi: 10.1016/j.apradiso.2011.03.020

PubMed Abstract | CrossRef Full Text | Google Scholar

9. IAEA. Current Status of Neutron Capture Therapy. Vienna (2001).

10. Chadwick MB, Herman M, Obložinský P, Dunn ME, Danon Y, Kahler AC, et al. ENDF/B-VII.1 nuclear data for science and technology: cross sections, covariances, fission product yields and decay data. Nucl Data Sheets. (2011) 112:2887–996. doi: 10.1016/j.nds.2011.11.002

CrossRef Full Text | Google Scholar

11. Zaidi L, Belgaid M, Taskaev S, Khelifi R. Beam shaping assembly design of 7Li(p,n)7Be neutron source for boron neutron capture therapy of deep-seated tumor. Appl Radiat Isot. (2018) 139:316–24. doi: 10.1016/j.apradiso.2018.05.029

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Ke G, Sun Z, Shen F, Liu T, Li Y, Zhou Y. The study of physics and thermal characteristics for in-hospital neutron irradiator (IHNI). Appl Radiat Isot. (2009) 67:S234–7. doi: 10.1016/j.apradiso.2009.03.117

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Sakurai Y, Kobayashi T. Characteristics of the KUR heavy water neutron irradiation facility as a neutron irradiation field with variable energy spectra. Nucl Instruments Methods Phys Res Sect A Accel Spectrometers Detect Assoc Equip. (2000) 453:569–96. doi: 10.1016/S0168-9002(00)00465-4

CrossRef Full Text | Google Scholar

14. Ceballos C, Esposito J, Agosteo S, Colautti P, Conte V, Moro D, et al. Towards the final BSA modeling for the accelerator-driven BNCT facility at INFN LNL. Appl Radiat Isot. (2011) 69:1660–3. doi: 10.1016/j.apradiso.2011.01.032

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Rasouli FS, Masoudi FS, Kasesaz Y. Design of a model for BSA to meet free beam parameters for BNCT based on multiplier system for D-T neutron source. Ann Nucl Energy. (2012) 39:18–25. doi: 10.1016/j.anucene.2011.08.025

CrossRef Full Text | Google Scholar

16. Rasouli FS, Masoudi SF. Simulation of the BNCT of brain tumors using MCNP code: Beam designing and dose evaluation. Iran J Med Phys. (2012) 9:183–92. doi: 10.22038/ijmp.2012.150

CrossRef Full Text | Google Scholar

17. McConn RJ, Gesh CJ, Pagh RT, Rucker RA, Williams R III. Compendium of Material Composition Data for Radiation Transport Modeling. Richland, WA (2011) doi: 10.2172/1023125

CrossRef Full Text | Google Scholar

18. Seki R, Wakisaka Y, Morimoto N, Takashina M, Koizumi M, Toki H, et al. Physics of epi-thermal boron neutron capture therapy (epi-thermal BNCT). Radiol Phys Technol. (2017) 10:387–408. doi: 10.1007/s12194-017-0430-5

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Montagnini B, Cerullo N, Esposito J, Giusti V, Mattioda F, Varone R. Spectrum shaping of accelerator-based neutron beams for BNCT. Nucl Instruments Methods Phys Res Sect A Accel Spectrometers Detect Assoc Equip. (2002) 476:90–8. doi: 10.1016/S0168-9002(01)01398-5

CrossRef Full Text | Google Scholar

20. Kim KO, Kim JK, Kim SY. Optimized therapeutic neutron beam for accelerator-based BNCT by analyzing the neutron angular distribution from 7Li(p,n)7Be reaction. Appl Radiat Isot. (2009) 67:1173–9. doi: 10.1016/j.apradiso.2009.02.004

CrossRef Full Text | Google Scholar

21. Fantidis JG. Beam shaping assembly study for BNCT facility based on a 2.5 MeV proton accelerator on Li target. J Theor Appl Phys. (2018) 12:249–56. doi: 10.1007/s40094-018-0312-1

CrossRef Full Text | Google Scholar

Keywords: accelerator-based BNCT, beam shaping assembly, thermal neutron, epithermal neutron, multi-terminal

Citation: Li G, Jiang W, Zhang L, Chen W and Li Q (2021) Design of Beam Shaping Assemblies for Accelerator-Based BNCT With Multi-Terminals. Front. Public Health 9:642561. doi: 10.3389/fpubh.2021.642561

Received: 16 December 2020; Accepted: 08 February 2021;
Published: 11 March 2021.

Edited by:

Lu Cai, University of Louisville, United States

Reviewed by:

Xingcai Guan, Lanzhou University, China
Jacob Fantidis, International Hellenic University, Greece

Copyright © 2021 Li, Jiang, Zhang, Chen 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: Qiang Li, liqiang@impcas.ac.cn

Disclaimer: 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.