Linan
Wang‡
ab,
Dongdong
Chu‡
ab,
Zhihua
Yang
ab,
Junjie
Li
*ab and
Shilie
Pan
*ab
aResearch Center for Crystal Materials, State Key Laboratory of Functional Materials and Devices for Special Environmental Conditions, Xinjiang Key Laboratory of Functional Crystal Materials, Xinjiang Technical Institute of Physics & Chemistry, CAS, 40-1 South Beijing Road, Urumqi 830011, China. E-mail: lijunjie@ms.xjb.ac.cn; slpan@ms.xjb.ac.cn
bCenter of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
First published on 2nd April 2024
Infrared (IR) nonlinear optical (NLO) materials with strong NLO response, wide band gap and high laser-induced damage threshold (LIDT) are highly expected in current laser technologies. Herein, by introducing double alkaline-earth metal (AEM) atoms, three wide band gap selenide IR NLO materials AIIMg6Ga6Se16 (AII = Ca, Sr, Ba) with excellent linear and NLO optical properties have been rationally designed and fabricated. AIIMg6Ga6Se16 (AII = Ca, Sr, Ba) are composed of unique [AIISe6] triangular prisms, [MgSe6] octahedra and [GaSe4] tetrahedra. The introduction of double AEMs effectively broadens the band gaps of selenide-based IR NLO materials. Among them, CaMg6Ga6Se16, achieving the best balance between the second-harmonic generation response (∼1.5 × AgGaS2), wide band gap (2.71 eV), high LIDT (∼9 × AgGaS2), and moderate birefringence of 0.052 @ 1064 nm, is a promising NLO candidate for high power IR laser. Theoretical calculations indicate that the NLO responses and band gaps among the three compounds are mainly determined by the NLO-active [GaSe4] units. The results enrich the chemical diversity of chalcogenides, and give some insight into the design of new functional materials based on the rare [AIISe6] prismatic units.
For an excellent IR NLO material, the following requirements should be satisfied:22 (1) noncentrosymmetric (NCS) crystal structure; (2) a large NLO response (≥0.5 × AGS, preferably ≥1.0 × AGS); (3) a wide IR transparent range that covers the two important IR atmospheric windows of 3–5 and 8–12 μm; (4) a high LIDT, which is usually related to the band gap and thermal conductivity; (5) suitable birefringence for phase-matching (PM) behavior; (6) good crystal growth habits. To achieve the large NLO response and wide IR transparent range, metal selenide has been demonstrated as a competitive system, due to the large hyperpolarizability of NLO-active [MSe4] (M = Ga, In, Ge, Sn, Zn, Cd, Hg, etc.) tetrahedral units, and red shift of the vibration peaks of M–Se compared to M–O, M–S or M–P bonding.23,24 However, for most selenides, their band gaps are narrow due to the presence of outer electrons on the 3d, 4s, 4p orbits in the Se atom.25,26
To enhance the band gap of chalcogenides, introducing an alkaline-earth metal (AEM) without d–d and f–f electron transitions into the crystal structures has been regarded as a feasible strategy.27,28 Recently, by coupling [MgQ6] (Q = S, Se) octahedral and [MQ4] (M = Al, Ga) tetrahedral units, a series of Mg-containing compounds AIB3IIC3IIIQ8VI (AI = Li, Na, Ag; BII = Mg; CIII = Al, Ga; QVI = S, Se) have been developed.29 In this work, by introducing double AEM atoms, three wide band gap selenide IR NLO materials AIIMg6Ga6Se16 (AII = Ca, Sr, Ba) have been synthesized by the high temperature solid-state reactions in sealed quartz tubes. AIIMg6Ga6Se16 (AII = Ca, Sr, Ba) crystallize in the NCS P space group, and are built by [AIISe6] triangular prisms, [MgSe6] octahedral and [GaSe4] tetrahedral units. The compounds show balanced optical properties including strong second-harmonic generation (SHG) responses (1.5 × AGS for CaMg6Ga6Se16, 1.0 × AGS for SrMg6Ga6Se16, 1.1 × AGS for BaMg6Ga6Se16), wide band gaps (in selenides) (2.69–2.71 eV), high LIDTs (8.0–9.0 × AGS), and moderate birefringences (0.044–0.052 @ 1064 nm). Density functional theory (DFT) calculations reveal that the SHG responses and band gaps of the compounds are mainly determined by the [GaSe4] tetrahedral units.
Since the three compounds show similar cell parameters and crystal structures, CaMg6Ga6Se16 is utilized as an example to illustrate their structures here. CaMg6Ga6Se16 crystallizes in the P space group with a = b = 17.5327(3) Å, c = 7.7603(2) Å. In the asymmetric unit of CaMg6Ga6Se16, there are three Ca atoms, six Mg atoms, three Ga atoms and eleven Se atoms. The Ca atom and Mg atom are coordinated with Se atoms to form [CaSe6] triangular prisms and [MgSe6] octahedral units, while Ga atoms are bonded to four Se atoms to construct [GaSe4] tetrahedra (Fig. 1a). [GaSe4] tetrahedra are connected with each other by vertex-sharing to form the [GaSe3] single-chains and [Ga2Se5] double-chains (Fig. 1a). [MgSe6] octahedra are linked with each other by edge-sharing to form [MgSe4] single-chains and [Mg2Se6] double-chains (Fig. 1b). The formed [GaSe3] and [MgSe4] single-chains are linked by vertex-sharing to build the windmill-like [Mg3Ga3Se24] channels (Fig. 1c), which are further grouped with [Ga2Se5] and [Mg2Se6] double-chains to form the channel-like [Mg9Ga9Se24] framework (Fig. 1d). The six-coordinated Ca ions are filled in the channels (Fig. 1e) to balance the charge, resulting in the final 3D structure of CaMg6Ga6Se16 (Fig. 1f). It is worth noting that, different from CaMg6Ga6Se16 and SrMg6Ga6Se16, the asymmetric unit of BaMg6Ga6Se16 contains two crystallographically independent Ba atoms at Wyckoff 1c, 1e positions, and two partially occupied Ba atoms with the site occupancy of 0.9 and 0.1 at Wyckoff 1a, 1b positions, respectively.
Compared to the crystal structures of AIBII3CIII3QVI8 family compounds, the title compounds can be regarded as the atomic aliovalent substitution of monovalent AI (AI = Li, Na, Ag) by divalent AEM AII (AII = Ca, Sr, Ba). The aliovalent substitution does not cause the change of space group, but induces a doubled cell parameter c in AIIMg6Ga6Se16 compared to AIMg3Ga3Se8, which can be attributed to the formed ordered vacancy in AIIMg6Ga6Se16 (AII = Ca, Sr), or the partially occupied Ba atoms in BaMg6Ga6Se16. In addition, accompanying the fluctuation of A-site cations (Fig. 2a), the channel size fluctuates from 2.49–2.56 Å to 2.39–2.65 Å (Fig. 2b). Most of them are smaller than the AI/AII–Se bond lengths in the compounds, except in the case of AgMg3Ga3Se8 (Fig. 2b). It is worth noting that Ag is three-coordinated with three Se atoms to form a planar [AgSe3] unit in AgMg3Ga3Se8, different from the six coordinated [Li/NaSe6] triangular prism units in AIMg3Ga3Se8 (AI = Li, Na) and [Ca/Sr/BaSe6] triangular prism units in AIIMg6Ga6Se16 (AII = Ca, Sr, Ba).
To evaluate the optical properties, the pure-phase powder samples of the title compounds were prepared and characterized. The experimental XRD patterns are matched well with the theoretical results derived from cif files (Fig. 3a and S3†). To obtain experimental band gaps, UV-vis-NIR diffuse reflectance spectra were measured. Based on the Kubelka–Munk function,30–32 the experimental band gaps of CaMg6Ga6Se16, SrMg6Ga6Se16 and BaMg6Ga6Se16 were measured to be 2.71, 2.71 and 2.69 eV, respectively (Fig. 3b). The large band gaps inherently contribute to high resistances to laser damage.33,34 The LIDTs of the three compounds were evaluated by the single-pulse LIDT method with AGS as the reference.35–37 The measured results show that the LIDTs of CaMg6Ga6Se16, SrMg6Ga6Se16 and BaMg6Ga6Se16 powder samples were ∼9.0, ∼9.0, and ∼8.0 × AGS, respectively.
Beyond the band gap and LIDT, SHG response is another important parameter for IR NLO materials. The SHG responses of the title compounds were evaluated by the Kurtz and Perry method with a 2.09 μm Q-switched laser.38–41 As shown in Fig. 3c, the SHG intensity of CaMg6Ga6Se16, SrMg6Ga6Se16 and BaMg6Ga6Se16 was tested to be ∼1.5, ∼1.0, and ∼1.1 × AGS, respectively, showing that the A-site AEM has an influence on the SHG responses. It can be attributed to the changes of distortion degree and hyperpolarizability of nonlinear-active [GaSe4] units induced by the different AII–Se chemical bonding (the bond length dCa–Se = 3.0609–3.0640 Å, dSr–Se = 3.1520–3.1590 Å and dBa–Se = 3.2450–3.2612 Å) in the structures. Specifically, the hyperpolarizabilities of [Ga(1–3)Se4] in CaMg6Ga6Se16 were calculated to 283.24, 640.08, and 665.56, while 253.14, 615.32, and 625.81 for SrMg6Ga6Se16, and 190.51, 575.27, and 591.85 for BaMg6Ga6Se16, thus resulting in a larger SHG response in CaMg6Ga6Se16. Moreover, the hyperpolarizabilities of [Ga(1–3)Se4] in NaMg3Ga3Se8 were calculated to be 278.91, 627.83, and 659.72, smaller than the ones in CaMg6Ga6Se16. The results indicate that the alkali metal (AM) and AEM can exert an influence on the SHG response in an indirect way by modifying the geometry of NLO tetrahedral units in the AIBII3CIII3QVI8 and AIIBII6CIII6QVI16 families. Notably, the experimental results imply that CaMg6Ga6Se16 achieves a good balance between wide band gap and large SHG response, which is comparable with the recently reported typical selenide IR NLO candidates like β-BaGa2Se4 (0.6 × AGS, 2.49 eV),25 DCL-MgGa2Se4 (0.9 × AGS, 2.96 eV),26 AgIn5Se8 (0.7 × AGS, 1.34 eV),42 Pb4SeBr6 (1.3 × AGS, 2.62 eV),43 KSmP2Se6 (1.08 × AGS, 1.92 eV)44 and Na2Ga3In3Se10 (1.5 × AGS, 2.47 eV)45 (Fig. 3d).
To understand the structure–performance relationship, the band structures, total and partial density of states (TDOS/PDOS) were computed by DFT calculations.46–48 As shown in Fig. 4a, S4a and c,† the title compounds could be direct band gap compounds with calculated GGA band gaps of 1.851 eV (CaMg6Ga6Se16), 1.866 eV (SrMg6Ga6Se16) and 1.853 eV (BaMg6Ga6Se16). Considering the underestimation of band gaps in the standard DFT calculations with GGA because of the discontinuity of the exchange–correlation energy functional,49,50 the HSE06 band gaps of the compounds were calculated to be 2.568 eV (CaMg6Ga6Se16), 2.564 eV for (SrMg6Ga6Se16) and 2.556 eV (BaMg6Ga6Se16) (Table S14†), which are close to the experimental results in Fig. 3b. From the results of TDOS and PDOS (Fig. 4b, S4b and d†), it can be seen that the tops of the valence band maximum (VBM) near the Fermi level are determined by Se-4p orbitals, and the bottoms of the conduction band minimum (CBM) are mainly occupied by Se-4s4p and Ga-4s orbitals, indicating that the band gaps of the three compounds are mainly determined by Ga–Se bonding in [GaSe4] tetrahedral units. Moreover, the calculated birefringence values are 0.052 for CaMg6Ga6Se16, 0.048 for SrMg6Ga6Se16, and 0.044 for BaMg6Ga6Se16 at 1064 nm (Fig. 4c), matched with the PM behaviours in Fig. 3c. Meanwhile, the bonding electron density difference (Δρb) and contributions (Δω) indicate that [MgSe6] and [GaSe4] units make the main contribution to the optical anisotropy in the compounds (Fig. 4d and S5†).51
To detect the origin of the NLO effect, the NLO coefficients and SHG density maps of the title compounds were investigated by DFT calculations.52–56 As shown in Table S14,† the calculated NLO coefficients are d11 = −11.05 and d22 = 3.92 pm V−1 for CaMg6Ga6Se16, d11 = −10.67 and d22 = 3.93 pm V−1 for SrMg6Ga6Se16, d11 = −10.58 and d22 = 3.71 pm V−1 for BaMg6Ga6Se16, which match the experimental results in Fig. 3c. Meanwhile, the SHG density maps of occupied and unoccupied orbitals in the Virtual-Electron (VE) process which occupied 93.74% in CaMg6Ga6Se16, 95.25% in SrMg6Ga6Se16, and 95.82% in BaMg6Ga6Se16 (the Virtual-Hole (VH) process is neglected here because it accounts for a small quantity of contributions) indicate that the NLO responses of AIIMg6Ga6Se16 (A = Ca, Sr, Ba) also mainly originate from the [GaSe4] tetrahedral units (Fig. 4e, f and S6†), matched with the anionic group theory promoted by Chen et al.57,58
Footnotes |
† Electronic supplementary information (ESI) available. CCDC 2311758–2311760. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4sc00334a |
‡ L. A. Wang and D. D. Chu contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2024 |