Preparation and characterisation of certain II-VI, I-III-VI2 semiconductor thin films and transparent conducting oxides

There is an increasing demand for renewable energies due to the limited availability of fossil and nuclear fuels and due to growing environmental problems. Photovoltaic (PV) energy conversion has the potential to contribute significantly to the electrical energy generation in the future. Currently, the cost for photovoltaic systems is one of the main obstacles preventing production and application on a large scale. The photovoltaic research is now focused on the development of materials that will allow mass production without compromising on the conversion efficiencies. Among important selection criteria of PV material and in particular for thin films, are a suitable band gap, high absorption coefficient and reproducible deposition processes capable of large-volume and low cost production. The chalcopyrite semiconductor thin films such as Copper indium selenide and Copper indium sulphide are the materials that are being intensively investigated for lowering the cost of solar cells. Conversion efficiencies of 19 % have been reported for laboratory scale solar cell based on CuInSe2 and its alloys. The main objective of this thesis work is to optimise the growth conditions of materials suitable for the fabrication of solar cell, employing cost effective techniques. A typical heterojunction thin film solar cell consists of an absorber layer, buffer layer and transparent conducting contacts. The most appropriate techniques have been used for depositing these different layers, viz; chemical bath deposition for the window layer, flash evaporation and two-stage process for the absorber layer, and RF magnetron sputtering for the transparent conducting layer. Low cost experimental setups were fabricated for selenisation and sulphurisation experiments, and the magnetron gun for the RF sputtering was indigenously fabricated. The films thus grown were characterised using different tools. A powder X-ray diffractometer was used to analyse the crystalline nature of the films. The energy dispersive X-ray analysis (EDX) and scanning electron microscopy i (SEM) were used for evaluating the composition and morphology of the films. Optical properties were investigated using the UV-Vis-NIR spectrophotometer by recording the transmission/absorption spectra. The electrical properties were studied using the two probe and four probe electrical measurements. Nature of conductivity of the films was determined by thermoprobe and thermopower measurements. The deposition conditions and the process parameters were optimised based on these characterisations.

Photoconductivity studies on some semiconducting thin films for photovoltaic applications

Photoconductivity (PC) processes may be the most suitable technique for obtaining information about the states in the gap. It finds applications in photovoItaics, photo detection and radiation measurements. The main task in the area of photovoltaics, is to increase the efficiency of the device and also to develop new materials with good optoelectronic properties useful for energy conversion, keeping the idea of cost effectiveness. Photoconduction includes generation and recombination of carriers and their transport to the electrodes. So thermal relaxation process, charge carrier statistics, effects of electrodes and several mechanisms of recombination are involved in photoconductivity.A major effect of trapping is to make the experimentally observed decay time of photocurrent, longer than carrier lifetime. If no trapping centers are present, then observed photocurrent will decay in the same way as the density of free carriers and the observed decay time will be equal to carrier lifetime. If the density of free carriers is much less than density of trapped carriers, the entire decay of photocurrent is effectively dominated by the rate of trap emptying rather than by the rate of recombination.In the present study, the decay time of carriers was measured using photoconductive decay (PCD) technique. For the measurements, the film was loaded in a liquid Helium cryostat and the temperature was controlled using Lakshore Auto tuning temperature controller (Model 321). White light was used to illuminate the required area of the sample. Heat radiation from the light source was avoided by passing the light beam through a water filter. The decay current. after switching off the illumination. was measured using a Kiethely 2000 multi meter. Sets of PCD measurements were taken varying sample temperature, sample preparation temperature, thickness of the film, partial pressure of Oxygen and concentration of a particular element in a compound. Decay times were calculated using the rate window technique, which is a decay sampling technique particularly suited to computerized analysis. For PCD curves with two well-defined regions, two windows were chosen, one at the fast decay region and the other at the slow decay region. The curves in a particular window were exponentially fitted using Microsoft Excel 2000 programme. These decay times were plotted against sample temperature and sample preparation temperature to study the effect of various defects in the film. These studies were done in order to optimize conditions of preparation technique so as to get good photosensitive samples. useful for photovoltaic applications.Materials selected for the study were CdS, In2Se3, CuIn2Se3 and CuInS2• Photoconductivity studies done on these samples are organised in six chapters including introduction and conclusion.

Photothermal Investigations on Certain Plasma Polymerized Thin Films and Near IR Overtone Studies of Some Organic Molecules

Photothermal effect refers to heating of a sample due to the absorption of electromagnetic radiation. Photothermal (PT) heat generation which is an example of energy conversion has in general three kinds of applications. 1. PT material probing 2. PT material processing and 3. PT material destruction. The temperatures involved increases from 1-. 3. Of the above three, PT material probing is the most important in making significant contribution to the field of science and technology. Photothermal material characterization relies on high sensitivity detection techniques to monitor the effects caused by PT material heating of a sample. Photothermal method is a powerful high sensitivity non-contact tool used for non-destructive thermal characterization of materials. The high sensitivity of the photothermal methods has led to its application for analysis of low absorbance samples. Laser calorimetry, photothermal radiometry, pyroelectric technique, photoacoustic technique, photothermal beam deflection technique, etc. come under the broad class ofphotothermal techniques. However the choice of a suitable technique depends upon the nature of the sample, purpose of measurement, nature of light source used, etc. The present investigations are done on polymer thin films employing photothermal beam deflection technique, for the successful determination of their thermal diffusivity. Here the sample is excited by a He-Ne laser (A = 6328...\ ) which acts as the pump beam. Due to the refractive index gradient established in the sample surface and in the adjacent coupling medium, another optical beam called probe beam (diode laser, A= 6500A ) when passed through this region experiences a deflection and is detected using a position sensitive detector and its output is fed to a lock-in amplifier from which the amplitude and phase of the deflection can be directly obtained. The amplitude and phase of the signal is suitably analysed for determining the thermal diffusivity.The production of polymer thin film samples has gained considerable attention for the past few years. Plasma polymerization is an inexpensive tool for fabricating organic thin films. It refers to formation of polymeric materials under the influence of plasma, which is generated by some kind of electric discharge. Here plasma of the monomer vapour is generated by employing radio frequency (MHz) techniques. Plasma polymerization technique results in homogeneous, highly adhesive, thermally stable, pinhole free, dielectric, highly branched and cross-linked polymer films. The possible linkage in the formation of the polymers is suggested by comparing the FTIR spectra of the monomer and the polymer.Near IR overtone investigations on some organic molecules using local mode model are also done. Higher vibrational overtones often provide spectral simplification and greater resolution of peaks corresponding to nonequivalent X-H bonds where X is typically C, N or O. Vibrational overtone spectroscopy of molecules containing X-H oscillators is now a well established tool for molecular investigations. Conformational and steric differences between bonds and structural inequivalence ofCH bonds (methyl, aryl, acetylenic, etc.) are resolvable in the higher overtone spectra. The local mode model in which the X-H oscillators are considered to be loosely coupled anharmonic oscillators has been widely used for the interpretation of overtone spectra. If we are exciting a single local oscillator from the vibrational ground state to the vibrational state v, then the transition energy of the local mode overtone is given by .:lE a......v = A v + B v2 • A plot of .:lE / v versus v will yield A, the local mode frequency as the intercept and B, the local mode diagonal anharmonicity as the slope. Here A - B gives the mechanical frequency XI of the oscillator and B = X2 is the anharmonicity of the bond. The local mode parameters XI and X2 vary for non-equivalent X-H bonds and are sensitive to the inter and intra molecular environment of the X-H oscillator.

Mass production of nitrifying bacterial consortia for the rapid establishment of nitrification in saline recirculating aquaculture systems

Two distinct nitrifying bacterial consortia, namely an ammonia oxidizing non-penaeid culture (AMO NPCU-1) and an ammonia oxidizing penaeid culture (AMOPCU-1), have been mass produced in a nitrifying bacterial consortia production unit (NBCPU). The consortia, maintained at 4 C were activated and cultured in a 2 l fermentor initially. At this stage the net biomass (0.105 and 0.112 g/l), maximum specific growth rate (0.112 and 0.105/h) and yield coefficients (1.315 and 2.08) were calculated respectively, for AMONPCU-1 and AMOPCU-1 on attaining stationary growth phase. Subsequently on mass production in a 200 l NBCPU under optimized culture conditions, the total amounts of NH4 ?–N removed by AMONPCU-1 and AMOPCU-1 were 1.948 and 1.242 g/l within 160 and 270 days, respectively. Total alkalinity reduction of 11.7–14.4 and 7.5–9.1 g/l were observed which led to the consumption of 78 and 62 g Na2CO3. The yield coefficient and biomass of AMONPCU-1 were 0.67 and 125.3 g/l and those of AMOPCU-1 were 1.23 and 165 g/l. The higher yield coefficient and growth rate of AMOPCU-1 suggest better energy conversion efficiency and higher CO2 fixation potential. Both of the consortia were dominated by Nitrosomonas-like organisms. The consortia may find application in the establishment of nitrification within marine and brackish water culture systems.

Dieselmotor-Kraft-Wärme-Kopplungsanlagen im Kontext der Integration Erneuerbarer Energien in die Energieversorgung

In dieser Arbeit werden die sich abzeichnenden zukünftigen Möglichkeiten, Stärken und Schwächen der Kraft-Wärme-Kopplung (KWK) untersucht. Dies geschieht vor dem Hintergrund des Klimawandels, der Integration steigender Anteile Erneuerbarer Energien in die Stromerzeugung und unter Berücksichtigung der sich damit ergebenden Herausforderungen, eine sichere und nachhaltige Stromversorgung zu gestalten. Der Fokus liegt auf der Dieselmotor-KWK und der Nutzung nachwachsender Kraftstoffe. Es wird davon ausgegangen, dass der Übergang zu einer reinen Stromerzeugung aus Erneuerbaren Energiequellen in Deutschland unter erheblicher Einbindung des hohen Potentials der kostengünstigen, umweltfreundlichen, aber in der Leistung extrem fluktuierenden Windenergie erfolgen wird. Als dezentrales Integrationswerkzeug wurde die Kraft-Wärme-Kopplung mit Dieselmotoren untersucht. Sie entspricht aufgrund ihrer großen Flexibilität und ihrer hohen Wirkungsgrade mit vergleichsweise kleinen Leistungen sehr gut den Anforderungen der gleichzeitigen dezentralen Wärmenutzung. In der Dissertation werden die Randbedingungen der Dieselmotor-KWK untersucht und beschrieben. Darauf aufbauend werden unterschiedliche Modelle der Windintegration durch KWK erarbeitet und in diversen Variationen wird der Ausgleich der Stromerzeugung aus Windenergie durch KWK simuliert. Darüber hinaus werden dezentrale KWK-Anlagen hinsichtlich eines koordinierten gemeinsamen Betriebs und hinsichtlich der optimalen Auslegung für den Windenergieausgleich betrachtet. Es wird für den beschriebenen Kontext der Erneuerbaren Energien und der Kraft-Wärme-Kopplung das Thema „Umweltwirkungen“ diskutiert. Es wird dargelegt, dass die heute verwendeten Ansätze zur Bewertung der KWK zu einer Verzerrung der Ergebnisse führen. Demgegenüber wurde mit der so genannten Outputmethode eine Methode der Ökobilanzierung vorgestellt, die, im Gegensatz zu den anderen Methoden, keine verzerrenden Annahmen in die Wirkungsabschätzung aufnimmt und somit eine eindeutige und rein wissenschaftliche Auswertung bleibt. Hiermit ist die Grundlage für die Bewertung der unterschiedlichen Technologien und Szenarien sowie für die Einordnung der KWK in den Kontext der Energieerzeugung gegeben. Mit der Outputmethode wird u.a. rechnerisch bewiesen, dass die gekoppelte Strom- und Wärmeerzeugung in KWK-Anlagen tatsächlich die optimale Nutzung der regenerativen Kraftstoffe „Biogas“ und „Pflanzenöl“ im Hinblick auf Ressourceneinsatz, Treibhausgaseinsparung und Exergieerzeugung ist. Es wurde darüber hinaus die Frage untersucht woher die für die Stromerzeugung durch Dieselmotor-KWK-Anlagen notwendige Bioenergie genommen werden kann. Es ist erwiesen, dass die in Deutschland nutzbare landwirtschaftliche Fläche nur zur Deckung eines Teils der Stromerzeugung ausreichen würde. Einheimisches Biogas und nachhaltiges importiertes Pflanzenöl, das in hohem Maße auf degradierten Böden angebaut werden sollte, können die notwendige Brennstoffenergie liefern. Um im Ausland ausreichend Pflanzenöl herstellen zu können, wird eine landwirtschaftliche Fläche von 6 bis 12 Mio. ha benötigt. Das Ergebnis ist, dass ein voller Ausgleich von Windenergie-Restlast durch KWK mit Erneuerbaren Energieträgern sinnvoll und machbar ist! Dieses Wind-KWK-DSM-System sollte durch ein Stromnetz ergänzt sein, das Wasserkraftstrom für den Großteil der Regelenergieaufgaben nutzt, und das den großräumigen Ausgleich Erneuerbarer Energien in Europa und den Nachbarregionen ermöglicht.

On the energetics of stratified turbulent mixing, irreversible thermodynamics, Boussinesq models and the ocean heat engine controversy

In this paper, the available potential energy (APE) framework of Winters et al. (J. Fluid Mech., vol. 289, 1995, p. 115) is extended to the fully compressible Navier– Stokes equations, with the aims of clarifying (i) the nature of the energy conversions taking place in turbulent thermally stratified fluids; and (ii) the role of surface buoyancy fluxes in the Munk & Wunsch (Deep-Sea Res., vol. 45, 1998, p. 1977) constraint on the mechanical energy sources of stirring required to maintain diapycnal mixing in the oceans. The new framework reveals that the observed turbulent rate of increase in the background gravitational potential energy GPEr , commonly thought to occur at the expense of the diffusively dissipated APE, actually occurs at the expense of internal energy, as in the laminar case. The APE dissipated by molecular diffusion, on the other hand, is found to be converted into internal energy (IE), similar to the viscously dissipated kinetic energy KE. Turbulent stirring, therefore, does not introduce a new APE/GPEr mechanical-to-mechanical energy conversion, but simply enhances the existing IE/GPEr conversion rate, in addition to enhancing the viscous dissipation and the entropy production rates. This, in turn, implies that molecular diffusion contributes to the dissipation of the available mechanical energy ME =APE +KE, along with viscous dissipation. This result has important implications for the interpretation of the concepts of mixing efficiency γmixing and flux Richardson number Rf , for which new physically based definitions are proposed and contrasted with previous definitions. The new framework allows for a more rigorous and general re-derivation from the first principles of Munk & Wunsch (1998, hereafter MW98)’s constraint, also valid for a non-Boussinesq ocean: G(KE) ≈ 1 − ξ Rf ξ Rf Wr, forcing = 1 + (1 − ξ )γmixing ξ γmixing Wr, forcing , where G(KE) is the work rate done by the mechanical forcing, Wr, forcing is the rate of loss of GPEr due to high-latitude cooling and ξ is a nonlinearity parameter such that ξ =1 for a linear equation of state (as considered by MW98), but ξ <1 otherwise. The most important result is that G(APE), the work rate done by the surface buoyancy fluxes, must be numerically as large as Wr, forcing and, therefore, as important as the mechanical forcing in stirring and driving the oceans. As a consequence, the overall mixing efficiency of the oceans is likely to be larger than the value γmixing =0.2 presently used, thereby possibly eliminating the apparent shortfall in mechanical stirring energy that results from using γmixing =0.2 in the above formula.

Understanding mixing efficiency in the oceans. Do the nonlinearities of the equation of state matter?

There exist two central measures of turbulent mixing in turbulent stratified fluids that are both caused by molecular diffusion: 1) the dissipation rate D(APE) of available potential energy APE; 2) the turbulent rate of change Wr, turbulent of background gravitational potential energy GPEr. So far, these two quantities have often been regarded as the same energy conversion, namely the irreversible conversion of APE into GPEr, owing to the well known exact equality D(APE)=Wr, turbulent for a Boussinesq fluid with a linear equation of state. Recently, however, Tailleux (2009) pointed out that the above equality no longer holds for a thermally-stratified compressible, with the ratio ξ=Wr, turbulent/D(APE) being generally lower than unity and sometimes even negative for water or seawater, and argued that D(APE) and Wr, turbulent actually represent two distinct types of energy conversion, respectively the dissipation of APE into one particular subcomponent of internal energy called the "dead" internal energy IE0, and the conversion between GPEr and a different subcomponent of internal energy called "exergy" IEexergy. In this paper, the behaviour of the ratio ξ is examined for different stratifications having all the same buoyancy frequency N vertical profile, but different vertical profiles of the parameter Υ=α P/(ρCp), where α is the thermal expansion coefficient, P the hydrostatic pressure, ρ the density, and Cp the specific heat capacity at constant pressure, the equation of state being that for seawater for different particular constant values of salinity. It is found that ξ and Wr, turbulent depend critically on the sign and magnitude of dΥ/dz, in contrast with D(APE), which appears largely unaffected by the latter. These results have important consequences for how the mixing efficiency should be defined and measured in practice, which are discussed.

Mass and size distribution measurement of particulates from a gas turbine combustor using modern mobility analyzer and particle sizer

SMPS and DMS500 analysers were used to measure particulate size distributions in the exhaust of a fully annular aero gas turbine engine at two operating conditions to compare and analyse sources of discrepancy. A number of different dilution ratio values were utilised for the comparative analysis, and a Dekati hot diluter operating at a temperature of 623°K was also utilised to remove volatile PM prior to measurements being made. Additional work focused on observing the effect of varying the sample line temperatures to ascertain the impact. Explanations are offered for most of the trends observed, although a new, repeatable event identified in the range from 417°K to 423°K – where there was a three order of magnitude increase in the nucleation mode of the sample – requires further study.

Convectively coupled equatorial waves in high resolution Hadley centre climate models

Current global atmospheric models fail to simulate well organised tropical phenomena in which convection interacts with dynamics and physics. A new methodology to identify convectively coupled equatorial waves, developed by NCAS-Climate, has been applied to output from the two latest models of the Met Office/Hadley Centre which have fundamental differences in dynamical formulation. Variability, horizontal and vertical structures, and propagation characteristics of tropical convection and equatorial waves, along with their coupled behaviour in the models are examined and evaluated against a previous comprehensive study of observations. It is shown that, in general, the models perform well for equatorial waves coupled with off-equatorial convection. However they perform poorly for waves coupled with equatorial convection. The vertical structure of the simulated wave is not conducive to energy conversion/growth and does not support the correct physical-dynamical coupling that occurs in the real world. The following figure shows an example of the Kelvin wave coupled with equatorial convection. It shows that the models fail to simulate a key feature of convectively coupled Kelvin wave in observations, namely near surface anomalous equatorial zonal winds together with intensified equatorial convection and westerly winds in phase with the convection. The models are also not able to capture the observed vertical tilt structure and the vertical propagation of the Kelvin wave into the lower stratosphere as well as the secondary peak in the mid-troposphere, particularly in HadAM3. These results can be used to provide a test-bed for experimentation to improve the coupling of physics and dynamics in climate and weather models.

Tracking and mean structure of easterly waves over the Intra-Americas Sea

Easterly waves (EWs) are prominent features of the intertropical convergence zone (ITCZ), found in both the Atlantic and Pacific during the Northern Hemisphere summer and fall, where they commonly serve as precursors to hurricanes over both basins.Alarge proportion of Atlantic EWs are known to form over Africa, but the origin of EWs over the Caribbean and east Pacific in particular has not been established in detail. In this study reanalyses are used to examine the coherence of the large-scale wave signatures and to obtain track statistics and energy conversion terms for EWs across this region. Regression analysis demonstrates that some EW kinematic structures readily propagate between the Atlantic and east Pacific, with the highest correlations observed across Costa Rica and Panama. Track statistics are consistent with this analysis and suggest that some individual waves are maintained as they pass from the Atlantic into the east Pacific, whereas others are generated locally in the Caribbean and east Pacific. Vortex anomalies associated with the waves are observed on the leeward side of the Sierra Madre, propagating northwestward along the coast, consistent with previous modeling studies of the interactions between zonal flow and EWs with model topography similar to the Sierra Madre. An energetics analysis additionally indicates that the Caribbean low-level jet and its extension into the east Pacific—known as the Papagayo jet—are a source of energy for EWs in the region. Two case studies support these statistics, as well as demonstrate the modulation of EW track and storm development location by the MJO.

Baroclinic stationary waves in aquaplanet models

An aquaplanet model is used to study the nature of the highly persistent low-frequency waves that have been observed in models forced by zonally symmetric boundary conditions. Using the Hayashi spectral analysis of the extratropical waves, the authors find that a quasi-stationary wave 5 belongs to a wave packet obeying a well-defined dispersion relation with eastward group velocity. The components of the dispersion relation with k ≥ 5 baroclinically convert eddy available potential energy into eddy kinetic energy, whereas those with k < 5 are baroclinically neutral. In agreement with Green’s model of baroclinic instability, wave 5 is weakly unstable, and the inverse energy cascade, which had been previously proposed as a main forcing for this type of wave, only acts as a positive feedback on its predominantly baroclinic energetics. The quasi-stationary wave is reinforced by a phase lock to an analogous pattern in the tropical convection, which provides further amplification to the wave. It is also found that the Pedlosky bounds on the phase speed of unstable waves provide guidance in explaining the latitudinal structure of the energy conversion, which is shown to be more enhanced where the zonal westerly surface wind is weaker. The wave’s energy is then trapped in the waveguide created by the upper tropospheric jet stream. In agreement with Green’s theory, as the equator-to-pole SST difference is reduced, the stationary marginally stable component shifts toward higher wavenumbers, while wave 5 becomes neutral and westward propagating. Some properties of the aquaplanet quasi-stationary waves are found to be in interesting agreement with a low frequency wave observed by Salby during December–February in the Southern Hemisphere so that this perspective on low frequency variability, apart from its value in terms of basic geophysical fluid dynamics, might be of specific interest for studying the earth’s atmosphere.

Irreversible compressible work and APE dissipation in turbulent stratified fluid

Although it plays a key role in the theory of stratified turbulence, the concept of available potential energy (APE) dissipation has remained until now a rather mysterious quantity, owing to the lack of rigorous result about its irreversible character or energy conversion type. Here, we show by using rigorous energetics considerations rooted in the analysis of the Navier-Stokes for a fully compressible fluid with a nonlinear equation of state that the APE dissipation is an irreversible energy conversion that dissipates kinetic energy into internal energy, exactly as viscous dissipation. These results are established by showing that APE dissipation contributes to the irreversible production of entropy, and by showing that it is a part of the work of expansion/contraction. Our results provide a new interpretation of the entropy budget, that leads to a new exact definition of turbulent effective diffusivity, which generalizes the Osborn-Cox model, as well as a rigorous decomposition of the work of expansion/contraction into reversible and irreversible components. In the context of turbulent mixing associated with parallel shear flow instability, our results suggests that there is no irreversible transfer of horizontal momentum into vertical momentum, as seems to be required when compressible effects are neglected, with potential consequences for the parameterisations of momentum dissipation in the coarse-grained Navier-Stokes equations.

Automated algorithms for multilayer thin film filter design using metamaterials

Currently, infrared filters for astronomical telescopes and satellite radiometers are based on multilayer thin film stacks of alternating high and low refractive index materials. However, the choice of suitable layer materials is limited and this places limitations on the filter performance that can be achieved. The ability to design materials with arbitrary refractive index allows for filter performance to be greatly increased but also increases the complexity of design. Here a differential algorithm was used as a method for optimised design of filters with arbitrary refractive indices, and then materials are designed to these specifications as mono-materials with sub wavelength structures using Bruggeman’s effective material approximation (EMA).

A comparison of greenhouse gas emissions in the residential sector of major Canadian cities

One of the most significant sources of greenhouse gas (GHG) emissions in Canada is the buildings sector, with over 30% of national energy end-use occurring in buildings. Energy use must be addressed to reduce emissions from the buildings sector, as nearly 70% of all Canada’s energy used in the residential sector comes from fossil sources. An analysis of GHG emissions from the existing residential building stock for the year 2010 has been conducted for six Canadian cities with different climates and development histories: Vancouver, Edmonton, Winnipeg, Toronto, Montreal, and Halifax. Variation across these cities is seen in their 2010 GHG emissions, due to climate, characteristics of the building stock, and energy conversion technologies, with Halifax having the highest per capita emissions at 5.55 tCO2e/capita and Montreal having the lowest at 0.32 tCO2e/capita. The importance of the provincial electricity grid’s carbon intensity is emphasized, along with era of construction, occupancy, floor area, and climate. Approaches to achieving deep emissions reductions include innovative retrofit financing and city level residential energy conservation by-laws; each region should seek location-appropriate measures to reduce energy demand within its residential housing stock, as well as associated GHG emissions.

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.