Open Access Article
Aleksandra
Khanina
*,
Tatyana
Sviridova
,
Alexandra
Ivanova
,
Andrey
Voronin
and
Vladimir
Khovaylo
*
Academic Research Center for Energy Efficiency, National University of Science and Technology MISIS, Leninsky Av. 4, Moscow, 119049, Russia. E-mail: khanina.as@misis.ru; khovaylo@misis.ru
First published on 29th October 2025
We report a novel rapid synthesis method for Ba-doped BiCuSeO using mechanochemical synthesis. The phase formation mechanism during high-energy milling as well as the combined effects of the powder-to-ball mass ratio and the heterovalent substitution of bismuth by barium on the thermoelectric properties was investigated. It was shown that the single-phase oxyselenide compounds are formed after 25 minutes of ball milling. The thermoelectric properties of these samples were found to be comparable with the properties of oxyselenides prepared via conventional solid-state reaction reported in the literature.
From this perspective, oxide-based thermoelectric materials, such as zinc oxides,13,14 perovskites,15 and oxychalcogenides,16 offer a significant advantage. Among bismuth-based oxyselenides, BiCuSeO is one of the most promising candidates. BiCuSeO-based compounds have been widely recognized for their exceptionally low thermal conductivity (∼0.5–0.9 W m−1 K−1 at 300 K)17–21, attributed to their layered structure, the presence of heavy elements, and weak interlayer bonding.22 However, despite their low thermal conductivity, these materials typically exhibit low charge carrier concentration and mobility. In this regard, the primary approaches towards enhancing thermoelectric properties involve lattice engineering and defect introduction to reduce thermal conductivity, texturing the material to control carrier scattering and optimizing synthesis methods and conditions.17,23–31 Notably, the intrinsic charge carrier concentration of undoped BiCuSeO is relatively low, approximately 1018 cm−3.19,21,32–37 Therefore, the most common strategy to improve electrical transport is doping at either Bi or Cu sites to increase carrier concentration. Several studies have reported partial substitution with alkaline earth metals such as Mg, Ca, Sr, and Ba,19,22,38,39 as well as other elements, including Ag, Na, La, Zn, Cd, and S,29,30,34,40,41,41–44 which have been shown to enhance the thermoelectric performance of oxyselenides in the mid-temperature range. Additionally, Pb doping, often in combination with other elements, has been used to optimize thermoelectric properties.45–48 It is worth noting that the highest reported zT value for BiCuSeO is 1.6, achieved primarily through shear exfoliation and Pb doping.49 However, the toxicity of lead has prompted scientists to focus on more environmentally friendly elements such as barium.19,50–55 Nevertheless, as mentioned earlier, from an industrial perspective, the thermoelectric figure of merit alone is not the sole determining factor for practical application.
Oxyselenides with composition Bi1−xBaxCuSeO (x = 0, 0.075, 0.125 and 0.175) were obtained by mechanochemical synthesis. To facilitate large-scale production of oxyselenide, this study utilizes cost-effective industrial precursors, specifically using copper oxide instead of bismuth oxide. Moreover, as previously reported, employing CuO is thermodynamically more favorable.79 Accordingly, the commercial high-purity powders of Bi, BaO, Cu, CuO and Se were used as the precursors. A stoichiometric mixture of the raw powders (20 g) was subjected to high-energy ball milling in an air atmosphere using a high-speed planetary ball mill (Activator 2S, Activator, Russia). All samples were milled in 0.25 l steel grinding jars with ∅ 5 mm diameter steel grinding balls. The rotation speed of the jars was 694 rpm, and the grinding time was up to 25 minutes. As previously reported, the phase composition of powders and consequently their transport properties vary with different milling rates or powder-to-ball mass ratios due to differences in the mechanical energy imparted during MS.80 In order to reveal the influence of mechanical energy applied during milling on the phase formation and resultant thermoelectric performance of BiCuSeO, the powder-to-ball mass ratio was varied from 1
:
10 to 1
:
20 and 1
:
40. The resulting powders were consolidated using the pressureless sintering (PLS) method (Fig. 1). This process involved cold pressing (10 t, no holding) into cylindrical tablets (∅ 10.15 mm, 1.5 mm in height), followed by annealing in a sealed quartz ampoule under a 40 kPa argon atmosphere. The annealing was performed at 973 K for 6 hours, with a heating and cooling rate of 1 K min−1 down to 473 K, after which the samples were allowed to cool naturally to room temperature. Phase composition analysis was carried out through powder X-ray diffraction (PXRD) using a DRON-4 diffractometer (IC Bourevestnik, Russia) with Co-Kα radiation (λ = 1.79021 Å) and a Difray-401k desktop diffractometer (SC Scientific Instruments, Russia) with Cr-Kα radiation (λ = 2.2909 Å). The PXRD patterns were refined by the Rietveld method with the self-developed software package.81 The uncertainty in the lattice parameters is 0.004 Å for parameter a and 0.009 Å for c, the uncertainty in volume fractions of phases is 5–10%. The fracture morphologies, distribution of the elements and actual chemical composition of all the samples were analyzed by scanning electron microscopy (SEM; Vega 3 SB, Tescan, Czech Republic) equipped with an energy dispersive X-ray spectrometer (EDS; x-act, Oxford Instruments, UK). Additionally, the average particle size of the powders was determined by dynamic light scattering (DLS; Zetasizer Nano ZS, Malvern Panalytical, UK). Consolidated specimens were cut into discs with ∅ 10 and 1.5 mm thickness for the thermal diffusivity measurements and bars with 10 mm length, 1.5 mm height and 2 mm width for the electrical transport measurements. The Seebeck coefficient α and the electrical conductivity σ were simultaneously measured at 300 K by using a homemade Seebeck coefficient/electrical resistivity measuring system (Cryotel Ltd, Russia). The thermal diffusivity χ was measured by using the laser flash diffusivity method (LFA 457, Netzsch, Germany). The thermal conductivity κtot was calculated using the equation κtot = χCpd where Cp represents the specific heat capacity, which was calculated using the Debye model. The density d was determined using the Archimedes method. The uncertainty in σ and α was 7%, along with 8% in κtot. The overall uncertainty in zT was estimated to be less than 20%.
![]() | ||
| Fig. 1 Schematic representation of the synthesis stages, corresponding reaction products, and consolidation process. | ||
:
40, and an air atmosphere. Fig. 2 shows the PXRD patterns of BiCuSeO samples milled for 15, 20, and 25 minutes. For reference, the manually mixed stoichiometric precursor mixture is designated as 0 minutes. As reported in the literature,73,75,77,79 the formation of the BiCuSeO phase during mechanical milling involves the formation and subsequent reaction of intermediate phases, including bismuth selenides and copper oxide. After 15 minutes of milling, the powder consists of an intermediate phase BiSe1−x (45 wt%), the target BiCuSeO phase (40 wt%) and unreacted CuO (15 wt%). Extending the milling time to 20 minutes significantly increases the content of the BiCuSeO phase to 85 wt%, and appearance of an intermediate phase, Cu2Se, whose content is 15 wt%. After 25 minutes of milling, the phase formation is complete, and reflections only from the BiCuSeO phase are observed on the PXRD patterns (Fig. 2). All major reflections can be attributed to the tetragonal BiCuSeO phase (PDF#01-076-6689) with a ZrCuSiAs structural type and a P4/mmm space group.
Based on the previously described results, a second series of experiments was conducted for 25 min in air to study the effect of the powder-to-ball mass ratio on the formation of the Ba-doped oxyselenide phase. Fig. 3 shows the PXRD patterns of the Bi1−xBaxCuSeO (x = 0, 0.075, 0.125, 0.175) samples. The undoped sample, milled at a powder-to-ball mass ratio of 1
:
10, retains significant amounts of unreacted precursors and intermediate phases, specifically BiSe1−x (35 wt%) and CuO (15 wt%). This indicates that the relatively mild milling conditions do not provide sufficient energy to complete the reaction. The results suggest that at this ratio, the amount of secondary phases decreases with increasing Ba concentration. For the sample with x = 0.075, only a small amount of unreacted CuO (5 wt%) was detected, while no secondary phases were observed in samples with x > 0.075. Similarly, when milling at a ratio of 1
:
40, the undoped sample still contains CuO (5 wt%). Based on these findings, a powder-to-ball mass ratio of 1
:
20 appears to be optimal for the BiCuSeO synthesis. After the PLS process, phase composition changes were observed in almost all samples, with the formation of secondary phases such as BiSe1−x, Bi2O3, and Cu2O (Fig. S1). Notably, in powders synthesized at a 1
:
10 powder-to-ball mass ratio, the BiCuSeO phase formation was completed during PLS. This suggests that to achieve a predominantly BiCuSeO phase, a minimum milling time of 25 minutes at a 1
:
10 ratio and a rotation speed of 694 rpm is required, followed by reaction sintering during the PLS process to ensure complete phase formation.
The refined lattice parameters of BiCuSeO are in good agreement with previously reported values for samples synthesized via SSR34 and MS79 (see Table S1 and Fig. 4a). When oxyselenide samples are prepared by MS, point defects such as copper vacancies (VCu) may form.77 It is important to note that copper vacancies can have a significant impact on the transport properties of oxyselenide.37,82,83 In this case, two competing factors can influence the evolution of lattice parameters in doped Bi1−xBaxCuSeO (x = 0, 0.075, 0.125, 0.175) oxyselenides. Copper deficiency tends to decrease the lattice parameters a and c,17 while doping with Ba2+ (1.42 Å), which has a larger ionic radius than bismuth Bi3+ (1.17 Å),55 should lead to an increase in the lattice parameters.19,51 Our experimental data indicate that doping is the dominant factor, resulting in a gradual linear increase in both lattice parameters a and c (Fig. 4 and Table S1). However, a slight decrease in the lattice parameters is observed with an increase in the powder-to-ball mass ratio, which is likely due to the formation of VCu.
![]() | ||
| Fig. 4 (a) Lattice parameters a and c (comparative data from Wu et al.:79 dashed line – parameter a and dotted line – parameter c), (b) crystallite size dCS and (c) microdeformation (md) as a function of nominal BaBi concentration. | ||
The Bi
:
Cu
:
Se elemental ratio for all compositions deviates slightly from the nominal 1
:
1
:
1 ratio, remaining within the error limits of the EDS method. In most cases, an excess of Bi and a deficiency of Cu are observed (Table S1). However, we propose that this non-stoichiometry is intrinsically linked to the reaction pathway, which proceeds through the formation of copper-deficient intermediate phases such as BiSe1−x and Cu2Se, as identified by XRD (Fig. 2). The kinetics of this diffusion-controlled solid-state process84 appear to favor the formation and persistence of these defective intermediates, ultimately leading to the incorporation of copper vacancies (VCu) into the final BiCuSeO lattice. The measured Ba/(Bi + Ba) ratio closely matches the nominal values, confirming the intended doping levels. The EDS data are also consistent with the PXRD results, confirming the presence of CuO and BiSe1−x inclusions in the samples milled at a powder-to-ball mass ratio of 1
:
10 (Fig. 5 and S3). The values of the crystallite size (dCS) systematically decrease with increasing powder-to-ball mass ratio, which is explained by the higher milling energies (Fig. 4b and Table S1). However, Ba substitution does not introduce significant changes, and the CS values remain in a similar range without exhibiting a clear monotonic trend. The impact of the powder-to-ball mass ratio on reaction kinetics can be explained in terms of stress-induced electron redistribution at interfaces of dissimilar solids. Our precursor mixture is a complex interface of metals, a semiconductor, and insulators. During high-energy ball milling, these interfaces experience extreme mechanical stress, which is quantitatively evidenced by the microdeformations (md). As shown in Fig. 4c, the systematic increase in microdeformation values with higher barium content is attributed to the larger ionic radius of Ba2+, which induces significant local lattice distortion and amplifies the internal stress fields. According to the model proposed by Sun et al.,85 the stress does not merely create fresh surfaces and defects; it directly modifies the electronic states at the interfaces, leading to significant charge redistribution (e.g., from metallic Bi/Cu to CuO). The powder-to-ball mass ratio is a direct control parameter for this process: a higher ratio increases the frequency and force of these triboelectrically active collisions, thereby amplifying the cumulative charge redistribution effect. At the same time, the high microstrain values for the 1
:
10 ratio, particularly in the doped samples, are interpreted as a consequence of incomplete reaction and the persistence of hard, unreacted precursor phases, which create localized stress and hinder homogeneous plastic deformation. Thus, increased mechanical energy input efficiently drives the diffusion-controlled synthesis of complex oxyselenides, far surpassing what thermal energy alone could achieve under similar conditions.
![]() | ||
Fig. 5 SEM micrograph of the polished surface of the BiCuSeO samples obtained with powder-to-ball mass ratios of (a) 1 : 10, (b) 1 : 20, and (c) 1 : 40. | ||
The microstructure of oxyselenides strongly depends on the synthesis method. MS produces fine powders with average particle sizes of approximately 189, 143, and 89 nm for the powders milled at 1
:
10, 1
:
20, and 1
:
40 ratios, respectively (Fig. S2). These values are significantly smaller than those observed in samples synthesized via conventional SSR methods.19,51,54,60–63,65,86 Moreover, MS promotes particle dispersion and introduces numerous defects and microstrains, leading to an increased grain boundary area. This, in turn, enhances phonon scattering and effectively reduces thermal conductivity. After 6 hours of PLS, the grain and agglomerate sizes consistently increase to 5–10 μm for all ratios. However, the density of the resulting volumetric images remains relatively low, ranging from 79% to 88%. Thus, all samples exhibit fine grains with a lath morphology, randomly oriented, which results from the high polycrystallinity of the material (Fig. S4).
For all samples, the Seebeck coefficient values are positive, indicating p-type conductivity. The absolute values range from 80 to 250 μV K−1 at room temperature (Fig. 6a and b). An increase in the BaBi concentration and the powder-to-ball mass ratio leads to a decrease in α at 300 K, likely due to a higher charge carrier concentration.1 However, for undoped compositions, the α values are slightly lower than those reported in the literature,19,51,52,75,87 which may be attributed to a higher charge carrier concentration resulting from copper vacancies formed during MS.
![]() | ||
| Fig. 6 Concentration and powder-to-ball mass ratio dependence of the (a and b) Seebeck coefficient α and (c and d) electrical conductivity σ for Bi1−xBaxCuSeO (x = 0, 0.075, 0.125, 0.175) samples consolidated by PLS. Literature data for other Bi1−xBaxCuSeO-based oxyselenides are also shown for comparison (Feng et al.,75 Li et al.,19 Park et al.,51 Pei et al.,52 Jiang et al.,88 Kim et al.,61 Lei et al.,89 Liang et al.,64 Liu et al.,90 Pan et al.,68 Park et al.,51 Wen et al.,53 Yang et al.,91 Zeng et al.,23 Li et al.,86 and Yang et al.54). All displayed data points correspond to values obtained at 300 K. | ||
The electrical conductivity σ exhibits an inverse dependence on concentration across all compositions, gradually increasing with the increasing concentration of barium substituting for bismuth (Fig. 6c and d). Notably, for each powder-to-ball mass ratio, the undoped samples demonstrate higher σ values than those reported in previous studies on oxyselenides synthesized via MS, self-propagating high-temperature synthesis, sol–gel, microwave, and SSR.19,23,51–53,61,64,68,75,88–91 This observation suggests a high carrier concentration, likely due to VCu. Increasing the powder-to-ball mass ratio to 1
:
20 naturally enhances σ, whereas a further increase to 1
:
40 results in a slight decrease in the absolute σ values, presumably due to an increased defect density and, consequently, enhanced charge carrier scattering. However, a different trend is observed for the sample with x = 0.075, where the highest σ value corresponds to the composition milled at a 1
:
40 ratio. For the remaining samples, the highest σ values are achieved at a powder-to-ball mass ratio of 1
:
20. The highest absolute σ values, approximately 350 Ω−1 cm−1 at 300 K, are obtained for the sample with the highest dopant concentration, Bi0.825Ba0.175CuSeO, milled at a 1
:
20 ratio. However, the absolute conductivity values for the Ba-doped samples remain lower than those reported for as-synthesized SSR composites in the literature,51,53,54,86 likely due to the MS process, which introduces numerous additional scattering centers, including grain boundaries, dislocations, and point defects.
The total thermal conductivity κtot for all samples is approximately 0.7–0.9 W m−1 K−1 at 300 K, which is in good agreement with Feng et al.'s results for MS oxyselenides (Fig. S5a and b).75 Thermal measurements indicate that barium doping leads to an increase in κtot, which can be understood by considering electron (κel) and lattice (κlat) components of κtot. The electron thermal conductivity κel is proportional to the electrical conductivity, σ, following the Wiedemann–Franz law: κel = σLT. Here, the Lorentz number, L, was temperature-dependent and estimated using experimental values of the Seebeck coefficient in the framework of the effective mass model. The primary reason for the observed increase in κtot is the rise in κel, which directly correlates with σ due to a higher charge carrier concentration (Fig. S5e and f). Comparison of the lattice component of thermal conductivity with literature data suggests that structural defects, including pores—whose formation was promoted by the MS process—introduce additional phonon scattering surfaces, thereby reducing their mean free path (Fig. S5c and d). Additionally, the low density of bulk samples contributes to a decrease in the speed of sound within the material, which is proportional to κlat.
Based on the measured electrical and thermal transport characteristics, the thermoelectric figure of merit, zT, was calculated and is presented as a function of temperature in Fig. 6g and h. Although the zT values exhibit only minor variations within the error range for different powder-to-ball mass ratios, it is noteworthy that the highest thermoelectric figure of merit for all compositions was achieved for the samples milled at a ratio of 1
:
10. While the obtained zT values are not the highest reported for Ba-doped oxyselenides, they fall within the range of 0.07 to 0.13 at T = 300 K, consistent with previously published results. Moreover, the zT values obtained for undoped compositions are higher (about 0.1 at 300 K) than those reported for finely dispersed BiCuSeO powders prepared via sol–gel or ball milling methods,23,75,90,92 as well as for BiCuSeO synthesized by more established techniques such as self-propagating high-temperature synthesis, microwave-assisted synthesis, and SSR.51–53,61,64,86,88,89,91 This confirms that mechanochemical synthesis, requiring only 25 minutes, is a promising approach for the large-scale production of Bi1−xBaxCuSeO-based oxyselenides.
| This journal is © The Royal Society of Chemistry 2026 |