Effective electrocatalytic hydrodechlorination of 2,4,6-trichlorophenol by a novel Pd/MnO2/Ni foam cathode

2022-09-16 05:24ZiMengZhngRuiChengJunNnXueQiChenCongHungDiCoCiHuBiJingLongHnBinLingZhiLingLiAiJieWng
Chinese Chemical Letters 2022年8期

Zi-Meng Zhng, Rui Cheng, Jun Nn, Xue-Qi Chen, Cong Hung, Di Co,Ci-Hu Bi, Jing-Long Hn, Bin Ling, Zhi-Ling Li,∗, Ai-Jie Wng,

a State Key Laboratory of Urban Water Resource and Environment, School of Environment, Harbin Institute of Technology, Harbin 150090, China

b National Technology Innovation Center of Synthetic Biology, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin 300308,China

c State Key Laboratory of Urban Water Resource and Environment, School of Civil & Environmental Engineering, Harbin Institute of Technology Shenzhen,Shenzhen 518055, China

ABSTRACT Pd modified electrodes possess problems such as easy agglomeration and low electrolytic ability, and the use of manganese dioxide (MnO2) to facilitate Pd reduction of organic pollutants is just started.However,there is still a limited understanding of how to match the Pd load and MnO2 to realize optimal dechlorination efficiency at minimum cost.Here, a Pd/MnO2/Ni foam cathode was successfully fabricated and applied for the efficient electrochemical dechlorination of 2,4,6-trichlorophenol (2,4,6-TCP).The optimal electrocatalytic hydrodechlorination (ECH) performance with 2,4,6-TCP dechlorination efficiency (92.58%in 180 min) was obtained when the concentration of PdCl2 precipitation was 1 mmol/L, the deposition time of MnO2 was 300 s and cathode potential was −0.8 V.Performance influenced by the exogenous factors (e.g., initial pH and coexisted ions) were further investigated.It was found that the neutral pH was the most favorable for ECH and a reduction in dechlorination efficiency (6%∼47.6%) was observed in presence of 5 mmol/L of NO−2, NO−3, S2−or SO2−3 .Cyclic voltammetry (CV) and quenching experiments verified the existence of three hydrogen species on Pd surface, including adsorbed atomic hydrogen (H∗ads), absorbed atomic hydrogen (H∗abs), and molecular hydrogen (H2).And the introduction of MnO2 promoted the generation of atomic H∗.Only adsorbed atomic hydrogen (H∗ads) was confirmed that it truly facilitated the ECH process.Besides H∗ads induced reduction, the direct reduction by cathode electrons also participated in the 2,4,6-TCP dechlorination process.Pd/MnO2/Ni foam cathode shows excellent dechlorination performance, fine stability and recyclable potential, which provides strategies for the effective degradation of persistent halogenated organic pollutants in groundwater.

Keywords:Pd/MnO2/Ni foam cathode 2,4,6-Trichlorophenol Electrocatalytic hydrodechlorination Dechlorination pathway Atomic H∗generation and utilization

2,4,6-Trichlorophenol (2,4,6-TCP) constitutes an important member of chlorinated organic compounds (COCs).2,4,6-TCP has been widely applied as a pesticide precursor or formed during drinking water disinfection [1,2].2,4,6-TCP is toxic, durable, and resistant due to the existence of chlorine atomic [3].Environmental risk of groundwater contamination by 2,4,6-TCP has raised great concerns due to its biotoxicity, persistence, and stability[4,5], urging the U.S.Environmental Protection Agency (EPA) to make it a priority pollutant [6].

Chemical remediation methods, especially chemical reduction,have been proposed as the most widely used method for detoxifying COCs due to their high efficiency [7].Among chemical reduction, electrocatalytic hydrodechlorination (ECH) gains ever-growing interest due to its outstanding performance, mild reaction conditions, and convenient operation conditions [8], which has been regarded as a promising alternative in contamination control of COCs.In general, the mechanism of ECH can be summarized in two pathways: i) Direct reduction by direct electron transfer between pollutants and cathode surface [9]; ii) indirect reduction by atomic H∗produced in process of water electrolysis [10].Various previous studies have demonstrated that the latter way, ECHviaatomic H∗,is the dominant pathway for contaminant removal [11].Details on the generation of atomic H∗in the ECH process are elucidated as following (M represents metal catalyst and adsorbed atomic hydrogen is denoted as H∗ads) (Eqs.1-4) [12]:

Atomic H∗plays a significant role in electrocatalytic hydrodechlorination [13], while hydrogen evolution reaction (HER),generating by both Heyrovsky step and Tafel step, consumes atomic H∗, and inhibits ECH efficiency to some extent [14].Considering this problem, many metals, such as Au [15], Pd [16], and Cu [17], serving as catalysts to modify the cathode surface, have been applied to strengthen the ECH process [18].Among them, Pd has been proved as the most excellent one due to its good performance in atomic H∗generation and storageviaboth H∗adsorption on Pd surface (denoted as H∗ads) and H∗absorption into Pd crystal lattice forming Pd hydride (denoted as H∗abs) [19].However, since Pd is precious, expensive and easy-to-aggregation [20], its commercial application is limited.

To solve the Pd aggregation problem, Ni foam is introduced to improve Pd utilization by dispersing Pd particles due to its desirable 3D reticular structure [21].Its high electrical conductivity accelerates electron transfer inside the Pd layer, which results in good reactivity [22].Meanwhile, considering that Pd is inefficient in water electrolysisviathe Volmer step [23], many transitional metal compounds have been induced to enhance atomic H∗generation/retainment, such as nickel oxides [24], manganese oxides[25] Among them, MnO2with low cost has been proved as the most promising one to provide more atomic H∗viawater dissociation, which is favorable to be stored by Pd [26].However, MnO2shows poor conductivity, and the excessive MnO2load inhibits both ECH and H2generation [23].

Recent studies have reported the introduction of MnO2to Pd/Ni foam cathode would enhance the production of atomic H∗[26,27].An appropriate amount of MnO2deposition is beneficial to enlarge the surface area of electrode and disperse Pd particles, while excessive MnO2can block the fold structure and reduce the conductivity, resulting in the inhabitation of ECH [28].Thus, how to pair the dose of Pd and MnO2to enhance ECH was one of the focuses.Also,cathode potential is pivotal that determines both the amounts of electrons and catalytic activity, and how cathode potential regulation influences ECH performance by Pd and MnO2loaded electrode is largely unknown.Moreover, the influence of exogenous factors in the groundwater environment on catalytic activity, such as the initial pH and coexisted ions, still needs to be further explored.

Herein, a novel electrode was fabricated by introducing MnO2and Pd to the surface of Ni foam, aiming at improving ECH performance and reducing the dose of precious Pd.As a typical COC in groundwater, 2,4,6-TCP was selected as a target contaminant.Electrode operations (Pd load, MnO2load, and cathode potential)were optimized to seek the supervisor conditions for 2,4,6-TCP hydrodechlorination.Meanwhile, 2,4,6-TCP removal efficiency and kinetics, the proposed pathway as well as the impacts of environmental exogenous factors were thoroughly investigated to assess ECH performance by using Pd/MnO2/Ni foam cathode (the detailed experimental content was given in the Supporting information).

Pd/MnO2/Ni foam cathodes with the different Pd and MnO2loads were prepared, and the electrode under the preparation conditions of 1 mmol/L deposition solution concentration and 600 s electrodeposition time was selected for micro-characterization(Figs.S2A-C in Supporting information).SEM showed the morphology of 3D reticular structure at different magnification and found both MnO2and Pd were uniformly attached to the surface of foamed nickel electrode with cracks and grooves, indicating specific surface area was increased after deposition.Meanwhile, the elemental mapping of Pd/MnO2/Ni foam confirmed the existence and uniform distribution of Ni, Pd, Mn and O (Figs.S2D and E in Supporting information).According to the elemental mapping, almost all MnO2were loaded on nickel foam skeleton, and a portion of Pd nano-particles covered on the MnO2layer, while the others dotted edges of the MnO2layer and the surface of Ni foam.

Subsequently, XRD pattern was performed to elucidate the formation of Pd phase (JCPDS No.05–0681) with the characteristic diffraction peaks at 2θ=39.888° and 46.408° for Pd (111) and(200) [29], respectively (Fig.S2F in Supporting information).Moreover, three typical diffraction peaks of Ni (JCPDS No.04–0850)were observed at 2θ=44.266°, 51.618° and 76.157° corresponding to Ni (111), (200) and (220) [30], respectively.However, no significant characteristic diffraction peaks of MnO2were detected due to its slight crystallographic form [31], which might be concealed by others.Another possible explanation came from the fact that the intensities of Mn peaks would be more week after Pd loaded on the MnO2layer, which suggested that XRD could hardly detect the typical diffraction peaks of MnO2.Nevertheless, the presence of Mn was proved by EDS and XPS, indicating the successful introduction of MnO2.

After characterization, electrode operations that influenced ECH performance were investigated.Pd is a desirable catalyst for its excellent performance in atomic H∗retainment, while SEM has shown that superfluous Pd nanoparticles lead to aggregation (Fig.S2 in Supporting information).Consequently, excessive Pd load resulted in a reduction in reactivity and ECH performance [23].The effect of Pd load was investigated to find a chemical solution concentration required for the optimal dose (Fig.1A).With the increase of Pd load, ECH efficiency of 2,4,6-TCP firstly increased from less than 5% to 89.21% and then decreased to 75.82%.The removal efficiency reached the optimal value of 89.21% when the concentration of sediment was 1 mmol/L.The electrocatalytic dechlorination followed pseudo-first-order kinetics.The reaction kinetic constant firstly increased to 11.37 × 10−3min−1, then decreased to 6.99 × 10−3min−1, and the reaction kinetic constant reached the maximum of 11.37 × 10−3min−1at 1 mmol/L (Fig.S4A, Table S1 in Supporting information).Moreover, cyclic voltammetry (CV) curves were conducted to illustrate the relationship between atomic H∗and Pd load at a sweep rate of 1 mV/s (Fig.1B).CV curves of different Pd loads showed two reduction peaks at −0.7 V and −0.2 V,which were probably assigned to H∗absand H∗ads[26].Pd load is positively correlated with atomic H production in a certain range.Moderate Pd load could store atomic H∗, while excessive Pd load would cause aggregation and accelerate the formation of H2bubblesviaa Heyrovsky or Tafel step [32].This could explain the phenomenon for the reduction of dechlorination efficiency and kinetic rate at 2 mmol/L.

MnO2was introduced to expedite the generation of atomic H∗,which could be absorbed into the Pd crystal lattice or adsorbed on the Pd surface.Since the incapability to store atomic H∗of MnO2,no significant dechlorination occurred without Pd load (Fig.1A).ECH efficiency of 2,4,6-TCP firstly increased from 41.49% to 89.21%and then decreased to less than 5% with the increase of electrodeposition time of MnO2(Fig.1C).The removal efficiency of 2,4,6-TCP of 300 s was 2.15 times higher than that of 0 s.Also, the reaction kinetic constantkof 300 s was 4.23 times higher than that of 0 s (Fig.S4B, Table S2 in Supporting information), indicating the improvement in generation of atomic H∗after the load of MnO2.Subsequently, electrochemical measurements including CV and LSV were applied to clarify the enhancement of atomic H∗production.On the one hand, introducing MnO2to Pd/Ni foam made its hydrogen evolution potential positively rise from −0.8 V to −0.6 V (Fig.S5 in Supporting information), verifying the promotion of atomic H∗generationviathe Volmer step.On the other hand, two oxidation peaks were observed at around −0.7 V and −0.2 V (Fig.S5),corresponding to the desorption of atomic H∗adsorbed on (111)and (200) planes of polycrystalline Pd [27].The peak current increased initially, and then reached its maximum (2.35 mA/cm2) at 300 s (Fig.1D), further verifying the promotion of water electrolysis by the addition of MnO2.However, 2,4,6-TCP dechlorination rate and efficiency exhibited a drop trend as the electrodeposition time of MnO2exceeded 300 s (Fig.1C).As previously reported, excessive MnO2would cause low electrical conductivity, which not only limited electron transfer but also inhibited the utilization of atomic H∗[26].This would be the major attribution for the decreased activity.

Fig.1.ECH performance of Pd/MnO2/Ni foam cathode at different (A) Pd loads, (C) MnO2 loads, and (E) cathode potentials in N2-saturated 0.1 mol/L Na2SO4 solution in the presence of 100 μmol/L of 2,4,6-TCP; Cyclic voltammograms recorded for Pd/MnO2/Ni foam cathode under different (B) Pd loads and (D) MnO2 loads at a sweep rate of 1 mV/s; and (F) the relationship between ln kobs and ln(−E) under different cathode potentials.

2,4,6-TCP dechlorination performance was investigated at cathode potentials from −0.4 V to −1.2 V by using the optimized Pd load (1 mmol/L) and MnO2load (300 s) (Fig.1E).Dechlorination efficiency of 2,4,6-TCP increased from 0.51% to 92.58% when the potential decreased from −0.4 V to −0.9 V.Nevertheless, the sharp drop of dechlorination efficiency occurred as the cathode potential continued to fall to −1.2 V.The most effective potential was found to be −0.9 V, resulting in the dechlorination effi-ciency of 92.58%.Furthermore, kinetics analysis also obtained similar trends: the reaction kinetic constant reached the maximum of 14.07±1.25×10−3min−1at −0.9 V (Fig.S4C, Table S3 in Supporting information), showing the optimized ECH performance at−0.9 V.Fig.1F showed the relationship between lnkobsand ln(−E), wherekobsindicated kinetics constant andEreferred to the cathode potential.A linear relationship (R2=0.978) was observed with potential ranging from −0.6 V to −0.9 V, and the apparent rate order was 8.40.The dechlorination efficiency was less than 10%with the potential of −0.4 V and −0.6 V because the non-Faraday current assisted adsorption process was the dominant process on the electrode surface rather than electron transfer.The improvement of lnkobsat the potentials from −0.6 V to −0.9 V contributed to the enhancement of atomic H∗generation, especially reactive H∗ads, at the more negative potentials.However, considering that the hydrogen evolution reaction began to dominate around −0.9 V,the decrease in lnkobswas attributed to the enhanced H2formation at the low potential, which competed with ECH and reduced utilization of atomic H∗.Moreover, it also inhibited the effective electron and mass transfer at the liquid-catalyst interface.

The effects of the initial pH and coexisted ions on ECH performance were performed.The pH of catholyte plays a vital role in dechlorination reaction since the pathway of atomic H∗formation depends on the pH of the catholyte.The reaction equations that relied on pH are demonstrated as follows (Eqs.5-8) [28]:

Fig.2.Effects of (A) the initial pH and (B) the coexisted ions on 2,4,6-TCP removal performance by Pd/MnO2/Ni foam cathode.

The initial pH influenced ECH performance was investigated as shown in Fig.2A.With the initial pH of 3, 5, 7, and 9, the dechlorination efficiency approached 43.61%, 62.23%, 92.59%, and 74.21%, respectively.Although previous studies have confirmed the lower initial pH could provide more H+that favored HER process[26], 2,4,6-TCP is insoluble within the pH range of acidic condition(pKa=−5.99) [33], which possibly causes the mass transfer limitation and decreases 2,4,6-TCP removal efficiency.With the increase of initial pH from 3 to 7, 48.98% higher dechlorination efficiency was achieved.On the one hand, HER was more inhabited by ECH,indicating the effective utilization of atomic H∗went up.On the other, more amount of 2,4,6-TCP was dissolved, which could relieve the limitation of mass transfer.Consequently, a higher dechlorination efficiency was achieved with the increase of initial pH from 3 to 7.However, a decrease in reaction efficiency occurred with the increase of initial pH from 7 to 9, which was probably attributed to the shift in the dominant process of atomic H∗production.Reactions of Eqs.7 and 8 were dominant throughout the electrolysis at the initial pH of 9.Since Eq.7 was much slower than Eq.5, the generation of atomic H∗was too slow to meet the normal need of ECH [34].Above all, the pH of catholyte influenced ECH performance via regulating the competition between ECH and HER, and the optimized dechlorination efficiency was obtained at the condition with initial pH of 7.

Fig.2B demonstrated the effects of coexisted ions (Cl−, S2−,,and NO−3) on ECH performance by Pd/MnO2/Ni foam cathode.The removal rate of 2,4,6-TCP with Cl−was 1.17 times higher than that of control, suggesting that Cl−worked as a redox medium that strengthened the electrochemical reduction process[35].In comparison, the presence of S2−andcompletely inhibited the dechlorination process.Chen [36] confirmed that S2−completely deactivated Pd catalytic activity even at a low concentration of 1.0 mmol/L.The inhibition of Pd catalytic activity by S2−andcould be attributed to the competitive adsorption of HS−, the change of electronic properties of Pd catalyst near sulfur atomic as well as the formation of Pd-S compound on the cathode surface [37].also possessed a significant inhibitory effect on ECH, and the removal efficiency of 2,4,6-TCP decreased to 47.60%.It was probably caused by the competition ofreduction for both electrons and atomic H∗[20].Moreover,showed a slight inhibition of dechlorination efficiency, and about 6% of decreased efficiency was obtained.This was becausecompeted for the available atomic H∗[38].

Electrocatalytic dechlorination mechanism and proposed pathway were determined by various methods.XPS spectra were conducted to reveal the physiochemical morphology of Pd and Mn in the Ni foam electrode before and after the dechlorination reaction(Fig.S6 in Supporting information).As shown in Fig.S6A, there were two pairs of spin-orbit components in both Pd 3d spectrum before and after 2,4,6-TCP dechlorination reaction, suggesting two valence states of Pd0and Pd2+[16].Characteristic peaks located at 337.39 eV and 342.69 eV were assigned to Pd03d5/2and Pd03d3/2, while the binding energy of 339.09 eV and 345.59 eV represented Pd2+3d5/2and Pd2+3d3/2, respectively.The binding energy of Pd03d3/2and Pd2+3d5/2showed a positive shift due to the changes in ligand environment around Pd atomics after the introduction of MnO2[39].Pd2+may come from the fact that oxygen vacancy cannot completely reduce Pd from divalent to zero or oxidation by O2[40].The proportion of Pd in different valence states was changed during the reaction.Of these, the proportion of Pd0reduced from 86.49% to 42.53%, suggesting that partial Pd in electrode surface converted from Pd0to Pd2+.A previous study confirmed that it was possibly caused by the formation of transitional Pd−Cl bond during dechlorination [13].In our trials, excessive NaCl was added to better dissolve PdCl2, resulting in the formation of Pd−Cl.Therefore, Pd2+was likely to be derived from Pd−Cl formed during dechlorination reaction, the active site of which was provided by Pd0.Note that Pd2+could be reduced to Pd0at a cathode potential of −0.85 V, which suggested that this Pd/MnO2/Ni foam cathode was renewable.

Fig.S6B (Supporting information) showed the XPS spectrum of Mn 2P before and after the reaction.Two characteristic peaks located at 642.01 eV and 653.80 eV were assigned to Mn 2p3/2and Mn 2p1/2, suggesting Mn4+acted as the dominating valence state of Mn [41].The binding energy of Mn4+2p3/2was shifted from 642.01 eV to 642.72 eV due to the fluctuation of Mn valence during the electrochemical dechlorination reaction, which was consistent with previous studies [23].The intensities of Mn peaks after reaction were almost the same as the condition before reaction, indicating the stability and reusability of Mn.Furthermore,the stability of Pd/MnO2/Ni foam cathode was evaluated by repeating the dechlorination process at −1.0 V for five cycles (Fig.S7 in Supporting information), which also proved the good stability and reusability.

Generally speaking, the electrocatalytic dechlorination mechanism included direct electron transfer and indirect reduction via atomic H∗, thus it was necessary to clarify their specific contribution during 2,4,6-TCP dechlorination process.t-BuOH was served as a radical scavenger in quenching experiments due to its high reactivity with atomic H∗[42].As shown in Fig.3A, the removal rate of TCP was inhibited from 47.70% to 23.30% when t-BuOH increased from 0 mmol/L to 500 mmol/L within 60 min, demonstrating the significant role of atomic H∗during dechlorination process.The inset of Fig.3A showed kinetic rate constants also decreased drastically with the increasing t-BuOH dose.Only 68.90% of kinetic rate constant without t-BuOH was obtained when t-BuOH was 500 mmol/L.The atomic H∗is usually formed on the Pd surface via the interaction between hydrogen and the d orbital of Pd atom, which was proved by the change of Pd states in XPS analysis(Fig.S6A).Subsequently, the atomic H∗attacked the C−Cl bond and made it cleavage.Nevertheless, the removal rate of 2,4,6-TCP with 500 mmol/L t-BuOH was still maintained at 87.09% after 180 min,suggesting the ECH performance was driven by both direct electron transfer and atomic H∗reduction.

CV curves of Pd/MnO2/Ni foam cathode with or without 2,4,6-TCP were conducted to identify the role of H∗adsand H∗absin the dechlorination process of 2,4,6-TCP.Fig.3B showed three oxidation peaks were located in potential ranges of −0.80 V to −0.60 V,−0.40 V to −0.20 V, and −0.20 V to −0.00 V assigning to H2, H∗absand H∗ads, respectively.The oxidation peak corresponding to H∗adsdisappeared after the addition of 2,4,6-TCP, indicating that H∗adswas the key active hydrogen species, while H2and H∗abswere negatively contributed to ECH.

To determine mechanistic insights into 2,4,6-TCP hydrodechlorination, we identified the electrolyzed intermediates by using GC–MS.Pd/MnO2/Ni foam cathode generated detectable concentration of ion cluster for chlorophenols, for instance, 2,4-DCP, 2,6-DCP (m/z=161.96), 2-CP, 4-CP (m/z=128.0) as well as phenol(m/z=94.04) (Fig.S8 in Supporting information).Based on GC–MS analysis results, the 2,4,6-TCP dechlorination pathway was proposed (Fig.3C).Firstly, the ortho and para-Cl of 2,4,6-TCP were attacked by atomic H∗to form 2,4-DCP and 2,6-DCP.Subsequently, as C−Cl was further attacked by atomic H∗, both 2-CP and 4-CP were detected in MS analysis.It was noteworthy that 2-CP was maintained at a relatively higher concentration all along compared to 4-CP.This was because chemical reactions with lower Gibbs free energy and redox potential tended to occur [43].As shown in Table S5 (Supporting information), when 2-CP was served as one of the products, reaction possessed the lower E and ΔG than 4-CP,while the opposite was true when it was applied as the reactant.It was likely to explain why 4-CP existed for a long time compared to 2-CP.Moreover, the dechlorination rates (R) of 2,4-DCP was in the order R(2-CP) > R(4-CP) due to the stereoscopic effect of orthochlorinated group [44].Finally, phenol serving as the major end metabolite revealed the ECH process.

To summarize, Fig.4 showed the schematic concept of the ECH reaction mechanism of 2,4,6-TCP by Pd/MnO2/Ni foam cathode.In brief, the presence of MnO2enhanced the generation of atomic H∗.The fate of atomic H∗mainly included H∗ads, H∗absand H2,in which only H∗adsparticipated in the hydrodechlorination process, while H∗abswere inert, and H2inhibited the hydrodechlorination process.The introduction of MnO2strengthened the HO−H bond cleavage, resulting in the enhanced dechlorination activity by atomic H∗.Both direct reduction and indirect reduction participated in the 2,4,6-TCP dechlorination process.

Fig.3.(A) Quenching trials as a function of various t-BuOH concentrations.The inset showed the removal rates of 2,4,6-TCP with different t-BuOH concentrations and corresponding pseudo-first-order rate constants.(B) CV curves of Pd/MnO2/Ni foam cathode with or without 2,4,6-TCP.(C) Time course of typical intermediates determined on Pd/MnO2/Ni foam cathode at the cathode potential of −1.0 V.

Fig.4.Schematic concept of ECH reaction mechanism of 2,4,6-TCP by Pd/MnO2/Ni foam cathode.

The Pd/MnO2/Ni foam was successfully fabricated by using chemical deposition and electrodeposition methods, and the superior dechlorination activity of 2,4,6-TCP was obtained by optimizing Pd load, MnO2load and cathode potentials in our study.The introduction of MnO2enhanced the generation of atomic H∗.Only H∗adswere confirmed that truly facilitated the ECH process,and besides the reduction by H∗ads, the direct reduction by cathode electrons also participated in the 2,4,6-TCP dechlorination process.As the major end metabolite, phenol revealed the ECH process.The highest ECH performance was obtained under the neutral condition and no significant inhibition was observed in presence of most common inorganic ions, representing the favorable application potential.This study gave strategies for the effective removal of persistent halogenated organic pollutants in groundwater.

Declaration of competing interest

We declare that we have no financial and personal relationships with other people or organizations that can inappropriately influence our work.

Acknowledgments

This study was supported by the NSFC-JSPS joint research program (No.51961145202) and the Natural Science Foundation of Heilongjiang Province, China (No.C2018035).

Supplementary materials

Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.cclet.2021.11.068.