Open Access Article
Khaoula Ben Brahim
a,
Arafat Ghoudi
a,
Ines Mbareka,
Noweir Ahmad Alghamdib,
Malika Ben Gzaiela and
Abderrazek Oueslati
*a
aLaboratory for Spectroscopic Characterization and Optics of Materials, Faculty of Sciences, University of Sfax, B. P. 1171, 3000 Sfax, Tunisia. E-mail: oueslatiabderrazek@yahoo.fr
bDepartement of Physics, Faculty of Science, Al-Baha University, Alaqiq, 65779, Saudi Arabia
First published on 11th November 2025
The present work focuses on elucidating the origin of the phase transition in [(CH3)4P]3Sb2Br9 and its associated electrical properties. The compound was synthesized by slow solvent evaporation at constant temperature and characterized using X-ray powder diffraction, FTIR spectroscopy, temperature-dependent Raman spectroscopy, and impedance spectroscopy. XRD confirmed the formation of a monoclinic structure with the polar P31c space group. Complementary FTIR and Raman analyses enabled the assignment of vibrational modes at room temperature, providing a reliable reference for monitoring temperature-induced phase transitions. Raman spectra collected between 294 and 563 K revealed specific molecular displacements and reorientations of both organic and inorganic components associated with the transition near 540 ± 7 K. Impedance spectroscopy in the 433–563 K range displayed a single semicircle in the Nyquist plots, evidencing a grain-dominated electrical response. The impedance data were successfully modeled by an equivalent circuit comprising a resistance in parallel with a constant phase element. Both ac and dc conductivities were analyzed to clarify the conduction mechanisms, with charge transport found to obey Jonscher's power law and consistent with the Correlated Barrier Hopping (CBH) and Non-overlapping Small Polaron Tunneling (NSPT) models. [(CH3)4P]3Sb2Br9 shows both phase transition and ionic conduction, highlighting its potential in memory devices, electrolytes, light-harvesting, and energy applications.
The synergy between organic and inorganic domains enables the engineering of multifunctional materials whose combined properties exceed those of their individual constituents.5 For instance, the structural versatility of hybrid materials allows for the enhancement of energy conversion efficiencies in third-generation solar cells and light-emitting diodes (LEDs), contributing to improved performance and longevity.6–8 In the realm of energy storage, hybrid materials have shown promise in reinforcing the capacity and durability of battery and supercapacitor electrodes.9–11 Furthermore, they are widely applied in smart coatings, offering attributes such as UV resistance, anti-corrosion, and anti-fouling properties.12,13 In catalysis, these materials act as highly efficient catalysts, particularly in photocatalytic processes for pollutant degradation and hydrogen generation.14,15
Various synthesis methods, including sol–gel processing, precipitation techniques, and slow evaporation, are employed to fabricate hybrid systems, each offering distinct advantages in controlling morphology, particle size, and component dispersion.16–18 Phase transitions in these materials are often governed by intricate mechanisms, such as atomic rearrangements, modifications in intermolecular forces, or changes in lattice structure.19 Due to the flexible nature of organic subunits and the nanoscale dimension of inorganic clusters, these transitions can exhibit behaviors not observed in homogeneous systems. Understanding and manipulating such transitions is crucial for tailoring material properties to suit specific functional applications.20 Despite notable progress, several challenges remain in the practical development of hybrid materials, including ensuring long-term stability, reducing production costs, and managing interfacial compatibility. Future studies are expected to focus on resolving these issues through advanced synthesis strategies and deeper insights into nanoscale interactions.
Hybrid materials incorporating antimony present particularly interesting opportunities across various high-tech domains. Thanks to their exceptional optical and electronic properties, Sb-based hybrids show great promise in applications such as waveguides, modulators, and energy storage systems. In batteries and supercapacitors, they contribute to enhanced electrode performance.21,22 They are also good candidates for catalytic and photocatalytic applications due to their potent absorption and effective charge separation capabilities.23 Antimony-containing perovskites, as (CH3NH3)3Sb2I9, are being investigated in solar technologies to increase thermal stability and conversion efficiency.24 Sb-based hybrids are positioned as essential materials for next-generation technologies because of their multifunctional properties. In this context, the present work investigates the antimony-based hybrid compound [(CH3)4P]3Sb2Br9, with a particular focus on elucidating the origin of its phase transitions through synthesis, temperature-dependent vibrational analysis, and electrical characterization by impedance spectroscopy.
To analyze the intermolecular interactions in the crystal, Hirshfeld surface analysis serves as an effective approach.27 This technique, applied using the Crystal Explorer software,28 separates the crystal's electron density into molecular fragments, allowing detailed investigation of how molecules are packed within the structure. The Hirshfeld surface itself is defined using two key metrics: de, denoting the distance from a given surface point to the closest atom outside the molecule, and di, representing the distance to the nearest atom inside the molecule.
A normalized contact measure, dnorm, is calculated from these distances combined with atomic van der Waals radii, ensuring accurate mapping of intermolecular contacts. For visualization, dnorm values are color-coded: red areas show close contacts (short intermolecular distances, negative dnorm), white areas represent contacts at van der Waals distances (zero dnorm), and blue areas indicate longer contacts (positive dnorm). The dnorm mapping, generated in Crystal Explorer, enables the creation of two-dimensional fingerprint plots that summarize and classify the types of intermolecular interactions present in the crystal. In this process, the van der Waals radii for atoms (rvdW) are accounted for, and the proximity of atoms internal (di) and external (de) to the surface is used to characterize specific crystal contact.
The infrared spectrum of [(CH3)4P]3Sb2Br9 was acquired at room temperature in the frequency range of 500–4000 cm−1 using a PerkinElmer FT-IR 1000 spectrophotometer.
The Raman spectra were acquired between 70 and 3100 cm−1 throughout the temperature range of 294–563 K using a T64000 (Jobin-Yvon, Horiba group, Kyoto, Japan) spectrometer. An argon–krypton ion laser (Innova, Coherent, France) provided the incident wavelength light, which was then focused at 647.1 nm using an MSPlan 50× (N.A. 0.55), on an Olympus BX41 microscope. The laser power used on the sample during the experiment was less than 2 mW, and two scans of the Raman spectra were produced with an integration time of 60 seconds per scan. Then, a nitrogen-cooled, front-illuminated charge-coupled device (CCD) detector connected to a spectrometer with a single monochromator (600 gratings per millimeter) was used to analyze the scattered light. The spectral profiles were deconvoluted using the LabSpec software with a combination of Gaussian and Lorentzian functions to enable a quantitative analysis of band evolution with temperature.
Impedance spectroscopy (IS) was employed to investigate the electrical properties of [(CH3)4P]3Sb2Br9 polycrystalline powder. This technique is widely recognized for its ability to provide detailed insights into the electrical behavior of materials and their interfaces across a broad range of frequencies and temperatures.29 In this study, the polycrystalline powder was compressed into disc-shaped pellets (8 mm diameter, 1.2 mm thickness) under uniaxial pressure of 3 t cm−2. To ensure optimal electrical contact and minimize interfacial resistance, silver conductive paste was applied to both flat surfaces of each pellet before impedance measurements. Electrodes were then affixed to the silvered faces, allowing for consistent electrical connections throughout the measurement range. The pellet and its electrodes were placed between the two contacts of the Solartron SI 1260 impedance analyzer inside a temperature-controlled chamber to allow for the acquisition of impedance spectra between 433 K and 563 K.
The compound crystallizes in the trigonal system and is made up of several [Sb2Br9]3− anions and tetramethylphosphonium cations (Fig. 2). Each anion is composed of two corner-sharing SbBr6 octahedra that form a bioctahedral [Sb2Br9]3− unit. With two various types of Sb–Br connections, the antimony atoms take on some form of distorted octahedral coordination: bridging bonds around 3.03 Å and terminal bonds between 2.64 and 2.80 Å (Table S1). With bond angles ranging from 91.9° to 93.2°, the Br–Sb–Br exhibits a minor departure from ideal octahedral symmetry. The geometry found in similar compounds like [(CH3)4N]3[Sb2Br9] and [(CH3)4N]3[Sb2Br3Cl6] is compatible with the Sb⋯Sb separation located within the anionic dimer, which is 3.933 Å.
![]() | ||
| Fig. 2 Crystal structure of [(CH3)4P]3Sb2Br9 viewed along the a-axis, generated using MERCURY software.31 | ||
Weak C–H⋯Br hydrogen bonds between the methyl groups of the cations and the terminal bromine atoms of the anions control the crystal packing, resulting in a three-dimensional supramolecular network. Via local polarization effects, these interactions maintain the crystal structure and might have a little impact on the electrical and dielectric responses.
Therefore, although the crystal structure has been previously described, this detailed recall is essential to understand how the local geometry and supramolecular interactions impact the macroscopic electrical behavior of the compound studied in this work.
In order to verify the formation and purity of the synthesized compound, X-ray powder diffraction (XRPD) measurements were carried out on the obtained powder sample.
Fig. 3 shows the experimental XRPD pattern and its refinement by the Lebai method at room temperature of the [(CH3)P]3Sb2Br9 sample. The product is extremely pure and crystalline, as evidenced by the absence of peaks from other contaminants. Every measured line has an exact index with the predicted space group (P31c). The calculated and experimental spectra exhibit good agreement, with χ2 = 17.1, Rp = 7.45, Rwp = 8.77, and Rexp = 2.12, indicating a high level of refinement. The obtained unit-cell parameters are reported a = b = 9.797(32) Å, c = 22.653(90) Å, α = β = 90° and ϒ = 120°. These parameters and the experimental powder diffraction pattern matches well with the simulated pattern derived from the single crystal data.32
To delineate the specific region occupied by each molecule within the crystal, the Hirshfeld surface approach was utilized, segmenting the electron density into molecular fragments. The three-dimensional surface was mapped across several shape metrics, including dnorm, shape index, di, and de. As revealed in Fig. 4, mapping over dnorm generates three main color-coded zones: white indicates contacts matching the sum of van der Waals radii, blue suggests longer contacts, and red corresponds to closer, higher-density regions. Red patches on the dnorm surface point to strong interactions, such as N–H⋯Cl hydrogen bonds, which mark the shortest intermolecular contacts in this structure. Additionally, visualization by shape index highlights donors and acceptors involved in bonding, showing up as red, orange, or blue features near specific atoms.
![]() | ||
| Fig. 4 Hirshfeld surfaces mapped with (a) dnorm, (b) shape index, (c) di, and (d) de surface for the title compound. | ||
Analysis of the two-dimensional fingerprint plots (see Fig. 5) provided a quantitative breakdown of surface contacts: Br⋯H interactions accounted for 62.2% of the Hirshfeld surface, making them the main contributor to molecular packing. H⋯H contacts followed at 33%, reflecting the high number of surface hydrogens present. Sb⋯Br (3.6%), Br⋯Br (0.6%), and Sb⋯H (0.4%) contacts made smaller but measurable contributions. These findings highlight that Br⋯H and H⋯H interactions are dominant in dictating the crystal packing arrangement.
| IR (cm−1) | Raman (cm−1) | Assignments | References | Compound |
|---|---|---|---|---|
| a W: weak; m: medium; s: strong. | ||||
| 91m | δ(Sb–Br) | 35 and 37 | Cs3Sb2Br9 (CH3NH3)3Sb2Br9 | |
| 133s | δ(Sb–Br) | |||
| 181m | δ(Br–Sb–Br) | |||
| 204s | ʋ(Sb–Br) | |||
| 240m | ʋ(Sb–Br) | |||
| 283s | ʋ(Sb–Br) | |||
| 648s | ʋs(PC4) | 33 and 36 | [(CH3)4P]3Bi2Cl9 [(CH3)4P]3Sb2Cl9 | |
| 775m | 778m | ʋas(PC4) | ||
| 973s | 943w, 981w | ρ(CH3) | ||
| 1300m | 1305w | δs(CH3) | ||
| 1419m | 1402m, 1408m, 1424w | δas(CH3) | ||
| 2787w | ʋ(CH) | 33 and 36 | [(CH3)4P]3Bi2Cl9 [(CH3)4P]3Sb2Cl9 | |
| 2840w | ʋ(CH) | |||
| 2907w | 2909s | ʋs(CH3) | ||
| 2985s | 2975m, 2986m | ʋas(CH3) | ||
In the low-frequency region below 400 cm−1, the spectra are primarily characterized by lattice vibrations and internal modes of the [Sb2Br9]3− anionic unit. Above 400 cm−1, the vibrational features are dominated by the [(CH3)4P]+ cations, including stretching and deformation modes of the PC4 and CH3 groups.36 The high-frequency region is associated with C–H stretching vibrations from the methyl groups. All specific mode assignments and corresponding frequencies are summarized in Table 1.
Differential scanning calorimetry (DSC) measurements reported in the literature30 reveal that [(CH3)4P]3Sb2Br9 undergoes two reversible phase transitions at 193 K and 540 K upon heating. To investigate the structural dynamics underlying the high-temperature transition, a series of fifty Raman spectra was recorded between 303 K and 563 K. This temperature-dependent study provides a direct probe of molecular reorientations in both the organic cations and inorganic anions. Careful analysis of the evolving vibrational bands allows the identification of specific lattice distortions and cationic motions that accompany the phase transition, offering a detailed understanding of the mechanisms driving the structural change. The temperature dependence of the main vibrational modes, including both peak positions and full widths at half maximum (FWHM), is presented in Fig. 8–12.
![]() | ||
| Fig. 8 Raman wavenumber temperature evolution of (a) the full width at half maximum (FWHM) of the Raman bands and (b) the position of the bands in the frequency range corresponding to δ(Sb–Br). | ||
![]() | ||
| Fig. 12 Raman wavenumber temperature evolution of (a) the full width at half maximum (FWHM) of the Raman bands and (b) the position of the bands in the frequency range corresponding to δ(CH3). | ||
Significant spectral changes are observed in the internal modes of the [Sb2Br9]3− anion. The Sb–Br bending vibration at 134 cm−1 exhibits a slight blue shift of less than 2 cm−1 accompanied by a noticeable broadening of approximately 4 cm−1 above 540 K. The symmetric Sb–Br stretching mode at 205 cm−1 shifts upward by 3 cm−1 and broadens to roughly 10 cm−1 near the transition, while the asymmetric Sb–Br stretching mode undergoes a similar broadening.
Vibrational modes associated with the tetramethylphosphonium cation are also affected by the transition. The symmetric and asymmetric PC4 stretching vibrations at 648 and 779 cm−1 shift by about 1 cm−1 toward higher frequencies, with the asymmetric mode showing a pronounced broadening of approximately 11 cm−1. The CH3 deformation modes at 1402 and 1409 cm−1 merge into a single, broader band with a widening of around 6 cm−1 as the temperature approaches 540 K. Similarly, the symmetric CH3 stretching mode between 2902 and 2909 cm−1 shifts by 4 cm−1, and the asymmetric CH3 stretching modes at 2974 and 2975 cm−1 coalesce into a single peak with an increase in FWHM of about 5 cm−1 above 540 K.
These observations highlight the active participation of both cationic and anionic substructures in the phase transition. The abrupt changes in peak positions and spectral widths, particularly beyond 540 K, indicate the onset of a disordered phase consistent with an order–disorder transition mechanism. Most vibrational modes display discontinuities at the transition, further supporting this interpretation.
The electric response of [(CH3)4P]3Sb2Br9 will be significantly impacted by the fact that temperature-driven structural dynamics, molecular vibrational, and lattice scatter occur at 540 K, according to Raman spectroscopy results. Impedance spectroscopy and dielectric measurements were carried out across a wide temperature range to investigate these effects. The rivalry between the structural dynamism processes and the charge transport in the compound is fully described, and correlations between vibrational changes and the evolution of electrical conduction and dielectric characteristics are made possible.
The time constants for grain conductivity and grain boundaries need to differ greatly for two distinct arcs to form. There are two situations in which a single arc may be seen in the Nyquist plot: (a) if the grain boundaries are highly conductive, meaning that they efficiently transmit charge without a significant impedance, they may not form a distinct arc; and (b) if the bulk conductivity is very high, the grain arc associated with grain conductivity is so small as to be imperceptible, and only the contribution from grain boundaries is observed.
At every temperature examined, the chemical under study exhibits a single depressed semicircular arc in its spectra. The observation of only one arc must be interpreted by hypothesis (b) above, because the impedance of the grain borders is crucial in this scenario, implying that the arcs seen in the Nyquist plots represent the behavior of the sample's grain boundaries. The hypothesized electrode-material interface contribution to the observed impedance is not present for the sample under study. Conversely, the fact that the measured arc radii decrease with temperature suggests that the material under study contains thermally active relaxation processes. Moreover, each arc's center is below the Z′-axis, indicating that the material does not behave in a Debye-type manner.38
The measured impedance data were fitted using the Z-view program to a typical equivalent circuit made up of two parallel elements: a fractal capacity element (also known as a constant phase element or CPE1) and a resistance element R1 (see insert in Fig. 13). In the equivalent circuit, the non-ideal capacitance behavior of the material is compensated for by the CPE element (ZCPE1 = 1/Q1(jω)α1, where j is the imaginary number, w the angular frequency, and α1 a constant whose value stays between 0 and 1, being zero for the pure resistive behavior and unity for the pure capacitive behavior39). The data in Fig. 13 makes it abundantly evident that the equivalent circuit model accurately reproduces the experimental data. Table 2 shows the corresponding circuit characteristics.
| T (K) | α1 | Q1 (×10−11 F) | R1 (×105 Ω) |
|---|---|---|---|
| 433 | 0.97722 | 4.71 | 578 |
| 443 | 0.97537 | 4.84 | 327 |
| 453 | 0.97377 | 4.98 | 193 |
| 463 | 0.97399 | 5.01 | 114 |
| 473 | 0.97299 | 5.10 | 70.3 |
| 483 | 0.97195 | 5.20 | 45.7 |
| 493 | 0.9714 | 5.27 | 30.2 |
| 503 | 0.97089 | 5.33 | 20.1 |
| 513 | 0.96992 | 5.52 | 13.2 |
| 523 | 0.96878 | 5.65 | 8.7287 |
| 533 | 0.96729 | 5.83 | 5.8442 |
| 543 | 0.96553 | 6.05 | 3.8893 |
| 553 | 0.96407 | 6.26 | 2.5199 |
| 563 | 0.96427 | 6.40 | 1.5026 |
Fig. 14(a and b) illustrates the behavior of the real (Z′) and imaginary (Z′′) components of the complex impedance as a function of frequency at various temperatures.
![]() | ||
| Fig. 14 (a) Temperature and frequency dependence of the real component of impedance variation. (b) Temperature and frequency dependence of the imaginary component of impedance variation. Fits shown as line plots represent the results of eqn (1) and (2) for Z′ and Z′′, respectively. | ||
The Z′ spectra (Fig. 14a) reveal two distinct regions: (i) a low-frequency plateau where Z′ remains nearly constant, attributed to the long-range motion of charge carriers traversing the material, and (ii) a frequency-dependent region at intermediate frequencies, associated with AC conductivity. In this regime, AC conduction arises from localized, short-range movements of charge carriers that occur after they overcome potential barriers within the material. At lower temperatures, Z′ is larger in the low-frequency range but decreases progressively with increasing temperature or frequency. This trend suggests a semiconductor-like behavior, indicating that [(CH3)4P]3Sb2Br9 possesses a negative temperature coefficient of resistance (NTCR).40 At higher frequencies, the Z′ curves converge, reflecting the presence of a space charge region. The NTCR behavior suggests that [(CH3)4P]3Sb2Br9 could serve as a thermistor material.41
The Z′′ spectra (Fig. 14b) exhibit an initial rise with frequency, reaching a peak at the material's characteristic relaxation frequency, ωmax (where resistance equals reactance), followed by a decline at higher frequencies. These relaxation peaks are broad, asymmetric, and indicative of non-Debye-type relaxation behavior.38,42 As the temperature increases, the peak height decreases while shifting towards higher frequencies, confirming thermally activated relaxation processes in [(CH3)4P]3Sb2Br9 and an increase in the electron hopping rate with temperature.
At higher frequencies, polarization decreases because the reduced relaxation time of hopping electrons prevents them from keeping pace with the alternating electric field. Moreover, space charge polarization is prominent in the low-frequency range but becomes negligible at higher frequencies, as evidenced by the complete overlap of Z′′ values in this region.
In this study, the experimental Z′ (real part) and Z′′ (imaginary part) spectra, presented in Fig. 14, were fitted using the theoretical models described by eqn (1) and (2), respectively. In Fig. 14a, the solid lines represent the fits to the Z′ data using eqn (1), while in Fig. 14b, the solid lines correspond to the fits of the Z′′ data using eqn (2).
![]() | (1) |
![]() | (2) |
These results confirm that the relaxation mechanism in [(CH3)4P]3Sb2Br9 is non-Debye in nature, indicating the presence of a single relaxation process over the studied frequency and temperature ranges.
Here, σdc & ZeroWidthSpace; represents the dc conductivity, which is independent of frequency, while Aωs denotes the dispersive component that varies with frequency. The parameter s is a dimensionless factor that depends on both frequency and temperature, while A is a temperature-dependent constant. According to Jonscher's theory, the relaxation dynamics of mobile charge carriers contribute significantly to the observed AC conductivity. When these carriers move between two potential energy minima, their motion can be categorized as translational (s < 1), indicative of extended charge transport, or localized (s > 1), where carriers remain confined within a restricted spatial domain.
| σac = σdc + Aωs | (3) |
At elevated temperatures, the ac conductivity spectra reveal two distinct regions. In the low-frequency domain, a plateau emerges where conductivity remains nearly independent of frequency. This region is associated with the long-range translational motion of charge carriers,46 as shown in Fig. 15a. In contrast, at higher frequencies, the conductivity increases following a power-law dependence on ωs. The transition frequency, marking the crossover between dc-dominated and dispersive regimes, shifts toward higher frequencies with increasing temperature. This behavior reflects the enhanced mobility of charge carriers at elevated thermal energies, facilitating more rapid hopping or tunneling processes over shorter length scales. Such temperature-dependent shifting of the crossover frequency indicates a thermally activated conduction mechanism, consistent with models based on correlated barrier hopping or quantum mechanical tunneling.
To ensure the reliability of the proposed conduction mechanisms, the experimental conductivity spectra were systematically fitted using Jonscher's universal power law. The fitting, presented in Fig. 15b, closely reproduces the experimental data, thereby confirming the adequacy of the model in describing the AC response of the studied material. The extracted parameters, including σdc, A, and s across different temperatures, are listed in Table 3.
| T (K) | A | σdc | s |
|---|---|---|---|
| 433 | 6.39 × 10−11 | 3.334 × 10−7 | 0.969 |
| 443 | 1.50 × 10−11 | 6.285 × 10−7 | 0.920 |
| 453 | 5.37 × 10−12 | 1.034 × 10−6 | 0.919 |
| 463 | 1.49 × 10−12 | 1.726 × 10−6 | 0.909 |
| 473 | 1.09 × 10−12 | 2.805 × 10−6 | 0.851 |
| 483 | 9.09 × 10−13 | 4.294 × 10−6 | 0.823 |
| 493 | 4.18 × 10−13 | 6.556 × 10−6 | 0.785 |
| 503 | 2.21 × 10−13 | 9.856 × 10−6 | 0.727 |
| 513 | 8.93 × 10−14 | 1.498 × 10−5 | 0.698 |
| 523 | 6.12 × 10−14 | 2.270 × 10−5 | 0.685 |
| 533 | 3.11 × 10−14 | 3.390 × 10−5 | 0.717 |
| 543 | 1.19 × 10−14 | 5.104 × 10−5 | 0.754 |
| 553 | 1.11 × 10−14 | 7.9 × 10−5 | 0.795 |
| 563 | 1.09 × 10−16 | 1.337 × 10−4 | 0.878 |
These results not only validate the theoretical framework but also provide quantitative evidence for the evolution of charge transport processes with temperature.
The temperature-dependent behavior of the parameter s provides critical insights into the nature of charge hopping and transport mechanisms in the material. As reported in the literature, s exhibit distinct variations depending on the underlying physical processes and charge carrier type (e.g., atoms, electrons, polarons). These mechanisms include:47–50
• Quantum mechanical tunneling (QMT): In this model, the exponent s is approximately 0.8 and either remains constant or shows a slight increase with temperature.
• Correlated barrier hopping (CBH): In this case, the exponent s decreases with rising temperature.
• Overlapping large-polaron tunneling (OLPT): where s decreases with increasing temperature to a minimum value, then rises again as the temperature increases.
• Non-overlapping small-polaron Tunneling (NSPT): In this model, s increases with temperature.
For the studied material, the graph in Fig. 15b reveals a discontinuity at 533 K, which corresponds to the phase transition observed in the DSC analysis.
• For the temperature range of 433 K to 533 K, the exponent s decreases as the temperature rises. This trend aligns closely with the CBH mechanism, suggesting that charge carriers hop between localized states separated by correlated potential barriers.
As the temperature increases, the thermal energy enables carriers to overcome these barriers more easily, resulting in enhanced hopping conductivity and a corresponding reduction in s. This behavior reflects a thermally activated transport process.
• For the temperature range of 533 K to 563 K, the (NSPT) model is shown to be appropriate since the exponents increase with increasing temperature.
Overall, the combined analysis of s across the phase transition provides valuable insight into the evolution of charge transport mechanisms: from hopping conduction governed by correlated barriers in the low-temperature phase, to polaronic tunneling influenced by lattice modifications in the high-temperature phase. This understanding is crucial for tailoring material properties to meet the requirements of applications that rely on temperature-dependent electrical conductivity.
![]() | ||
| Fig. 16 (a) Frequency dependence of M′ at different temperatures of [(CH3)4P]3Sb2Br9. (b) Frequency dependence of M′′ at different temperatures of [(CH3)4P]3Sb2Br9. The peak is fitted using the model described by eqn (4). | ||
Fig. 16b illustrates the frequency-dependent behavior of the imaginary part of the complex modulus (M′′). In the low-frequency region, no M′′ peaks are observed, indicating that charge carriers can hop effectively between sites without significant energy barriers. This is consistent with long-range charge transport, where carriers experience minimal impedance. However, at higher frequencies, a pronounced M′′ peak emerges. This peak corresponds to the transition of charge carriers from long-range translational motion to localized motion within their potential energy wells. At these frequencies, charge carriers become confined and can only oscillate locally within the bounds of their energy barriers, reflecting a short-range transport mechanism.52
The position of the M′′ peak shifts to higher frequencies with increasing temperature, signifying thermally activated dielectric relaxation. This behavior indicates that as the temperature rises, the charge carriers gain sufficient thermal energy to overcome energy barriers more rapidly, thereby increasing the relaxation frequency. Moreover, the asymmetric shape of the observed M′′ peaks support the presence of non-Debye relaxation, which deviates from the ideal Debye model. This is a typical characteristic of materials like [(CH3)4P]3Sb2Br9, where complex structural and dynamic interactions among charge carriers and the material lattice lead to a distribution of relaxation times.53
To further analyze the relaxation behavior, numerical simulations of the modulus spectra were performed. These simulations enable the extraction of important parameters, such as activation energy and relaxation frequency, which provide insights into the dynamics of charge carriers. Bergman's frequency-dependent representation of the Kohlrausch–Williams–Watts (KWW) function is given as:54
![]() | (4) |
Here,
represents the maximum value of the imaginary component at the corresponding angular frequency ωmax, while β is the Kohlrausch parameter that quantifies the degree of non-Debye behavior. The parameter β ranges between 0 and 1. In this study, the β < 1 value (Fig. 17) confirms the predominance of non-Debye relaxation.
![]() | (5) |
![]() | ||
| Fig. 18 Arrhenius plots constructed using (a) the σdc conductivity values, (b) the measured Z′′ peak values, (c) the measured M′′ peak values, at different temperatures. | ||
Additional evaluations of Ea were conducted for the studied material using alternative approaches. These involve analyzing the peak values of Z′′ and M′′ at various temperatures. These methods yielded activation energy values of Ea1 = 0.91 eV, Ea1 = 0.92 eV, and Ea2 = 1.05 eV and Ea2 = 1.15 eV, respectively. These results closely align with the value obtained from the dc conductivity and confirm that the same underlying mechanism governs both electric conduction and dielectric polarization in the material.38
Several factors can affect the activation energy, with the ionic radius playing a crucial role in ionic conduction. Ions with smaller ionic radii typically move more easily compared to larger ions. Among the ions present, Sb3+ has the smallest ionic radius, followed by Br−, while the [(CH3)4P]3+ cation has a much larger ionic radius.55,56
Considering the ionic radii and electronegativities of the constituent species, the [Sb2Br9]3− anion framework provides a pathway for limited anionic motion, while the high electronegativity of Br− restricts its mobility. Combined with the impedance and σdc results indicating ionic conduction, this suggests that H+ is the dominant charge carrier.38 Such protonic conduction, facilitated by hydrogen-bond networks, is well-documented in hybrid halometallates.57,58
These electrical results are consistent with the Raman analysis as a function of temperature. In both cases, the modifications observed near the phase transition reveal significant changes in the vibrational dynamics of the lattice, involving both the organic and inorganic sublattices, with more pronounced variations in the inorganic framework. Such coupling between structural vibrations and electrical response highlights the key role of the inorganic moieties in governing the overall behavior of the hybrid system. Such combined structural and electrical features open interesting perspectives for potential applications of this material, whether as a hybrid compound exhibiting phase transitions, as an ionic conductor, or as a functional material for energy-related technologies.
Supplementary information is available. See DOI: https://doi.org/10.1039/d5ra06870c.
| This journal is © The Royal Society of Chemistry 2025 |