
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39127766
69430
10.1038/s41598-024-69430-3
Article
Comprehensive study on structural, electronic, optical, elastic, and transport properties of natural mercury sulphohalides via DFT computation
Hariharan M.
Eithiraj R. D. eithiraj.rd@vit.ac.in

grid.412813.d 0000 0001 0687 4946 Department of Physics, School of Advanced Sciences, Vellore Institute of Technology (VIT), Chennai, Tamil Nadu 600127 India
10 8 2024
10 8 2024
2024
14 1859314 5 2024
5 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
The Mercury Sulphohalides have attracted significant attention in the fields of solar cells and thermoelectric applications. This study delves into the fundamental characteristics, including structural, elasticity, electronic behavior, phonon stability, optical properties, and transport features of AgHgSZ (Z = Br, I) through computational simulations based on Density Functional Theory (DFT) using WIEN2k software. Meticulous calculations of the phonon band structure ensure dynamic stability. The semiconductor nature with indirect band gaps (1.833 eV and 1.832 eV) for Mercury Sulphohalides (Br, I), as revealed by their band structures, suggests diverse photovoltaic and transport applications. Mechanical assessments show stable ductility for AgHgSBr and brittleness for AgHgSI, along with anisotropy and resistance to scratching. Optical properties exhibit anisotropy and significant UV absorption. Analysis of effective masses, exciton binding energy, and exciton Bohr radius suggests low exciton binding energy and classification under Mott-Wannier excitons. Positive thermopower results indicate holes as the predominant charge carriers in AgHgSBr and AgHgSI materials. Moreover, essential thermoelectric factors are examined, revealing the compounds’ potential for thermoelectric applications. Notably, the figure of merit (ZT) at 300 K for AgHgSBr and AgHgSI are calculated to be 0.41 and 0.13, respectively. While these values are low at 300 K, they indicate promising potential for thermoelectric applications at higher temperatures. In summary, this investigation provides valuable understanding into the photovoltaic and thermoelectric properties of AgHgSZ (Z = Br, I) materials, potentially paving the way for further exploration in this domain.

Keywords

Mercury sulphohalides
DFT
Indirect band gap
Elastic anisotropy
Lattice thermal conductivity
ZT
Subject terms

Condensed-matter physics
Theory and computation
Vellore Institute of Technology, ChennaiOpen access funding provided by Vellore Institute of Technology.

issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Natural mercury sulphohalides constitute a distinct group of minerals characterized by their specific geochemistry and crystal structure. Typically found in the oxidized zones of ore deposits containing primary mercury-bearing minerals like cinnabar or fahlore enriched with mercury, these minerals have the capacity to concentrate other elements such as bromine and iodine, which are often sparsely distributed in nature. This mineral family comprises ten distinct species, further categorized into two subgroups. The first subgroup consists around six minerals that have a general formula of Hg3S2Hal2 (where Hal represents halogen elements such as chlorine, bromine, or iodine), including three polymorphs of Hg3S2Cl2 (known as corderoite1, lavrentievite2,3, and kenhsuite4), along with dimorphous arzakite2,3 and grechishchevite5–7 via the ideal composition of Hg3S2Br2, and radtkeite with the formula Hg3S2ICl8–10. The second group is made up of four sulphohalide minerals that include both mercury and silver as metal cations: dimorphs of AgHgSCl (capgaronnite11 and iltisite12), perroudite with the composition Ag4Hg5S5(I,Br)2Cl213–15, and the newly found hanauerite with the formula AgHgSI16. The crystal structure of hanauerite closely resembles that of its synthetic counterpart AgHgSI, as well as similar compounds such as CuHgSeBr and AgHgSBr17. Additionally, high-temperature phases like CuHgSCl18 and CuHgSBr19 share similar orthorhombic crystal structures in space group Pmma, with unit cell dimensions close to those of hanauerite.

Numerous experiments have been conducted on “coinage metal mercury chalcogenide halides”20. Beck et al. investigated the crystal structures of CuHgSBr20, CuHgSCl20, CuHgSeBr, AgHgSBr, and AgHgSI17 using single crystal specimens synthesised hydrothermally in hydrohalic acids at distinctive temperatures. Additionally, we already conducted first principles calculations on CuHgSBr and CuHgSeBr21,22 with a focus on their thermoelectric potential. However, theoretical calculations are not available for AgHgSZ (Z = Br, I); only experimental techniques are available17. Therefore, we conducted first principles calculations on the compounds AgHgSBr and AgHgSI (Mercury Sulphohalides) for thermoelectric applications.

In different energy utilization scenarios such as fossil fuel combustion, electricity generation, and solar energy harvesting, the production of waste heat is an unavoidable byproduct23. Frequently disregarded, this excess heat can pose risks, including the potential for damage to computational and data storage systems due to heat buildup. Thermoelectric (TE) devices provide a solution by converting a portion of this waste heat directly into electrical energy to operate the system, while the remaining electrical energy is used at a specific level to maintain system stability. It can be used to cool areas24–27. Recent research reveals that the ZT value of n-type Mg3Bi2-based materials is about ~ 0.9 at 350 K, and TE coolers can generate temperature differences of about ~ 91 K28, which is in contrast to previous data. This is higher than the reported 81 K29. With no moving parts, noise, or greenhouse gas emissions, solid-state thermoelectric (TE) devices are ideal for miniaturization. The effectiveness of TE devices is dependent on the performance of TE materials, which enable effective heat conversion into electricity and refrigeration, as indicated by high figure-of-merit (ZT) values. At the moment, most TE devices and materials have energy conversion efficiencies of less than 10% and ZT values of about 1.030. However, achieving a ZT value of four could potentially lead to an energy conversion efficiency of 30%, presenting an enticing prospect for practical systems31–33. Furthermore, obtaining a ZT value of three may make thermoelectric (TE) household refrigerators more economically competitive than traditional compressor-based ones. Slack was the first to introduce the concept of a “phonon glass electron crystal” (PGEC), referring to a material with ideal ZT features including high mobility charge carriers and glass-like thermal conductivity. In contrast, fine-tuning specific characteristics in bulk materials to achieve this goal is a difficult process34. In 2012, researchers expanded on the PGEC concept by introducing a fresh technique called as “phonon-liquid electron-crystal” (PLEC). This method leverages superionic conductors' liquid-like behaviour to reduce thermal conductivity (κL) to glass levels35.

For instance, many of the materials with exceptionally high ZT values considered for power generation, such as Bi2Te329,36, PbTe37,38, and GeTe39,40, rely on tellurium, a rare and consequently expensive element. SnSe demonstrates impressive ZT values attributed to its exceptionally low thermal conductivity and advantageous electronic characteristics upon doping41. However, its high-temperature stability is constrained42. Although Half-Heusler materials demonstrate substantial power factors (S2σ), their overall ZT remains constrained by their inherently high thermal conductivity43,44. In the present study, we performed Density Fuctional Theory (DFT) computation on Mercury Sulphohalides AgHgSZ (Z = Br, I) to estimate the various physical properties.

Method of calculations

In this current investigation, computations were carried out on Mercury Sulphohalides represented by the chemical formula AgHgSZ (Z = Br, I) using the WIEN2k code45,46 (Version: 18.2, http://susi.theochem.tuwien.ac.at/). Using Murnaghan's equations, the PBE-GGA method was carefully used to compute a number of parameters, including the bulk modulus, ground state characteristics, and lattice constant. Although this approximation method is effective in determining structural parameters, it has limitations in accurately estimating electronic properties, particularly the band gap. To do exact band gap computations, the TB-mBJ functional was used, which takes less computing time than the HSE06 function47. For inorganic materials, the outcomes of the TB-mBJ and HSE06 functionals are essentially the same. As for organic-based materials with doping compositions that fall below 0.25, the HSE06 potential accurately represents them, whereas the TB-mBJ functional tends to underestimate the band gap values. Despite the fact that the HSE06 potential requires substantial computational resources and time, the TB-mBJ functional is significantly less computationally expensive and time-consuming. In the case of simple inorganic materials, we propose using TB-mBJ because it provides a good mix of reliability and computational capacity48.

Hence, the TB-mBJ functional is utilized in conjuction with PBE-GGA for precise band gap calculations. All ground state properties discussed in this paper are determined using PBE-GGA, based on the volume optimization and total energy. Furthermore, the table includes the reported band gap values obtained with TB-mBJ. For precise computations, essential parameters including “Lmax = 10, Gmax = 12, and RMT x Kmax = 8” are specified as inputs in the WIEN2k package. The muffin tin radii values for Silver (Ag), Mercury (Hg), Sulfur (S), Bromine (Br), and Iodine (I) are 2.5, 2.5, 1.95, 2.45, and 2.5 atomic units (a.u.), respectively. Additonally, a K-point of 100 has been employed to achieve converged parameters, with a merging threshold set at 10–4 Ry. The optical properties were calculated using OPTIC, as implemented in WIEN2k. The effective masses of electrons and holes were estimated from the band structure diagram. Thermoelectric properties are assessed using BoltzTrap code, whereas the phonon thermal conductivity (κL) is elucidated manually by utilizing the “Slack equation”49.

Results and discussion

Crystal structure and structural optimization

AgHgSZ (Z = Br, I) exhibit iso-typical crystallization. The unit cell structure and crystallographic arrangement of AgHgSBr and AgHgSI reveals a planar arrangement comprising Hg and S atoms forming chains, accompanied by markedly distorted planar M2X2 rings, depicted in Figs. 1 and 2. The mercury atoms, denoted as Hg(1) and Hg(2), along with the halogen atoms, occupy specific twofold lattice positions characterized by the symmetry mm2 or 2/m. Conversely, atoms of type M (where M = Cu, Ag) and the chalcogen atoms occupy distinct fourfold lattice positions under the symmetry m. Each elementary unit consists of four formula units. These structures consist of two primary units: a flat chain made up of Hg and S atoms, displaying considerable distortion. The construction of these HgY chains involves two distinct crystallographic configurations of mercury atoms, Hg(1) and Hg(2), along with a chalcogen atom S. Hg(1) and Hg(2) have strikingly comparable coordination environments, with both being essentially linearly coordinated by two chalcogen atoms S. The Hg–S distances in the sulfur-containing compounds range from 2.36 (Å) to 2.38 (Å) are represented in the Table 1, consistent with experimental observations in cinnabar (Hg–S: 2.36 (Å)17).Figure 1 Crystal structure of Mercury Sulphohalides AgHgSZ (Z = Br, I).

Figure 2 Unit cell of Hg-S chains for Mercury Sulphohalides AgHgSZ (Z = Br, I).

Table 1 The computed bond length in the crystal structure for AgHgSBr and AgHgSI compounds.

Material	Bond	Bond distance (Å)
(Theoretical)	Bond distance (Å)
(Experimental17)	
AgHgSBr	Ag–Br	3.1095	3.1090	
Hg–S	2.3626	2.3540	
Ag–S	2.4980	2.4980	
Hg–Hg	3.4957	–	
AgHgSI	Ag–I	3.2093	3.2090	
Hg–S	2.3758	2.3760	
Ag–S	2.5349	2.5350	
Hg–Hg	3.5975	–	

As seen in Fig. 3a,b, we evaluated energy vs volume and used Murnaghan's equation of states to analyse the results in order to find the lattice parameter for AgHgSBr and AgHgSI as well as the best configuration for the unit cell. As compared to the TB-mBJ functional, the equilibrium volume predicted using GGA-PBE closely matches experimental results, demonstrating the superiority of the former. Consequently, we determine all physical characteristics using the GGA approximation. The bandgap and total energy computed using mBJ and GGA-PBE functionals are summarised in Table 2. Furthermore, Table 3 shows the volume variation, equilibrium volume and optimised lattice constants of the GGA-PBE functional, which demonstrates strong agreement with the experimental work.Figure 3 Geometry optimization (volume vs energy) for (a) AgHgSBr and (b) AgHgSI.

Table 2 Computed band gaps and total energy of AgHgSBr and AgHgSI compounds using different functionals.

Material	XC-Functionals	Eg (eV)	E0 (Ry)	
AgHgSBr	GGA-PBE	1.833	−223,880.236	
mBJ	2.717	−223,716.697	
AgHgSI	GGA-PBE	1.832	−259,979.517	
mBJ	2.660	−259,796.289	

Table 3 Elucidated the optimized lattice constants, unit cell volume and volume variation for AgHgSBr and AgHgSI compounds using different functionals.

Material (GGA)	a (Å)	b (Å)	c (Å)	V (Å3)	Volume variation, (%)	
AgHgSBr Experimental17	9.648	4.661	9.426	423.88	−5.766	
Theoretical	9.841	4.754	9.615	449.82	
AgHgSI Experimental17	10.159	4.647	9.849	464.99	−5.721	
Theoretical	10.361	4.739	10.045	493.21	

Electronic properties

Semiconductors, insulators, metals, and semimetals can each be distinguished by their unique electronic properties, which play a crucial role in understanding carrier transport mechanisms. Consequently, electronic band structure calculations can provide valuable insights into the suitability of materials for specific commercial applications. The band structures (BS) of AgHgSZ (Z = Br, I) were analyzed using both the TB-mBJ and PBE-GGA functionals. According to Table 2, the computed indirect band gaps (BS), which are shown in Fig. 4a,b, for AgHgSBr and AgHgSI using PBE-GGA are 1.833 eV and 1.832 eV, respectively. In AgHgSBr, the valence band lies within the asymmetry point, bridging the M and X directions, while conduction band occupies the Γ direction. For AgHgSI, valence band extends between the Γ and X directions, whereas the conduction band remains in the Γ direction. This arrangement indicates an indirect transition that benefits thermoelectric qualities but may have a little impact on optical properties. The bandgaps investigated make these materials promising candidates for solar cells and photovoltaic applications. AgHgSI has a lower band gap than AgHgSBr due to iodine's (I) higher ionic radius. This is explained by a downward shift of the atomic orbitals and a decrease in energy levels inside the conduction band. Moreover, the valence band (VB) shows a relatively flat distribution with fewer states and suggests that holes are the dominant carriers, whereas the conduction band (CB) has a tighter packing of states. This discovery also implies that holes have greater effective masses than electrons, which impacts thermoelectric properties.Figure 4 Energy band structure of (a) AgHgSBr and (b) AgHgSI.

The optical performance of the device is significantly influenced by electron transitions between VB and CB, which are the basis for understanding optoelectronic attributes. Additionally, Fig.s 5a,b, 6a,b demonstrate the density of states for AgHgSBr, AgHgSI, Ag, Hg, S, and Br/I. The CB is dominated by Br d states with energies ranging from 5 to 17.31 eV for AgHgSBr and 2.42–17.31 eV for AgHgSI. The VB band is predominantly made up of Ag d states with energies ranging from 0 to − 5.0 eV, as well as Br p and Hg d states. Consequently, our analysis demonstrates considerable electron transfers from Ag d states to Br d states, altering the optical characteristics of solar cells and other optical consequences through the recombination process. Furthermore, we removed certain orbitals due to their low energy.Figure 5 TDOS of (a) AgHgSBr and (b) AgHgSI using GGA-PBE.

Figure 6 PDOS of (a) AgHgSBr and (b) AgHgSI using GGA-PBE.

Dynamical stability

To confirm the dynamic stability, we carried phonon band structure computations with the PHONOPY algorithm50 combined with the WIEN2k code. The dimensions of the supercell were set as “1 1 1”, while the q-point mesh was specified as 11 × 11 × 11, resulting in 1331 points. Figure 7a,b present the phonon dispersion spectra for AgHgSBr and AgHgSI, showcasing their distinctive behaviors along the k-points direction R-Γ-X-M-Γ. The presence of imaginary (negative) frequency components within these spectra indicates dynamic instability, hinting at the potential for a structural phase transition. As a result, the current configuration is regarded as metastable. The Ag, Br, and Hg atoms are thought to be the origin of this instability as they produce negative frequencies and a significant imaginary mode amplitude. Moving forward, adjustments to the lattice parameter are required, along with a thorough examination of the phonon band structure to assess structural stability. Future investigations could explore the nature of the imaginary mode to understand its implications and investigate potential properties of alternative structures, determining whether the current configuration represents an averaging of two local minima.Figure 7 Phonon band structure (a) AgHgSBr and (b) AgHgSI.

Elastic properties

Studying the mechanical properties of materials is crucial for assessing their potential applications within specific technical domains. These properties, intimately linked with the material's elastic constants, offer insights into various characteristics such as bonding nature, ductility, hardness, and brittleness. In an orthorhombic system, nine elastic constants (C11, C22, C33, C44, C55, C66, C12, C13, and C23) play critical roles in determining elastic behaviour. C11, C22, and C33 represent the material’s compression resistance with axes a, b, and c, whilst the remaining constants (C44, C55, C66, C12, C13, and C23) represent the crystal’s geometric properties. Furthermore, for an orthorhombic system to achieve mechanical stability, it must satisfy the subsequent Born–Huang criteria51 as indicated in Eqs. (1–7).1 C11>0

2 C11C22>C122

3 C11C22C33+2C12C13C23

4 -C11C232-C22C132-C33C122>0

5 C44>0

6 C55>0

7 C66>0

Elastic constant (Cij) for AgHgSZ (Z = Br, I) have been estimated and represented in the Table 4. Materials having larger elastic constants Cij are more resistant to compression, making them less prone to deformation under external pressures. In contrast, materials with smaller elastic constants flex more easily. The C11 and C33 indicate “the resistance to uniform compression in the [100] and [001] directions”. Analysis with GGA-PBE functional demonstrates that the estimated elastic constant C33 for AgHgSBr and AgHgSI crystals exceeds C11, suggesting stronger bonding along the [001] directions than the [100] directions. The elastic constant C44 measures a material's capacity to endure shear deformation generated by tangential stress implemented to the [100] plane in the [010] direction. The estimated C44 value of AgHgSBr is lower than that of C11 and C33, implying that compression along any of its three primary crystallographic directions could lead to material deformation.Table 4 The estimated elastic constants for orthorhombic structure AgHgSBr and AgHgSI using GGA-PBE functional.

Material	C11	C22	C33	C44	C55	C66	C12	C13	C23	
AgHgSBr	40.83	82.48	65.23	25.21	21.88	18.55	28.61	34.14	26.46	
AgHgSI	33.04	103.03	42.51	130.05	42.00	25.99	38.53	32.77	22.76	

Conversely, for AgHgSI, C44 exceeds C11 and C33, suggesting that the material is less susceptible to shear deformation along these principal directions. To determine internal stresses, the Kleinman parameter (ς) was calculated using Eq. (9)52. This measure offers understanding regarding a compound's ability to withstand stresses resulting from stretching and bending forces. The (ς) values for AgHgSBr and AgHgSI suggest a limited impact of bond stretching/contracting in resisting externally applied stress. In AgHgSBr and AgHgSI, mechanical strength is primarily determined by contributions from bond stretching/contracting, with (ς) values of 0.78 and 1.10.8 ς=C11+8C127C11+2C12

The elastic constants Cij were used to calculate a variety of other elastic characteristics important for engineering applications, including as “bulk modulus (B), shear modulus (G), Young's modulus (E), and Poisson's ratio (n)”21, as described in Eqs. 9–12. Table 5 summarises the calculated “bulk moduli (B), shear moduli (G), and Young's moduli (E)” for AgHgSBr and AgHgSI crystals. Remarkably, the Young’s moduli (E) for both materials were notably higher than their bulk moduli (B) and shear moduli (G), indicating increased stiffness. Furthermore, several thermo-physical properties, such as “Young’s modulus (E) and lattice thermal conductivity (κL)”, are interrelated, as expressed by κL ~ E53. The shear moduli (G) for AgHgSBr and AgHgSI were found to be lower than their bulk moduli (B), as shown in Table 5, indicating that shear deformation has the greatest effect on the mechanical strength of the materials examined here.9 B=Bv+BR2

10 G=Gv+GR2

11 E=9GB3B+G

12 n=3B-2G6B+2G

Table 5 The calculated Bulk modulus (GPa), Young’s modulus (GPa), Shear modulus (GPa), Poisson’s ratio, and Pugh’s ratio of AgHgSBr and AgHgSI.

Material	Bv	BR	B	Ev	ER	E	Gv	GR	G	nv	nR	n	B/G.	
AgHgSBr	40.77	37.66	39.21	51.01	42.49	46.75	19.75	16.19	17.97	0.29	0.31	0.30	2.18	
AgHgSI	40.74	32.31	36.52	99.06	12.71	55.88	45.24	4.43	24.83	0.09	0.43	0.22	1.47	

To distinguish between a material's brittle and ductile qualities, one can utilise Poisson's ratio (n) and Pugh's ratio (B/G). A (B/G) ratio greater than 1.75 implies ductility, whereas values less than 1.75 show brittleness. Pure covalent materials typically have a Poisson's ratio of 0.1–0.25, whereas ionic materials have a range of 0.25 to 0.5. As a result, AgHgSBr has ionic properties and ductility, whereas AgHgSI is covalent and brittle.

Furthermore, we conducted an examination of elastic anisotropy, a significant property influencing microcrack formation. For example, Eqs. (13–15) outline the shear anisotropy ratio for the (100), (010), (001) directions. Table 6 summarises the estimated shear anisotropic factors for AgHgSBr and AgHgSI using GGA-PBE. The shear anisotropic ratio values for AgHgSBr, labelled as A010, are less than one, but A100 and A001 are more than one, demonstrating the anisotropy of the AgHgSBr structure. In contrast, for AgHgSI, A100, A010, and A001 are all above one, suggesting the anisotropic nature of the AgHgSI structure. Moreover, the “elastic anisotropy of orthorhombic structures originates from the combination of shear anisotropy and the anisotropy of the uniform bulk modulus”.13 A100=4C44C11+C33-2C13

14 A010=4C55C22+C33-2C23

15 A001=4C66C11+C22-2C12

Table 6 Elucidated shear anisotropy factors (A100, A010, A001) and percent (AB, AG) for AgHgSBr and AgHgSI.

Material	A100	A010	A001	AB	AG	
AgHgSBr	2.66	0.92	1.12	0.03	0.09	
AgHgSI	51.96	1.67	1.76	0.11	0.82	

The variations in bulk modulus with respect to the b-axis along the a- and c-axes can be expressed using Eqs. (16 and 17). In this case, the B and G from the Voigt approximation are represented by BV and GV, respectively, and the B and G from the Reuss approximation by BR and GR.16 AB=BV-BRBV+BR

17 AG=GV-GRGV+GR

When AB = AG = 0 (elastically isotropic), and when AB = AG = 1 (elastic anisotropy). The computed values of AB and AG for AgHgSBr and AgHgSI, as depicted in Table 6, are observed to exceed 0, signifying the anisotropic nature of both AgHgSBr and AgHgSI.

Optical properties

The optical properties of Mercury Sulphohalides AgHgSZ (Z = Br, I) were examined. These properties are contingent upon the interaction between electrons and photons, which governs the transition from the VB to CB. Moreover, the interband transition within semiconductor materials holds considerable importance for applications in optoelectronics54. Furthermore, the optical properties for AgHgSBr and AgHgSI are computed as follows: ɛ(ω) = ɛ1(ω) + iɛ2(ω), where ɛ1(ω) represents the real component and ɛ2(ω) represents the imaginary component of the dielectric function. Several optical parameters including the refractive index n(ω), absorption coefficient α(ω), reflectivity R(ω), extinction coefficient κ(ω) and energy loss L(ω) are determined as function of energy within the range 0–20 eV. These parameters are depicted in Figs. 8a–g and 9a–g. All the three polarization directions (xx, yy, zz) are taken into account. The calculation of ɛ1(ω) relies on the polarization of light intensity using Kramer’s-Kronig relationship.Figure 8 Elucidated optical properties (a) Absorption coefficient, (b) Energy loss, (c) and (d) real and imaginary dielectric constants, (e) Refractive index, (f) Extinction coefficient, and (g) Reflectivity for AgHgSBr.

Figure 9 Elucidated optical properties (a) Absorption coefficient, (b) Energy loss, (c) and (d) real and imaginary dielectric constants, (e) Refractive index, (f) Extinction coefficient, and (g) Reflectivity for AgHgSI.

A material’s absorption coefficient α(ω) reveals its capacity to absorb light at a certain wavelength. Figures 8a and 9a shows the computed energy-absorption relationships for AgHgSBr and AgHgSI, respectively. As energy exceeds the threshold, the absorption coefficient, α(ω), rises. The major absorption peaks for AgHgSBr are roughly 9.70, 11.52, and 6.78 eV, whereas the important peaks for AgHgSI are around 9.02, 10.40, and 6.87 eV, respectively. These materials absorb a lot of UV light, suggesting that they are not very transparent. Based on the absorption coefficient study, the examined compounds show potential for photovoltaic applications.

Figures 8b and 9b illustrate the computed electron energy loss L(ω) for the examined mercury sulphohalides complexes. Plasmon oscillation energy is predicted to occur when L(ω) reaches its peak, which corresponds to the point when the real component hits zero. Peak values represent the energy of plasmon oscillations and transitions between bands. The most notable peak in the L(ω) spectra of AgHgSBr and AgHgSI corresponds to the characteristic plasmon resonance, positioned at 19.95 eV.

The understanding of a material's electronic polarizability can be attained through the real dielectric function, as shown in Figs. 8c and 9c. The static dielectric constants ɛ1(0) along the three crystallographic directions are determined to be 6.14 (xx), 5.50 (yy), and 6.16 (zz) for AgHgSBr, and 7.07 (xx), 6.05 (yy), and 7.06 (zz) for AgHgSI, respectively. The average values of the zero-frequency dielectric constants ɛ1(0) for AgHgSBr and AgHgSI are 5.93 and 6.72, respectively. Nevertheless, there is no experimental polarized zero-frequency dielectric constant available for comparison. These findings distinctly indicate anisotropy in the optical characteristics of orthorhombic AgHgSZ (Z = Br, I). The ratio ε1yy(0)ε1zz(0) is calculated to be 0.892 for AgHgSBr and 0.856 for AgHgSI to assess the degree of anisotropy. Commencing from the zero-frequency limit, in AgHgSBr, they steadily increase until they reach peak values of 8.97 at 2.97 eV (xx), 8.32 at 3.06 eV (yy), and 9.33 at 3.03 eV (zz), respectively. Similarly, in AgHgSI, they begin to ascend and peak at maximum values of 11.26 at 2.73 eV (xx), 9.50 at 2.73 eV (yy), and 11.66 at 2.81 eV (zz). These peaks represent the distribution of light intensity polarisation across different frequencies inside the materials. The polarisation curve reaches its lowest point when there is a little shift in the resonance frequency value. The materials under investigation exhibits semiconductor-like response to incident light. In the case of AgHgSBr, the energy values are 6.87 eV to 20.36 for xx, 8.80 eV to 20.36 eV for yy, and 6.76 eV to 20.36 eV for zz indicates negative values. Similarly, for AgHgSI the energy values are 6.38 eV to 20.36 eV for xx, 8.06 eV to 20.36 eV for yy, 6.21 eV to 20.36 eV exhibits negative values.

Both of the halides (AgHgSBr and AgHgSI) under study undergo a transition where their real dielectric values become negative as the energy level rises over 6 eV. The electromagnetic waves are hindered in their propagation by this negative dispersion feature. The material's characteristics shift from dielectric to metallic as a result. Furthermore, as can be shown in Figs. 8g and 9g, reflectivity must be maximum in the negative regions of epsilon (real) due to the inverse relationship between reflectivity and real epsilon. The energy-wise absorption sections of ɛ2(ω) are depicted in Figs. 8d and 9d. There are direct transitions between bands indicated by peaks in the imaginary components of ɛ2(ω). AgHgSBr and AgHgSI threshold energies are determined to be 1.833 and 1.832 eV, respectively. Both of ɛ2(ω)'s components improve beyond these points. Figures 8e and 9e illustrate a gradual rise in the parallel components of the static refractive index n(0), for each of the materials under investigation as energy escalates, ultimately reaching peak values. However, when energy levels climb, both components of the refractive index decline significantly, occasionally going below unity at energies of 9 eV and 7 eV for AgHgSBr and AgHgSI, respectively. The calculated values for both components of the static refractive index, n(0), which are 2.47 (xx), 2.35 (yy), and 2.48 (zz) for AgHgSBr and 2.65 (xx), 2.46 (yy), and 2.65 (zz) for AgHgSI. Both compounds' plots show a minor anisotropy in the xx, yy, and zz directions.

The imaginary component of refraction, κ(ω), also known as the extinction coefficient, represents the loss in light energy as it travels through materials. Figures 8f and 9f show determined values for κ(ω). These plots closely resemble those of ε2(ω), ensuring the precision of calculations. Whether measured through dielectric constants or the refraction process, the findings yield identical results in terms of the decay of light energy. κ(ω) are depicted in Figs. 8f and 9f. These plots closely resemble those of ε2(ω), ensuring the precision of calculations. Whether measured through dielectric constants or the refraction process, the findings yield identical results in terms of the decay of light energy. The estimated reflectivity R(ω) spectra, shown in Figs. 8g and 9g shows that the compounds under investigation found near-UV reflectivity values ranging from around 30% to 40%. Reflectivity increases with energy levels, peaking at 50% at 20.36 eV. This discovery shows that these materials might be good protective coatings in the energy spectrum with the maximum reflecting characteristics. AgHgSBr has zero-frequency reflectivity values of 18%, 16%, and 18% at xx, yy, and zz. Similarly, AgHgSI possesses zero-frequency reflectivity values of 20%, 17%, and 20% for xx, yy, and zz, respectively.

Transport properties

Band structures made simpler by comparing their effective mass to the behaviour of a free element of that mass. We study the effective mass of AgHgSZ (Z = Br, I) compounds to examine the influence of substitution on conductivity, using the theory of effective mass as a direct tool for understanding the electrical characteristics of materials. In order to calculate the effective band mass, one must look at the band's curvature at point Γ, which is the extreme point in k-space. To put it simply, parabolic curves given in Eq. (18) match certain band mass components, including the valence band maximum (VBM) and conduction band minimum (CBM) seen in Fig. 10a,b:18 1meff=1ħ2∂2E∂k2

Figure 10 The estimated VBM and CBM for (a) AgHgSBr and (b) AgHgSI.

The Table 7 shows effective masses (meff) that in AgHgSBr and AgHgSI, electrons have lower effective masses than holes, implying that electrons may have more mobility. The band structures shown in Fig. 4 for AgHgSBr and AgHgSI indicate a p-type semiconductor, implying that holes are the favoured charge carriers. As a result, holes have a higher effective mass than electrons do. Equation (19) makes it easier to calculate exciton binding energy, with values for AgHgSBr and AgHgSI estimated at 36.3, 45.2, and 36.1 meV, as well as 24.2, 32.9, and 24.3 meV across the x, y, and z directions of the ɛ1(0). Using Eq. (20)55, the exciton Bohr radius may be computed, obtaining values for AgHgSBr and AgHgSI of 32.24 Å, 28.87 Å, and 32.34 Å, 42.05 Å, 36.03 Å, and 41.99 Å, in the x, y, and z directions. Interestingly, AgHgSBr and AgHgSI indicates a lower exciton binding energy throughout, with its a* value exceeding the lattice parameters exhibits as a weak exciton of the Mott-Wannier type.19 Eb=13.6mμm01ε102

20 a∗=ε10m0mμa0

Table 7 The estimated effective masses for AgHgSBr and AgHgSI compounds.

Material	Electron	Hole	
AgHgSBr	0.1210 me	0.6051 me	
AgHgSI	0.1008 me	0.7563 me	

The exploration of the electronic behavior and optical properties suggests that AgHgSBr and AgHgSI hold potential for photovoltaic applications. Harnessing waste heat to generate electrical energy and prevent electronic devices from overheating by converting this thermal energy into electricity represents an emerging domain within renewable energy. The effectiveness of AgHgSBr and AgHgSI are evaluated using figure of merit (ZT), calculated as (S2σ/κ). Therefore, achieving high ZT values relies on maximizing the potential gradient between the metal contacts, which is computed by the Seebeck coefficient (S) and electrical conductivity (σ) at room temperature. The thermal conductivity, comprising electronic (κe) and phononic contributions (κL), is also addressed in this context.

Thermoelectric evaluations for AgHgSBr and AgHgSI are done using Boltzmann transport theory, which is based on frequency band dispersion. This study examines σ, κe, and κL in relation to relaxation time (t = 10–14 s), which is factored into the computation method. The thermoelectric parameters are illustrated in Figs. 11a–d, 12a,b, thoroughly discussed below. In this research, the investigators examined the fluctuations in thermoelectric (TE) properties across a temperature range spanning from 50 to 1000 K.Figure 11 The estimated thermoelectric properties (a) Seebeck coefficient, (b) Electrical conductivity, (c) Electronic thermal conductivity, and (d) Lattice thermal conductivity for AgHgSBr and AgHgSI.

Figure 12 The estimated (a) Power factor, and (b) ZT for AgHgSBr and AgHgSI.

The Seebeck coefficient (S), also known as thermopower, is a thermoelectric material's capacity to create an electric current in reaction to a change in temperature. This factor is crucial in assessing the effectiveness of thermocouples and other thermoelectric apparatus. In Fig. 11a, the Seebeck coefficient for AgHgSZ (Z = Br, I) compounds predominantly carry holes as charge carriers. This implies a surplus of positively charged carriers compared to electrons within these materials. When electrons went from VB to CB, holes are created. The Seebeck coefficient's positive value denotes an excess of holes relative to electrons since it is proportional to the number of charge carriers. Because of their weaker atomic affinity, holes are more mobile than electrons, which increases how useful they are for producing energy. This positive thermopower (S) raises the prospect of thermoelectric applications, wherein devices that generate power from temperature changes employ the Seebeck effect. Higher S denote increased thermoelectric device efficiency. Therefore, the positive Seebeck coefficient observed in AgHgSZ (Z = Br, I) compounds underscores their promise for such applications. Additionally, the Fig. 11a illustrates substantial variations in the Seebeck coefficient of AgHgSBr and AgHgSI with rising temperature. At ambient temperature, AgHgSBr and AgHgSI exhibit Seebeck coefficients of 176 μV/K and 167 μV/K, respectively. However, at a significantly elevated temperature of 1000 K, the S values decrease to 145 μV/K for AgHgSBr and 156 μV/K for AgHgSI. The Seebeck coefficient of AgHgSBr demonstrates a gradual rise up to 300 K, followed by a slow decline until reaching 1000 K. Similarly, for AgHgSI, the S exhibits an increase up to 350 K, with a subsequent gradual decrease observed up to 1000 K.

The electrical conductivity (σ) of thermoelectric materials pertains to their ability to transport electrical charge. Typically, semiconductor materials can have their electrical conductivity adjusted through a process known as doping, which entails introducing impurities to modify their electrical properties. The level of electrical conductivity holds significance in thermoelectric materials as it directly impacts the efficiency of the thermoelectric phenomenon. In our investigation, the electrical conductivity of AgHgSZ (Z = Br, I) compounds under scrutiny is illustrated in Fig. 11b. The graph depicts a linear augmentation in electrical conductivity for both materials as temperature increases, attributed to the escalating clusters of charge carriers. Notably, AgHgSI demonstrates higher electrical conductivity compared to AgHgSBr, a characteristic contrasting with the temperature-dependent alterations observed in thermopower a common trait of semiconducting materials. At 1000 K, both materials exhibit a maximum electrical conductivity nearing 7.890 × 10–18 1/(Ω.cm.s) and 8.940 × 10–18 1/(Ω.cm.s), respectively.

Thermal conductivity refers to a material’s capacity to transmit heat by measuring the rate of heat transfer across a given area per unit temperature differential. Watts per metre Kelvin (W/mK) is a popular unit of measurement. Elevated thermal conductivity indicates enhanced efficiency in heat propagation. Consequently, materials exhibiting high thermal conductivity find utility in applications necessitating rapid heat dissipation, such as heat sinks. Conversely, materials characterized by low thermal conductivity serve as effective insulators, impeding the transfer of heat between distinct locations. The lattice vibrations associated with atomic mobility in semiconductor materials increase as temperatures rise. Furthermore, the movement of heat energy by convection from locations with greater to lower potential accentuates this impact. The interplay between phonons and electrons estimates semiconductor thermal conductivity. The temperature dependent electronic thermal conductivity (κe) as shown in Fig. 11c, demonstrating that the κe of AgHgSZ (Z = Br, I) materials escalates with increasing temperature. The value κe at room temperature for AgHgSBr and AgHgSI are 5.91 × 10–13 (W/mKs) and 6.28 × 10–13 (W/mKs).

The lattice thermal conductivity stands as a fundamental property of solids. In this study, we have examined the fluctuation in the lattice thermal conductivity (κL) of our compounds across varying temperatures as shown in Fig. 11d, employing the Slack model49 in Eq. (21).21 κl=AMaθ3δγ2Tn2/3

In the aforementioned equation, the symbols “A, Ma, θ, δ, γ, and n stand for respective physical parameters: average atomic mass, Debye temperature, cube root of the average atomic volume, Grüneisen parameter, and the number of atoms in a primitive unit cell”21. The determination of the physical quantity A involves a calculation based on the following relationship in Eq. (22).22 A=2.43×10-81-0.514γ+0.228γ2

For AgHgSBr and AgHgSI, we determined the Grüneisen parameter (γ) and Debye temperature (θ) using the WIEN2k-interfaced IRelast package45. For AgHgSBr, we determined θ to be 184.28 K, accompanied by a γ value of 1.78. Similarly, for AgHgSI, we calculated θ as 206.71 K, with a corresponding γ of 1.37. As temperature rises, the lattice contribution to thermal conductivity (κL) decreases in both compounds. Specifically, at 300 K, κL is approximately 1.01 (W/mK) for AgHgSBr, declining to 0.30 (W/mK) at 1000 K. Similarly, for AgHgSBr, κL is around 2.86 (W/mK) at 300 K, dropping to 0.85 (W/mK) at 1000 K. Presently, there are no experimental values available for comparison. When compared to other quaternary coinage compounds that have low lattice thermal conductivity, such as the newly found CuHgSeBr (with a value of 0.69 W/mK at 300 K), it appears to be lower than AgHgSBr and AgHgSI. The increase in electronic thermal conductivity is caused by the formation of electron–hole pairs as a result of heat, resulting in a decrease in the mean free path of phonons and a considerable drop in κL.

Furthermore, Fig. 12a indicates the computed relationship of the power factor (PF) with temperature, it is calculated as PF = S2σ. As seen in Fig. 12a, the PF increases significantly as temperature rises, culminating at high temperatures. This tendency results from the simultaneous increase in carrier concentration and growing temperature. Peak power factor is associated with a lower S and higher σ. This behaviour highlights how important it is to strike a balance between these two elements in order to improve thermoelectric material power production efficiency.

Without a doubt, the most crucial metric for determining the thermoelectric characteristics of materials is the ZT. Using formula ZT = S2σ/κ. Figure 12b shows how the figure of merit (ZT) fluctuates in response to temperature. Notably, it is discernible that both the AgHgSBr and AgHgSI compounds exhibit an upward trend in ZT with increasing temperature. This observed augmentation in ZT with temperature for both AgHgSBr and AgHgSI suggests an enhancement in their thermoelectric efficiency as temperature escalates. The recorded ZT of 0.41 and 0.13 at room temperature for AgHgSBr and AgHgSI, respectively, underscore the significance of these studied materials. Furthermore, juxtaposing with CuHgSeBr22, the ZT values are notably higher for AgHgSBr and comparatively lower for AgHgSI. Based on this ZT analysis, our investigated compounds demonstrate suitability for thermoelectric applications, particularly at elevated temperatures.

Conclusion

In this comprehensive investigation, we delve deeply into the optical, electronic, elastic, and thermoelectric properties of AgHgSBr and AgHgSI materials through meticulous density functional theory (DFT) computations. Employing the GGA-PBE method, we establish band gap values of 1.833 eV for AgHgSBr and 1.832 eV for AgHgSI. Furthermore, we evaluate the dynamic and structural stability of these compounds. Our analysis of their elastic properties confirms the mechanical resilience and ductility of AgHgSBr, in stark contrast to the brittle nature of AgHgSI, indicative of a prevalent combination of ionic and covalent bonding traits. Additionally, considering their anticipated band gap, high refractive index, and significant absorption coefficient, particularly in the UV spectrum, our findings unveil compelling prospects for photovoltaic applications. Moreover, we meticulously investigate the thermoelectric characteristics of AgHgSZ (Z = Br, I) using the semi-classical BoltzTrap method. Through a detailed examination of the band structure, we derive crucial parameters such as effective mass and exciton binding energy. Depending on the exciton Bohr radius, our results suggest the presence of a Mott-Wannier type excitons. Notably, both materials show considerable potential as thermoelectric devices, as evidenced by power factor and ZT measurements. This extensive study significantly advances the field of materials science by enriching our understanding and application of AgHgSBr and AgHgSI materials in photovoltaic and thermoelectric domains. Moreover, our theoretical findings provide a foundation for experimental researchers to further explore the thermoelectric and optical properties of these materials.

Author contributions

R.E. provided ideas and drafted the manuscript. M. H. calculated the structural and electronic properties in detail.

Funding

Open access funding provided by Vellore Institute of Technology.

Data availability

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Foord EE Berendsen P Storey LO Corderoite, first natural occurrence of α-Hg3S2CI2, from the Cordero mercury deposit, Humboldt County, Nevada1 Am. Mineral. J. Earth Planet. Mater. 1974 59 7–8 652 655
Foord, E. E., Berendsen, P. & Storey, L. O. Corderoite, first natural occurrence of α-Hg3S2CI2, from the Cordero mercury deposit, Humboldt County, Nevada1. Am. Mineral. J. Earth Planet. Mater. 59(7–8), 652–655 (1974).
2. Vasil’ev VL Pal’chik NA Grechishev OK Lavrentievite and arzakite, new natural sulfohalogenides of mercury Geol. Geofizika 1984 7 54 63
Vasil’ev, V. L., Pal’chik, N. A. & Grechishev, O. K. Lavrentievite and arzakite, new natural sulfohalogenides of mercury. Geol. Geofizika 7, 54–63 (1984).
3. Vasil’ev VI Lavrent’ev YG Pal’chik NA New data on arzakite and lavrentievite Dokl. Akad. Nauk. SSSR 1986 290 948 951
Vasil’ev, V. I., Lavrent’ev, Y. G. & Pal’chik, N. A. New data on arzakite and lavrentievite. Dokl. Akad. Nauk. SSSR 290, 948–951 (1986).
4. McCormack JK Dickson FW Kenhsuite, gamma-Hg3S2Cl2, a new mineral species from the McDermitt mercury deposit, Humboldt County, Nevada Can. Mineral. 1998 36 1 201 206
McCormack, J. K. & Dickson, F. W. Kenhsuite, gamma-Hg3S2Cl2, a new mineral species from the McDermitt mercury deposit, Humboldt County, Nevada. Can. Mineral. 36(1), 201–206 (1998).
5. Vasil’ev VI Usova LV Pal’chik NA Grechishchevite—Hg3S2(Br, Cl, I)2—A new supergene mercury sulfohalide Geol. Geofiz 1989 30 61 69
Vasil’ev, V. I., Usova, L. V. & Pal’chik, N. A. Grechishchevite—Hg3S2(Br, Cl, I)2—A new supergene mercury sulfohalide. Geol. Geofiz 30, 61–69 (1989).
6. Pervukhina, N.V., Borisov, S.V., Magarill, S.A., Naumov, D.Y. & Vasil'ev, V.I. The crystal structure of a polymorph of Hg2+ 3S2Br1.0Cl0.5I0.5 (2003).
7. Pervukhina NV Borisov SV Magarill SA Naumov DY Vasiliev VI Nenashev BG Crystal structure of the synthetic analog of grechishchevite and Hg-X (X= S, Se, Te) bond topology in mercury chalcogen halides J. Struct. Chem. 2004 45 437 445 10.1007/s10947-005-0011-z
Pervukhina, N. V. et al. Crystal structure of the synthetic analog of grechishchevite and Hg-X (X= S, Se, Te) bond topology in mercury chalcogen halides. J. Struct. Chem. 45, 437–445 (2004).10.1007/s10947-005-0011-z
8. McCormack JK Dickson FW Leshendok MP Radtkeite, Hg3S2CII, a new mineral from the McDermitt mercury deposit, Humboldt County, Nevada Am. Mineral. 1991 76 9–10 1715 1721
McCormack, J. K., Dickson, F. W. & Leshendok, M. P. Radtkeite, Hg3S2CII, a new mineral from the McDermitt mercury deposit, Humboldt County, Nevada. Am. Mineral. 76(9–10), 1715–1721 (1991).
9. Pervukhina NV Vasil’ev VI Naumov DY Borisov SV Magarill SA The crystal structure of synthetic radtkeite, Hg3S2ClI Can. Mineral. 2004 42 1 87 94 10.2113/gscanmin.42.1.87
Pervukhina, N. V., Vasil’ev, V. I., Naumov, D. Y., Borisov, S. V. & Magarill, S. A. The crystal structure of synthetic radtkeite, Hg3S2ClI. Can. Mineral. 42(1), 87–94 (2004).10.2113/gscanmin.42.1.87
10. Pervukhina NV Borisov SV Magarill SA Naumov DY Vasiliev VI Nenashev BG Crystal structure of the synthetic analog of radtkeite Hg3S2Cl 1.00 I 1.00 J. Struct. Chem. 2004 45 720 723 10.1007/s10947-005-0050-5
Pervukhina, N. V. et al. Crystal structure of the synthetic analog of radtkeite Hg3S2Cl 1.00 I 1.00. J. Struct. Chem. 45, 720–723 (2004).10.1007/s10947-005-0050-5
11. Mason B Mumme WG Sarp H Capgaronnite, HgS·Ag (Cl, Br, I), a new sulfide-halide mineral from Var, France Am. Mineral. 1992 77 1–2 197 200
Mason, B., Mumme, W. G. & Sarp, H. Capgaronnite, HgS·Ag (Cl, Br, I), a new sulfide-halide mineral from Var, France. Am. Mineral. 77(1–2), 197–200 (1992).
12. Sarp H Sanz-Gysler J Perroud P L'iltisite, HgSAg (Cl, Br), un nouveau minéral de la mine de Cap Garonne, Var (France), polymorphe de la capgaronnite Archiv. Sci. 1997 50 1 5
Sarp, H., Sanz-Gysler, J. & Perroud, P. L’iltisite, HgSAg (Cl, Br), un nouveau minéral de la mine de Cap Garonne, Var (France), polymorphe de la capgaronnite. Archiv. Sci. 50, 1–5 (1997).
13. Sarp H Birch WD Hlava PF Pring A Sewell DK Perroudite, a new sulfide-halide of Hg and Ag from Cap-Garonne, Var, France, and from Broken Hill, New South Wales, and Coppin Pool, Western Australia Am. Mineral. 1987 72 11–12 1251 1256
Sarp, H., Birch, W. D., Hlava, P. F., Pring, A. & Sewell, D. K. Perroudite, a new sulfide-halide of Hg and Ag from Cap-Garonne, Var, France, and from Broken Hill, New South Wales, and Coppin Pool, Western Australia. Am. Mineral. 72(11–12), 1251–1256 (1987).
14. Mumme WG Nickel EH Crystal structure and crystal chemistry of perroudite; a mineral from Coppin Pool, Western Australia Am. Mineral. 1987 72 11–12 1257 1262
Mumme, W. G. & Nickel, E. H. Crystal structure and crystal chemistry of perroudite; a mineral from Coppin Pool, Western Australia. Am. Mineral. 72(11–12), 1257–1262 (1987).
15. Keller P Lissner F Schleid T Single-crystal structure determination of perroudite, Hg5Ag4S5(I, Br)2Cl2, from Tsumeb (Namibia), and its structural relationships to other sulfide halides of mercury and cinnabar Neues Jahrbuch fur Mineralogie-Abhandlungen 2005 181 1 1 10 10.1127/0077-7757/2005/0181-0004
Keller, P., Lissner, F. & Schleid, T. Single-crystal structure determination of perroudite, Hg5Ag4S5(I, Br)2Cl2, from Tsumeb (Namibia), and its structural relationships to other sulfide halides of mercury and cinnabar. Neues Jahrbuch fur Mineralogie-Abhandlungen 181(1), 1–10 (2005).10.1127/0077-7757/2005/0181-0004
16. Pekov IV Zubkova NV Britvin SN Agakhanov AA Polekhovsky YS Pushcharovsky DY Möhn G Desor J Blass G A new mineral hanauerite, AgHgSI, and common crystal chemical features of natural mercury sulphohalides Crystals 2023 13 8 1218 10.3390/cryst13081218
Pekov, I. V. et al. A new mineral hanauerite, AgHgSI, and common crystal chemical features of natural mercury sulphohalides. Crystals 13(8), 1218 (2023).10.3390/cryst13081218
17. Beck J Keller HL Rompel M Wimbert L Hydrothermalsynthese und Kristallstruktur der Münzmetall-Quecksilber-Chalkogenidhalogenide CuHgSeBr, AgHgSBr und AgHgSI Zeitschrift für anorganische und allgemeine Chemie 2001 627 9 2289 2294 10.1002/1521-3749(200109)627:9<2289::AID-ZAAC2289>3.0.CO;2-E
Beck, J., Keller, H. L., Rompel, M. & Wimbert, L. Hydrothermalsynthese und Kristallstruktur der Münzmetall-Quecksilber-Chalkogenidhalogenide CuHgSeBr, AgHgSBr und AgHgSI. Zeitschrift für anorganische und allgemeine Chemie 627(9), 2289–2294 (2001).10.1002/1521-3749(200109)627:9<2289::AID-ZAAC2289>3.0.CO;2-E
18. Moro’oka M Ohki H Yamada K Okuda T Crystal structure of high temperature phase and ionic conductivity mechanism of CuHgSX (X= Cl, Br) Bull. Chem. Soc. Jpn. 2003 76 11 2111 2115 10.1246/bcsj.76.2111
Moro’oka, M., Ohki, H., Yamada, K. & Okuda, T. Crystal structure of high temperature phase and ionic conductivity mechanism of CuHgSX (X= Cl, Br). Bull. Chem. Soc. Jpn. 76(11), 2111–2115 (2003).10.1246/bcsj.76.2111
19. Beck J Rompel M Über Münzmetall-Quecksilber-Chalkogenidhalogenide. II Hydrothermalsynthese, Kristallstruktur und Phasenumwandlung von CuHgSCl und CuHgSBr Zeitschrift für anorganische und allgemeine Chemie 2003 629 3 421 428 10.1002/zaac.200390070
Beck, J. & Rompel, M. Über Münzmetall-Quecksilber-Chalkogenidhalogenide. II Hydrothermalsynthese, Kristallstruktur und Phasenumwandlung von CuHgSCl und CuHgSBr. Zeitschrift für anorganische und allgemeine Chemie 629(3), 421–428 (2003).10.1002/zaac.200390070
20. Beck J Rompel M On coinage metal mercury chalcogenide halides II: Hydrothermal synthesis, crystal structure, and solid state phase transition of CuHgSCl and CuHgSBr Zeitschrift fuer Anorganische und Allgemeine Chemie 2003 629 421 428 10.1002/zaac.200390070
Beck, J. & Rompel, M. On coinage metal mercury chalcogenide halides II: Hydrothermal synthesis, crystal structure, and solid state phase transition of CuHgSCl and CuHgSBr. Zeitschrift fuer Anorganische und Allgemeine Chemie 629, 421–428 (2003).10.1002/zaac.200390070
21. Hariharan M Eithiraj RD Investigations on the thermoelectric and thermodynamic properties of quaternary coinage metal HgSBr Heliyon 2023 9 9 e14938 37095997
Hariharan, M. & Eithiraj, R. D. Investigations on the thermoelectric and thermodynamic properties of quaternary coinage metal HgSBr. Heliyon 9(9), e14938 (2023).37095997
22. Hariharan M Investigations on the electronic, optical, thermoelectric and thermodynamic properties of quaternary coinage metal HgSeBr: A DFT study Phys. Scripta 2024 99 035937 10.1088/1402-4896/ad26c0
Hariharan, M. Investigations on the electronic, optical, thermoelectric and thermodynamic properties of quaternary coinage metal HgSeBr: A DFT study. Phys. Scripta 99, 035937 (2024).10.1088/1402-4896/ad26c0
23. Bell LE Cooling, heating, generating power, and recovering waste heat with thermoelectric systems Science 2008 321 5895 1457 1461 10.1126/science.1158899 18787160
Bell, L. E. Cooling, heating, generating power, and recovering waste heat with thermoelectric systems. Science 321(5895), 1457–1461 (2008).18787160 10.1126/science.1158899
24. DiSalvo FJ Thermoelectric cooling and power generation Science 1999 285 5428 703 706 10.1126/science.285.5428.703 10426986
DiSalvo, F. J. Thermoelectric cooling and power generation. Science 285(5428), 703–706 (1999).10426986 10.1126/science.285.5428.703
25. Liu W Jie Q Kim HS Ren Z Current progress and future challenges in thermoelectric power generation: From materials to devices Acta Mater. 2015 87 357 376 10.1016/j.actamat.2014.12.042
Liu, W., Jie, Q., Kim, H. S. & Ren, Z. Current progress and future challenges in thermoelectric power generation: From materials to devices. Acta Mater. 87, 357–376 (2015).10.1016/j.actamat.2014.12.042
26. He W Zhang G Zhang X Ji J Li G Zhao X Recent development and application of thermoelectric generator and cooler Appl. Energy 2015 143 1 25 10.1016/j.apenergy.2014.12.075
He, W. et al. Recent development and application of thermoelectric generator and cooler. Appl. Energy 143, 1–25 (2015).10.1016/j.apenergy.2014.12.075
27. Chowdhury I Prasher R Lofgreen K Chrysler G Narasimhan S Mahajan R Koester D Alley R Venkatasubramanian R On-chip cooling by superlattice-based thin-film thermoelectrics Nat. Nanotechnol. 2009 4 4 235 238 10.1038/nnano.2008.417 19350033
Chowdhury, I. et al. On-chip cooling by superlattice-based thin-film thermoelectrics. Nat. Nanotechnol. 4(4), 235–238 (2009).19350033 10.1038/nnano.2008.417
28. Mao J Zhu H Ding Z Liu Z Gamage GA Chen G Ren Z High thermoelectric cooling performance of n-type Mg3Bi2-based materials Science 2019 365 6452 495 498 10.1126/science.aax7792 31320557
Mao, J. et al. High thermoelectric cooling performance of n-type Mg3Bi2-based materials. Science 365(6452), 495–498 (2019).31320557 10.1126/science.aax7792
29. Kim SI Lee KH Mun HA Kim HS Hwang SW Roh JW Yang DJ Shin WH Li XS Lee YH Snyder GJ Dense dislocation arrays embedded in grain boundaries for high-performance bulk thermoelectrics Science 2015 348 6230 109 114 10.1126/science.aaa4166 25838382
Kim, S. I. et al. Dense dislocation arrays embedded in grain boundaries for high-performance bulk thermoelectrics. Science 348(6230), 109–114 (2015).25838382 10.1126/science.aaa4166
30. Vining CB An inconvenient truth about thermoelectrics Nat. Mater. 2009 8 2 83 85 10.1038/nmat2361 19165205
Vining, C. B. An inconvenient truth about thermoelectrics. Nat. Mater. 8(2), 83–85 (2009).19165205 10.1038/nmat2361
31. Sales BC Mandrus DWRK Williams RK Filled skutterudite antimonides: A new class of thermoelectric materials Science 1996 272 5266 1325 1328 10.1126/science.272.5266.1325 8662465
Sales, B. C., Mandrus, D. W. R. K. & Williams, R. K. Filled skutterudite antimonides: A new class of thermoelectric materials. Science 272(5266), 1325–1328 (1996).8662465 10.1126/science.272.5266.1325
32. Kim HS Liu W Chen G Chu CW Ren Z Relationship between thermoelectric figure of merit and energy conversion efficiency Proc. Natl. Acad. Sci. 2015 112 27 8205 8210 10.1073/pnas.1510231112 26100905
Kim, H. S., Liu, W., Chen, G., Chu, C. W. & Ren, Z. Relationship between thermoelectric figure of merit and energy conversion efficiency. Proc. Natl. Acad. Sci. 112(27), 8205–8210 (2015).26100905 10.1073/pnas.1510231112
33. Goldsmid HJ Douglas RW The use of semiconductors in thermoelectric refrigeration Br. J. Appl. Phys. 1954 5 11 386 10.1088/0508-3443/5/11/303
Goldsmid, H. J. & Douglas, R. W. The use of semiconductors in thermoelectric refrigeration. Br. J. Appl. Phys. 5(11), 386 (1954).10.1088/0508-3443/5/11/303
34. Rowe DM CRC Handbook of Thermoelectrics 2018 CRC Press
Rowe, D. M. (ed.) CRC Handbook of Thermoelectrics (CRC Press, 2018).
35. Zhao K Qiu P Shi X Chen L Recent advances in liquid-like thermoelectric materials Adv. Funct. Mater. 2020 30 8 1903867 10.1002/adfm.201903867
Zhao, K., Qiu, P., Shi, X. & Chen, L. Recent advances in liquid-like thermoelectric materials. Adv. Funct. Mater. 30(8), 1903867 (2020).10.1002/adfm.201903867
36. Lv HY Liu HJ Shi J Tang XF Uher C Optimized thermoelectric performance of Bi2Te3 nanowires J. Mater. Chem. A 2013 1 23 6831 6838 10.1039/c3ta10804j
Lv, H. Y., Liu, H. J., Shi, J., Tang, X. F. & Uher, C. Optimized thermoelectric performance of Bi2Te3 nanowires. J. Mater. Chem. A 1(23), 6831–6838 (2013).10.1039/c3ta10804j
37. Wu Y Nan P Chen Z Zeng Z Liu R Dong H Xie L Xiao Y Chen Z Gu H Li W Thermoelectric enhancements in PbTe alloys due to dislocation-induced strains and converged bands Adv. Sci. 2020 7 12 1902628 10.1002/advs.201902628
Wu, Y. et al. Thermoelectric enhancements in PbTe alloys due to dislocation-induced strains and converged bands. Adv. Sci. 7(12), 1902628 (2020).10.1002/advs.201902628
38. Biswas K He J Blum ID Wu CI Hogan TP Seidman DN Dravid VP Kanatzidis MG High-performance bulk thermoelectrics with all-scale hierarchical architectures Nature 2012 489 7416 414 418 10.1038/nature11439 22996556
Biswas, K. et al. High-performance bulk thermoelectrics with all-scale hierarchical architectures. Nature 489(7416), 414–418 (2012).22996556 10.1038/nature11439
39. Suwardi A Cao J Hu L Wei F Wu J Zhao Y Lim SH Yang L Tan XY Chien SW Yin Y Tailoring the phase transition temperature to achieve high-performance cubic GeTe-based thermoelectrics J. Mater. Chem. A 2020 8 36 18880 18890 10.1039/D0TA06013E
Suwardi, A. et al. Tailoring the phase transition temperature to achieve high-performance cubic GeTe-based thermoelectrics. J. Mater. Chem. A 8(36), 18880–18890 (2020).10.1039/D0TA06013E
40. Zhang X Bu Z Lin S Chen Z Li W Pei Y GeTe thermoelectrics Joule 2020 4 5 986 1003 10.1016/j.joule.2020.03.004
Zhang, X. et al. GeTe thermoelectrics. Joule 4(5), 986–1003 (2020).10.1016/j.joule.2020.03.004
41. Zhao LD Lo SH Zhang Y Sun H Tan G Uher C Wolverton C Dravid VP Kanatzidis MG Ultralow thermal conductivity and high thermoelectric figure of merit in SnSe crystals Nature 2014 508 7496 373 377 10.1038/nature13184 24740068
Zhao, L. D. et al. Ultralow thermal conductivity and high thermoelectric figure of merit in SnSe crystals. Nature 508(7496), 373–377 (2014).24740068 10.1038/nature13184
42. Popuri SR Pollet M Decourt R Morrison FD Bennett NS Bos JWG Large thermoelectric power factors and impact of texturing on the thermal conductivity in polycrystalline SnSe J. Mater. Chem. C 2016 4 8 1685 1691 10.1039/C6TC00204H
Popuri, S. R. et al. Large thermoelectric power factors and impact of texturing on the thermal conductivity in polycrystalline SnSe. J. Mater. Chem. C 4(8), 1685–1691 (2016).10.1039/C6TC00204H
43. Quinn RJ Bos JWG Recent progress in phosphide materials for thermoelectric conversion J. Mater. Chem. A 2023 11 16 8453 8469 10.1039/D3TA00620D
Quinn, R. J. & Bos, J. W. G. Recent progress in phosphide materials for thermoelectric conversion. J. Mater. Chem. A 11(16), 8453–8469 (2023).10.1039/D3TA00620D
44. Quinn RJ Stenning GB Bos JWG Electronic scattering in half-Heusler thermoelectrics from resistivity data J. Phys. Energy 2022 4 2 024005 10.1088/2515-7655/ac5f37
Quinn, R. J., Stenning, G. B. & Bos, J. W. G. Electronic scattering in half-Heusler thermoelectrics from resistivity data. J. Phys. Energy 4(2), 024005 (2022).10.1088/2515-7655/ac5f37
45. Madsen GK Blaha P Schwarz K Sjöstedt E Nordström L Efficient linearization of the augmented plane-wave method Phys. Rev. B 2001 64 19 195134 10.1103/PhysRevB.64.195134
Madsen, G. K., Blaha, P., Schwarz, K., Sjöstedt, E. & Nordström, L. Efficient linearization of the augmented plane-wave method. Phys. Rev. B 64(19), 195134 (2001).10.1103/PhysRevB.64.195134
46. Schwarz K Blaha P Madsen GK Electronic structure calculations of solids using the WIEN2k package for material sciences Comput. Phys. Commun. 2002 147 1–2 71 76 10.1016/S0010-4655(02)00206-0
Schwarz, K., Blaha, P. & Madsen, G. K. Electronic structure calculations of solids using the WIEN2k package for material sciences. Comput. Phys. Commun. 147(1–2), 71–76 (2002).10.1016/S0010-4655(02)00206-0
47. Perdew JP Burke K Ernzerhof M Generalized gradient approximation made simple Phys. Rev. Lett. 1996 77 18 3865 10.1103/PhysRevLett.77.3865 10062328
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77(18), 3865 (1996).10062328 10.1103/PhysRevLett.77.3865
48. Ashiq MGB Albalawi H Study of electronic, mechanical, thermoelectric and optical properties of K2ScAuZ6 (Z= Br, I) for energy harvesting Eur. Phys. J. Plus 2023 138 6 1 12 10.1140/epjp/s13360-023-04136-5
Ashiq, M. G. B. & Albalawi, H. Study of electronic, mechanical, thermoelectric and optical properties of K2ScAuZ6 (Z= Br, I) for energy harvesting. Eur. Phys. J. Plus 138(6), 1–12 (2023).10.1140/epjp/s13360-023-04136-5
49. Morelli, D.T. & Slack, G.A. High lattice thermal conductivity solids. In High thermal conductivity materials (pp. 37–68). Springer (2006).
50. Togo A First-principles phonon calculations with phonopy and phono3py J. Phys. Soc. Jpn. 2023 92 1 012001 10.7566/JPSJ.92.012001
Togo, A. First-principles phonon calculations with phonopy and phono3py. J. Phys. Soc. Jpn. 92(1), 012001 (2023).10.7566/JPSJ.92.012001
51. Fan Q Wei Q Yan H Zhang M Zhang D Zhang J A new potential superhard phase of OsN2 Acta Phys. Polon. A 2014 126 3 740 746 10.12693/APhysPolA.126.740
Fan, Q. et al. A new potential superhard phase of OsN2. Acta Phys. Polon. A 126(3), 740–746 (2014).10.12693/APhysPolA.126.740
52. Kleinman L Deformation potentials in silicon. I. Uniaxial strain Phys. Rev. 1962 128 6 2614 10.1103/PhysRev.128.2614
Kleinman, L. Deformation potentials in silicon. I. Uniaxial strain. Phys. Rev. 128(6), 2614 (1962).10.1103/PhysRev.128.2614
53. Kim W Strategies for engineering phonon transport in thermoelectrics J. Mater. Chem. C 2015 3 40 10336 10348 10.1039/C5TC01670C
Kim, W. Strategies for engineering phonon transport in thermoelectrics. J. Mater. Chem. C 3(40), 10336–10348 (2015).10.1039/C5TC01670C
54. Kattan NA Rouf SA Sfina N mana Al-Anazy M Ullah H Hakamy A Mera A Mahmood Q Amin MA Tuning of band gap by anion variation of double perovskites K2AgInX6 (X= Cl, Br) for solar cells and thermoelectric applications J. Solid State Chem. 2023 319 123820 10.1016/j.jssc.2022.123820
Kattan, N. A. et al. Tuning of band gap by anion variation of double perovskites K2AgInX6 (X= Cl, Br) for solar cells and thermoelectric applications. J. Solid State Chem. 319, 123820 (2023).10.1016/j.jssc.2022.123820
55. Yu CJ Advances in modelling and simulation of halide perovskites for solar cell applications J. Phys. Energy 2019 1 2 022001 10.1088/2515-7655/aaf143
Yu, C. J. Advances in modelling and simulation of halide perovskites for solar cell applications. J. Phys. Energy 1(2), 022001 (2019).10.1088/2515-7655/aaf143
