The role of copper in the thermal conductivity of thermoelectric oxychalcogenides: do lone pairs matter?

Understanding the underlying mechanisms that suppress thermal conduction in solids is of paramount importance for the targeted design of materials for thermal management and thermoelectric energy conversion applications. Bismuth copper oxychalcogenides, BiOCuQ (Q = Se, Te), are highly crystalline thermoelectric materials with an unusually low lattice thermal conductivity of approx. 0.5 Wm-1K-1, a value normally found in amorphous materials. Here we unveil the origin of the unusual thermal transport properties of these phases. First principles calculations of the vibrational properties combined with analysis of in-situ neutron diffraction data, demonstrate that weak bonding of copper atoms within the structure leads to an unexpected vibrational mode at low frequencies, which is likely to be a major contributor to the low thermal conductivity of these materials. In addition, we show that anharmonicity and the large Grüneisen parameter in these oxychalcogenides are mainly related to the low frequency copper vibrations, rather than to the Bi3+ lone pairs.

Synthesis, characterisation and thermoelectric properties of the oxytelluride Bi2O2Te

Bi2O2Te was synthesised from a stoichiometric mixture of Bi, Bi2O3 and Te by a solid state reaction. Analysis of powder X-ray diffraction data indicates that this material crystallises in the anti-ThCr2Si2 structure type (space group I4/mmm), with lattice parameters a = 3.98025(4) and c = 12.70391(16) Å. The electrical and thermal transport properties of Bi2O2Te were investigated as a function of temperature over the temperature range 300 ≤ T/K ≤ 665. These measurements indicate that Bi2O2Te is an n-type semiconductor, with a band gap of 0.23 eV. The thermal conductivity of Bi2O2Te is remarkably low for a crystalline material, with a value of only 0.91 W m-1 K-1 at room temperature.

The interplay of metal-atom ordering, Fermi leveltuning and thermoelectric properties in cobalt shandites Co3M2S2 (M = Sn, In)

A combination of structural, physical and computational techniques including powder X-ray and neutron diffraction, SQUID magnetometry, electrical and thermal transport measurements, DFT calculations and 119Sn Mössbauer and X-ray photoelec-tron spectroscopies has been applied to Co3Sn2-xInxS2 (0 ≤ x ≤ 2) in an effort to understand the relationship between metal-atom ordering and physical properties as the Fermi level is systematically varied. Whilst solid solution behavior is found throughout the composition region, powder neutron diffraction reveals that indium preferentially occupies an inter-layer site over an alternative kagome-like intra-layer site. DFT calculations indicate that this ordering, which leads to a lowering of energy, is related to the dif-fering bonding properties of tin and indium. Spectroscopic data suggest that throughout the composition range 0 ≤ x ≤ 2, all ele-ments adopt oxidation states that are significantly reduced from expectations based on formal charges. Chemical substitution ena-bles the electrical transport properties to be controlled through tuning of the Fermi level within a region of the density of states, which comprises narrow bands of predominantly Co d-character. This leads to a compositionally-induced double metal-to-semiconductor-to-metal transition. The marked increase in the Seebeck coefficient as the semiconducting region is approached leads to a substantial improvement in the thermoelectric figure of merit, ZT, which exhibits a maximum of ZT = 0.32 at 673 K. At 425 K, the figure of merit for phases in the region 0.8 ≤ x ≤ 0.85 is amongst the highest reported for sulphide phases, suggesting these materials may have applications in low-grade waste heat recovery.

Ball milling as an effective route for the preparation of doped bornite: synthesis, stability and thermoelectric properties

Bornite, Cu5FeS4, is a naturally-occuring mineral with an ultralow thermal conductivity and potential for thermoelectric power generation. We describe here a new, easy and scalable route to synthesise bornite, together with the thermoelectric behaviour of manganese-substituted derivatives, Cu5Fe1-xMnxS4 (0 ≤ x ≤ 0.10). The electrical and thermal transport properties of Cu5Fe1-xMnxS4 (0 ≤ x ≤ 0.10), which are p-type semiconductors, were measured from room temperature to 573 K. The stability of bornite was investigated by thermogravimetric analysis under inert and oxidising atmospheres. Repeated measurements of the electrical transport properties confirm that bornite is stable up to 580 K under an inert atmosphere, while heating to 890 K results in rapid degradation. Ball milling leads to a substantial improvement in the thermoelectric figure of merit of unsusbtituted bornite (ZT = 0.55 at 543 K), when compared to bornite prepared by conventional high-temperature synthesis (ZT < 0.3 at 543 K). Manganese-substituted samples have a ZT comparable to that of unsubstituted bornite.

Up-scaled synthesis process of sulphur-based thermoelectric materials

The scale up of Spark Plasma Sintering (SPS) for the consolidation of large square monoliths (50 × 50 × 3 mm3) of thermoelectric material is demonstrated and the properties of the fabricated samples compared with those from laboratory scale SPS. The SPS processing of n-type TiS2 and p-type Cu10.4Ni1.6Sb4S13 produces highly dense compacts of phase pure material. Electrical and thermal transport property measurements reveal that the thermoelectric performance of the consolidated n- and p-type materials is comparable with that of material processed using laboratory scale SPS, with ZT values that approach 0.8 and 0.35 at 700 K for Cu10.4Ni1.6Sb4S13 and TiS2, respectively. Mechanical properties of the consolidated materials shows that large-scale SPS processing produces highly homogeneous materials with hardness and elastic moduli that deviate little from values obtained on materials processed on the laboratory scale.

The impact of charge transfer and structural disorder on the thermoelectric properties of cobalt intercalated TiS2

A family of phases, CoxTiS2 (0 ≤ x ≤ 0.75) has been prepared and characterised by powder X-ray and neutron diffraction, electrical and thermal transport property measurements, thermal analysis and SQUID magnetometry. With increasing cobalt content, the structure evolves from a disordered arrangement of cobalt ions in octahedral sites located in the van der Waals’ gap (x ≤ 0.2), through three different ordered vacancy phases, to a second disordered phase at x ≥ 0.67. Powder neutron diffraction reveals that both octahedral and tetrahedral inter-layer sites are occupied in Co0.67TiS2. Charge transfer from the cobalt guest to the TiS2 host affords a systematic tuning of the electrical and thermal transport properties. At low levels of cobalt intercalation (x < 0.1), the charge transfer increases the electrical conductivity sufficiently to offset the concomitant reduction in |S|. This, together with a reduction in the overall thermal conductivity leads to thermoelectric figures of merit that are 25 % higher than that of TiS2, ZT reaching 0.30 at 573 K for CoxTiS2 with 0.04 ≤ x ≤ 0.08. Whilst the electrical conductivity is further increased at higher cobalt contents, the reduction in |S| is more marked due to the higher charge carrier concentration. Furthermore both the charge carrier and lattice contributions to the thermal conductivity are increased in the electrically conductive ordered-vacancy phases, with the result that the thermoelectric performance is significantly degraded. These results illustrate the competition between the effects of charge transfer from guest to host and the disorder generated when cobalt cations are incorporated in the inter-layer space.

Kinetic multi-layer model of aerosol surface and bulk chemistry (KM-SUB): the influence of interfacial transport and bulk diffusion on the oxidation of oleic acid by ozone

We present a novel kinetic multi-layer model that explicitly resolves mass transport and chemical reaction at the surface and in the bulk of aerosol particles (KM-SUB). The model is based on the PRA framework of gas-particle interactions (Poschl-Rudich-Ammann, 2007), and it includes reversible adsorption, surface reactions and surface-bulk exchange as well as bulk diffusion and reaction. Unlike earlier models, KM-SUB does not require simplifying assumptions about steady-state conditions and radial mixing. The temporal evolution and concentration profiles of volatile and non-volatile species at the gas-particle interface and in the particle bulk can be modeled along with surface concentrations and gas uptake coefficients. In this study we explore and exemplify the effects of bulk diffusion on the rate of reactive gas uptake for a simple reference system, the ozonolysis of oleic acid particles, in comparison to experimental data and earlier model studies. We demonstrate how KM-SUB can be used to interpret and analyze experimental data from laboratory studies, and how the results can be extrapolated to atmospheric conditions. In particular, we show how interfacial and bulk transport, i.e., surface accommodation, bulk accommodation and bulk diffusion, influence the kinetics of the chemical reaction. Sensitivity studies suggest that in fine air particulate matter oleic acid and compounds with similar reactivity against ozone (carbon-carbon double bonds) can reach chemical lifetimes of many hours only if they are embedded in a (semi-)solid matrix with very low diffusion coefficients (< 10(-10) cm(2) s(-1)). Depending on the complexity of the investigated system, unlimited numbers of volatile and non-volatile species and chemical reactions can be flexibly added and treated with KM-SUB. We propose and intend to pursue the application of KM-SUB as a basis for the development of a detailed master mechanism of aerosol chemistry as well as for the derivation of simplified but realistic parameterizations for large-scale atmospheric and climate models.

Kinetic multi-layer model of aerosol surface and bulk chemistry (KM-SUB): the influence of interfacial transport and bulk diffusion on the oxidation of oleic acid by ozone

We present a novel kinetic multi-layer model that explicitly resolves mass transport and chemical reaction at the surface and in the bulk of aerosol particles (KM-SUB). The model is based on the PRA framework of gas–particle interactions (P¨oschl et al., 5 2007), and it includes reversible adsorption, surface reactions and surface-bulk exchange as well as bulk diffusion and reaction. Unlike earlier models, KM-SUB does not require simplifying assumptions about steady-state conditions and radial mixing. The temporal evolution and concentration profiles of volatile and non-volatile species at the gas-particle interface and in the particle bulk can be modeled along with surface 10 concentrations and gas uptake coefficients. In this study we explore and exemplify the effects of bulk diffusion on the rate of reactive gas uptake for a simple reference system, the ozonolysis of oleic acid particles, in comparison to experimental data and earlier model studies. We demonstrate how KM-SUB can be used to interpret and analyze experimental data from laboratory stud15 ies, and how the results can be extrapolated to atmospheric conditions. In particular, we show how interfacial transport and bulk transport, i.e., surface accommodation, bulk accommodation and bulk diffusion, influence the kinetics of the chemical reaction. Sensitivity studies suggest that in fine air particulate matter oleic acid and compounds with similar reactivity against ozone (C=C double bonds) can reach chemical lifetimes of 20 multiple hours only if they are embedded in a (semi-)solid matrix with very low diffusion coefficients (10−10 cm2 s−1). Depending on the complexity of the investigated system, unlimited numbers of volatile and non-volatile species and chemical reactions can be flexibly added and treated with KM-SUB. We propose and intend to pursue the application of KM-SUB 25 as a basis for the development of a detailed master mechanism of aerosol chemistry as well as for the derivation of simplified but realistic parameterizations for large-scale atmospheric and climate models.

The impact of oceanic heat transport on the atmospheric circulation

A general circulation model of intermediate complexity with an idealized Earth-like aquaplanet setup is used to study the impact of changes in the oceanic heat transport on the global atmospheric circulation. Focus is on the atmospheric mean meridional circulation and global thermodynamic properties. The atmosphere counterbalances to a large extent the imposed changes in the oceanic heat transport, but, nonetheless, significant modifications to the atmospheric general circulation are found. Increasing the strength of the oceanic heat transport up to 2.5 PW leads to an increase in the global mean near-surface temperature and to a decrease in its equator-to-pole gradient. For stronger transports, the gradient is reduced further, but the global mean remains approximately constant. This is linked to a cooling and a reversal of the temperature gradient in the tropics. Additionally, a stronger oceanic heat transport leads to a decline in the intensity and a poleward shift of the maxima of both the Hadley and Ferrel cells. Changes in zonal mean diabatic heating and friction impact the properties of the Hadley cell, while the behavior of the Ferrel cell is mostly controlled by friction. The efficiency of the climate machine, the intensity of the Lorenz energy cycle and the material entropy production of the system decline with increased oceanic heat transport. This suggests that the climate system becomes less efficient and turns into a state of reduced entropy production as the enhanced oceanic transport performs a stronger large-scale mixing between geophysical fluids with different temperatures, thus reducing the available energy in the climate system and bringing it closer to a state of thermal equilibrium.