Because of the great importance in nuclear fuel cycle and environmental decontamination, actinide chemistry has gained increasing attentions from chemists, geologists and materials specialists[
Journal of Inorganic Materials, Volume. 35, Issue 2, 243(2020)
In this work, we report a novel octa-nuclear uranyl (U8) motif [(UO2)8O4(μ3-OH)2(μ2-OH)2] 4+ embedded in a uranyl-oxalate coordination polymer (compound 1) based on a U-shaped linker with extra-long xylylene chains for stabilizing the resulting high-nuclear motif through additional cross-linking connectivity. A comparison with dimeric and monomeric uranyl compounds obtained at different pH value from the same hydrothermal system reveals that, solution pH plays a vital role in formation of this octa-nuclear uranyl motif by promoting hydrolysis of uranyl source. Since high similarity of eight uranium centers in this nearly planar U8 motif here, overlapping and broadening of signals in fluorescence, infra-red (IR) and Raman spectra can be found.
Because of the great importance in nuclear fuel cycle and environmental decontamination, actinide chemistry has gained increasing attentions from chemists, geologists and materials specialists[
While many oligomeric uranyl units as well as infinite chain-like polynuclear uranyl motifs have been characterized in uranyl-organic compounds, few cases of uranyl-organic compounds with oligomeric octa-nuclear uranyl motifs[
Herein, a uranyl-oxalate coordination polymer containing a new type of oligomeric octa-nuclear uranyl motif, [(UO2)8O4-(μ3-OH)2(μ2-OH)2]4+, was isolated successfully and identified with the aid of a flexible U-shaped dicarboxylate linker (m-Xyl-BPy4CA) with an extra-long xylylene chain.
1 Experimental
1.1 Synthesis
(UO2)8O4(μ3-OH)2(μ2-OH)2(C2O4)2(C20H16N2O4)2 (1) UO2(NO3)2·6H2O (200 μL, 0.10 mmol), [m-Xyl- BPy4CEt]Br2 (28.0 mg, 0.05 mmol), NaOH (3.2 mg, 0.08 mmol), ultrapure water (1.0 mL) was loaded into a 15 mL autoclave (the initial pH is 2.72), sealed and heated to 150 ℃ in an oven for 2 d, then automatically cooled to ambient temperature. Dark yellow crystals of 1 (Fig. S1) were produced, filtered off, and rinsed with ultrapure water and subjected to air-drying at room temperature to give 1 in pure phase in spite of preferred alignment of some crystal indices (Fig. S2). Yield: 14 mg (35% based on uranium).
[(UO2)2(μ2-OH)2(C20H16N2O4)2(H2O)2](NO3)2 (2) compound 2 was synthesized using the same protocol except a reducing amount of NaOH (2.4 mg, 0.06 mmol). Prism yellow-green crystals were filtered off, and rinsed with ultrapure water and subjected to air-drying at room temperature to give 2 in pure phase (Fig. S1 and S3). Yield: 22 mg (31% based on uranium).
(UO2)(C2O4)(C6H5NO2) (3) compound 3 was synthesized using the same protocol like 1 and 2 but without NaOH added. Stick-like yellow-green crystals of 3 can be obtained as accompanied by considerable amount of 2 and other unidentified impurities (Fig. S1 and S4).
1.2 Characterization
1H-NMR spectra were recorded on a Bruker AVANCE III (500 MHz, Bruker, Switzerland) with deuterium oxide (D2O) as solvent. ESI-MS spectra were obtained with a Bruker AmaZon SL ion trap mass spectrometer (Bruker, USA). Powder X-ray diffraction (PXRD) measurements were recorded on a Bruker D8 Advance diffractometer with Cu Kα radiation (λ=0.15406 nm) in the range 5°-60° (step size: 0.02°). Thermogravimetric analysis (TGA) was performed on a TA Q500 analyzer over the temperature range of 25-800 ℃ in air with a heating rate of 10 ℃/min. Solid-state fluorescence spectra were measured on a Hitachi F-4600 fluorescence spectrophotometer. The Fourier transform infrared (IR) spectra were recorded from KBr pellets in the range of 4000-400 cm-1 at Bruker Tensor 27 spectrometer. Raman spectra were recorded on a Thermo Fisher Scientific DXRxi Micro Raman imaging spectrometer excited at 780 nm.
Uranyl-oxalate coordination polymer (namely compound 1, Fig. 1(a)) containing a new type of oligomeric octa- nuclear uranyl motif, [(UO2)8O4-(μ3-OH)2(μ2-OH)2]4+ (Fig. 1(b)), was isolated successfully and identified with the aid of a flexible U-shaped dicarboxylate linker (m-Xyl-BPy4CA) with an extra-long xylylene chain (Fig. 1(c-d)). By a direct comparison with another uranyl compounds (namely compound 2 and compound 3) with monomeric or dimeric uranyl nodes from the same organic ligand, factors affecting the formation of compound 1 such as aqueous pH and hydrothermal stability of organic ligand have been discussed. The physicochemical properties of compound 1 were also characterized and compared in detail.
Figure .Octa-nuclear uranyl-oxalate network reinforced by U-shaped zwitterionic dicarboxylate linkers
2 Results and discussion
2.1 Structural description
Crystallographic analysis shows that compound 1 crystallizes in the monoclinic crystal system with a space group of P21/c, and its asymmetric unit contains four unique uranium centers, a linker of m-Xyl- BPy4CA ([C20H16N2O4]) and an in situ-formed oxalate ligand ([C2O4]2-) (Fig. 2(a)). The formation of oxalate ion is more or less a little surprising in spite of several similar cases found in other heterocycle-based or labile organic ligands[
Figure .Crystal structure of compound 1
Each U8 motif is eight-connected with four oxalate and four m-Xyl-BPy4CA moieties surrounded (Fig S6(a)). The oxalate ligands, always go together with a U-shaped bidentate m-Xyl-BPy4CA, promote the extension of U8 motif from four directions through connecting four adjacent ones in a bridging mode (one side is in μ2-η2 : η1 mode and the other side is in η2 mode) (Fig. S6(b)), and subsequently lead to a 2D network with the minimum rhombic loop in size of 1.193 nm×1.077 nm (Fig. S6(c-d)). Detailed analysis reveals that, every U8 motif displays a different overall orientation from that of its adjacent U8 with an angle of inclination of 36.6(4)°, which is contributed to the distortion of the above rhombic loop (Fig. S7).
When the amount of sodium hydroxide added to the hydrothermal system of uranyl and [m-Xyl- BPy4CEt]Br2 decreased gradually, another two uranyl compounds 2 and 3 emerge. Compound 2 is a molecular compound with a dimeric uranyl motif containing two equivalent uranyl centers (Fig. 3(a-b)). Each uranyl is in a pentagonal bipyramid geometry with equatorial U-O bond distances in the range of 0.2325(2)-0.2576(2) nm (Table S1). Besides two bridging hydroxyl groups and two water molecules, the uranyl dimer is coordinated by two bidentate η1-mode m-Xyl-BPy4CA linkers both in a head-to-tail way, and finally gives a pattern of twinned double loops (Fig. 3(b)). The simple double loops interact with two nitrate anions by hydrogen bonds (Fig. S8 and Table S2) and stack with each other to form the final crystal lattice of molecular compound (Fig. 3(c-d)).
Figure .Crystal structure of compound 2
Another uranyl compound 3 with a one-dimensional (1D) chain-like topology was produced under the more acidic condition. As the case of compound 1, the in situ-formed oxalate ligand ([C2O4]2-) is present in compound 3, indicating the frequent occurrence of oxalate ion for the hydrothermal systems of pyridine derivatives[
Figure .Crystal structure of compound 3
2.2 pH-dependent regulation on the formation of U8 motif
The distinct structural difference among compounds 1, 2 and 3 suggests the significant impact of pH on the uranyl hydrothermal systems, especially the uranyl speciation and crystallization process (Fig. 5). The strong dependence on pH of uranyl speciation is largely attributed to hydrolysis of uranyl ion in aqueous solution. Generally speaking, uranyl monomer is likely to be stable at low pH, whereas high pH of aqueous solution always induces uranyl hydrolysis and promotes its oligomerization through the olation or oxolation process. With the aid of certain organic ligands, these oligomeric uranyl species could further transfer to the solid state as poly-nuclear uranyl compounds. Resembling most uranyl systems in aqueous solution, the formation of U8 motif in compound 1 also follows a pH-regulated mechanism. As revealed by the synthesis of compound 1-3, the aqueous pH used for compound 1 (pH ~2.72) with octa-nuclear uranyl motifs is higher than that for compounds 2 (pH ~2.29) and 3 (pH ~1.35) with only dimeric or monomeric uranyl motifs. Interestingly, pH also exerts a non-negligible influence on the stabilization of organic ligands under hydrothermal conditions. In contrast to the large degree of decomposition of [m-Xyl-BPy4CEt]Br2 to simple isonicotinate linkers found in compound 3, the original zwitterionic dicarboxylate linkers still retains its main molecular skeleton with controlled hydrolysis to its acid form of m-Xyl-BPy4CA during hydrothermal reaction with uranyl under less acidic conditions for the synthesis of compound 1. It can be speculated that the preservation of total skeleton of m-Xyl-BPy4CA with strong coordination capability toward uranyl units as demonstrated in compound 1 also contributes to the formation of U8 motifs by affording additional stabilization through coordination and cross- linkage. The role of U-shaped m-Xyl-BPy4CA linkers in the stabilization of high-nuclear U8 motifs will be discussed as followed.
Figure .pH-dependent regulation of hydrothermal reactions of
2.3 Structural stabilization of high-nuclear U8 motifs
It is interesting to find that most of uranyl-oxalate compounds[
Moreover, a comprehensive survey on the relationship between high-nuclear uranyl motifs and organic ligands used reveals that, multi-topic organic ligand with a nonlinear configuration is more likely to fit for capture, fixation and stabilization of high-nuclear uranyl motifs by coordination to uranyl units in a compact way (Fig. S11). Another important feature of m-Xyl-BPy4CA linker in molecular structure is its flexible conformation (Fig. S12). This U-shaped ligand could display varying molecular conformations in different compounds so as to be adapted to different coordination environment of uranyl motif. For example, the m-Xyl-BPy4CA linker gives a more open conformation when located at the cavity of uranyl-oxalate network in compound 1 (angle between two pyridium groups is ~47.9°, left in Fig. S12), while it becomes more compact with two dangled pyridium groups nearly paralleled to each other (angle between two pyridium groups is ~13.8°, right in Fig. S12). As shown above, the flexibility of molecular skeleton for m-Xyl-BPy4CA linker ensures it structural adaptivity and facilitates the formation of new U8 motif.
2.4 Physicochemical Properties
Characterizations of properties including thermogravimetric analysis (TGA), Fourier transform infrared (IR) spectra, Raman spectra, and solid-state fluorescence spectra for compounds 1 and 2 in pure phase were conducted. In terms of thermal stability, compound 1 does not undergo thermal decomposition of oxalate group on the backbone of 2D network until the temperature increases up to ~295 ℃ (Fig. S13), which is in sharp contrast to compound 2 (Fig. S14). The good stability of compound 1 is similar to the thermal behavior of uranyl oxalate[
As shown in Fig. S17, solid-state fluorescence spectra of compounds 1 and 2 are more or less different as compared to that of uranyl nitrate (UO2(NO3)2). Compound 2 gives the typical vibronic progression of uranium (VI) with the five main emission bands located at 499, 520, 543, 568 and 596 nm corresponding to the S11 → S00 and S10 → S0ν (ν=0-4) electronic transitions[
3 Conclusions
In summary, we have successfully synthesized a uranyl- oxalate coordination polymer with a new type of octa- nuclear uranyl motif, [(UO2)8O4-(μ3-OH)2(μ2-OH)2]4+ by utilizing U-shaped pyridinium-based organic ligand [m-Xyl-BPy4CEt]Br2 via a pH-dependent regulation. The mixed ligand system of oxalate and m-Xyl-BPy4CA formed in situ plays important roles for stabilizing high-nuclear uranyl species by bridging coordination linkage and superimposed structural connectivity, respectively. The U8-bearing uranyl compound found here provides more information on uranyl hydrolysis and speciation in aqueous solution, and also facilitates the exploration of new actinide materials by employing well- designed organic linkers.
Supporting Materials:
Supporting Materials related to this article can be found at
S1. General Methods
S1. 1 Synthesis
1,1°-(1,3-phenylenebis(methylene))bis(4-(ethoxycarbonyl)pyridin-1-ium) bromide ([m-Xyl-BPy4CEt]Br2). [m-Xyl-BPy4CEt]Br2 was synthesized according to the reported procedure[
S1. 2 X-ray single crystal structure determination
X-ray diffraction data of compound 1 and 2 was performed on Bruker D8 VENTURE X-ray CMOS diffractometer with a Mo Kα X-ray source (λ=0.071073 nm) at 296 K. X-ray diffraction data of compound 3 was collected on a Agilent SuperNova X-ray CCD diffractometer with a Cu Kα X-ray source (λ=0.154184 nm) with higher diffraction capability at 293 K. Standard Agilent Crysalis software was used for the determination of the unit cells and data collection control. All the crystal structures were solved by means of direct methods and refined with full-matrix least squares on SHELXL- 2014. SIMU were used to constrain the displacement parameters of the phenyl and pyridyl groups and ISOR were used to even out the electron density associated with disordered portions of the moieties for both 1 and 3. OMIT were used to eliminate bad reflections obscured by the beamstop for all compounds. Since there is high disorder of pyridyl groups dangling aside, the pyridyl ligand was forced to be half occupied to create a chemically sensible model for 3. Solvent molecules (water) in the structure are highly disordered and impossible to be modelled as discrete atomic sites. To resolve this issue, the contribution of solvent-electron density was removed using the SQUEEZE/PLATON procedure, thereby producing a set of solvent-free diffraction intensities used for improving the structure refinements. The crystal data of both compounds are given in Table S3. Crystallographic data for the structures in this study were deposited with the Cambridge Crystallographic Data Centre as supplementary publication nos. CCDC- 1510791 (1), CCDC-1898268 (2) and CCDC-1510792 (3).
S2. Figures
Figure .Different optical morphologies of 1 with octa-nuclear uranyl (U8) motifs, 2 with binuclear uranyl (U2) motifs and 3 with monomeric uranyl (U1) motifs
Figure .Experimental and simulated patterns of powder X-ray diffraction (PXRD) of compound 1
Figure .Experimental and simulated patterns of powder X-ray diffraction (PXRD) of compound 2
Figure .Experimental and simulated patterns of powder X-ray diffraction (PXRD) of compound 3
Figure .(a) A nearly planar geometry of U8 motif found in this work; (b) a non-planar U8 motif with cation-cation interactions (CCIs) reported by Loiseau,
Figure .(a-b) Eight-connected U8 motif with four oxalate (Ox) and four
Figure .Each U8 motif displays a different overall orientation from that of its adjacent U8 with an angle of inclination of 36.6(4)° (a), resulting in a distortion of the rhombic loop (b)
Figure .Hydrogen bonds between double loops and two nitrate anions
Figure .Two ‘U’-shaped bidentate
Figure .Hydrogen bonds between adjacent layers of 2D sheets through U8 motifs that interact with neighboured
Figure .Some examples of high-nuclear uranyl motif based on nonlinear multi-topic organic ligands, as suggested by the cases of pentanuclear (U5), hexanuclear (U6) and octanuclear (U8) uranyl motifs derived from sulfobenzoate precursors[2],
Figure .Different molecular conformation of
Figure .Thermogravimetric analysis (TGA) of compounds 1, where 1 starts to decompose at ~295 ℃, and finally transforms to U3O8 with residual weight of 69.31% (theoretical value: 70.25%)
Figure .Thermogravimetric analysis (TGA) of compounds 2, where 2 starts to decompose at ~233 ℃, and finally transforms to U3O8 with residual weight of 40.95% (theoretical value: 40.20%)
Figure .Fourier transform infrared (IR) spectra of compounds 1 (U8 motif, blue line), 2 (U2 motif, red line) and 3 (U1 motif, black line) with characteristic symmetric
Figure .The Raman spectra of compounds s 1 (U8 motif) and 3 (U1 motif) with characteristic asymmetric ν3 vibrations (1: 833 and 863 cm-1; 3: 829 and 860 cm-1)
Figure .Solid-state fluorescence spectra of compound 1 and 2 as compared to that of uranyl nitrate (UO2(NO3)2): 1, a broad peak ranging from 530 to 550 nm; 2, five main emission bands located at 499, 520, 543, 568 and 596 nm; UO2(NO3)2, 488, 511, 534, 561 and 589 nm
Figure .1H NMR of [
S3. Tables
Selected bond distances related to uranyl centers in compounds 1, 2 and 3
Selected bond distances related to uranyl centers in compounds 1, 2 and 3
|
Distances and angles for hydrogen bonds observed in compounds 1 and 2
Distances and angles for hydrogen bonds observed in compounds 1 and 2
|
Crystal data and structure refinement for compounds 1, 2 and 3
Crystal data and structure refinement for compounds 1, 2 and 3
|
References:
[1] MIHALCEA I, HENRY N, CLAVIER N, et al. Occurence of an octanuclear motif of uranyl isophthalate with cation-cation interactions through edge-sharing connection mode. Inorganic Chemistry, 2011, 50(13): 6243-6249.
[2] THU RY P. Sulfonate complexes of actinide ions: structural diversity in uranyl complexes with 2-sulfobenzoate. Inorganic Chemistry, 2013, 52(1): 435-447.
[3] ZHENG Y Z, TONG M L, CHEN X M, et al. Synthesis, structure and photoluminescent studies of two novel layered uranium coordination polymers constructed from UO (OH) polyhedra and pyridinedicarboxylates. European Journal of Inorganic Chemistry, 2005, (20): 4109-4117.
[4] THU RY P, NIERLICH M, SOULEY B, et al. Complexation of a hexameric uranium (VI) cluster by p-benzylcalix [7] arene. Journal of the Chemical Society, Dalton Transactions, 1999, (15): 2589- 2594.
[5] SINDELAR V, MOON K, KAIFER A E, et al. Binding selectivity of cucurbit[7]uril: bis(pyridinium)-1,4-xylylene versus 4,4°- bipyridinium guest sites. Organic Letters, 2004, 6(16): 2665- 2668.
[6] HUANG F, SLEBODNICK C, MAHAN E J, et al. [3]Pseudorotaxanes based on the cryptand/monopyridinium salt recognition motif. Tetrahedron, 2007, 63(13): 2875-2881.
[7] MEI L, WANG L, YUAN L Y, et al. Supramolecular inclusion- based molecular integral rigidity: a feasible strategy for controlling the structural connectivity of uranyl polyrotaxane networks. Chemical Communications, 2015, 51(60): 11990-11993.
[1] M ALTMAIER, X GAONA, T FANGHANEL et al. Recent advances in aqueous actinide chemistry and thermodynamics. Chemical Reviews, 113, 901-943(2013).
[2] B JONES M, J GAUNT A. Recent developments in synthesis and structural chemistry of nonaqueous actinide complexes. Chemical Reviews, 113, 1137-1198(2013).
[3] X WANG K, S CHEN J. Extended structures and physicochemical properties of uranyl-organic compounds. Accounts of Chemical Research, 44, 531-540(2011).
[4] B ANDREWS M, L CAHILL C. Uranyl bearing hybrid materials: synthesis, speciation, and solid-state structures. Chemical Reviews, 113, 1121-1136(2013).
[5] T YANG W, G PARKER T, M SUN Z et al. Structural chemistry of uranium phosphonates. Coordination Chemistry Reviews, 303, 86-109(2015).
[6] T LOISEAU, I MIHALCEA, N HENRY et al. The crystal chemistry of uranium carboxylates. Coordination Chemistry Reviews, 266, 69-109(2014).
[7] D RAI, R FELMY A, L RYAN J et al. Uranium (IV) hydrolysis constants and solubility product of UO2·
[8] S AHRLAND. On the complex chemistry of the uranyl ion Ι. The hydrolysis of the 6-valent uranium in aqueous solutions. Acta Chemica Scandinavica, 3, 374-400(1949).
[9] P ZANONATO, P DI BERNARDO, A BISMONDO et al. Hydrolysis of uranium (VI) at variable temperatures (10-85 ℃). Journal of the American Chemical Society, 126, 5515-5522(2004).
[10] L SALMON, P THUERY, M EPHRITIKHINE et al. Crystal structure of the first octanuclear uranium (IV) complex with compartmental schiff base ligands. Polyhedron, 23, 623-627(2004).
[11] I MIHALCEA, N HENRY, N CLAVIER et al. Occurence of an octanuclear motif of uranyl isophthalate with cation-cation interactions through edge-sharing connection mode. Inorganic Chemistry, 50, 6243-6249(2011).
[12] S PASQUALE, S SATTIN, E C ESCUDERO-ADAN et al. Giant regular polyhedra from calixarene carboxylates and uranyl. Nature Communications, 3, 785(2012).
[13] P THUERY. A highly adjustable coordination system: nanotubular and molecular cage species in uranyl ion complexes with kemp's triacid. Crystal Growth & Design, 14, 901-904(2014).
[14] H WANG L, R SHANG, Z ZHENG et al. Two systems of [DabcoH2]2+/[PipH2]2+-uranyl-oxalate showing reversible crystal-to- crystal transformations controlled by the diammonium/uranyl/oxalate ratios in aqueous solutions ([DabcoH2]2+=1,4-diazabicyclo- [2.2.2]-octaneH2 and [PipH2]2+ = PiperazineH2). Crystal Growth & Design, 13, 2597-2606(2013).
[15] B CHAPELET-ARAB, G NOWOGROCKI, E ABRAHAM et al. Crystal structure of new uranyl oxalates (NH4)2[UO2(C2O4)·2H2O] and (NH4)2-
[16] A GIESTING P, J PORTER N, C BURNS P et al. A series of sheet-structured alkali metal uranyl oxalate hydrates: structures and IR spectra. Zeitschrift für Kristallographie, 221, 589-599(2006).
[17] A GIESTING P, J PORTER N, C BURNS P et al. Uranyl oxalate hydrates: structures and IR spectra. Zeitschrift für Kristallographie, 221, 252-259(2006).
[18] L DUVIEUBOURG, G NOWOGROCKI, F ABRAHAM et al. Hydrothermal synthesis and crystal structures of new uranyl oxalate hydroxides:
[19] P THUERY. Reaction of uranyl nitrate with carboxylic diacids under hydrothermal conditions. Crystal structure of complexes with L(+)-tartaric and oxalic acids. Polyhedron, 26, 101-106(2007).
[20] V VOLOGZHANINA A, B SEREZHKINA L, A NEKLYUDOVA N et al. Synthesis and characterisation of a trinuclear uranyl complex: crystal structure of (CN3H6)5[(UO2)3O(OH)2(CH3COO)(C2O4)3]. Inorganica Chimica Acta, 362, 4921-4925(2009).
[21] A CHUGH C, A SHARMA, A SHARMA et al. Kinetics and mechanism of thermal decomposition of uranyl oxalate. Asian Journal of Chemistry, 23, 1865-1866(2011).
[22] R BARTLETT J, P COONEY R et al. On the determination of uranium oxygen bond lengths in dioxouranium (VI) compounds by raman-spectroscopy. Journal of Molecular Structure, 193, 295-300(1989).
[23] A BRACHMANN, G GEIPEL, G BERNHARD et al. Study of uranyl (VI) malonate complexation by time resolved laser-induced fluorescence spectroscopy (TRLFS). Radiochimica Acta, 90, 147-153(2002).
[24] L MEI, Z WANG C, Z ZHU L et al. Exploring new assembly modes of uranyl terephthalate: templated syntheses and structural regulation of a series of rare 2d→3d polycatenated frameworks. Inorganic Chemistry, 56, 7694-7706(2017).
[25] S NATRAJAN L. Developments in the photophysics and photochemistry of actinide ions and their coordination compounds. Coordination Chemistry Reviews, 256, 1583-1603(2012).
[26] P THUERY, J HARROWFIELD. Solvent effects in solvo-hydrothermal synthesis of uranyl ion complexes with 1,3-adamantanediacetate. CrystEngComm, 17, 4006-4018(2015).
[27] P THUERY, J HARROWFIELD. Structural variations in the uranyl/4,4'-biphenyldicarboxylate system. rare examples of 2d→3d polycatenated uranyl-organic networks. Inorganic Chemistry, 54, 8093-8102(2015).
[28] P THUERY, E RIVIERE, J HARROWFIELD et al. Uranyl and uranyl-3d block cation complexes with 1,3-adamantanedicarboxylate: crystal structures, luminescence, and magnetic properties. Inorganic Chemistry, 54, 2838-2850(2015).
Get Citation
Copy Citation Text
Si WU, Lei MEI, Kong-Qiu HU, Zhi-Fang CHAI, Chang-Ming NIE, Wei-Qun SHI.
Category: RESEARCH LETTERS
Received: Mar. 21, 2019
Accepted: --
Published Online: Jan. 27, 2021
The Author Email: NIE Chang-Ming (niecm196132@163.com), SHI Wei-Qun (shiwq@ihep.ac.cn)