Introduction
Ferroelectric material can reveal the switching of electric polarization, and exhibit many excellent characteristics, such as pyroelectricity, piezoelectricity, optical nonlinearities, and the electro-optic effect. Since ferroelectricity was discovered in Rochelle salt (NaKC4H4O6∙4H2O) (Valasek Reference Valasek1921), ferroelectric material has always been an important topic in condensed matter science, and the application is important as a basic element of devices and memories (Lines & Glass Reference Lines and Glass2001). Recently, much attention has been focused on developing inorganic-organic hybrid ferroelectric materials (Guloy et al. Reference Guloy, Tang, Miranda and Srdanov2001; Zhao et al. Reference Zhao, Li, Ren, Song and Jin2009; Chen et al. Reference Chen, Shi, Zheng, Zhou, Zhu and Priya2015; Huang et al. Reference Huang, Sun, Wang, Zhang, Ji, Luo, Zhang and Chen2017; Xu et al. Reference Xu, Zhang, Gao, Gan, Sun, Cui, Zepeng, Cai and Wu2017; Ji et al. Reference Ji, Liu, Han, Tao, Sun and Luo2018) because of the advantages of combining an inorganic unit (rigid, dimensional stability, and thermal stability) with organic components (straightforward synthesis approach, easily tailored molecular structure, and functional properties). The technologically useful ferroelectricity requires that the compound crystallizes only in a polar space group. Over the past decades, supramolecular chemistry and crystal engineering strategies have exhibited remarkable capabilities in self-assembly of designed functional materials (Lehn Reference Lehn1985; Desiraju & Parshall Reference Desiraju and Parshall1989; Fabbrizzi & Poggi Reference Fabbrizzi and Poggi1995; Braga & Grepioni Reference Braga and Grepioni1999; Guloy et al. Reference Guloy, Tang, Miranda and Srdanov2001; Desiraju Reference Desiraju2007; Radhakrishnan Reference Radhakrishnan2008).
As a natural clay mineral, kaolinite (K, chemical formula Al2Si2O5(OH)4) with a typical 1:1 dioctahedral layered structure (Murray & Keller Reference Murray, Keller, Murray, Bundy and Harvey1993) forms an excellent host material. The layers are stacked along the c axis via hydrogen bonds and dipolar interactions, and due to attractive Van der Waals forces. The polar interlayer space can accommodate guest molecules to form intercalation compounds. Due to the tight bonding between the kaolinite layers, only a few small molecules with high polarity can be intercalated directly into kaolinite. These include dimethyl sulfoxide (DMSO) (Fafard et al. Reference Fafard, Terskikh and Detellier2017), potassium acetate (KAc) (Kristof et al. Reference Kristof, Mink, Horvath and Gábor1993; Yariv et al. Reference Yariv, Lapides, Nasser, Lahav, Brodsky and Michaelian2000; Frost et al. Reference Frost, Kristof, Schmidt and Kloprogge2001), hydrazine (Barrios Reference Barrios1977; Frost et al. Reference Frost, Kloprogge, Kristof and Horvath1999; Cruz & Franco Reference Cruz and Franco2000), formamide (Adams et al. Reference Adams, Reid, Thomas and Walters1976; Theng Reference Theng1984), phenylphosphonic acid (Breen et al. Reference Breen, D'Mello and Yarwood2002), and urea (Seifi et al. Reference Seifi, Diatta-Dieme, Blanchart, Lecomte-Nana, Kobor and Petit2016; Zhang et al. Reference Zhang, Liu, Gao, Li, Liu, Li, Boyd, Johnston and Teppen2016), whereas most other functional guests are usually indirectly intercalated by replacing previously intercalated molecules, which limits the application of kaolinite intercalation compounds in environmental protection (He et al. Reference He, Guo, Zhu and Yang2001), polymer-clay nanocomposites (Liu Reference Liu2007; Letaief & Detellier Reference Letaief and Detellier2009; Matusik et al. Reference Matusik, Stodolak and Bahranowski2011), electrochemistry (Braggs et al. Reference Braggs, Fornasiero, Ralston and Smart1994; Letaief et al. Reference Letaief, Tonle, Diaco and Detellier2008; Dedzo et al. Reference Dedzo, Letaief and Detellier2012), catalysis (Nakagaki et al. Reference Nakagaki, Benedito and Wypych2004; Ayodele & Hameed Reference Ayodele and Hameed2013), pharmaceuticals (Carretero & Pozo Reference Carretero and Pozo2009), ionic liquids (Letaief & Detellier Reference Letaief and Detellier2009; Dedzo & Detellier Reference Dedzo and Detellier2013), and dielectrics (Orzechowski et al. Reference Orzechowski, Słonka and Głowinski2006; Wang et al. Reference Wang, Zhou and Du2018), etc. The ferroelectricity of kaolinite intercalation compounds, however, has received much less attention. Using Grand canonical Monte Carlo simulations, Croteau et al. (Reference Croteau, Bertram and Patey2010) investigated the structure of water adsorbed on kaolinite surfaces with trenchlike defects, which indicated that kaolinite has unique advantages in the construction of inorganic-organic hybrid ferroelectrics. Others (Zhao et al. Reference Zhao, Gao, Ren, Yuan and Lu2012; Qiao et al. Reference Qiao, Ding, Zhao, Li, Liu and Ren2017) reported that kaolinite intercalated with selected small molecules (DMSO, urea, and ethylene glycol) exhibited ferroelectricity. Consequently, insofar as crystal-oriented synthesis of kaolinite-based inorganic-organic hybrid ferroelectrics is concerned, a feasible strategy for achieving polar hybrid crystals is to select or design suitable organic guest molecules to intercalate kaolinite.
In the present study, PABA was selected as the candidate for intercalation for two reasons: (1) PABA with an acceptor amide group on the first carbon atom of the benzene and a donor amino group on the fourth carbon atom (Fig. 1) has high polarity, so the formation of an intercalation compound is possible through hydrogen bonds interacting with SiO4 or hydroxyl groups (Al-OH) on the inner surface of kaolinite layers. (2) PABA can provide maximum acentricity for the molecule and give rise to a non-centrosymmetric crystal structure (Almeida et al. Reference Almeida, Monte, Matos and Morais2013) due to the aforementioned special structure, thus the intercalation compound is more likely to crystallize in a polar space group, which is favorable for the construction of hybrid ferroelectrics with large polarity. However, negative factors still exist: firstly, PABA has a high melting point (182–184°C) and a larger molecular size than DMSO and urea, so kaolinite intercalation is difficult to achieve using traditional methods (melting or mechanical stirring for a long time); secondly, direct observation of the binding mode between PABA and kaolinite is difficult, which means the mechanism by which ferroelectricity is achieved is unclear. In contrast, the autoclave method is suitable for higher-temperature synthesis reaction, because of the advantages of the lack of requirement for stirring and a definite pressure being exerted to facilitate the intercalation, so the above problems of synthesis can be solved easily. At the same time, theoretical calculations can be used in the design of new materials by simulating the bonding mode between guest molecules and the kaolinite layer, by supplementing the information obtained from chemical experiments, and by predicting ferroelectricity of the intercalation compound that have not been observed so far. The purpose of the present study was to test the autoclave method for preparing K-PABA and to characterize the resulting intercalated structure and explain a possible mechanism for its ferroelectric properties.
Experimental
Materials
The kaolinite sample (KGa-1b) was obtained from the Source Clays Repository of The Clay Minerals Society. It has an average particle size of 1250 mesh, contains some Fe2O3 (≤0.03 wt.%), Al2O3 (36±0.5 wt.%), and SiO2 (47±0.5 wt.%). Analytical grade DMSO and PABA purchased from Aladdin Reagent Co., Ltd. (Shanghai, China) were used without further purification, and all operations were carried out in air.
Preparation of K-PABA Intercalation Compound
The K-PABA intercalation compound was prepared in an autoclave by replacing a DMSO precursor. The K-DMSO precursor was prepared first according to the disclosed procedure (Takenawa et al. Reference Takenawa, Komori, Hayashi, Kawamata and Kuroda2001). Secondly, PABA (1.0 g) and K-DMSO (1.0 g) were sealed in a 25 mL autoclave for heating at 180°C for 12 h and then cooled to ambient temperature. 30 mL of hot ethanol was added and then the sample was stirred for 10 min, separated by a Buchner funnel suction filter (100 mm diameter membrane) to remove adsorbed guest molecules, and then the solid powder was dried at 60°C for >24 h.
In addition, to achieve the optimum preparation conditions of K-PABA, a series of K-PABA samples was synthesized at various preparation temperatures (from 150 to 200°C) and times (from 7 to 17 h at 180oC), and the relationship between the rate of intercalation (I.R.) and preparation conditions was investigated in detail.
Characterization
X-ray diffraction (XRD) patterns of kaolinite and the intercalated samples were measured using a Bruker D8 Advance Diffractometer (Bruker AXS, Rheinstetten, Baden-Wuerttemberg, Germany) with CuKα radiation (λ = 1.5418 Å) at ambient temperature. Fourier-transform Infrared spectroscopy (FTIR) was carried out using a Nicolet AVATAR-360 spectrometer (Nicolet Instrument Corp., Madison, Wisconsin, USA) (KBr disks, 4000–400 cm–1). Thermogravimetric (TG) analysis was performed using a simultaneous SDT 2960 DTA-TG (TA instrument, Inc., New Castle, Delaware, USA); the sample was held in a Pt pan under a N2 flow rate of 100 mL min–1 and heated at a ramping rate of 10°C·min–1 from ambient temperature to 700°C. The morphology of the samples was investigated using an Hitachi S-4800 scanning electron microscope (SEM) (Hitachi High-Technologies Corp., Tokyo, Japan). The polarization–electric field (P-E) hysteresis loops for the powdered samples were performed using a Precision Multi-ferroeletric tester (Radiant Technologies, Inc. Albuquerque, New Mexico, USA) in an AC electric field, the powdered samples were pressed into pellets at 10 MPa for 1 min with a thickness of ~0.359 mm and an area of 132.67 mm2, and each pellet was sandwiched by copper electrodes and immersed in insulating oil while being measured.
Computational Details
The primitive unit cell of kaolinite is known as a = 5.15 Å, b = 8.94 Å, c = 7.39 Å, and α = 91.93°, β = 105.04°, γ = 89.79°, then one supercell was built with the basal spacing of layers set to 14.5 Å obtained from the XRD measurements. The possible structural models of K-PABA were calculated using the quantum mechanical program CASTEP (Segall et al. Reference Segall, Lindan, Probert, Pickard, Hasnip, Clark and Payne2002; Clark et al. Reference Clark, Segall, Pickard, Hasnip, Probert, Refson and Payne2005) with the GGA-PW91 (generalized gradient approximation-Perdew-Wang 91) function (Perdew et al. Reference Perdew, Burke and Ernzerhof1996), which is appropriate for the relatively weak interactions (Martorell et al. Reference Martorell, Kremleva, Krüger and Rösch2010). The cutoff energy based on the plane wave was set at 300 eV for K-PABA. The convergence parameters were set as follows: self-consistent field (SCF) tolerance, 1×10–6 eV per atom; total energy tolerance, 1×10–5 eV per atom; maximum force tolerance, 0.03 eV Å–1; maximum stress component, 0.05 GPa; and displacement of convergence tolerance, 0.001 Å. The other calculation parameters were set at the default values in the CASTEP code. All the aforementioned limits were able to constrain accurately the simulative process to end with a stable construction of K-PABA.
Results and Discussion
XRD Analysis
The XRD pattern (Fig. 2a) displayed the 001 reflection at 7.2 Å (12.12°2θ) of pure kaolinite. A new, sharp, high-intensity peak appeared in K-DMSO with a d 001 value of 11.2 Å (7.94°2θ), which represents a lattice expansion of 4.0 Å relative to pure kaolinite. The PABA-intercalated kaolinite obtained from K-DMSO was a white powder with a basal lattice spacing of 14.5 Å (6.12°2θ), which represents a lattice expansion of 7.3 Å in relation to pure kaolinite. In this case, the I.R. can be calculated directly as 78.4% based on the change of diffraction peak intensity before and after intercalation (Ledoux & White Reference Ledoux and White1964). Meanwhile, the XRD pattern (Fig. 2a) showed no evidence of a DMSO residue in the K-PABA, implying that the substitution of PABA was completed and that the only crystalline materials observed within the matrix were kaolinite and the expected K-PABA intercalated composite. In addition, based on the plots of intercalation rate versus temperature and time (Fig. 2b and c), the ideal condition for K-PABA synthesis was 180oC for 12 h.
FTIR Analysis
FTIR measurements (Fig. 3) revealed three characteristic bands for kaolinite at 3695, 3662, and 3652 cm–1, being ascribed to weak, acid-active, inner-surface hydroxyl groups (inner-surface OH) and a sharp peak at 3620 cm–1, being assigned to an inner hydroxyl band (inner OH) (Balan et al. Reference Balan, Delattre, Guillaumet and Salje2010). K-DMSO produced two bands at 3695 and 3620 cm–1 similar to pure kaolinite, while the band at 3652 cm–1 disappeared and that at 3662 cm–1 became sharper, suggesting the formation of hydrogen bonds between inner-surface OH and the Si=O group of the DMSO. Slightly stronger new bands were observed at 3536 and 3500 cm–1, suggesting that the inner surface OH bands were disturbed by intercalation molecules. In addition, K-DMSO also had two weak, broad bands at 3023 and 2937 cm–1; while liquid DMSO appeared at 2994 and 2911 cm–1, which is related to the symmetric stretching of –CH3 in DMSO.
In contrast, only two sharper characteristic OH bands of kaolinite were observed at 3695 and 3620 cm–1 in K-PABA, as well as a band at 3500 cm–1 and a new shoulder peak at 3635 cm–1, indicating that PABA successfully replaced DMSO and formed hydrogen bonds with kaolinite layers. Bands at 3481 cm–1 and 3468 cm–1 correspond to the stretching vibration of N–H of PABA. A strong, sharp peak at 3408 cm–1 and two broad bands at 3357 cm–1 and 3224 cm–1 indicate the associations existing between guest molecules, similar to that observed for pure PABA molecules at 3326 cm–1 and 3323 cm–1. In addition, the stretching vibration bands of the aromatic ring skeleton (Gardolinski et al. Reference Gardolinski, Ramos, de Souza and Wypych2000) were observed in the range 1606–1443 cm–1, while a band at 840 cm–1 confirmed the guest molecule was a para-disubstituted compound on the benzene ring. The band at ~1630 cm–1 was attributed to the overlap of N–H bending and C=O stretching vibrations, which moved ~20 cm–1 toward the low-wavenumber region compared with C=O of the free amide group, indicating that C=O of the amide group in the guest molecule also formed hydrogen bonds. Therefore, some conclusions can be drawn: K-PABA was prepared successfully, which was consistent with the XRD results, and both C=O and N–H from PABA molecules participated in the formation of hydrogen bonds. However, whether N–H comes from the amino or amide group needs further research.
TG Analysis
The TG curves for the pure kaolinite (Fig. 4) showed a mass-loss step (at ~500oC) related to the thermal dehydroxylation of kaolinite into meta-kaolinite. For K-DMSO, the initial decrease in the mass was probably related to the removal of adsorbed water, but the intense decrease in mass (~13%) between 100 and 207oC was associated with the removal of DMSO molecules. The jump at ~500oC was similar to that observed in pure kaolinite. The TG curve for K-PABA proved that heating at 207oC decreased the mass of the sample by 21.19% and then, after heating at 344oC, to a relative decrease of 11.46%. The 21.19% decrease was assigned to deintercalation and decomposition of PABA molecules, indicating that the kaolinite inorganic layers improved the thermal stability of PABA (Prasad et al. Reference Prasad, Krishnakumar, Nagalakshmi and Manohar2011). The second decrease (at ~344–620oC) was interpreted as a dehydroxylation of kaolinite similar to pure kaolinite. This process generated a mass loss of 11.46%, similar to the 13.9% value predicted from the theoretical formula of pure kaolinite (Al2Si2O5(OH)4) (Kristof et al. Reference Kristof, Frost, Kloprogge, Horváth and Gábor1999). Finally, the chemical formula of K-PABA can be expressed as Al2Si2O5(OH)4∙(PABA)0.7 on the basis of mass loss (21.19%) and intercalation rate (78.4%).
SEM Analysis
The SEM image of pure kaolinite morphology (Fig. 5a) revealed a well ordered stack with a tightly layered structure. The morphology of K-DMSO is almost the same as pure kaolinite (Fig. 5b), indicating that K-DMSO still has a long-range ordered structure. In contrast, the images for K-PABA (Fig. 5c, d) showed small stacks and some individual, thin flakes related to exfoliation or a departure from a layered structure, confirming that the long-range order of kaolinite was destroyed but the short-range order remained. The SEM results are direct evidence that the kaolinite layered structure was changed by PABA intercalation.
Structural Optimization
From the calculated chemical formula Al2Si2O5(OH)4∙(PABA)0.7 of K-PABA, each unit cell contains a PABA molecule. The molecular dimension of PABA, measured between the p-substituted amino group hydrogen and the hydrogen atom of the amide group was determined to be ~7.9 Å (Fig. S3) using Material Studio 6.0 software. Therefore, based on the variations in the 7.3 Å interlayer spacing of the kaolinite after intercalation, each intercalated molecule was oriented at an angle of ~70° to the plane of the kaolinite layer. In fact, this conclusion agrees with the 50–75o orientation angle that is normally assumed by hydroxyl groups at the surface of the kaolinite layer (Olejnik et al. Reference Olejnik, Aylmore, Posner and Quirk1968; Gardolinski et al. Reference Gardolinski, Ramos, de Souza and Wypych2000; Kloprogge Reference Kloprogge2019). Moreover, two orientations of PABA are possible between kaolinite layers because amide and amino groups on benzene rings have the ability to form hydrogen bonds with the kaolinite layer. In the structural optimizations (Fig. 6), the amino group close to the Si-sheet was labeled model A and that near the Al-sheet was marked as model B (Fig. S2). In the process of geometric optimization of K-PABA performed using the CASTEP module, all the atoms and unit-cell parameters (a, b, α, β, and γ) were relaxed, but the c-axis length was fixed at 14.5 Å.
Both optimization models (A and B) (Table S1) yielded similar crystal structure parameters. The remarkable difference was that the calculated final energy of model A was 166.6 kJ/mol less than model B, which suggested that model A was more reasonable. The rationality of model A can be confirmed further by simulating XRD (Fig. S1), where the position of the 001 peak in the simulated (6.12°2θ) patterns deviated from the experimental ones (6.12°2θ) by <1%.
The structural parameters of optimized and crystalline PABA are given in Table S2. Some optimized bonds (C7–O1, C7–N1, and C4–N2) were longer than those in the crystalline state. In addition, angle 1 became larger and angle 2 became smaller, indicating that the guest molecule was subjected to a force in the direction of the c axis of kaolinite. This force was ascribed to the hydrogen-bond interactions between kaolinite and PABA. On the other hand, torsion angles observed in both states indicated that the benzene ring of PABA was not coplanar with two substituent groups, suggesting amino and amide groups were affected by different forces. The resulting differences of torsion angles provide strong evidence that PABA in K-PABA was also affected by inorganic laminates, except for the intermolecular interactions of organic parts similar to those in the crystal state; in other words, strong hydrogen bond interactions exist between kaolinite and PABA, as well as between PABA molecules, indicating that both amino and amide groups participate in the interaction process.
With respect to H···Y hydrogen-bond studies, the bond distance was limited to <3.2 Å and a minimum value for the X-H···Y angle was set as 90° (Schuster et al. Reference Schuster, Zundel and Sandorfy1976; Jeffrey Reference Jeffrey1997; Steiner Reference Steiner2002; Desiraju & Steiner Reference Desiraju2001; Benco et al. Reference Benco, Tunega, Hafner and Lischka2001). Notable hydrogen bonds were observed in model A in three places: (1) inside kaolinite laminates, including the Al-sheet and Si-sheet; (2) between PABA and kaolinite layers (PABA-kaolinite), such as N2–H···O(Al-sheet), N1···H–O(Si-sheet), and O1···H–O(Si-sheet); and (3) between PABA molecules (PABA-PABA), i.e. O1···H–N1 (3a = 2.185 Å, 3b = 2.996 Å) of amide groups formed an infinite one-dimensional chain along the a axis. Apparently, the hydrogen bonds of (1) were responsible for the stability of the inorganic layer. The latter two in (2) and (3) played an important role in the formation of K-PABA. As a consequence, the guest molecules were actually inclined between kaolinite layers, and formed an infinite two-dimensional network on the ac plane through the aforementioned two intermolecular hydrogen bonds, further revealing that an amide group is involved in the formation of hydrogen bonds of PABA-kaolinite and PABA-PABA, while the amino group only participates in that of PABA-kaolinite.
Ferroelectricity
A hysteresis loop is one of the prominent characteristics of ferroelectric materials, which usually means that the sample contains some macroscopic polarizations, and the macroscopic polarization axis can be reversed under the action of an external electric field. The polarization–electric field (P-E) measurements (Fig. 7) displayed a hysteresis loop of K-PABA powder at room temperature, 14 kV/cm electric field, and 10 Hz frequency, at which the saturation polarization strength (P s), residual polarization strength (P r), and coercive field (E c) were 0.06 μC/cm2, 0.0352 μC/cm2, and 7.1 kV/cm, respectively. In contrast, pure kaolinite exhibits a linear hysteresis loop suggesting no ferroelectricity; while K-DMSO has a hysteresis loop (Zhao et al. Reference Zhao, Gao, Ren, Yuan and Lu2012), but the E c of K-PABA is significantly larger, probably due to the larger size of the PABA molecule. In addition, the P s of K-PABA is smaller than that of a typical ferroelectric crystal BaTiO3 (Xiong et al. Reference Xiong, Pernice, Ngai, Reiner, Kumah, Walker and Tang2014).
Free PABA with high polarity has a permanent dipole along the axis from the amino group to the amide group, and this dipole becomes tilted at an angle of ~70o relative to the plane of the kaolinite layer in K-PABA, as shown by arrow 1 in Fig. 8. Arrow 2 refers to the dipole of intermolecular hydrogen bonds between PABA molecules on the a axis, so the macroscopic polarization axis of the whole K-PABA should be almost parallel to the c axis of kaolinite as indicated by arrow 3. In fact, the intermolecular hydrogen bond energy is usually ~10–65 kJ/mol (Horiuchi & Tokura Reference Horiuchi and Tokura2008; Szafrański & Katrusiak Reference Szafrański and Katrusiak2008), which is relatively easily destroyed by an external electrics (Williams & Moore Reference Williams and Moore1972; Ye et al. Reference Ye, Fu, Zhang, Zhang, Xiong and Huang2008; Cornelissen et al. Reference Cornelissen, Biler, Urbanaviciute, Norman, Linares and Kemerink2019). On the other hand, the unit-cell volumes of K-PABA and kaolinite are known to be 756.19 Å3 and 327.84 Å3, respectively; thus the interlayer space of the intercalated compound is 428.35 Å3 (∆V = 751.19 – 327.84 Å3), which is greater than the 167.54 Å3 of each PABA volume (molecular unit-cell volume is 335.08 Å3, Z = 2) (Alléaume Reference Alléaume1967; Berkovitch-Yellin et al. Reference Berkovitch-Yellin, Van Mil, Addadi, Idelson, Lahav and Leiserowitz1985). This implied that the kaolinite interlayer can provide a space for PABA to move, although the guest molecule cannot flip between kaolinite layers like small molecules (DMSO and urea) can, due to their large molecular size (~7.9 Å). For these reasons, when an external AC electric field was applied to K-PABA along the c axis, the intermolecular hydrogen bonds of PABA-kaolinite are reversed (arrow 1), while the original direction of the a axis was maintained (arrow 2) because it was almost perpendicular to the direction of the external electric field. As a result, the dipole of the entire K-PABA (arrow 3) was reversed under the action of the external AC electric field. Obviously, the real ferroelectric mechanism of K-PABA is more complicated than the above analysis, and further research is necessary.
Conclusions
In summary, autoclaving is an effective method to intercalate kaolinite with PABA at 180°C over 12 h. The acquisition of a K-PABA polar hybrid crystal further proved that the strategy of intercalating kaolinite with a polar guest is successful. The interlayer spacing, intercalation ratio, and mass-loss rate of K-PABA reached 14.5 Å, 78.4%, and 21.19%, respectively. Theoretical calculation disclosed correctly the arrangement of PABA between kaolinite layers. The intercalation of PABA not only changed the morphology and structure of kaolinite, but also generated the intermolecular hydrogen bonds in kaolinite-PABA and PABA-PABA, which lead to ferroelectricity of K-PABA. This work provided a useful exploration for the design and application of kaolinite hybrid ferroelectrics.
Electronic supplementary material
The online version of this article (https://doi.org/10.1007/s42860-019-00036-x) contains supplementary material, which is available to authorized users.
Acknowledgements
The authors gratefully acknowledge the financial support provided by the Natural Science Foundation of Anhui Education Department (No. KJ2018A0381) and the Anhui Provincial Key Laboratory of Optoelectronic and Magnetism Functional Material (ZD2017003).