• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Complex coordinated functional groups:A great genes for nonlinear optical materials

    2022-06-20 07:59:10WeikngWngDjingMeiShoguoWenJinWngYundongWu
    Chinese Chemical Letters 2022年5期

    Weikng Wng, Djing Mei,c,*, Shoguo Wen, Jin Wng, Yundong Wu

    a College of Chemistry and Chemical Engineering, Shanghai University of Engineering Science, Shanghai 201620, China

    b Department of Chemistry and Biochemistry, Wichita State University, Wichita, KS 67260, United States

    c State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, China

    Keywords:Complex coordinated Nonlinear optical Second harmonic generation Optical band gap Structure-performance relationship

    ABSTRACT Complex coordinated functional groups [MAxBy] (M = Central coordination element; A, B = P, O, S, Se,F, Cl, Br or I) are composed of different types of anions A, B jointly linked to the same central cation M,which are in high potential to tune the physical properties of materials, e.g., second-order susceptibility,energy gaps and birefringence.Recently, Compound containing complex coordinated functional groups have attracted great attention in the nonlinear optical (NLO) field and a large number of this type crystals exhibit promising NLO performance.However, the inherent relationship between ionic group structure and optical properties of complex coordinated NLO materials have not been systematically studied.This article systematically summarizes complex coordinated NLO materials in recent five years from the perspective of the internal relationship between crystal structure and optical properties.In addition, we propose the ideal combination and arrangement modes for structural building units, and also reveal the influence of complex coordinated functional groups [MAxBy] toward the NLO response, optical band gap and phase matching ability of complex coordinated NLO materials.

    1.Introduction

    Nonlinear optical (NLO) crystals can generate tunable laser beams from deep ultraviolet (DUV, wavelength less than 200 nm)to mid IR (MIR, wavelength between 3 μm and 30 μm)viafrequency conversion technology, thus plays an important role in optoelectronic and laser technology [1–14].Researchers have been dedicated to explore NLO crystals with good performance since Frankenet al.discovered second harmonic generation phenomenon in quartz crystal in 1961 [15].In the past few decades, many NLO crystals with comprehensively excellent properties have emerged,such asβ-BaB2O4(β-BBO) [16], LiB3O5(LBO) [17], LiNbO3(LN)[18], KH2PO4(KDP) [19] and KTiOPO4(KTP) [20].With the extensive application of laser from the 21stcentury, the demand for DUV and MIR laser keeps increasing [21–25].However, the commercialized crystals in the DUV and MIR region cannot meet the practical needs due to the lack of comprehensive performances [26–29].The prerequisites of the ideal NLO crystal in these two important wavelength regions are large second harmonic generation (SHG) coefficientdij, suitable birefringence valueΔnand wide optical band gapEg[30,31].Large SHG coefficientdijcorresponds to high frequency conversion efficiency, and suitable birefringence valueΔncan ensure the phase matching output [32–34].Wide optical band gapEgcorresponds to short ultraviolet cut-off edge and large laser damage threshold [35].

    According to the anion group and functional moieties theory,the optical properties of the compounds are mainly affected by the anion functional groups in the structure [36–38].Keeping a good balance amongdij,ΔnandEgis not easy for NLO crystals with only one anionic functional group [39].Rational combination of various anionic functional groups is critical to obtain the ideal NLO crystals.In the last few years, different research groups have systematically summarized the NLO materials with mixed anion.For example, De Yoreoet al.comprehensively reviewed the crystal chemistry and optical properties of metal halide borate NLO materials, according to the arrangement of anionic groups (0D clusters, 1D chain, 2D layer, 3D network) [19].Panet al.specifically analyzed the crystal chemistry and optical properties of the mixed anion NLO materials [39].Liet al.summarized different types of IR mixed anion NLO inorganic compounds from the perspective of structure-property relationship [40].

    Inorganic salts containing two or more anions or anionic groups are normally defined as mixed anion compounds.For example,two type ions Cl-and S2-exist in the structure of salt sulfide halide compounds [K3Cl][Ga3PS8] [41].There are two single anion functional groups [IO3] and [NO3] in iodate nitrates Sc(IO3)2(NO3)[42].It should be noted that the compounds above are composed of different types of anions or anionic groups.Different from the traditional mixed anion compounds, the different anions A and B of [MAxBy] groups in complex coordinated crystals jointly link to the same central atom M [43].Complex coordinated functional groups exhibit superior propertiesthan the single-anion group from the following aspects:(1) The structure of mixed anion group[MAxBy] has a stronger distortion and polarity, which is conducive to a strong SHG effect.(2) Due to different size and electronic properties of A and B anions, the mixed anion group [MAxBy]has stronger anisotropy and large birefringence values, which is conducive to the appropriate phase matching ability.(3) In the DUV nonlinear optics field, the compounds containing mixed anion group produced by the introduction of halogen element with strong electronegativity,e.g., F and oxygen in [B/PO4] have a wider DUV transparent region and also strong SHG effects and a suitable birefringence value.(4) In the field of IR NLO region, the compounds containing [MOxSy] gene containing O and S can have large band gap and strong NLO effect [43–46].The mixed anion functional group [MAxBy] has been proven to be superior NLO active units, and the compounds containing this group are more likely to keep a good balance among three key property parameters (dij,ΔnandEg) [27,40].For example, the well-performed DUV NLO crystal KBe2BO3F2(KBBF) contains the mixed anion group [BeO3F] [47,48].Moreover, the strongly polarized mixed anion group [SbOS4] has been recognized as the origin of the strong SHG effect in the new MIR NLO crystal Sr6Cd2Sb6O7S10[46].Since 2017, complex coordinated NLO materials have attracted great attention in the NLO field.Different research groups have designed and synthesized a large number of NLO optical crystals with mixed anion groups, which have potential to become candidates for the next generation of DUV and MIR NLO crystals.For example, AB4O6F(A = Rb, Cs, K, Na, NH4) series [42,49–51], A(B5O7)F3(A = Sr,Ca, Pb) series [52–54], Ba3Mg3(BO3)3F3[55], CsAlB3O6F [56],γ-Be2BO3F [57], C(NH2)3SO3F [58], Ba2NaP2O7Cl [59], CsSiP2O7F [60],NaNH4PO3F·H2O [61], have short ultraviolet cut-off edges and strong SHG effect, thus are regarded as potential source of DUV NLO crystals.K5(W3O9F4)(IO3) [62], BiIO3F [63], (NH4)Bi2(IO3)2F5[64], CsVO2F(IO3) [65], K2Bi2(SeO3)3F2[66], Pb2GaF2(SeO3)2Cl [67],Pb13O6Cl9Br5[68] and MII3PnI3(MII= Zn, Cd; Pn = P, As) [69],etc.,have a wide MIR transparent area, which can achieve a good balance between a wideEgand a largedij, regarded as a potential candidate of MIR NLO crystals.The difference between traditional complex coordinated compounds and those with complex coordinated groups as structural building units (SBUs) are still not clear,especially systematic summary of how complex coordinated functional groups [MAxBy] influence on the optical properties of NLO crystal materials.

    This article focuses on complex coordinated NLO materials in recent five years.The preparation methods, crystal structure, and optical properties of this type NLO crystals are comprehensively reviewed.The regulation law of multiple critical properties (SHG effect, phasemaching capability and band gap) by the permutation pattern of the complex coordinated functional group [MAxBy]was highlighted.Viathe structure-performance analysis of the NLO complex coordinated materials, the ideal combination and arrangement modes of the anion group of UV/DUV and MIR NLO complex coordinated crystals are proposed, respectively.We hope that this article can help future research minimize time-consuming "trial and error" and provide some insights on the design and synthesis of next-generation complex coordinated NLO crystals.

    2.IR complex coordinated NLO materials

    In the field of infrared (IR) nonlinear optics, the common singleanion NLO active units (NAUs) consist of tri-coordinated and tetracoordinated groups.Tri-coordinated groups play a key role in NLO response, and exhibit large microscopic SHG effects and good optical anisotropy [70–72].Tetrahedron-based sulfides have been widely used in the research on IR NLO crystals due to their spacious IR transparent region and good balance betweenEganddij[73].In the last few years, a number of research groups have combined these single-anion units with mixed anion groups to design and synthesize a series of new IR NLO materials.In this part, some representative complex coordinated compounds are presented and the relationship between crystal structure and optical properties of this type compounds are summarized in detail (Table 1).According to the types of anions, the IR NLO crystals containing single anions and mixed anions as SBUs can be divided into:(1) [MAxBy] and[IO3] as SBUs; (2) [MAxBy] and [Se/TeO3] as SBUs; (3) [MAxBy] and[MS4] as SBUs.In addition, another type IR complex coordinated NLO materials contain only mixed-anion [MAxBy].In the field of IR nonlinear optics, this type of compound mainly includes the following three groups:(1) PbO-PbCl2-PbBr2, PbO-PbCl2-PbI2series of IR NLO crystals containing O-Pb-X complex coordinated groups;(2) MII3PnI3containing [MPn3I] mixed anion (MII= Zn, Cd; Pn = P,As) series pnictides; (3) oxychalcogenides.

    Table 1 Optical properties of IR complex coordinated NLO crystals.

    2.1.[MAxBy] and [IO3] as NLO active units

    The iodates with the lone pair electron benefit the formation of the NCS structure required by the compound for the SHG response.[IO3] group has received extensive attention in NLO region due to their wide IR transparent region [74].Recently, Mao, Ye, Poeppelmeier and other research groups combined the [IO3] group with the mixed anion group [MAxBy] to design and synthesize a series of IR NLO crystals with excellent performance, such as Ba2[GaF4(IO3)2](IO3) [75], (NH4)Bi2(IO3)2F5[64], Ba2[VO2F2(IO3)2]IO3[76].This section summarizes this type of NLO crystal from the perspective of the structure-performance relationship and anionic group structure dimensions from 0D to 3D.

    2.1.1.0D cluster formed by [MAxBy] and [IO3] anionic groups

    Ba2[GaF4(IO3)2](IO3) polymorphs are representative compounds of this type iodates.Maoet al.obtained colorless and transparentα-Ba2[GaF4(IO3)2](IO3) andβ-Ba2[GaF4(IO3)2](IO3) [75] on the basis of Ba2[VO2F2(IO3)2](IO3) [76] by chemical substitution.They have a similar crystal structure.Takingα-Ba2[GaF4(IO3)2](IO3) as an example, its structure is composed of one mix-anion [GaF4O2]group and two [IO3]-groups to form OD isolation molecular structure (Fig.1a).Ba2+is trapped in molecular voids to maintain charge balance.Both of them exhibit wide transparent areas of 0.26–12.5 μm and 0.28–12.9 μm, respectively.The powder test shows that both of them have phase matching ability and strong SHG response of 6 times that of KDP (d36= 0.38 pm/V) under 1064 nm laser irradiation.Optical experiments revealed that they have a largeEg:4.61 and 4.35 eV.The theoretical analysis of partial densities of states (PDOS) shows that the valence bands (VBs) maximum is occupied by O 2p orbitals, while the conduction bands(CBs) minimum is contributed by O 2p and I 5p orbitals, which indicates that its optical properties are mainly dominated by [IO3]groups.The 2p orbital of F is above top of the VBs, which indicates that the mix-anion [GaF4O2] group contribute the large energy band gap.

    Fig.1.Crystal structure of NLO materials with [MAxBy] and [IO3] as NLO active units.

    2.1.2.1D chain formed by [MAxBy] and [IO3] anionic groups

    α,β-Ba2[VO2F2(IO3)2]IO3:Poeppelmeieretal.synthesizedβ-Ba[VFO2(IO3)2],α,β-Ba2[VO2F2(IO3)2]IO3series[76] fluoroiodate compounds by hydrothermal method.β-Ba[VFO2(IO3)2] has a 0D isolated molecular structure which composed of [VO4F] mixed-anion polyhedra and two [IO3] units(Fig.1c).α,β-Ba2[VO2F2(IO3)2]IO3polymorphs have similar crystal structure.Both of their crystal structures are composed of[VO4F2] polyhedra and two [IO3] units to form a 1D chain, and Ba2+are incorporated between the chains.However, in theαphase of polymorphic Ba2[VO2F2(IO3)2]IO3, the order of [VO4F]mixed-anion polyhedra and [IO3] units are better than that ofβphase (Fig.1b).Optical experiments show that they have a wide transparent region (~0.5–10.5 μm) in the IR band.Polymorphic Ba2[VO2F2(IO3)2]IO3have type I phase matching ability, and can maintain a good balance between the wideEg(>2.5 eV) and strong SHG effect (~9 × KDP).

    K5(W3O9F4)(IO3):Wuet al.obtained the colorless and transparent K5(W3O9F4)(IO3) crystal by hydrothermal method [62].As shown in Fig.1d, in the crystal structure of K5(W3O9F4)(IO3),[WO5F] and [WO4F2] mixed-anion tetrahedrons form [W3O12F4]∞1D chain through common vertex connection and [IO3] groups are isolated and scattered in the chain.Optical experiments show that the transparent region of K5(W3O9F4)(IO3) is 0.32–10.5 μm,and the indirect band gap is 3.83 eV.K5(W3O9F4)(IO3) had phase matching (PM) ability, and the SHG response was 11 × KDP and 0.5 × AgGaS2(AGS) (d36= 23.6 ± 2.4 pm/V) under 1064 nm and 2100 nm laser irradiation, respectively.DFT calculations show that VBs maximum and CBs minimum of K5(W3O9F4)(IO3) are mainly determined by I 5p, W 5d, O 2p, and F 2p orbitals, which means that the optical properties of this compound are determined by mixed-anion group [WO3F] and single anion groups [WO4], [IO3].

    Sn(IO3)2F2:Luoet al.obtained a colorless and transparent fluorocompound iodate Sn(IO3)2F2by heating and cooling in a sealed high-pressure reactor fitted with the gold liner [77].As shown in Fig.1e, the 1D chain structure of Sn(IO3)2F2is connected by two[IO3] groups and a low symmetry mixed-anion [SnO4F2] polyhedral.Optical experiments show that Sn(IO3)2F2has a wide IR band transparent region and anEgof 4.08 eV.In addition, the SHG powder test show that Sn(IO3)2F2exhibited type-I PM behavior, and the SHG response was 3 × KDP@1064 nm.DFT calculation shows that the introduction of highly polar halogen F based mixed-anion group in iodate is not only beneficial to SHG response, but also beneficial to increase the optical band gap.

    2.1.3.2D layered formed by [MAxBy] and [IO3] anionic groups

    (NH4)Bi2(IO3)2F5is representative compound of this type iodates.In 2019, Fanet al.prepared the first ammonium containing multianion fluoride iodate (NH4)Bi2(IO3)2F5by hydrothermal method with KBi2(IO3)2F5as template [64].As shown in the Fig.1f, mixed-anion groups [Bi(1)O2F5] and [Bi(2)O4F4] form a 2D layered structure through co-edge connection and co-vertex connection with [IO3] group.It is worth noting that the NH4+and F distributed between layers form N–H–F hydrogen bonds, which can strengthen the interaction between layers.Optical experimental results show that (NH4)Bi2(IO3)2F5has a wideEg(3.88 eV)and a high transparency in the range of 2.5–6.8 μm.Kurtz powder test show that it had type-I PM ability, and SHG response was 9.2 × KDP@1064 nm.DFT calculations show that the VBs maximum and CBs minimum of (NH4)Bi2(IO3)2F5are mainly determined by the atomic orbitals of Bi, I, O and F atoms which indicates that optical properties of (NH4)Bi2(IO3)2F5are dominated by the [Bi(1)O2F5], [Bi(2)O4F4] and [IO3] groups.

    2.1.4.3D framework formed by [MAxBy] and [IO3] anionic groups

    BiIO3F, Bi3OF3(IO3)4:Mao and Panet al.designed and synthesized Bi-based fluorinated iodate BiIO3F, Bi3OF3(IO3)4by hydrothermal method [63,78].As shown in Fig.1g, the crystal structure of BiIO3F is built by [BiO4F4] polyhedrons and interconnection of [IO3] groups to build a 3D framework structure.The network tunnel structure of Bi3OF3(IO3)4is composed of [BiO7F2] polyhedrons connected with [IO3] groups, and part of [IO3] groups are scattered in the tunnel holes (Fig.1h).The experimental results show that BiIO3F and Bi3OF3(IO3)4both have wide MIR transparent regions, withEgof 3.97 eV and 3.70 eV, respectively.Kurtz powder test showed that they all had PM ability, and SHG strength was 11.5 and 6.0 × KDP@1064 nm, respectively.

    Ce(IO3)2F2·H2O:Abudouwufuet al.synthesized Ce(IO3)2F2·H2O by hydrothermal method [79].As shown in Fig.1i, its 3D framework structure formed by mixed-anionic [CeO5F4] groups and [IO3]groups.In addition, two H atoms are connected with O atoms to form an H2O unit interwoven with the tunnel holes.The experimental results show that theEgof Ce(IO3)2F2·H2O is 2.6 eV, and SHG response is 3 × KDP@1064 nm with PM behavior.Theoretical analysis shows that VBs maximum of Ce(IO3)2F2·H2O is composed of I and O atomic orbitals, and the CBs minimum is composed of atomic orbitals of Ce.Therefore, its optical properties are mainly determined by the mixed-anionic [CeO5F4] groups and[IO3] groups.

    Y(IO3)2F, Cd(IO3)F:Caoet al.synthesized rare earth fluorocompound iodate Y(IO3)2F and transition metal fluorocompound iodate Cd(IO3)F by hydrothermal method [80,81].As shown in the Fig.1j, the crystal structure of Y(IO3)2F is connected by complex coordinated group [YO6F2] polyhedrons through co-edge to form a curved 1D chain, which is further connected by [IO3] group and finally forms a 3D skeleton structure.The 3D skeleton structure of Cd(IO3)F is composed of the mixed anion group [CdO5F2]polyhedrons and [IO3] groups (Fig.1k).Optical experiments show that Y(IO3)2F and Cd(IO3)F have wide transparent region in the IR band, withEgof 3.91 eV and 4.22 eV, respectively.Kurtz powder test show that they all had type I PM ability, anddijwas 2 and 6.2 × KDP@1064 nm, respectively.

    CsVO2F(IO3):Chenet al.synthesized the mixed-anion fluoroiodate CSVO2F(IO3) by hydrothermal method [65].CSVO2F(IO3) is the first alkali metal vanadifluoroioioate compound with 3D anion skeleton.It is formed by the interconnection of [IO3] and [VO5F]polyhedral groups (Fig.1l).CSVO2F(IO3) has a wide MIR transparent region (0.5–10.5 μm) andEgis 2.39 eV.Kurtz powder test show that it has type I PM behavior, and SHG response is 1.1 × KTP(d33= 17.4 pm/V).DFT calculations show that the VBs maximum and the CBs minimum of CSVO2F(IO3) are mainly occupied by V,F and O atomic orbitals, which reveals that the optical properties of CSVO2F(IO3) are determined by the [IO3] and [VO5F] polyhedral groups.

    2.2.[MAxBy] and [Se/TeO3] as NLO active units

    In addition to iodate, selenite and tellurite containing three coordinated [SeO3] and [TeO3] groups have also received more and more attention in the nonlinear optics field [67].Similar to [IO3] groups, the combination of planarπ-conjugated groups[SeO3] and [TeO3] with complex coordinated groups [MAxBy] favors to the generation of IR NLO crystals with excellent comprehensive performance.Recently, Mao, Ye, OK and other research groups have found a series of IR NLO materials of this type,such as Ba(MoO2F)2(QO3)2(Q = Se, Te) [82], K2Bi2(SeO3)3F2[66],Bi3(SeO3)3(Se2O5)F [83].In this part, based on the anionic crystal structure dimension from 1D to 3D, this type of NLO crystals were comprehensively summarized from the perspective of structureproperty relationship.

    2.2.1.1D chain formed by [MAxBy] and [SeO3] anionic groups

    Lin and Youet al.synthesized isomorphically Pb2GaF2(SeO3)2X(X = Cl, Br) by hydrothermal method [67,84].They have similar crystal structure.Taking Pb2GaF2(SeO3)2Cl as an example, its crystal structure is connected by [GaO3F3] polyhedron and [SeO3]group to form a 1D chain structure, and each Cl atom connects two Pb atoms along theabdirection (Fig.2a).The results of optical experiment show that both of them have wide MIR transparent region, and theEgare 4.32 eV and 3.73 eV for Pb2GaF2(SeO3)2Cl and Pb2GaF2(SeO3)2Br, respectively.Under 1064 nm laser irradiation, their SHG response is 4.5 × KDP, and they have type I PM capability.Theoretical analysis shows that [GaO3F3] and [SeO3] units play a major role in their largedij.

    2.2.2.2D layered formed by [MAxBy] and [Se/TeO3] anionic groups

    Ba(MoO2F)2(QO3)2(Q= Se, Te):Lianget al.synthesized the colorless and transparent alkaline earth metal fluoroselenite and fluorotellurite compounds Ba(MoO2F)2(QO3)2(Q = Se,Te) by hydrothermal method [82].They belong to NCS orthogonal space groupAba2 (No.41).They have similar crystal structures.Taking Ba(MoO2F)2(SeO3)2as an example, there are two NAUs of mixed anion [MoO5F] and [SeO3] in the structure.They form a 2D layer by sharing vertexs, and Ba2+is scattered between the layers (Fig.2b).It is worth noting that the SHG intensity (7.8 × KDP@1064 nm) of Ba(MoO2F)2(TeO3)2with PM ability is significantly greater than that of Ba(MoO2F)2(SeO3)2(2.8 × KDP@1064 nm).This indicates that the SHG activity of[TeO3] group is stronger than that of [SeO3] group.TheirEgare 2.96 eV and 3.23 eV, respectively.DFT calculations show that the VBs maximum of Ba(MoO2F)2(TeO3)2is mainly occupied by atomic orbitals of O 2p and F 2p, while CBs minimum is mainly occupied by atomic orbitals of O 2p, Te 5p and Mo 4d, which indicates that its optical properties are mainly determined by [MoO5F] and[TeO3] groups.

    Fig.2.Crystal structure and NLO performance of (a) Pb2GaF2(SeO3)2Cl, (b) Ba(MoO2F)2(TeO3)2 and (c) RbGa3F6(SeO3)2.

    RbGa3F6(SeO3)2:Wuet al.designed and synthesized the alkali metal fluoride selenite compounds RbGa3F6(SeO3)2by hydrothermal method [85].In the crystal structure of RbGa3F6(SeO3)2, [SeO3]groups connect the mixed-anion groups [GaO2F4] (Fig.2c).Optical experiments show that RbGa3F6(SeO3)2have a wide transparent region in IR wavelengths.Kurtz powder test showed that RbGa3F6(SeO3)2have type I phase matching ability, and SHG strength was 5.6 × KDP@1064 nm andEgis 3.57 eV.According to the calculated PDOS diagram, the VBs maximum and CBs minimum are mainly occupied by the atomic orbitals of Ga, F,Se and O atoms, which indicates that the optical properties of RbGa3F6(SeO3)2are mainly determined by the [SeO3] and [GaO2F4]units.

    2.2.3.3D framework formed by [MAxBy] and [SeO3] anionic groups

    Lu3F(SeO3)4:Wuet al.synthesized the colorless and transparent rare earth metal fluoroselenite compound Lu3F(SeO3)4by hydrothermal method [86].The 3D frame structure of Lu3F(SeO3)4is formed by twisted [LuO7F] polyhedron and planar p-conjugated[SeO3] groups (Fig.3a).Lu3F(SeO3)4has a wide transparent region in IR wavelengths andEgis 3.57 eV.Kurtz powder test show that Lu3F(SeO3)4exhibits PM behavior, and the SHG strength of Lu3F(SeO3)4is 2.5 × KDP@1064 nm.DFT calculation revealed that the optical properties of Lu3F(SeO3)4are mainly determined by[LuO7F] and [SeO3] groups.

    Fig.3.Crystal structure and NLO performance of 5 representative 3D framework selenities.

    K2Bi2(SeO3)3F2, Bi3(SeO3)3(Se2O5)F:Shiet al.synthesized a colorless and transparent Bi-based alkali metal selenite K2Bi2(SeO3)3F2by hydrothermal method [66].In the crystal structure of K2Bi2(SeO3)3F2, Bi(1) atoms connect six O atoms and two F atoms to form a [BiO6F2] group, and Bi(2) atoms link to five O atoms and two F atoms to form [BiO5F2] group.[BiO6F2] group and [BiO5F2] group share one O and one F atom to form a 1D chain structure, and [IO3] connects these 1D chains to finally form a 3D framework structure (Fig.3b).The transparent region of the growing block K2Bi2(SeO3)3F2crystal covers the important band 3–5 μm in the IR region.K2Bi2(SeO3)3F2with phase matching ability can maintain a balance between the wideEg(3.72 eV) and the strong SHG response (15 × KDP @1064 nm).DFT calculation indicated that the optical properties are mainly determined by [BiOxFy]and [SeO3] groups.

    Chunget al.synthesized colorless and transparent Bi-based selenite fluoride Bi3(SeO3)3(Se2O5)F by hydrothermal method [83].There are three active NLO groups [BiO7], [BiO6F] and [SeO3] in the crystal structure, in which [BiO7] and [BiO6F] polyhedra form a 2D layer through co-edge connection, and [SeO3] further connects the 2D layer to form a 3D frame structure (Fig.3c).The optical results show that Bi3(SeO3)3(Se2O5)F has a wide IR transparent window.Furthermore, Bi3(SeO3)3(Se2O5)F exhibits wide band gap(3.80 eV) and strong SHG response (8.0 × KDP@1064 nm) with type-I PM behavior.Theoretical analysis shows that the strongly twisted mixed-anion [BiO6F] polyhedron play an important role in its good comprehensive performance.

    Pb2Cd(SeO3)2Cl2, Pb3(SeO3)Br4:Chen and Maet al.synthesized lead halogen selenite Pb2Cd(SeO3)2Cl2and Pb3(SeO3)Br4by hydrothermal method [87,88].As shown in Fig.3d, the crystal structure of Pb2Cd(SeO3)2Cl2is composed of [CdO6Br], [PbOxBry]and [SeO3] groups.They interpenetrate and connect to form a 3D framework structure.The crystal structure of Pb3(SeO3)Br4has three kinds of NAUs,e.g., [PbBr5O2], [PbBr4O3] and [SeO3], in which[PbBr5O2] and [PbBr4O3] form 3D tunnel structure by sharing Br atoms (Fig.3e).Under 1064 nm laser irradiation, Pb2Cd(SeO3)2Cl2and Pb3(SeO3)Br4with type I phase matching ability exhibit 1.8 and 1.0 times SHG effects that of KDP, respectively.TheirEgare 4.10 eV for Pb2Cd(SeO3)2Cl2and 3.35 eV for Pb3(SeO3)Br4, respectively.DFT calculations suggest the VBs maximum and CBs minimum of Pb2Cd(SeO3)2Cl2mainly occupy by Pb, Se, Cl, O atomic orbitals, which suggests that its optical properties are mainly depends on [PbOxBry] and [SeO3] groups.The VBs maximum and CBs minimum of Pb3(SeO3)Br4mainly occupy by the orbitals all constituent atoms, suggesting that its optical properties are mainly determinated by [PbBr5O2], [PbBr4O3] and [SeO3] groups.

    In summary, most of the IR NLO crystals composed of the complex coordinated group [MAxBy] and tri-coordinated group[Se/TeO3] have a 3D framework structure.Similar to iodates, they all have good optical anisotropy and all have phase matching capabilities.According to theoretical analysis, the introduction of highly electronegative halogens, such as F, Cl, Br, is in favor of enhancing theEgof iodate.Most compounds with complex coordinated groups [MAxBy] and [Se/TeO3] SBUs have ideal optical band gaps(larger than 3.5 eV).It is worth noting that the compounds with[BiOxFy] and [SeO3] as SBUs have good overall performance (good phase matching ability, large optical band gap, strong SHG effect).

    2.3.[MAxBy] and [MS4] as NLO active units

    Compared with oxygen-based tri-coordinated groups, few IR NLO compounds are composed of mixed anions [MAxBy] and [MS4]as SBUs.Recently, Wanget al.designed and synthesized two compounds Sr6Cd2Sb6O7S10[46] and Sr5Ga8O3S14[89] by combining the complex coordinated group [MOxSy] with the tetra-coordinated group [MS4].Herein, these two compounds are selected as representatives to discuss the influence of SBUs functional groups on optical properties of this type compound.

    Sr6Cd2Sb6O7S10:Wanget al.obtained the red Sr6Cd2Sb6O7S10crystal by solid-phase synthesis [46].As shown in Fig.4a, its 2D layered structure is formed by the mixed-anion SBUs [SbOS4],the single anion poly-coordinated group [SbS5] and the tetracoordinated group [CdS4] connected to each other on thebcplane.The powder frequency doubling test shows that Sr6Cd2Sb6O7S10has phase matching ability.When the particle size is 150–210 μm,the SHG effect is 4.0 × AGS in 2090 nm.The bandgap of Sr6Cd2Sb6O7S10is 1.89 eV.Atomistic first principles calculations show that the 5S2-electron [SbOxS5-x]7-(x = 0, 1) groups make a major contribution to its strong SHG effect.In addition, further theoretical analysis shows that the mixed-anionic group [SbOS4]has a greater contribution to thedijthan the single-anionic unit[SbS5].

    Fig.4.Crystal structure and NLO performance of Sr6Cd2Sb6O7S10 and Sr5Ga8O3S14.

    Sr5Ga8O3S14:Wanget al.obtained colorless and transparent melilite-derived oxysulfide Sr5Ga8O3S14crystals by solid-phase synthesis [89].As shown in Fig.4b, the crystal structure consists of a mixed polyanion group [GaOS3] tetrahedron and [GaS4]tetrahedron group connected to form a 3D framework structure,with Sr2+scattered in the tunnel holes.Optical research shows that Sr5Ga8O3S14with NPM behavior exhibits a largeEg(3.9 eV)and a strong SHG effect (0.8 × AGS @2090 nm).It has a wide IR transparent area of 0.3–13.4 μm.Theoretical calculations show that the mixed polyanionic group [GaOS3] tetrahedron and [GaS4]tetrahedron group play the main contribution to the SHG effect of Sr5Ga8O3S14.

    To sum up, the reasonable combination of complex coordinated groups [MAxBy] and other single anionic groups, such as [IO3],[SeO3], [TeO3], [MS4], is a very effective strategy to design and synthesize great-performance IR NLO crystals.At present, this type of NLO material, complex coordinated groups [MAxBy] and tricoordinated groups [IO3], [SeO3], [TeO3] as building units, has been widely explored.They all have good phase matching capabilities,and most of them have a wide optical band gap.The key issue is how to effectively improve the NLO effect.Interestingly, [MOxFy](M = V,Mo, W, Bi,etc.) composed of mixed anion groups centered on transition metal or metal with lone pair electrons not only provide strong SHG contribution, but also help guide tri-coordinated groups [IO3], [SeO3], [TeO3] to arrange orderly.Therefore, most compounds with this group have excellent comprehensive properties, such as K5(W3O9F4)(IO3), BiIO3F, Ba2[VO2F2(IO3)2]IO3,(NH4)Bi2(IO3)2F5, CsVO2F(IO3), Ba(MoO2F)2(TeO3)2, K2Bi2(SeO3)3F2and Bi3(SeO3)3(Se2O5)F.All of compounds above can keep a good balance between strong SHG reponse with phase matching ability and wideEg.The combination of complex coordinated group [MAxBy] and tetra-coordinated group [MS4] has the potential to design NLO crystals with strong SHG effect.For example,Sr6Cd2Sb6O7S10has an extremely strong SHG intensity of 4 times that of AGS.

    2.4.Halides-based pnictides

    Chenet al.solved the shortcomings of the band gap of phosphorus compounds by introducing strong electronegative heavy halogen "I" into phosphorus compounds [69].They obtain four kinds of iodinated phosphorus NLO crystals by heterovalent anion substitution strategy, namely MII3PnI3(MII= Zn, Cd; Pn = P,As).It has a defective diamond-like structure, and the mixed anion [MIIPnI3] tetrahedral groups arrange toward the same direction.Taking ZnPI3as an example, as shown in Fig.5a, its structure stacked by [ZnPI3] tetrahedrons along [111] direction by sharing vertexes.The optical experiment results show that theirEgare 2.85 eV and 2.38 eV, respectively.Their SHG effect are 2.7 and 4.3 × AGS, respectively, without PM capability.Taking CdPI3as an example, as shown in figure, its structure stacked by [CdPI3] tetrahedrons along [001] direction by sharing vertexes (Fig.5b).The optical experiment results show that theirEgare 2.44 eV and 2.05 eV,respectively.Their SHG effect are 2.7 and 4.3 × AGS, respectively,with type-I PM behavior.

    Fig.5.Crystal structure of ZnPI3 and CdPI3.

    2.5.Oxychalcogenides

    Compared with traditional oxides and chalcogenides, oxychalcogenides have received few attentions in NLO region.In recent years, oxychalcogenides has attracted people’s attention.Different research groups have designed and synthesized several IR NLO crystals with mixed anions [MOxQy] (Q= S,Se) as SBUs [90–92].Herein, the crystal structure and optical properties of this type compounds are analyzed in details.

    AEGeOSe2(AE = Ba, Sr):Liuet al.synthesized oxychalcogenides BaGeOSe2and SrGeOSe2by solid-phase method [44,45].They have similar crystal structures.As shown in Fig.6a, their crystal structures are composed of mixed-anion groups [GeO2Se2]connected by common edges to form a 1D chain structure, with Ba2+and Sr2+scattered among the chains.The optical band gap of BaGeOSe2and SrGeOSe2are 3.20 eV and 3.14 eV, respectively.Under 1400 and 2050 nm laser irradiation, both of them exhibit type-I PM behavior and SHG signal output intensity are 1.1 and 1.3 times that of benchmark AGS.Theoretical analysis shows that complex coordinated groups [GeO2Se2] make major contributions to their strong SHG effect.

    Fig.6.Crystal structure and NLO performance of (a) AEGeOSe2(AE = Ba, Sr), (b) SrZn2S2O and (c) Sr3Ge2O4Se3.

    SrZn2S2O:Loyeet al.obtained NCS polar oxysulfide SrZn2S2O by solid-phase synthesis.As shown in Fig.6b, its 2D layered structure consists of complex coordinated tetrahedral groups [ZnS3O] connected to each other along theacplane.The powder frequency doubling absorption spectrum curve shows that SrZn2S2O has a largeEgof 3.86 eV.DFT calculations show that the VBs maximum of SrZn2S2O is mainly occupied by atomic orbitals of O 2p and S 3p, while the CBs minimum are mainly occupied by atomic orbitals of Zn 4s and Zn 4p, which indicates that its optical properties are mainly determined by mixed-anion [ZnS3O] groups.

    Sr3Ge2O4Se3:Xinget al.obtained colorless and transparent Sr3Ge2O4Se3crystals through solid-phase synthesis [93].As shown in Fig.6c, the crystal structure of Sr3Ge2O4Se3contains [GeOSe3]mixed anion group and [GeO4] tetra-coordinated group, which are connected to form a unique dumbbell-like dimer [Ge2O4Se3]6-.The energy band gap of Sr3Ge2O4Se3is 2.96 eV.Under 1064 nm laser irradiation, Sr3Ge2O4Se3exhibits phase matching behavior,and the NLO response is 0.8 times that of AGS.DFT calculation indicates that the VBs maximum of Sr3Ge2O4Se3are mainly determined by the atomic orbits of O 2p and Se 4s, and the CBs minimum are predominantly composed of Ge 4s and 4p orbits.It is worth noting that the mixed anion group [GeOSe3] ’s contribution to SHG effect is significantly stronger than that of the single anionic group [GeO4].

    In conclusion, the comprehensive performance of IR NLO crystals with complex coordinated groups as SBUs exhibit superior performance to that of traditional single-anionic compounds.In comparison to traditional phosphides, MII3PnI3halides-based pnictides inherit the large effect of phosphides, and also have an increased optical band gap.The oxychalcogenides series also exhibit superior optical properties to conventional oxides and chalcogenides.

    3.UV/DUV complex coordinated NLO materials

    In the field of UV and DUV nonlinear optics, common singleanion NAUs are composed of oxygen-based tri-coordinated groups and tetra-coordinated groups,e.g., [BO3], [NO3], [CO3], [PO4], [SO4][20,90–96].In the last few years, a large number of UV/DUV NLO crystals with good performance have been designed and synthesized by different research groups in combination of the above groups with the mixed anion groups [MAxBy].Pan, Ye, Halasyamani and other teams applied KBBF as a template to design and synthesize many DUV NLO crystals with excellent performance.Many of these compounds not only have better comprehensive properties than KBBF, but also are expected to overcome the limitations of the layered growth habit.In addition, Mao, Ye, Zou and Halasyamani have synthesized a series of novel UV and DUV NLO materials by combining tri-coordinated ([NO3], [CO3]) and tetra-coordinated groups ([PO4], [SO4]) with mixed-anion group[MOxFy].According to the types and structural characteristics of the anions, the crystals can be divided into:(1) KBBF-derived crystals; (2) nitrate and carbonate; (3) phosphates; (4) sulfate.This part specifically summarizes UV/DUV complex coordinated NLO materials in recent five years on the relationship between crystal structure and optical properties (Table 2) in combination to theoretical analysis results.

    Table 2 Optical properties of UV/DUV complex coordinated NLO crystals.

    3.1.KBBF-derived crystals

    KBBF has been regarded as the benchmark in this field as one of the most commercialized DUV NLO materials [47,48].The discussion on the structure-property analysis of KBBF crystals has been reported in other literatures will not be summarized here.Thispart focuses on the analysis of a series of new DUV NLO materials derived from KBBF template and taking representative compounds as an example.Figs.7–10 show their crystal structure.Performance enhancement method originated from ntrinsic property is analyzed by relationship between their crystal structure and optical properties of KBBF.Depending on the type of compound, it is divided into six sections according to the type of compounds.

    3.1.1.AB4O6F (A=Rb, Cs, K, Na, NH4) family

    Wang and their team designed and synthesized a series of colorless and transparent KBBF-derived fluorooxoborate, including CsB4O6F, RbB4O6F, CsKB8O12F2, NaB4O6F, NH4B4O6F [32,49–51].These crystals were prepared by solid phase reaction method.CsKB8O12F2and NaB4O6F belong to NCS space groupsP321andC2,respectively.The crystal structure of these compounds is similar to KBBF, in which the cation K+is replaced by Cs+, Rb+, Cs+and K+mixed cations, Na+and NH4+, respectively.As shown in Fig.7,the interlayer spacing of CsB4O6F, CsKB8O12F2and NH4B4O6F is significantly smaller than that of KBBF (6.25 ?), indicating that their interlayer forces are stronger than KBBF’s.In addition, the substitution of [BO3F] for [BeO3F] can eliminate the toxicity of Be element.Optical experiments show that they all have a wide transparent region in the UV/DUV band, and these compounds have short ultraviolet cut-off edge:CsB4O6F (155 nm), RbB4O6F(<200 nm), CsKB8O12F2(<200 nm), NaB4O6F (<180 nm),NH4B4O6F (156 nm).In addition, all of these compounds have phase matching capability.As shown in Table 2, the SHG effect of CsB4O6F, CsKB8O12F2and NH4B4O6F is stronger than that of KBBF(1.2 × KDP@1064 nm) under the irradiation of 1064 nm laser.More importantly, the theoretical analysis shows that their shortest second harmonic phasematching wavelength are less than 200 nm.All these properties indicate that these crystals are promising to be the next generation DUV materials.DFT calculation shows that the optical properties of AB4O6F (A = Rb, Cs, K, Na, NH4) family are mainly determined by the anion layer [B4O6F]-∞.

    Fig.7.Crystal structure of AB4O6F (A = Rb, Cs, K, Na, NH4) series compounds.

    3.1.2.M(B5O7)F3 (M=Sr, Ca, Pb)

    Pan‘s team designed and synthesized KBBF-type SrB5O7F3and CaB5O7F3crystals [52,54].SrB5O7F3and CaB5O7F3have a 2D layered structure.The structural units [B5O7F3] are composed of three[BO3F] units and two [BO3] units.These [B5O7F3] units form 2D layer in theacplane by co-apex connection (Fig.8a).The NLO response of SrB5O7F3and CaB5O7F3are 1.6 and 2.0 × KDP@1064 nm, respectively.Their ultraviolet cut-off edges are less than 180 nm.In addition, theoretical calculations show that their shortest second harmonic phasematching wavelengths are 180 nm and 183 nm, respectively.All these properties indicate that SrB5O7F3and CaB5O7F3have the potential to be the next generation DUV materials.According to DFT calculation, the VBs maximum of CaB5O7F3is mainly determined by O 2p and F 2p atomic orbitals,while CBs minimum is mainly determined by Ca 3p and B 2p atomic orbitals.Therefore, the optical properties of CaB5O7F3are determined by the Ca2+and [B5O7F3] units.This is different from AB4O6F (A = Rb, Cs, K, Na, NH4) family compounds.

    Panet al.designed and synthesized a Pb-based PbB5O7F3crystal that is isomorphic to SrB5O7F3[53].Its anionic layered framework is similar to SrB5O7F3, and stereochemically active lonepair (SCALP) cations are connected by Pb-F bonds to form a 3D framework structure (Fig.8b).The distance between adjacent layers (4.38 ?) is significantly smaller than that of KBBF (6.25 ?).PbB5O7F3has a strong SHG effect (6.0 × KDP@1064 nm), showing a PM capability.In addition, the ultraviolet cut-off edge of PbB5O7F3is less than 225 nm, which is significantly less than other Pb-based NLO compounds.According to theoretical calculations, the VBs maximum of PbB5O7F3are predominately determined by O 2p, F 2p, B 2p and Pb 6s, p orbitals.The CBs minimum is mainly determined by O 2p, B 2p, Pb 6s, p atomic orbitals.Obviously, due to the influence of Pb, the band gap of PbB5O7F3is smaller than that of Sr-based compounds.Further theoretical analysis shows that due to the high electronegativity of F and the O-B-F bands, the band gap of PbB5O7F3is larger than other Pbcontaining compounds.The SHG-density method proves that the strong SHG response of PbB5O7F3is mainly derived from the cooperation of [BO3], [BO3F], [PbO6F3] groups.

    Fig.8.Crystal structure and NLO performance of SrB5O7F3 and PbB5O7F3.

    3.1.3.Ba3Mg3(BO3)3F3, CsAlB3O6F

    Mutailipuet al.designed and synthesized polymorphs Ba3Mg3(BO3)3F3by solid-phase synthesis method using SBBO as a template [55].α-andβ-Ba3Mg3(BO3)3F3have similar 2D layered structure.Their crystal structures are connected by[BO3] groups and [MgO4F2] polyhedrons to form 2D layers inbcandabplanes, respectively, and further connected by Mg-F bond to form 3D framework structure (Fig.9a).It is worth noting thatα-Ba3Mg3(BO3)3F3has grown high-quality large-size(16 × 14 × 8 mm3) crystals.Optical experiment tests show that phasematchingα-Ba3Mg3(BO3)3F3has a strong SHG response (1.8 × KDP@1064 nm) and a short ultraviolet cut-off edge(184 nm).All these properties indicate thatα-Ba3Mg3(BO3)3F3is a potential DUV material.

    Liuet al.designed and synthesized KBBF-type CsAlB3O6F crystal through flux method [56].CsAlB3O6F also has a KBBF-type 2D layered structure.Cs+is replaced by K+, and [BeO3F] is replaced by [AlO3F].[AlO3F] and [BO3] groups are connected in thebcplane to form a 2D layer, and Cs+is located in the 18 membered ring tunnel which composed of [AlO3F] and [BO3] units (Fig.9b).Notably, the interlaminar distance of CsAlB3O6F (4.03 ?) is significantly smaller than that of KBBF (6.25 ?), indicating that the interlaminar force of CsAlB3O6F is stronger than that of KBBF.Consequently, the growth characteristic of CsAlB3O6F would be better than that of KBBF.Optical experiments show that phasematching CsAlB3O6F has a strong SHG response (2.0 × KDP@1064 nm) and a short ultraviolet cut-off edge (<190 nm).In addition, theoretical calculation shows that the shortest second harmonic phasematching wavelength is 182 nm.All these properties indicate that CsAlB3O6F is a potential DUV material.According to calculation of density of states, the VBs maximum of CsAlB3O6F is mainly occupied by O 2p and F 2p orbitals, while the CBs minimum is mainly occupied by atomic orbitals of Ca 3p and B 2p.

    Fig.9.Crystal structures and NLO performance of α-Ba3Mg3(BO3)3F3, CsAlB3O6F, LaB4O6(OH)2Cl and Pb2BO3I.

    3.1.4.LaB4O6(OH)2Cl

    Guoet al.designed and synthesized colorless transparent KBBF type LaB4O6(OH)2Cl crystal by hydrothermal method [97].The crystal structure of LaB4O6(OH)2Cl is composed of 18-membered ring [B4O8(OH)2] structural motifs which formed by two [BO3(OH)]tetrahedrons and three [BO3] groups.These structural motifs are connected to each other to form a 2D layer (Fig.9c).Optical experiments show that phasematching LaB4O6(OH)2Cl has a strong frequency doubling effect (2.3 × KDP@1064 nm) and short UV cutoff edge (<180 nm).It can be seen from DFT calculation that the VBs maximum of LaB4O6(OH)2Cl are mainly determined by the atomic orbits of O 2p and Cl 3p, and the CBs minimum are predominantly composed of La 5p and B 2p orbits.

    3.1.5.Pb2BO3I

    Halasyamani and his team designed and synthesized Pb-based KBBF-type Pb2BO3I crystal by hydrothermal method [98].The 3D framework crystal structure of Pb2BO3I is formed by [PbO3I2]and [BO3] groups inabplane.Compared with KBBF, the [PbO3I2]groups in Pb2BO3I replaces [BeO3F] groups, and the K-F bond is replaced by Pb-I bond (Fig.9d).The calculation results show that the interlayer force of Pb2BO3I is stronger than that of KBBF.It should be noted that Pb2BO3I has the strongest SHG efficiency(10.0 × KDP@1064 nm) of KBBF family NLO compounds.It shows type-I phase matching ability.Th optical experiment test proves that the ultraviolet cut-off edge of Pb2BO3I is 330 nm.Pb2BO3I has the potential to become the ultraviolet NLO materials.

    3.1.6.NH4Be2BO3F2, γ-Be2BO3F, C(NH2)3SO3F

    Yeet al.designed and synthesized colorless and transparent KBBF-type NH4Be2BO3F2andγ-Be2BO3F crystals by high temperature hydrothermal method [57].NH4Be2BO3F2andγ-Be2BO3F have similar structure with KBBF, but their interlayer force is stronger than KBBF.Especially forγ-Be2BO3F, the cation at A site is directly removed from its crystal structure, and the layers are directly connected by Be-F-Be bonds, which eliminates the disadvantage of KBBF layer-by-layer growth habit (Figs.10c and d).Under 1064 nm laser irradiation, phase matching NH4Be2BO3F2andγ-Be2BO3F have strong SHG response (1.2 and 2.3 × KDP, respectively.The ultraviolet cut-off edges of these two compounds are less than 190 nm.DFT calculation shows that their shortest second harmonic phase matching wavelengths are 173.9 nm and 146 nm,respectively.All these optical properties indicate that they have potential to become DUV NLO materials.

    Recently, Yeet al.have grown colorless plate-like KBBF type C(NH2)3SO3F crystal in solution [77].[SO3] units replace [BO3]units, and [SO3F] units replace [BeO3F] units from KBBF to C(NH2)3SO3F.The difference is that the [BeO3F] tetrahedrons in KBBF are arranged in opposite intervals, while the [SO3F] groups in C(NH2)3SO3F are arranged in the same direction (Fig.10b), which is in favor of the superposition ofdij.C(NH2)3SO3F exhibits strong SHG response (5 × KDP) and moderate birefringence (0.133).The ultraviolet cut-off edge of C(NH2)3SO3F is 200 nm.Theoretical calculation shows that its shortest second harmonic phase matching wavelength is less than 200 nm.

    Fig.10.Crystal structure of NH4Be2BO3F2, γ-Be2BO3F and C(NH2)3SO3F.

    In general, it is a very effective method to design DUV NLO materials using KBBF as a template.Many crystals in this part inherit the excellent properties of KBBF and overcome some intrinsic performance limitations.In the AB4O6F (A = Rb, Cs, Na, NH4),A(B5O7)F3(A = Sr, Ca, Pb) series compounds, K+is replaced by other alkali metal/alkaline earth metal cations or NH4+, and mixed anion SBUs [BeO3F] are replaced by [BO3F].These compounds exhibit a good balance among short ultraviolet cut-off edge, strong SHG response and suitable phase matching ability.In addition,the substitution of the non-toxic element B with Be and much stronger interlayer forces make these compounds more applicable in practical use.γ-Be2BO3F and C(NH2)3SO3F are potential candidate as DUV NLO materials with strong interlayer forces and excellent comprehensive properties.Pb(B5O7)F3and Pb2BO3I showed strong SHG response with 6.0 and 10.0 × KDP, respectively.

    3.2.Nitrate and carbonate

    In this part, the NLO crystals containing tri-coordinated groups[NO3], [CO3] and complex coordinated groups [MAxBy] are analyzed.The four compounds Pb2(NO3)2(H2O)F2[99], Rb3SbF3(NO3)3[100], KMgCO3F [101] and NaZnCO3F [102] emerged in recent years are represented.

    3.2.1.Pb2(NO3)2(H2O)F2

    Ye and his team designed and synthesized the first fluoronitrate crystal of Pb2(NO3)2(H2O)F2by hydrothermal method [103].The 3D framework crystal structure of Pb2(NO3)2(H2O)F2consists of mixed-anion [PbO9F2] and [NO3] groups (Fig.11a).Optical experimental tests show that phasematching Pb2(NO3)2(H2O)F2has a very strong frequency doubling effect (12.0 × KDP@1064 nm)and a short ultraviolet cut-off edge (300 nm).The DFT calculation shows that the VBs maximum and CBs minimum of Pb2(NO3)2(H2O)F2are predominately determined by the atomic orbitals of O, N, Pb and F.Therefore, the linear and NLO properties of Pb2(NO3)2(H2O)F2are mainly affected by the [NO3] groups, the Pb-O and Pb-F bonds.

    Fig.11.Crystal structures and NLO performance of (a) Pb2(NO3)2(H2O)F2, (b) Rb3SbF3(NO3)3 and (c) KMgCO3F.

    3.2.2.Rb3SbF3(NO3)3

    Wanget al.designed and synthesized Sb-based Rb3SbF3(NO3)3crystal [100].Rb3SbF3(NO3)3has a 0D isolated molecular structure, which consists of two [SbO3F3] polyhedrons and three [NO3]groups with holes hosting Rb+cations (Fig.11b).The NLO response of Rb3SbF3(NO3)3is 2.2 × KDP@1064 nm, and it shows PM behavior.The bandgap of Rb3SbF3(NO3)3is 3.75 eV (corresponding to 331 nm of ultraviolet cut-off edge).Rb3SbF3(NO3)3has the potential to become the next generation UV NLO material.DFT calculation shows that the VBs maximum and CBs minimum of Rb3SbF3(NO3)3are determined by the atomic orbitals of O, N, Sb and F, thus its optical properties are mainly affected by [NO3] and[SbO3F3] groups.

    3.2.3.KMgCO3F, NaZnCO3F

    Halasyamani and Tranet al.designed and synthesized alkali metal fluorocarbonates KMgCO3F [101] and NaZnCO3F [102], respectively.They have a similar crystal structure.The crystal structure of KMgCO3F is composed of mixed anion [MgO3F2] group and[CO3] group to form a 3D framework structure (Fig.11c).The optical experiment test shows that both KMgCO3F and NaZnCO3F have phase matching ability and show a strong frequency doubling effect (3.0, 2.75 × KDP@1064 nm).In addition, they have short ultraviolet cut-off edges (<180 nm, 230 nm).All these properties indicate that KMgCO3F and NaZnCO3F have the potential to become NLO materials for the next generation of UV applications.

    In a word, compounds obtained by the combination of tri-coordinated groups [NO3], [CO3] and complex coordinated groups [MAxBy] are critical for UV NLO materials exploration.All crystals of this type have good phase matching capability.Pb2(NO3)2(H2O)F2had the largest SHG effect, which was 12 × KDP.In addition, Rb3SbF3(NO3)3, KMgCO3F and NaZnCO3F can keep a balance between the strong SHG response and the short ultraviolet cut-off edge, and they are potential candidates for UV NLO crystals.

    3.3.Phosphates

    This part mainly studies the new type of NLO crystals containingp-block coordination tetrahedral group [PO4] and mixed anion group [MAxBy].These compounds can be divided into the following three parts according to the combined tetrahedral units, such as [PO4], [P2O7].The comprehensive performance, including optical band gap, SHG effect and phase matching ability, are evaluated in combination with crystal structure and complex coordinated groups arrangement ([MAxBy] and [PO4] groups).

    3.3.1.BaZn(PO4)F, RbBPO4F, (NH4)2BPO4F2

    Tanget al.synthesized the fluorophosphate BaZn(PO4)F crystal by hydrothermal method [103].The crystal structure of BaZn(PO4)F consists of the interconnection of [PO4] groups and [ZnO4F] mixedanion polyhedrons to form a 3D framework structure (Fig.12a).Optical experiments show that the ultraviolet cut-off edge of BaZn(PO4)F is less than 190 nm, and thedijis 0.26 × KDP, showing type-I PM capability.DFT calculation indicates that the VBs maximum and CBs minimum of BaZn(PO4)F are mainly determined by the atomic orbitals of Zn, P, O and F, which indicates that the effect of alkaline earth metal element Ba on its optical properties is negligible.

    Subsequently, Wuet al.synthesized fluoroborophosphates RbBPO4F and (NH4)2BPO4F2crystals by hydrothermal method[104].In the crystal structure of RbBPO4F, the mixed-anion [BO3F]and [PO4] groups are connected to each other to form a 3D tunnel structure.The alkali metal Rb+is located in the tunnel (Fig.12b).As shown in Fig.12c, the 1D chain structure of (NH4)2BPO4F2is composed of [PO4] tetrahedrons and complex coordinated [BO2F2]groups.Optical experimental tests show that the ultraviolet cut-off edges of RbBPO4F and (NH4)2BPO4F2are both less than 200 nm,and their SHG effects are 0.4, 0.8 × KDP@1064 nm, respectively.RbBPO4F exhibits non-phasematching behavior due to its cubic structure, while (NH4)2BPO4F exhibits type-I phase matching ability.Theoretical calculations reveal that the NLO response of these two compounds mainly comes from the [PO4]+groups.

    Fig.12.Crystal structures and NLO performance of (a) BaZn(PO4)F, (b) RbBPO4F and (c) (NH4)2BPO4F2.

    3.3.2.Ba2NaP2O7Cl, CsSiP2O7F, K2Sb(P2O7)F

    Chenet al.synthesized colorless and transparent fresnoitelike structure Ba2NaP2O7Cl crystal by a high temperature solidstate reaction technology [59].Ba2NaP2O7Cl is a 3D tunnel structure formed by the connection of [P2O7] dimer groups and[NaO4Cl2] polyhedra, and Ba2+is located in the center of the pore (Fig.13a).Ba2NaP2O7Cl can maintain a good balance among short ultraviolet cut-off edge (<176 nm), strong SHG response(0.9 × KDP@1064 nm), and sufficient birefringence value (0.017).All these properties indicate that Ba2NaP2O7Cl has the potential to be DUV NLO materials.

    Dinget al.synthesized CsSiP2O7F crystal by high temperature solid-state reaction technology [60].As shown in Fig.13b, the 3D frame structure of CsSiP2O7F is composed of the [P2O7] dimer unit and [SiO5F] polyhedrons, and Cs+are scattered in the tunnel holes.CsSiP2O7F can also maintain a good balance between short cut-off edge (<190 nm), strong SHG effect (0.7 × KDP@1064 nm), and good phase matching ability.DFT calculation shows that the VBs maximum of CsSiP2O7F is mainly occupied by the atomic orbitals of Cs, P and Si, and the CBs minimum is occupied by the atomic orbitals of O and F, which indicates that all atoms have an effect on the linear optical properties.The SHG-density method proves that the anions [P2O7] dimer units and [SiO5F] groups make major contributions to whole SHG effect.

    Fig.13.Crystal structures and NLO performance of (a) Ba2NaP2O7Cl, (b) CsSiP2O7F and (c) K2Sb(P2O7)F.

    Denget al.synthesized colorless and transparent K2Sb(P2O7)F crystal by means of solvent-free synthesis technology [105].K2Sb(P2O7)F has a 2D layered structure, which is built by [P2O7]dimer and [SbO4F] polyhedra along theabplane (Fig.13c).Optical experiments show that K2Sb(P2O7)F has a strong frequency doubling effect (4 × KDP@1064 nm), showing the PM behavior.The ultraviolet cut-off edge of K2Sb(P2O7)F is 261 nm.DFT calculations show that VBs maximum of K2Sb(P2O7)F is mainly occupied by atomic orbitals of F, P and O, and the CBs minimum is mainly occupied by atomic orbitals of Sb and F, which indicates that the effect of alkali metal K on its optical properties is minimal.

    3.3.3.Na3Sc2(PO4)2F3

    Xuet al.synthesized rare earth fluorophosphate Na3Sc2(PO4)2F3under hydrothermal conditions [106].The 3D framework structure of Na3Sc2(PO4)2F3is composed of [PO4] tetrahedrons and[ScO4F2] mixed-anion polyhedrons.Na+is scattered in the tunnel holes (Fig.14).The SHG response of Na3Sc2(PO4)2F3is 0.26 × KDP under 1064 nm laser irradiation, and it has type-I phase matching ability.Na3Sc2(PO4)2F3exhibits short cutoff edge of 255 nm and a large birefringence value (0.0978).All these properties indicate that Na3Sc2(PO4)2F3has the potential to be a UV NLO material.

    Fig.14.Crystal structure and NLO performance of Na3Sc2(PO4)2F3.

    To summarize, compounds containing mixed anion groups[MAxBy] and [PO4] are important candidates for UV/DUV NLO materials.This type of compound has strong optical anisotropy and exhibits good phase matching ability.It is worth noting that when the mixed anion groups [MAxBy] and [PO4] are arranged in 2D layers, it is easier to maintain a balance between strong SHG response and the short ultraviolet cut-off edge.For example, Ba2NaP2O7Cl,CsSiP2O7F, K2Sb(P2O7)F all show excellent comprehensive performance.

    3.4.Sulfate

    In addition to phosphates containing [PO4] tetrahedrons, metal sulfates are also well-known UV/DUV NLO materials.In this part,CsSbF2SO4[107], NH4SbF2SO4[108] and RbSbSO4Cl2[109] are taken as representative compounds to analyze the influence of the binding and arrangement manner of [SO4] tetra-coordination groups and mixed-anion group [MAxBy] on optical band gap, SHG effect and phase matching ability.

    3.4.1.CsSbF2SO4, NH4SbF2SO4

    Donget al.designed and synthesized CsSbF2SO4[107] by electrothermal synthesis.The 1D chain structure of CsSbF2SO4is composed of twisted [SbO2F2] groups and [SO4] tetrahedral groups.Cs+is distributed among the chains (Fig.15a).The ultraviolet cutoff edge of CsSbF2SO4is 240 nm.Under 1064 nm laser irradiation, it exhibits strong frequency doubling effect (3 × KDP) with phase matching behavior.The theoretical calculation indicates that the NLO response of CsSbF2SO4mainly comes from the highly distorted [SbO2F2] units, while the [SO4] group is beneficial to the parallel arrangement of [SbO2F2] units.

    Fig.15.Crystal structures and NLO performance of CsSbF2SO4, NH4SbF2SO4 and RbSbSO4Cl2.

    Jantzet al.designed and synthesized NCS ammonia fluorine chalcogenide salt NH4SbF2SO4[108].The 1D chain structure of NH4SbF2SO4is formed by connecting twisted [SbO2F2] groups and [SO4] tetrahedral groups, with NH4+distributed between the chains (Fig.15b).The optical experiment test shows that the ultraviolet cut-off edge of NH4SbF2SO4is 266 nm.It has PM capability,and the SHG effect is about 0.7 × KDP.

    3.4.2.RbSbSO4Cl2

    Heet al.synthesized RbSbSO4Cl2by hydrothermal method[109].RbSbSO4Cl2has a 1D chain structure, which is composed of twisted [SbO3Cl2] polyhedrons and [SO4] units, while Rb+is scattered among the chains (Fig.15c).The ultraviolet cut-off edge of RbSbSO4Cl2is 347 nm.Under 1064 nm laser irradiation, the frequency doubling effect is 2.7 × KDP with type-I PM capability.RbSbSO4Cl2has the potential to become UV NLO materials.Theoretical calculations show that the VBs maximum and CBs minimum of RbSbSO4Cl2are determined by the atomic orbitals of Rb, O, Sb,S and Cl.Therefore, the optical properties of RbSbSO4Cl2mainly depend on [RbO5Cl4], [SbO3Cl2], [SO4] groups.

    3.5.Complex coordinated group [BO3F] as NLO active units

    Sn[B2O3F2]:H?PPE and his team synthesized colorless massive Sn[B2O3F2] crystals [109].The power X-ray diffraction test shows that Sn[B2O3F2] crystallizes in space groupP31m (No.157).The crystal structure of Sn[B2O3F2] is composed of two-dimensional layers connected by the [BO3F] cells along theabplane, and Sn2+is distributed among the layers to maintain the charge balance(Fig.16a).The UV cutoff edge of Sn[B2O3F2] is 243 nm.Under 1054 nm laser irradiation, its SHG effect is even stronger that of KDP, and it has phase matching ability.Sn[B2O3F2] has the potential to become the next generation of UV NLO materials.Theoretical calculations show that the VBs maximum and CBs minimum of Sn[B2O3F2] are mainly determined by the atomic orbitals of Rb,O, Sb, S and Cl.Its optical properties mainly depend on [RbO5Cl4],[SbO3Cl2] and [SO4] groups.

    3.6.Mixed-anion group [TeO2F2] as NLO active units

    BaF2TeF2(OH)2:In 2020, Heet al.designed and synthesized colorless and transparent BaF2TeF2(OH)2crystal by hydrothermal method [110].Power X-ray diffraction results show that BaF2TeF2(OH)2belongs to NCS space groupPmn21(No.31).BaF2TeF2(OH)2has a typical aurivillius structure.The [TeF2(OH)2]groups are connected with each other in theacplane to form a two-dimensional layer, and the [BaF4] tetrahedrons are also connected by common edge to form 2D layers.The crystal structure of BaF2TeF2(OH)2is composed of the two layers which are mutually spaced along theb-axis (Fig.16b).The UV cut-off edge of BaF2TeF2(OH)2is 205 nm.Under 1064 nm laser irradiation, it shows strong SHG effect (3 × KDP) and good phase matching ability.All these properties indicate that BaF2TeF2(OH)2has the potential to become the next generation of UV NLO materials.

    3.7.Mixed-anion group [PO3F] as NLO active units

    NaNH4PO3F·H2O:Luet al.designed and synthesized colorless and transparent NaNH4PO3F·H2O crystal [61].Power x-ray diffraction results show that NaNH4PO3F·H2O belongs to monoclinic space groupPc(No.7).The crystal structure of NaNH4PO3F·H2O is composed of [PO3F] groups, which are connected with [NH4] and H2O molecules through hydrogen bonds, and finally form a twodimensional layered structure, with Na+scattered among the layers (Fig.16c).NaNH4PO3F·H2O exhibits a strong SHG effect of 1.1 times that of KDP under 1064 nm laser irradiation.Theoretical calculation shows that its shortest second harmonic phase matching wavelength is 194 nm.There is no doubt that NaNH4PO3F·H2O has the potential to apply in the next generation of DUV NLO materials.

    Fig.16.Crystal structures and NLO performance of Sn[B2O3F2], NaNH4PO3F·H2O and BaF2TeF2(OH)2.

    4.Summary and outlook

    To sum up, complex coordinated functional groups [MAxBy] are excellent NLO active units, and compounds containing this group are more likely to guarantee excellent comprehensive performance.This group is combined with traditional NLO active group ([MQ4](M = metal cationic, Q = S,Se, Te, P), [IO3], [SeO3], [TeO3], [CO3],[NO3], [SO4], [PO4], B-O groups,etc.) is an important direction to explore new high-performance nonlinear optical materials.In this paper, a comprehensive review of NLO crystal materials with mixed anion functional groups [MAxBy] as the basic structural unit is made based on theoretical analysis results, anion combination methods, and available experimental performance test results.Through systematic summary, it can be concluded that:(1) The combination of complex coordinated groups [MAxBy] and [C(1)O3](C(1) = I, Se, Te) oxygen-containing groups is conducive to the synthesis of large band gap mid-IR NLO crystals with good phase matching ability.In order to solve the shortcoming of small NLO effect of this kind of compound, it is necessary to introduce d0-transition metal element,e.g., W6+, V5+, or metal element containing stereochemically active lone pairs, which is easy to produce strong distorted multi-coordination groupse.g., Bi and Pb.(2)The combination of the complex coordinated groups [MAxBy] and[MQ4] (Q = S,P) tetra-coordinated units to form a 3D frame structure is conducive to the synthesis of ideal MIR NLO crystals with phase matching ability, large band gap and strong SHG effect.(3)The comprehensive performance of IR NLO crystals with complex coordinated groups [MAxBy] as SBUs is improved compared with that of traditional single anionic compounds.(4) When the mixed anion functional groups [MOxF4-x] and [C(2)O3] (C(2) = B, N, C) tricoordinated oxygen-containing groups are arranged in a KBBF-like structure model (three coordination groups and four coordination groups containing halogen F are connected to form a 2D layer),it is beneficial to the formation of excellent deep ultraviolet NLO crystals (short ultraviolet cut-off edge, strong SHG effect, small second harmonic phasematching wavelength).(5) When the complex coordinated functional groups [MAxBy] and [C(3)O4] (C(3) = P,S)tetra-coordination groups are arranged in 2D layers, and there are scattered alkali/alkaline earth metal cations in the structure, it is more likely to achieve a balance between short cut-off edge and strong SHG effect.(6) UV/DUV NLO crystals composed of a 2D layered arrangement of complex coordinated groups [MAxBy] exhibit excellent comprehensive properties.It is believed that the combination mode of the anion NLO active units and the regulation of the arrangement rules on the optical properties above will be of guiding significance for the future design of MIR/DUV NLO materials.

    For the future design and synthesis of complex coordinated NLO crystals, we should focus on the following aspects:(1) For the field of IR nonlinear optics, the research of tetra-coordinated compounds is relatively less, and we should pay more attention to the diamond-like compounds composed of [MAxB4-x] (A, B is one of the elements of P, O, S, Se, F, Cl, Br, I) with four coordination groups arranged in parallel.(2) For the field of UV/deep UV nonlinear optics, more alternative structural templates similar to KBBF should be explored.In addition, at present, there are only a few compounds in type IV NLO materials with only [MAxB4-x] as SBUs,which has a lot of exploration potential in the future.(3) The practical application of NLO materials needs large single crystals with good quality, so the best growth conditions of NLO crystals with mixed anion group [MAxBy] as SBUs should be explored in the future.(4) The design and synthesis of [MAxB4-x] as SBUs NLO crystal is a time-consuming and laborious work.In order to improve the efficiency, it is necessary to combine theoretical analysis with experimental exploration.

    Declaration of competing interest

    The authors declare no conflict of interest.

    Acknowledgment

    This work was supported by National Natural Science Foundation of China (No.51972208).

    精品高清国产在线一区| 午夜福利视频精品| 国产精品99久久99久久久不卡| 18禁观看日本| 久久久欧美国产精品| 国产三级黄色录像| 亚洲人成电影免费在线| av天堂久久9| 一本久久精品| 欧美av亚洲av综合av国产av| 女性生殖器流出的白浆| 欧美亚洲 丝袜 人妻 在线| 精品久久久久久电影网| 18禁黄网站禁片午夜丰满| 1024视频免费在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品一二三区在线看| 又紧又爽又黄一区二区| 狂野欧美激情性bbbbbb| 国产亚洲av高清不卡| 一区二区三区激情视频| 欧美+亚洲+日韩+国产| 国产在线观看jvid| 宅男免费午夜| 欧美激情极品国产一区二区三区| 99精国产麻豆久久婷婷| 母亲3免费完整高清在线观看| 日本a在线网址| 黄片小视频在线播放| videosex国产| 麻豆av在线久日| 欧美精品一区二区免费开放| 日韩一卡2卡3卡4卡2021年| 美女扒开内裤让男人捅视频| 黄色怎么调成土黄色| 欧美日韩av久久| 亚洲精品一区蜜桃| 男女无遮挡免费网站观看| 在线精品无人区一区二区三| 99国产精品免费福利视频| 国产伦理片在线播放av一区| 欧美日韩av久久| 丰满饥渴人妻一区二区三| 麻豆乱淫一区二区| 黑人巨大精品欧美一区二区蜜桃| videos熟女内射| av欧美777| 黄色视频不卡| 亚洲七黄色美女视频| 亚洲美女黄色视频免费看| 最新在线观看一区二区三区 | 一级毛片 在线播放| 丰满人妻熟妇乱又伦精品不卡| 久久人人爽av亚洲精品天堂| 亚洲欧美清纯卡通| 午夜激情av网站| 国产亚洲av高清不卡| 热re99久久国产66热| 亚洲自偷自拍图片 自拍| 看免费av毛片| 成年人午夜在线观看视频| 久久久久久人人人人人| 亚洲国产毛片av蜜桃av| 亚洲国产看品久久| 色94色欧美一区二区| 亚洲av日韩在线播放| 国产精品人妻久久久影院| 亚洲欧美清纯卡通| 国产av一区二区精品久久| 欧美精品人与动牲交sv欧美| 人人妻人人添人人爽欧美一区卜| 嫩草影视91久久| 夫妻性生交免费视频一级片| 国产精品亚洲av一区麻豆| 真人做人爱边吃奶动态| 波多野结衣av一区二区av| 欧美日韩亚洲综合一区二区三区_| 国产精品秋霞免费鲁丝片| 亚洲,一卡二卡三卡| 欧美日韩亚洲高清精品| 黄色视频不卡| 久久久久国产一级毛片高清牌| 又黄又粗又硬又大视频| 少妇被粗大的猛进出69影院| 18在线观看网站| 中国国产av一级| a级毛片黄视频| 精品一品国产午夜福利视频| 91国产中文字幕| 欧美日韩视频高清一区二区三区二| 免费看av在线观看网站| 19禁男女啪啪无遮挡网站| 欧美日韩一级在线毛片| 美女福利国产在线| 久久狼人影院| 香蕉国产在线看| 脱女人内裤的视频| 日本五十路高清| 精品国产乱码久久久久久男人| 亚洲精品乱久久久久久| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品 国内视频| 午夜福利视频在线观看免费| 免费少妇av软件| 国产成人一区二区在线| 99精国产麻豆久久婷婷| 亚洲人成电影观看| 国产精品 欧美亚洲| 日本欧美国产在线视频| 中文字幕制服av| 欧美在线黄色| 19禁男女啪啪无遮挡网站| 欧美精品一区二区大全| 成人黄色视频免费在线看| av网站在线播放免费| 色视频在线一区二区三区| 婷婷色综合www| 日韩中文字幕欧美一区二区 | 免费看十八禁软件| 午夜激情久久久久久久| 亚洲成人手机| 亚洲伊人色综图| 51午夜福利影视在线观看| 亚洲国产中文字幕在线视频| 99九九在线精品视频| 亚洲专区中文字幕在线| 欧美+亚洲+日韩+国产| 亚洲精品国产色婷婷电影| 又大又黄又爽视频免费| 成人国产一区最新在线观看 | 久久精品人人爽人人爽视色| 国产伦人伦偷精品视频| 国产有黄有色有爽视频| 一区二区三区激情视频| 女人被躁到高潮嗷嗷叫费观| 99久久人妻综合| 十分钟在线观看高清视频www| 老熟女久久久| 美女午夜性视频免费| 亚洲精品美女久久av网站| 狂野欧美激情性bbbbbb| 国产av国产精品国产| 婷婷色综合www| 建设人人有责人人尽责人人享有的| 亚洲欧美清纯卡通| 最新在线观看一区二区三区 | 亚洲精品国产一区二区精华液| 日日爽夜夜爽网站| 手机成人av网站| av线在线观看网站| 欧美激情极品国产一区二区三区| 亚洲天堂av无毛| 久久午夜综合久久蜜桃| 亚洲成人手机| 日韩大片免费观看网站| 50天的宝宝边吃奶边哭怎么回事| 99久久综合免费| 久久精品国产a三级三级三级| 性高湖久久久久久久久免费观看| 亚洲av成人精品一二三区| 夫妻午夜视频| 9热在线视频观看99| av欧美777| 少妇裸体淫交视频免费看高清 | 日韩av不卡免费在线播放| 国产精品免费大片| 欧美国产精品va在线观看不卡| 女人高潮潮喷娇喘18禁视频| 美女高潮到喷水免费观看| 国产精品一国产av| 国产亚洲午夜精品一区二区久久| a级毛片黄视频| 欧美日韩av久久| 欧美变态另类bdsm刘玥| 精品国产一区二区三区久久久樱花| 国产97色在线日韩免费| 免费日韩欧美在线观看| 新久久久久国产一级毛片| 一级毛片黄色毛片免费观看视频| 91国产中文字幕| 亚洲国产精品一区三区| 亚洲欧美中文字幕日韩二区| 一区二区日韩欧美中文字幕| 久久免费观看电影| 国产精品三级大全| 久久久久国产精品人妻一区二区| 十八禁高潮呻吟视频| av有码第一页| 日韩大片免费观看网站| 91老司机精品| 精品少妇一区二区三区视频日本电影| 无遮挡黄片免费观看| 少妇裸体淫交视频免费看高清 | 亚洲中文日韩欧美视频| 啦啦啦啦在线视频资源| 国产在线视频一区二区| 国产成人影院久久av| 熟女少妇亚洲综合色aaa.| 少妇猛男粗大的猛烈进出视频| av不卡在线播放| av天堂在线播放| 亚洲成色77777| 日本vs欧美在线观看视频| 在线观看免费高清a一片| 日韩 亚洲 欧美在线| 99久久人妻综合| 国产欧美亚洲国产| 天天添夜夜摸| 久久人妻熟女aⅴ| 亚洲欧洲国产日韩| 一级毛片女人18水好多 | 90打野战视频偷拍视频| cao死你这个sao货| e午夜精品久久久久久久| 脱女人内裤的视频| 又黄又粗又硬又大视频| 亚洲欧美日韩另类电影网站| av在线app专区| 亚洲人成77777在线视频| 大陆偷拍与自拍| 国产成人精品在线电影| 国产精品久久久久久人妻精品电影 | 久久人妻熟女aⅴ| 亚洲综合色网址| 亚洲欧美清纯卡通| 精品国产一区二区三区久久久樱花| 亚洲精品第二区| 亚洲av在线观看美女高潮| 丁香六月天网| 美女福利国产在线| 亚洲欧美清纯卡通| 黑人巨大精品欧美一区二区蜜桃| 亚洲av美国av| 成在线人永久免费视频| av有码第一页| 国产精品偷伦视频观看了| 久久99精品国语久久久| 国语对白做爰xxxⅹ性视频网站| av网站免费在线观看视频| 日韩大片免费观看网站| 在线 av 中文字幕| 巨乳人妻的诱惑在线观看| 国产精品二区激情视频| 69精品国产乱码久久久| 99国产精品99久久久久| 老熟女久久久| 免费女性裸体啪啪无遮挡网站| 日韩一卡2卡3卡4卡2021年| 欧美乱码精品一区二区三区| 精品人妻一区二区三区麻豆| 大型av网站在线播放| 国产精品三级大全| 日韩一区二区三区影片| 在线天堂中文资源库| 老司机午夜十八禁免费视频| av电影中文网址| av在线app专区| 天堂中文最新版在线下载| 中文精品一卡2卡3卡4更新| 成年人午夜在线观看视频| 一区二区日韩欧美中文字幕| 超色免费av| 日韩中文字幕视频在线看片| 亚洲国产欧美一区二区综合| 自拍欧美九色日韩亚洲蝌蚪91| 热99国产精品久久久久久7| 国产激情久久老熟女| 国产一区二区在线观看av| 国产精品熟女久久久久浪| 午夜两性在线视频| 九草在线视频观看| 首页视频小说图片口味搜索 | 考比视频在线观看| 在线观看免费日韩欧美大片| 啦啦啦在线免费观看视频4| 婷婷色av中文字幕| 国产亚洲欧美在线一区二区| 大片电影免费在线观看免费| 精品国产超薄肉色丝袜足j| 精品久久久久久久毛片微露脸 | 国产精品偷伦视频观看了| 久久国产精品人妻蜜桃| 国产无遮挡羞羞视频在线观看| 女人久久www免费人成看片| 亚洲精品美女久久av网站| 人人妻,人人澡人人爽秒播 | 狂野欧美激情性bbbbbb| avwww免费| 成人手机av| 亚洲 欧美一区二区三区| 少妇粗大呻吟视频| 国产欧美亚洲国产| 欧美日韩成人在线一区二区| 人人妻人人澡人人看| 老司机在亚洲福利影院| 免费在线观看完整版高清| 一级毛片 在线播放| 久久久久国产精品人妻一区二区| 女人久久www免费人成看片| 中文欧美无线码| 中文字幕精品免费在线观看视频| 少妇被粗大的猛进出69影院| 国产麻豆69| 一区二区三区激情视频| 大陆偷拍与自拍| 色视频在线一区二区三区| 精品一区二区三卡| 国产一区二区三区av在线| 国产成人精品久久二区二区免费| 两个人免费观看高清视频| 男女下面插进去视频免费观看| 这个男人来自地球电影免费观看| 黑人欧美特级aaaaaa片| 亚洲一区中文字幕在线| 亚洲第一青青草原| 18禁裸乳无遮挡动漫免费视频| 一级黄片播放器| 国产一区二区激情短视频 | 十八禁人妻一区二区| 欧美大码av| 男人舔女人的私密视频| 国产精品久久久久久人妻精品电影 | 丝瓜视频免费看黄片| 亚洲激情五月婷婷啪啪| 中文字幕高清在线视频| 成人午夜精彩视频在线观看| av片东京热男人的天堂| 亚洲熟女毛片儿| 少妇裸体淫交视频免费看高清 | 亚洲男人天堂网一区| 成人国语在线视频| 亚洲精品在线美女| 中文字幕人妻丝袜制服| 伦理电影免费视频| 亚洲第一av免费看| 久久性视频一级片| 午夜福利视频精品| 久久久国产欧美日韩av| 男女下面插进去视频免费观看| 国产亚洲午夜精品一区二区久久| 一级黄片播放器| 女人久久www免费人成看片| 国产97色在线日韩免费| 亚洲自偷自拍图片 自拍| 男人舔女人的私密视频| 考比视频在线观看| 美女脱内裤让男人舔精品视频| 又大又黄又爽视频免费| 国产99久久九九免费精品| 久久精品久久久久久久性| 日韩一区二区三区影片| 性色av一级| 久久久久久免费高清国产稀缺| 亚洲欧美色中文字幕在线| 国产在线视频一区二区| 母亲3免费完整高清在线观看| 亚洲欧美日韩高清在线视频 | 丁香六月欧美| 丝袜喷水一区| 日韩一本色道免费dvd| 久久青草综合色| 国产真人三级小视频在线观看| 最新的欧美精品一区二区| 欧美久久黑人一区二区| 久久久精品区二区三区| 男女下面插进去视频免费观看| 视频区图区小说| 99久久综合免费| 一区福利在线观看| 欧美日韩一级在线毛片| 不卡av一区二区三区| 高清视频免费观看一区二区| 中文字幕色久视频| 久热爱精品视频在线9| 中文字幕人妻丝袜一区二区| 只有这里有精品99| 亚洲精品日本国产第一区| 国产精品一区二区在线观看99| 少妇的丰满在线观看| 丝袜在线中文字幕| 岛国毛片在线播放| 日韩av免费高清视频| 久久久久久人人人人人| 日韩人妻精品一区2区三区| 又大又黄又爽视频免费| 啦啦啦中文免费视频观看日本| 一边摸一边抽搐一进一出视频| 狠狠精品人妻久久久久久综合| h视频一区二区三区| 亚洲中文av在线| 国产黄色免费在线视频| 女人精品久久久久毛片| 黄色怎么调成土黄色| 亚洲欧美中文字幕日韩二区| 老熟女久久久| 欧美黑人欧美精品刺激| 51午夜福利影视在线观看| 美女主播在线视频| 男女高潮啪啪啪动态图| 最近手机中文字幕大全| 欧美日韩亚洲高清精品| 久久久欧美国产精品| 黄频高清免费视频| 国产成人av教育| 制服人妻中文乱码| 美女国产高潮福利片在线看| av在线app专区| 黄色视频不卡| 免费人妻精品一区二区三区视频| av又黄又爽大尺度在线免费看| 高清不卡的av网站| 亚洲三区欧美一区| 久久亚洲精品不卡| 亚洲九九香蕉| 啦啦啦在线观看免费高清www| 91字幕亚洲| 久久天躁狠狠躁夜夜2o2o | 免费少妇av软件| 国产男女超爽视频在线观看| 久久九九热精品免费| 欧美日韩亚洲国产一区二区在线观看 | 亚洲三区欧美一区| 一区二区日韩欧美中文字幕| 精品免费久久久久久久清纯 | 性色av乱码一区二区三区2| 最黄视频免费看| 我的亚洲天堂| 亚洲国产看品久久| 天天添夜夜摸| 亚洲av成人精品一二三区| 日本av手机在线免费观看| 99热全是精品| 在线观看免费视频网站a站| 久久精品国产亚洲av高清一级| 免费不卡黄色视频| 精品高清国产在线一区| 久久久久久久大尺度免费视频| 欧美黄色片欧美黄色片| 中文欧美无线码| 亚洲精品久久午夜乱码| 另类精品久久| 国产在线视频一区二区| 成年人黄色毛片网站| 丰满迷人的少妇在线观看| 精品少妇一区二区三区视频日本电影| 在线观看免费视频网站a站| 久久99一区二区三区| 亚洲成国产人片在线观看| 精品国产国语对白av| www.自偷自拍.com| 成人亚洲精品一区在线观看| bbb黄色大片| 中文字幕高清在线视频| www.自偷自拍.com| 午夜福利免费观看在线| 日韩伦理黄色片| 国产99久久九九免费精品| 每晚都被弄得嗷嗷叫到高潮| 日本黄色日本黄色录像| 亚洲美女黄色视频免费看| 久久人人爽av亚洲精品天堂| 国产精品国产三级国产专区5o| 久久天堂一区二区三区四区| 丰满人妻熟妇乱又伦精品不卡| 欧美精品亚洲一区二区| 午夜福利乱码中文字幕| 1024视频免费在线观看| 制服人妻中文乱码| 免费看av在线观看网站| 捣出白浆h1v1| 欧美日韩一级在线毛片| 中文字幕高清在线视频| 涩涩av久久男人的天堂| 欧美xxⅹ黑人| 国产亚洲午夜精品一区二区久久| 国产精品免费视频内射| 亚洲人成网站在线观看播放| 久久99精品国语久久久| 日本av手机在线免费观看| 男人舔女人的私密视频| 成人亚洲欧美一区二区av| 国产精品一区二区在线不卡| 在线观看免费视频网站a站| 国产成人一区二区三区免费视频网站 | 精品国产一区二区久久| 丝袜脚勾引网站| 99国产综合亚洲精品| 精品卡一卡二卡四卡免费| 19禁男女啪啪无遮挡网站| 久久ye,这里只有精品| 国产成人系列免费观看| 女人精品久久久久毛片| 无遮挡黄片免费观看| 首页视频小说图片口味搜索 | 亚洲精品乱久久久久久| 一二三四在线观看免费中文在| 交换朋友夫妻互换小说| 在现免费观看毛片| 欧美激情高清一区二区三区| 亚洲一区中文字幕在线| 看十八女毛片水多多多| 色视频在线一区二区三区| 亚洲国产日韩一区二区| 丁香六月欧美| 国产男女超爽视频在线观看| 精品卡一卡二卡四卡免费| 波野结衣二区三区在线| 日本av手机在线免费观看| 国产精品香港三级国产av潘金莲 | 飞空精品影院首页| 久久精品国产亚洲av高清一级| 亚洲成人手机| 国产三级黄色录像| 男女国产视频网站| 久久亚洲国产成人精品v| 黑人欧美特级aaaaaa片| 久久av网站| 天堂8中文在线网| 午夜福利,免费看| 精品一区二区三区av网在线观看 | 性少妇av在线| 91字幕亚洲| bbb黄色大片| 少妇粗大呻吟视频| 最黄视频免费看| 老熟女久久久| 亚洲国产欧美在线一区| 亚洲视频免费观看视频| 亚洲国产精品一区三区| 1024香蕉在线观看| 欧美另类一区| 国产午夜精品一二区理论片| 女人精品久久久久毛片| 国产精品麻豆人妻色哟哟久久| 在线av久久热| 在线观看一区二区三区激情| 国产精品亚洲av一区麻豆| 亚洲国产毛片av蜜桃av| 各种免费的搞黄视频| 午夜福利视频在线观看免费| 女人精品久久久久毛片| 国产亚洲av片在线观看秒播厂| 七月丁香在线播放| 美女午夜性视频免费| 日韩av免费高清视频| 国产午夜精品一二区理论片| 两个人看的免费小视频| 深夜精品福利| 黄片播放在线免费| 久久久久精品人妻al黑| 美女扒开内裤让男人捅视频| 精品国产一区二区三区四区第35| 亚洲欧洲日产国产| 午夜福利乱码中文字幕| 午夜久久久在线观看| 国产主播在线观看一区二区 | 欧美日韩黄片免| 丝袜人妻中文字幕| 午夜激情av网站| 亚洲一区二区三区欧美精品| 一边摸一边做爽爽视频免费| 中文欧美无线码| 少妇粗大呻吟视频| 一区二区av电影网| 国产野战对白在线观看| 爱豆传媒免费全集在线观看| 精品国产一区二区三区四区第35| 美女扒开内裤让男人捅视频| 精品人妻在线不人妻| 亚洲av电影在线观看一区二区三区| 1024香蕉在线观看| 国产99久久九九免费精品| cao死你这个sao货| 老熟女久久久| 丰满迷人的少妇在线观看| 天堂中文最新版在线下载| 久久久久久久大尺度免费视频| 国产亚洲午夜精品一区二区久久| 老司机影院毛片| 国产日韩欧美视频二区| 日韩中文字幕视频在线看片| 岛国毛片在线播放| 亚洲精品成人av观看孕妇| e午夜精品久久久久久久| 中文字幕人妻丝袜一区二区| 亚洲午夜精品一区,二区,三区| 99热国产这里只有精品6| 免费看不卡的av| 午夜福利影视在线免费观看| 欧美人与性动交α欧美精品济南到| 青春草视频在线免费观看| 99久久综合免费| 青草久久国产| 久久狼人影院| netflix在线观看网站| 只有这里有精品99| 自线自在国产av| 日韩中文字幕视频在线看片| 国产日韩欧美亚洲二区| 欧美激情极品国产一区二区三区| 久久精品国产综合久久久| 天堂中文最新版在线下载| 一边摸一边做爽爽视频免费| 视频区图区小说| 国产精品久久久av美女十八| 欧美精品av麻豆av| 女人久久www免费人成看片| 一本久久精品| 国产精品 国内视频| 悠悠久久av| 少妇人妻久久综合中文| 精品一品国产午夜福利视频| 久久久久久久久免费视频了| 大陆偷拍与自拍| 亚洲专区国产一区二区|