Chunlan Tangab,
Aoge Yaoc,
Wenhao Xingd,
Wenwen Jiangb,
Jian Tangb,
Lei Kang
c,
Jieyun Wu
*a,
Wenlong Yin
*b and
Bin Kangb
aSchool of Optoelectronic Science and Engineering, University of Electronic Science and Technology of China, Chengdu 611731, P.R. China. E-mail: jieyunwu@uestc.edu.cn
bInstitute of Chemical Materials, China Academy of Engineering Physics, Mianyang 621900, P.R. China. E-mail: wlyin@caep.cn
cBeijing Center for Crystal Research and Development, Key Laboratory of Functional Crystals and Laser Technology, Technical Institute of Physics and Chemistry, Chinese Academy of Sciences, Beijing 100190, P.R. China
dSchool of Mechanical and Electrical Engineering, Henan Institute of Science and Technology, Xinxiang 453003, P. R. China
First published on 9th June 2025
Temperature-induced phase transitions offer a promising route to engineer nonlinear optical materials, particularly for infrared applications where conventional design approaches face fundamental limitations. Herein, a temperature-induced centrosymmetric (CS) to noncentrosymmetric (NCS) irreversible phase transition strategy was employed to successfully prepare a novel NCS diamond-like (DL) chalcogenide, β-Ag4P2S7, which was derived from the CS phase α-Ag4P2S7 transformation. Structural analysis reveals that this transformation involves a reorganization of [Ag2PS4]2− layers from an AA′AA′ to AA′A′′AA′A′′ stacking pattern, facilitated by bridging [Ag1S4]7− tetrahedra. β-Ag4P2S7 has exceptional IR NLO properties, including a strong phase-matchable second-harmonic generation (SHG) response (1.02 × AgGaS2) and a wide band gap of 2.90 eV (the largest one in ternary Ag-based DL chalcogenides), which balances excellent NLO response with a wide band gap. Further structure–property relationship analyses show that superior NLO properties and band gap broadening of β-Ag4P2S7 mainly originate from the alteration of the [Ag2PS4]2− layer stacking configuration, which is driven by a temperature-induced irreversible phase transition. This work not only presents a new paradigm for designing high-performance NLO materials through phase transitions but also significantly advances the potential of temperature-mediated crystal engineering for optical applications.
The design of IR NLO crystals with outstanding NLO performance hinges on achieving a noncentrosymmetric (NCS) crystal structure, an essential yet highly challenging prerequisite.11–15 Molecular engineering strategies, which focus on the selection of functional motifs and spatial optimization of NLO-active chromophores, have been extensively employed to develop novel NCS crystals. Traditional approaches for constructing NCS structures often rely on the partial substitution of elements or functional motifs within the CS structure16–20 (e.g., CS: RbGaS2 to NCS: [ABa2Cl][Ga4S8] (A = Rb, Cs);21 CS: Rb4P2S6 to NCS: RbBiP2S6;22 CS: ZnS to NCS: M[M4Cl][Ga11S20](M = A/Ba, A = K, Rb)).23 However, these methods frequently introduce structural and functional complexity that hinders the simultaneous optimization of the band gap and NLO properties. Notably, strong second-harmonic generation (SHG) effect is often accompanied by a narrowing of the band gap, presenting a fundamental challenge in material design.
Overcoming this limitation requires innovative strategies to break centrosymmetry while maintaining a wide band gap. One of the promising approaches is phase transition engineering, which enables the formation of NCS structures through gradual structural modifications driven by high energy barriers. This method stabilizes metastable phases and generates diverse coordination environments for structural ligands under external stimuli such as temperature, pressure, or atmospheric changes.24,25 By leveraging phase transitions rather than chemical substitution, this strategy provides a more straightforward and controllable route to designing NCS materials.
Phase transitions can be classified as reversible and irreversible. Reversible phase transitions are predominantly explored in organic material systems, where they enable cyclic regulation of NLO switching function.26,27 In contrast, inorganic compounds tend to undergo irreversible phase transitions due to their high thermodynamic stability (ΔG < 0) and significant phase transition barrier (>2 eV).28,29 Such irreversible transitions offer key advantages: they not only facilitate the growth of large-size single crystals via the Bridgman–Stockbarger method, but also enhance the structural and optoelectronic stability of devices under extreme operating conditions, such as temperature fluctuations, mechanical stress, and radiation environment.
In chalcogenides, even minor rearrangements of the tetrahedral ligands in diamond-like (DL) structures can induce phase transitions, leading to significant modification in crystal structures. For instance, Cu2ZnSiS4 undergoes a transition between the α-phase (Pmn21) and the β-phase (Pn),30 and Ga2Se3 exhibits multiple transitions (α: F3m; β: Cc; γ: F
3m; σ: C2221).31 Additionally, Ag+ cations, possessing a d10 closed-shell configuration, can couple filled d orbitals and empty s orbitals of similar energy due to the second-order Jahn–Teller (SOJT) effect, lowering the energy barriers, which is easily overwhelmed by thermal perturbations.32–34
Building upon those insights, we propose a high-temperature-induced irreversible phase transition strategy to construct an NCS chalcogenide based on the ternary Ag/P/S system. This approach successfully led to the transformation of CS α-Ag4P2S7 into a novel NCS phase, β-Ag4P2S7, both of which adopt DL structures. The synthesis, crystal structures, phase transitions, NLO properties, and structure–property relationship were systematically investigated. The resulting β-Ag4P2S7 exhibits a strong and phase-matching SHG response (1.02 × AgGaS2), a wide band gap (2.9 eV), large birefringence (Δn = 0.045), and excellent thermal/chemical stability. Moreover, the ternary Ag/P/S system offers significant advantages over quaternary or quinary systems, particularly in terms of large-scale crystal growth. Its relatively simple composition facilitates better control over crystallization processes, making it a highly attractive platform for developing next-generation IR NLO materials.
The thermal behavior of α-Ag4P2S7 was studied by differential scanning calorimetry (DSC) measurements (Fig. 2a). It showed three absorption peaks at 421 °C, 463 °C, and 558 °C, and only one exothermic peak at 490 °C during the entire heating and cooling process. In addition, to gain a comprehensive insight into the relationship between the phase transitions and temperatures, a series of solid-state reactions at different calcination temperatures was performed in sealed silica tubes (Fig. 2b) based on the DSC results. For α-Ag4P2S7, it is stable until 420 °C. Phase transition begins at 420 °C. The β-phase content increases progressively between 420–570 °C. Notably, Ag7P3S11 emerges as a byproduct through peritectic reaction with the melt (Fig. S2†).39 When the temperature exceeds 570 °C, a pure phase of β-Ag4P2S7 is obtained. When the temperature begins to decrease, the β-Ag4P2S7 fails to transform into the α-Ag4P2S7, indicating the irreversible phase transition. Once the irreversible phase transition occurs, the system will move towards more stable thermal properties, leading to a decrease in the free energy, implying that the new phase β-Ag4P2S7 is resistant to thermal perturbations and can maintain the stability of the phases and properties over a certain temperature range. In addition, irreversible phase transitions may lead to significant changes in the energy band structure of the material, such as an increase or decrease in the band gap. It is therefore necessary to probe the band gap of α-Ag4P2S7, which has an absorption edge at 430 nm (Fig. S3†) and a band gap of about 2.65 eV (Fig. 2c), which fits well with the yellow color of the crystals.
![]() | ||
Fig. 2 (a) DSC curve of the α-Ag4P2S7; (b) the variable-temperature powder X-ray diffraction patterns of the Ag4P2S7; (c) Eg of α-Ag4P2S7. |
Single-crystal XRD analysis indicates that β-Ag4P2S7 crystallizes in the trigonal chiral space group P3221 (no. 154) (cell parameters: a = b = 6.3563(5), c = 23.3477(18) Å and Z = 3). The energy-dispersive spectroscopy (EDS) measurement corroborated the homogenous distribution of Ag/P/S with an approximate molar ratio of 4:
2
:
7 (Fig. S4†). Its asymmetric unit also consists of two crystallographically independent Ag atoms, one P atom, and four S atoms, all at the Wyckoff position 6c except S2 at 3a (Table S4†). It is found that the structure of β-Ag4P2S7 is composed of the typical [Ag1S4]7− tetrahedra with d(Ag–S) = 2.51(3)–2.882(3) Å, [Ag2S4]7− tetrahedra with d(Ag–S) = 2.526(3)–2.812(3) Å, and the [PS4]3− tetrahedra with d(P–S) = 1.996(3)–2.123 (3) Å (Table S5 and Fig. S1†), which are in accordance with those of related compounds.36–38 The crystal structure of β-Ag4P2S7 is shown in Fig. 3b and can be regarded as 3D anionic networks consisting of 2D [Ag2PS4]2− layers and [Ag1S4]7− tetrahedra, which are interconnected by sharing S atoms. Similar to α-Ag4P2S7, the [Ag2PS4]2− layers of β-Ag4P2S7 are bridged by 1D [Ag2S3]5− chains of [Ag2S4]7− tetrahedra linked by corner-sharing and isolated [P2S7]4−dimers constructed from two [PS4]3− tetrahedra via corner-sharing (Fig. 3d). Different from α-Ag4P2S7, however, [Ag2PS4]2− layers adopt a more complex alternating stacking manner, i.e., AA′A′′AA′A′′…, leading to a superposition of microscopic SHG coefficients and macroscopically enhancing the NLO effect. So, a larger unit cell is necessary for β-Ag4P2S7. In addition, the bond valence calculations resulted in values of 0.96–1.05 for the Ag+ cation, 5.26 for the P5+ cation, and 1.97–2.13 for the S2− anion (Table S4†). All of these are consistent with their respective ideal oxidation states.
Further structural comparisons of α- and β-Ag4P2S7 were performed to gain insight into the phase transition behavior (Fig. 3). The density of high-temperature phase β-Ag4P2S7 (4.377 g cm−3) is lower than that of low-temperature phase α-Ag4P2S7 (4.381 g cm−3), which is consistent with the principle of thermal expansion and contraction. Moreover, the unit cells gradually decrease from α- to β-Ag4P2S7, that is, a(α) = 10.7634(8) Å ≈ 1.69 × b(β), b(α) = 6.5280(4) Å ≈ b(β), and c(α) = 23.3477(18) Å ≈ 0.69 × c(β). As shown in Fig. 3a and b, along the a-axis, the a-Ag4P2S7 unit cell can accommodate two [P2S7]4− dimers, while that of β-Ag4P2S7 can only accommodate one [P2S7]4− dimer. Along the c-axis, the unit cell of a-Ag4P2S7 can accommodate two [P2S7]4− dimers, whereas that of β-Ag4P2S7 can accommodate three [P2S7]4− dimers. Those indicate that the phase transition from α- to β-Ag4P2S7 has an abnormal anisotropic thermal expansion, such as contraction along the a direction and expansion along the c direction. In addition, the Ag1–S distances of in the [Ag1S4]7− bridges of β-Ag4P2S7 (2.551–2.882 Å) are longer than the Ag1–S distances of a-Ag4P2S7 (2.538–2.876 Å) (Fig. S5†), facilitating the rotation of the [Ag2PS4]2− layer and constructing to NCS structure accompanied by a significant SHG response. Intriguingly, the [Ag2PS4]2− layer shows a higher symmetry of the 32 screw axis in β-Ag4P2S7 compared to the 2-fold axis in α-Ag4P2S7 (Fig. 3c and d). As can be seen in Fig. 3e and f, the spatial symmetry operation increases from four (in α-Ag4P2S7) to six (in β-Ag4P2S7).
The microstructural characteristics of a material, such as atomic arrangement, crystal structure, phase composition, etc., fundamentally determine the material's various physical and chemical properties. We evaluated the NLO properties of the β-Ag4P2S7 by using the Kurtz–Perry method under the Q-switched 2100 nm laser.40 As illustrated in Fig. 4a, the SHG intensity tended to attain saturation with an increase in the particle size of the sample, which indicated that β-Ag4P2S7 was phase-matchable. β-Ag4P2S7 demonstrates a strong SHG intensity of about 1.02 times that of AgGaS2 at 150–200 μm (Fig. S6†), and is comparable to previously reported IR NLO materials [RbBa2Cl][Ga4S8] (1.0 × AgGaS2),19 NaAgPS4 (1.05 × AgGaS2),41 NaMg3Ga3Se8 (1.0 × AgGaS2),42 and [K4Cl][CdGa9S16] (0.9 × AgGaS2).43 The large SHG response of β-Ag4P2S7 will be favorable for generating high application conversion efficiencies. Note that the key scientific issue in the design of promising IR NLO crystals, especially for high-power laser applications, is how to strike a balance between strong SHG effects and high LDT, which have opposite structural requirements. In general, the barrier of photon absorption increases with a larger optical band gap, and materials tend to have a higher LDT. However, Ag-based chalcogenides typically possess a shallow LDT further dramatically shortening the device longevity which stems from the 4d, 5s and 5p orbitals of Ag+ always occupying the bottom of the conduction band (CB) and the top of the valence band (VB) leading to the narrow band gap and the photodarkening effect of silver. The absorption edge of β-Ag4P2S7 is at 414 nm (Fig. S7†) and the band gap of β-Ag4P2S7 is deduced to be 2.90 eV based on the Kubelka–Munk equation, which is in good agreement with the buff color of the crystals (Fig. 4b) and, to our knowledge, is the largest band gap in the ternary Ag-based DL chalcogenides (Fig. 4e and Table S6†).44–51 We hypothesize that the strong covalent P–S anionic frameworks promote the expansion of the band gap and drive the blue shift of the shortwave absorption edge, thus achieving high LDT, which is verified in subsequent theoretical calculations.
The thermal behavior of β-Ag4P2S7 was investigated by DSC, as illustrated in Fig. 4c. There are two endothermic peaks at 460 °C and 558 °C, corresponding to the peritectic reaction of Ag7P3S11 with the melt and melting point, respectively, and one exothermic peak at 489 °C during one heating/cooling cycle. In addition, the match of the powder X-ray diffraction (PXRD) patterns between the DSC measurement of β-Ag4P2S7 and the β-Ag4P2S7 sample further illustrates the melts congruence and the irreversible behavior of the phase transition (Fig. 4d). Thereby, β-Ag4P2S7 is thermodynamically stable and can be grown as bulk single crystals by the Bridgeman–Stockbarger method. With its stable DL crystal structure, β-Ag4P2S7 possesses desirable water resistance. After being immersed in water for over two weeks, PXRD analysis showed that the crystalline structure of β-Ag4P2S7 remained stable for further practical applications (Fig. 4d). The Raman spectrum recorded on the β-Ag4P2S7 powder sample is displayed in Fig. S8.† The peaks located at 450–600 cm−1 correspond to vibration modes of the P–S bond, while peaks existing between 200 and 450 cm−1 can be designated to vibrations of the Ag–S bond.52,53
Electronic band structure calculations of α-Ag4P2S7 and β-Ag4P2S7 have been carried out to comprehensively understand the relationship between the structure and optical properties based on first principles. The α-Ag4P2S7 undergoes an indirect transition from the M point to the other position points, corresponding to an indirect band gap value of 2.16 eV (Fig. 5a). The β-Ag4P2S7 transitions directly from the M point to the M point, corresponding to a direct band gap of 2.42 eV (Fig. 5b). The calculated band gap value is significantly underestimated compared to the experimental result (2.65 eV and 2.90 eV, respectively) due to the discontinuity in the exchange–correlation potential.54 It is generally accepted that optical properties are mainly determined by electronic transitions close to the band gap. To profoundly understand their band structure characteristics, partial densities of states (PDOS) are presented in Fig. 5a and b. Both structures have similar PDOS; the top of the valence band (VB) is predominantly dominated by Ag 4d and S 3p orbitals, while the bottom of the conduction band (CB) is primarily contributed by P 3p and S 3p orbitals. In addition, orbital analyses are carried out to understand the structure-performance relationship more intuitively. The partial charge density distributions of the highest occupied atomic orbitals (HOMO) and the lowest unoccupied atomic orbitals (LUMO) of α-Ag4P2S7 and β-Ag4P2S7 are plotted in Fig. 5c and d, respectively. It is observed that the HOMO states of α-Ag4P2S7 are mainly composed of Ag 4d and S 3p orbitals, and the LUMO states of α-Ag4P2S7 consist of unoccupied P 3p, and S 3p orbitals. The HOMO state of β-Ag4P2S7 consists mainly of Ag2 4d and S 3p orbitals and the LUMO state consists mainly of unoccupied P 3p and S 3p orbitals. The stacking mode of the [Ag2PS4]2− layer change selectively strengthens the covalency of the anionic framework (e.g., P–S) and weakens the Ag1 4d–S 3p orbital hybridization. Therefore, the band gap increase induced by the CS to NCS structural phase transition is mainly altered by the stacking mode of the [Ag2PS4]2− layer. In addition, the charge transitions responsible for the SHG efficiency in β-Ag4P2S7 are mainly dominated by the charge transfer from Ag2 4d and S 3p to P 3p and S 3p orbitals in the [Ag2PS4]2− layer. The calculated SHG coefficients of d16 = d22 = 6.7 pm V−1 are much smaller than second-order nonlinear optical coefficient (d36 = 13.7 pm V−1) of the reference material AgGaS255 and the experimental value (deff = 13.8 pm V−1) of β-Ag4P2S7, and a similar underestimate of the second-order susceptibilities has also been observed in other reports.56,57 Furthermore, the birefringence of α-Ag4P2S7 is slightly larger than that of β-Ag4P2S7 at 2100 nm, which is consistent with the general rule that low crystal symmetry usually leads to higher birefringence (Fig. 5e).58 The birefringence of β-Ag4P2S7 is suitable for a mid-infrared material to achieve phase-matching conditions. To understand the contribution of NLO active electron states and units, the SHG-weighted electron density analysis of β-Ag4P2S7 was adopted with electron-including occupied and unoccupied states (Fig. 5f).59 The result shows that SHG efficiency in β-Ag4P2S7 is mainly determined by the synergistic effect of the anionic groups [Ag2S4]7−, [PS4]3−, and the small contribution of [Ag1S4]7− tetrahedra. Furthermore, the dipole moments calculations for [Ag1S4]7−, [Ag2S4]7−, and [PS4]3− tetrahedra in the unit cell were calculated based on the bond valence method60 to analyze the origin of the strong SHG effect, as displayed in Table 1. For per cell, the z components of the [Ag1S4]7−, [Ag2S4]7−, and [PS4]3− tetrahedral polarizations almost cancel out. The [PS4]3− tetrahedra possess a very large net dipole moment of 45.43 Debye (D), followed by [Ag2S4]7− tetrahedra with a net dipole moment of 18.78 D, whereas [Ag1S4]7− tetrahedra has the smallest net dipole moment of 4.35 D. Overall, the synergistic action of the [Ag2S4]7− and [PS4]3− tetrahedra in [Ag2PS4]2− layer contribute mainly to the net dipole moment, which facilitates the superposition of their polarizations. In other words, the temperature-induced changes in the alternating stacking pattern of the [Ag2PS4]2−layers contribute to the superposition of microscopic SHG coefficients and macroscopically enhance the NLO effect.
Units | Dipole moment | ||||
---|---|---|---|---|---|
x(a) | y(b) | z(c) | Magnitude | ||
Debye | esu cm Å−3 | ||||
[Ag1S4]7− | −2.1760 | −3.7696 | 0 | 4.35 | 0.016 |
[Ag2S4]7− | 9.3872 | 16.2600 | 0 | 18.78 | 0.069 |
[PS4]3− | 22.7225 | 39.3409 | 0 | 45.43 | 0.167 |
Per cell | 29.9337 | 51.8312 | 0 | 59.85 | 0.220 |
Footnote |
† Electronic supplementary information (ESI) available: Crystallographic data. CCDC 2425474 for (α-Ag4P2S7) and 2425475 for (β-Ag4P2S7). For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d5qi00969c |
This journal is © the Partner Organisations 2025 |