Heat transport in harmonic oscillator systems with thermal baths: application to optomechanical arrays

We investigate the transport of phonons between N harmonic oscillators in contact with independent thermal baths and coupled to a common oscillator, and derive an expression for the steady state heat flow between the oscillators in the weak coupling limit. We apply these results to an optomechanical array consisting of a pair of mechanical resonators coupled to a single quantized electromagnetic field mode by radiation pressure as well as to thermal baths with different temperatures. In the weak coupling limit this system is shown to be equivalent to two mutually-coupled harmonic oscillators in contact with an effective common thermal bath in addition to their independent baths. The steady state occupation numbers and heat flows are derived and discussed in various regimes of interest.

Nonconservative current-driven dynamics: beyond the nanoscale

Long metallic nanowires combine crucial factors for nonconservative current-driven atomic motion. These systems have degenerate vibrational frequencies, clustered about a Kohn anomaly in the dispersion relation, that can couple under current to form nonequilibrium modes of motion growing exponentially in time. Such motion is made possible by nonconservative current-induced forces on atoms, and we refer to it generically as the waterwheel effect. Here the connection between the waterwheel effect and the stimulated directional emission of phonons propagating along the electron flow is discussed in an intuitive manner. Nonadiabatic molecular dynamics show that waterwheel modes self-regulate by reducing the current and by populating modes in nearby frequency, leading to a dynamical steady state in which nonconservative forces are counter-balanced by the electronic friction. The waterwheel effect can be described by an appropriate effective nonequilibrium dynamical response matrix. We show that the current-induced parts of this matrix in metallic systems are long-ranged, especially at low bias. This nonlocality is essential for the characterisation of nonconservative atomic dynamics under current beyond the nanoscale.

Electromechanical properties of engineered lead free potassium sodium niobate based materials

K0.5Na0.5NbO3 (KNN), is the most promising lead free material for substituting lead zirconate titanate (PZT) which is still the market leader used for sensors and actuators. To make KNN a real competitor, it is necessary to understand and to improve its properties. This goal is pursued in the present work via different approaches aiming to study KNN intrinsic properties and then to identify appropriate strategies like doping and texturing for designing better KNN materials for an intended application. Hence, polycrystalline KNN ceramics (undoped, non-stoichiometric; NST and doped), high-quality KNN single crystals and textured KNN based ceramics were successfully synthesized and characterized in this work. Polycrystalline undoped, non-stoichiometric (NST) and Mn doped KNN ceramics were prepared by conventional ceramic processing. Structure, microstructure and electrical properties were measured. It was observed that the window for mono-phasic compositions was very narrow for both NST ceramics and Mn doped ceramics. For NST ceramics the variation of A/B ratio influenced the polarization (P-E) hysteresis loop and better piezoelectric and dielectric responses could be found for small stoichiometry deviations (A/B = 0.97). Regarding Mn doping, as compared to undoped KNN which showed leaky polarization (P-E) hysteresis loops, B-site Mn doped ceramics showed a well saturated, less-leaky hysteresis loop and a significant properties improvement. Impedance spectroscopy was used to assess the role of Mn and a relation between charge transport – defects and ferroelectric response in K0.5Na0.5NbO3 (KNN) and Mn doped KNN ceramics could be established. At room temperature the conduction in KNN which is associated with holes transport is suppressed by Mn doping. Hence Mn addition increases the resistivity of the ceramic, which proved to be very helpful for improving the saturation of the P-E loop. At high temperatures the conduction is dominated by the motion of ionized oxygen vacancies whose concentration increases with Mn doping. Single crystals of potassium sodium niobate (KNN) were grown by a modified high temperature flux method. A boron-modified flux was used to obtain the crystals at a relatively low temperature. XRD, EDS and ICP analysis proved the chemical and crystallographic quality of the crystals. The grown KNN crystals exhibit higher dielectric permittivity (29,100) at the tetragonal-to-cubic phase transition temperature, higher remnant polarization (19.4 μC/cm2) and piezoelectric coefficient (160 pC/N) when compared with the standard KNN ceramics. KNN single crystals domain structure was characterized for the first time by piezoforce response microscopy. It could be observed that <001> - oriented potassium sodium niobate (KNN) single crystals reveal a long range ordered domain pattern of parallel 180° domains with zig-zag 90° domains. From the comparison of KNN Single crystals to ceramics, It is argued that the presence in KNN single crystal (and absence in KNN ceramics) of such a long range order specific domain pattern that is its fingerprint accounts for the improved properties of single crystals. These results have broad implications for the expanded use of KNN materials, by establishing a relation between the domain patterns and the dielectric and ferroelectric response of single crystals and ceramics and by indicating ways of achieving maximised properties in KNN materials. Polarized Raman analysis of ferroelectric potassium sodium niobate (K0.5Na0.5)NbO3 (KNN) single crystals was performed. For the first time, an evidence is provided that supports the assignment of KNN single crystals structure to the monoclinic symmetry at room temperature. Intensities of A′, A″ and mixed A′+A″ phonons have been theoretically calculated and compared with the experimental data in dependence of crystal rotation, which allowed the precise determination of the Raman tensor coefficients for (non-leaking) modes in monoclinic KNN. In relation to the previous literature, this study clarifies that assigning monoclinic phase is more suitable than the orthorhombic one. In addition, this study is the basis for non-destructive assessments of domain distribution by Raman spectroscopy in KNN-based lead-free ferroelectrics with complex structures. Searching a deeper understanding of the electrical behaviour of both KNN single crystal and polycrystalline materials for the sake of designing optimized KNN materials, a comparative study at the level of charge transport and point defects was carried out by impedance spectroscopy. KNN single crystals showed lower conductivity than polycrystals from room temperature up to 200 ºC, but above this temperature polycrystalline KNN displays lower conductivity. The low temperature (T < 200 ºC) behaviour reflects the different processing conditions of both ceramics and single crystals, which account for less defects prone to charge transport in the case of single crystals. As temperature increases (T > 200 ºC) single crystals become more conductive than polycrystalline samples, in which grain boundaries act as barriers to charge transport. For even higher temperatures the conductivity difference between both is increased due to the contribution of ionic conduction in single crystals. Indeed the values of activation energy calculated to the high temperature range (T > 300 ºC) were 1.60 and 0.97 eV, confirming the charge transport due to ionic conduction and ionized oxygen vacancies in single crystals and polycrystalline KNN, respectively. It is suggested that single crystals with low defects content and improved electromechanical properties could be a better choice for room temperature applications, though at high temperatures less conductive ceramics may be the choice, depending on the targeted use. Aiming at engineering the properties of KNN polycrystals towards the performance of single crystals, the preparation and properties study of (001) – oriented (K0.5Na0.5)0.98Li0.02NbO3 (KNNL) ceramics obtained by templated grain growth (TGG) using KNN single crystals as templates was undertaken. The choice of KNN single crystals templates is related with their better properties and to their unique domain structure which were envisaged as a tool for templating better properties in KNN ceramics too. X-ray diffraction analysis revealed for the templated ceramics a monoclinic structure at room temperature and a Lotgering factor (f) of 40% which confirmed texture development. These textured ceramics exhibit a long range ordered domain pattern consisting of 90º and 180º domains, similar to the one observed in the single crystals. Enhanced dielectric (13017 at TC), ferroelectric (2Pr = 42.8 μC/cm2) and piezoelectric (d33 = 280 pC/N) properties are observed for textured KNNL ceramics as compared to the randomly oriented ones. This behaviour is suggested to be due to the long range ordered domain patterns observed in the textured ceramics. The obtained results as compared with the data previously reported on texture KNN based ceramics confirm that superior properties were found due to ordered repeated domain pattern. This study provides an useful approach towards properties improvement of KNN-based piezoelectric ceramics. Overall, the present results bring a significant contribution to the pool of knowledge on the properties of sodium potassium niobate materials: a relation between the domain patterns and di-, ferro-, and piezo-electric response of single crystals and ceramics was demonstrated and ways of engineering maximised properties in KNN materials, for example by texturing were established. This contribution is envisaged to have broad implications for the expanded use of KNN over the alternative lead-based materials.

Phonon spectra and thermodynamic properties of rare gas solids based on empirical and semi-empirical (ab initio) two-body potentials : a comparative study /

We study the phonon dispersion, cohesive and thermal properties of raxe gas solids Ne, Ar, Kr, and Xe, using a variety of potentials obtained from different approaches; such as, fitting to crystal properties, purely ab initio calculations for molecules and dimers or ab initio calculations for solid crystalline phase, a combination of ab initio calculations and fitting to either gas phase data or sohd state properties. We explore whether potentials derived with a certain approaxih have any obvious benefit over the others in reproducing the solid state properties. In particular, we study phonon dispersion, isothermal ajid adiabatic bulk moduli, thermal expansion, and elastic (shear) constants as a function of temperatiue. Anharmonic effects on thermal expansion, specific heat, and bulk moduli have been studied using A^ perturbation theory in the high temperature limit using the neaxest-neighbor central force (nncf) model as developed by Shukla and MacDonald [4]. In our study, we find that potentials based on fitting to the crystal properties have some advantage, particularly for Kr and Xe, in terms of reproducing the thermodynamic properties over an extended range of temperatiures, but agreement with the phonon frequencies with the measured values is not guaranteed. For the lighter element Ne, the LJ potential which is based on fitting to the gas phase data produces best results for the thermodynamic properties; however, the Eggenberger potential for Ne, where the potential is based on combining ab initio quantum chemical calculations and molecular dynamics simulations, produces results that have better agreement with the measured dispersion, and elastic (shear) values. For At, the Morse-type potential, which is based on M0ller-Plesset perturbation theory to fourth order (MP4) ab initio calculations, yields the best results for the thermodynamic properties, elastic (shear) constants, and the phonon dispersion curves.

Phonon spectra and thermal properties of some fcc metals using embedded-atom potentials /

By employing the embedded-atom potentials of Mei et ai.[l], we have calculated the dynamical matrices and phonon dispersion curves for six fee metals (Cu,Ag,Au,Ni,Pd and Pt). We have also investigated, within the quasiharmonic approximation, some other thermal properties of these metals which depend on the phonon density of states, such as the temperature dependence of lattice constant, coefficient of linear thermal expansion, isothermal and adiabatic bulk moduli, heat capacities at constant volume and constant pressure, Griineisen parameter and Debye temperature. The computed results are compared with the experimental findings wherever possible. The comparison shows a generally good agreement between the theoretical values and experimental data for all properties except the discrepancies of phonon frequencies and Debye temperature for Pd, Pt and Au. Further, we modify the parameters of this model for Pd and Pt and obtain the phonon dispersion curves which is in good agreement with experimental data.

An analysis of decoupling procedures in generating thermal Green's functions of O([lambda]â ´) by the Zubarev equation of motion method

The Zubarev equation of motion method has been applied to an anharmonic crystal of O( ,,4). All possible decoupling schemes have been interpreted in order to determine finite temperature expressions for the one phonon Green's function (and self energy) to 0()\4) for a crystal in which every atom is on a site of inversion symmetry. In order to provide a check of these results, the Helmholtz free energy expressions derived from the self energy expressions, have been shown to agree in the high temperature limit with the results obtained from the diagrammatic method. Expressions for the correlation functions that are related to the mean square displacement have been derived to 0(1\4) in the high temperature limit.

On the anharmonic, multi-phonon, Debye-Waller contributions to the phonon-limited resistivity of metals : applications to Na and K

The anharmonic, multi-phonon (MP), and Oebye-Waller factor (OW) contributions to the phonon limited resistivity (;0) of metals derived by Shukla and Muller (1979) by the doubletime temperature dependent Green function method have been numerically evaluated for Na and K in the high temperature limit. The anharmonic contributions arise from the cubic and quartic shift of phonons (CS, QS), and phonon width (W) and the interference term (1). The QS, MP and OW contributions to I' are also derived by the matrix element method and the results are in agreement with those of Shukla and Muller (1979). In the high temperature limit, the contributions to;O from each of the above mentioned terms are of the type BT2 For numerical calculations suitable expressions are derived for the anharmonic contributions to ~ in terms of the third and fourth rank tensors obtained by the Ewald procedure. The numerical calculation of the contributions to;O from the OW, MP term and the QS have been done exactly and from the CS, Wand I terms only approximately in the partial and total Einstein approximations (PEA, TEA), using a first principle approach (Shukla and Taylor (1976)). The results obtained indicate that there is a strong pairwise cancellation between the: OW and MP terms, the QS and CS and the Wand I terms. The sum total of these contributions to;O for Na and K amounts to 4 to 11% and 2 to 7%, respectively, in the PEA while in the TEA they amount to 3 to 7% and 1 to 4%, respectively, in the temperature range.

The finite bandwidth model for spin fluctuations in Pd

The frequency dependence of the electron-spin fluctuation spectrum, P(Q), is calculated in the finite bandwidth model. We find that for Pd, which has a nearly full d-band, the magnitude, the range, and the peak frequency of P(Q) are greatly reduced from those in the standard spin fluctuation theory. The electron self-energy due to spin fluctuations is calculated within the finite bandwidth model. Vertex corrections are examined, and we find that Migdal's theorem is valid for spin fluctuations in the nearly full band. The conductance of a normal metal-insulator-normal metal tunnel junction is examined when spin fluctuations are present in one electrode. We find that for the nearly full band, the momentum independent self-energy due to spin fluctuations enters the expression for the tunneling conductance with approximately the same weight as the self-energy due to phonons. The effect of spin fluctuations on the tunneling conductance is slight within the finite bandwidth model for Pd. The effect of spin fluctuations on the tunneling conductance of a metal with a less full d-band than Pd may be more pronounced. However, in this case the tunneling conductance is not simply proportional to the self-energy.

Thermodynamic properties and Debye-Waller factor of fee materials

We have calculated the equation of state and the various thermodynamic properties of monatomic fcc crystals by minimizing the Helmholtz free energy derived in the high temperature limit for the quasiharmonic theory, QH, and the lowest-order (cubic and quartic), 'A2, anharmonic terms of the perturbation theory, PT. The total energy in each case is obtained by adding the static energy. The calculation of the thermal properties was carried out for a nearest-neighbour central-force model of the fcc lattice by means of the appropriate thermodynamic relations. We have calculated the lattice constant, the thermal expansion, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the isothermal and adiabatic bulk moduli, and the Griineisen parameter, for the rare-gas solids Kr and Xe, and gold. Morse potential and modified Morse potential were each used to represent the atomic interaction for the three fcc materials. For most of the calculated thermodynamic properties from the QH theory, the results for Kr and Xe with the modified Morse potential show an improvement over the results for the Morse potential when compared with the experimental data. However, the results of the 'A 2 equation of state with the modified Morse potential are in good agreement with experiment only in the case of the specific heat at constant volume and at constant pressure. For Au we have calculated the lattice contribution from the QH and 'A 2 PT and the electronic contribution to the thermal properties. The electronic contribution was taken into account by using the free electron model. The results of the thermodynamic properties calculated with the modified Morse potential were similar to those obtained with the Morse potential. U sing the minimized equation of state we also calculated the Mossbauer recoilless fraction for Kr and Xe and the Debye-Waller factor (DWF) for Pb, AI, eu, Ag, and Au. The Mossbauer recoilless fraction was obtained for the above two potentials and Lennard-Jones potential. The L-J potential gives the best agreement with experiment for Kr. No experimental data exists for Xe. At low temperature the calculated DWF results for Pb, AI, and eu show a good agreement with experimental values, but at high temperature the experimental DWF results increase very rapidly. For Ag the computed values were below the expected results at all temperatures. The DWF results of the modified Morse potential for Pb, AI, eu and Ag were slightly better than those of the Morse potential. In the case of Au the calculated values were in poor agreement with experimental results. We have calculated the quasiharmonic phonon dispersion curves for Kr, Xe, eu, Ag, and Au. The calculated and experimental results of the frequencies agree quite well for all the materials except for Au where the longitudinal modes show serious discrepancies with the experimental results. In addition, the two lowest-order anharmonic contributions to the phonon frequency were derived using the Green's function method. The A 2 phonon dispersion curves have been calculated only for eu, and the results were similar to those of the QH dispersion curves. Finally, an expression for the Griineisen parameter "( has been derived from the anharmonic frequencies, and calculated for these materials. The "( results are comparable with those obtained from the thermodynamic definition.

Etude ab initio du pnicture de fer supraconducteur LaOFeAs

Le présent mémoire traite de la description du LaOFeAs, le premier matériau découvert de la famille des pnictures de fer, par la théorie de la fonctionnelle de la densité (DFT). Plus particulièrement, nous allons exposer l’état actuel de la recherche concernant ce matériau avant d’introduire rapidement la DFT. Ensuite, nous allons regarder comment se comparent les paramètres structuraux que nous allons calculer sous différentes phases par rapport aux résultats expérimentaux et avec les autres calculs DFT dans la littérature. Nous allons aussi étudier en détails la structure électronique du matériau sous ses différentes phases magnétiques et structurales. Nous emploierons donc les outils normalement utilisés pour mieux comprendre la structure électronique : structures de bandes, densités d’états, surfaces de Fermi, nesting au niveau de Fermi. Nous tirerons profit de la théorie des groupes afin de trouver les modes phononiques permis par la symétrie de notre cristal. De plus, nous étudierons le couplage électrons-phonons pour quelques modes. Enfin, nous regarderons l’effet de différentes fonctionnelles sur nos résultats pour voir à quel point ceux-ci sont sensibles à ce choix. Ainsi, nous utiliserons la LDA et la PBE, mais aussi la LDA+U et la PBE+U.

Propriétés optiques dans l'infrarouge des nanotubes de carbone et du graphène

Les nanotubes de carbone et le graphène sont des nanostructures de carbone hybridé en sp2 dont les propriétés électriques et optiques soulèvent un intérêt considérable pour la conception d’une nouvelle génération de dispositifs électroniques et de matériaux actifs optiquement. Or, de nombreux défis demeurent avant leur mise en œuvre dans des procédés industriels à grande échelle. La chimie des matériaux, et spécialement la fonctionnalisation covalente, est une avenue privilégiée afin de résoudre les difficultés reliées à la mise en œuvre de ces nanostructures. La fonctionnalisation covalente a néanmoins pour effet de perturber la structure cristalline des nanostructures de carbone sp2 et, par conséquent, d’affecter non seulement lesdites propriétés électriques, mais aussi les propriétés optiques en émanant. Il est donc primordial de caractériser les effets des défauts et du désordre dans le but d’en comprendre les conséquences, mais aussi potentiellement d’en exploiter les retombées. Cette thèse traite des propriétés optiques dans l’infrarouge des nanotubes de carbone et du graphène, avec pour but de comprendre et d’expliquer les mécanismes fondamentaux à l’origine de la réponse optique dans l’infrarouge des nanostructures de carbone sp2. Soumise à des règles de sélection strictes, la spectroscopie infrarouge permet de mesurer la conductivité en courant alternatif à haute fréquence des matériaux, dans une gamme d’énergie correspondant aux vibrations moléculaires, aux modes de phonons et aux excitations électroniques de faible énergie. Notre méthode expérimentale consiste donc à explorer un espace de paramètres défini par les trois axes que sont i. la dimensionnalité du matériau, ii. le potentiel chimique et iii. le niveau de désordre, ce qui nous permet de dégager les diverses contributions aux propriétés optiques dans l’infrarouge des nanostructures de carbone sp2. Dans un premier temps, nous nous intéressons à la spectroscopie infrarouge des nanotubes de carbone monoparois sous l’effet tout d’abord du dopage et ensuite du niveau de désordre. Premièrement, nous amendons l’origine couramment acceptée du spectre vibrationnel des nanotubes de carbone monoparois. Par des expériences de dopage chimique contrôlé, nous démontrons en effet que les anomalies dans lespectre apparaissent grâce à des interactions électron-phonon. Le modèle de la résonance de Fano procure une explication phénoménologique aux observations. Ensuite, nous établissons l’existence d’états localisés induits par la fonctionnalisation covalente, ce qui se traduit optiquement par l’apparition d’une bande de résonance de polaritons plasmons de surface (nanoantenne) participant au pic de conductivité dans le térahertz. Le dosage du désordre dans des films de nanotubes de carbone permet d’observer l’évolution de la résonance des nanoantennes. Nous concluons donc à une segmentation effective des nanotubes par les greffons. Enfin, nous montrons que le désordre active des modes de phonons normalement interdits par les règles de sélection de la spectroscopie infrarouge. Les collisions élastiques sur les défauts donnent ainsi accès à des modes ayant des vecteurs d’onde non nuls. Dans une deuxième partie, nous focalisons sur les propriétés du graphène. Tout d’abord, nous démontrons une méthode d’électrogreffage qui permet de fonctionnaliser rapidement et à haute densité le graphène sans égard au substrat. Par la suite, nous utilisons l’électrogreffage pour faire la preuve que le désordre active aussi des anomalies dépendantes du potentiel chimique dans le spectre vibrationnel du graphène monocouche, des attributs absents du spectre d’un échantillon non fonctionnalisé. Afin d’expliquer le phénomène, nous présentons une théorie basée sur l’interaction de transitions optiques intrabandes, de modes de phonons et de collisions élastiques. Nous terminons par l’étude du spectre infrarouge du graphène comportant des îlots de bicouches, pour lequel nous proposons de revoir la nature du mécanisme de couplage à l’œuvre à la lumière de nos découvertes concernant le graphène monocouche.

Photothermal deflection measurement on heat transport in GaAs epitaxial layers

In this paper, we report the in-plane and cross-plane measurements of the thermal diffusivity of double epitaxial layers of n-type GaAs doped with various concentrations of Si and a p-type Be-doped GaAs layer grown on a GaAs substrate by the molecular beam epitaxial method, using the laser-induced nondestructive photothermal deflection technique. The thermal diffusivity value is evaluated from the slope of the graph of the phase of the photothermal deflection signal as a function of pump-probe offset. Analysis of the data shows that the cross-plane thermal diffusivity is less than that of the in-plane thermal diffusivity. It is also seen that the doping concentration has a great influence on the thermal diffusivity value. Measurement of p-type Be-doped samples shows that the nature of the dopant also influences the effective thermal diffusivity value. The results are interpreted in terms of a phonon-assisted heat transfer mechanism and the various scattering process involved in the propagation of phonons.

Photothermal deflection measurement on heat transport in GaAs epitaxial layers

In this paper, we report the in-plane and cross-plane measurements of the thermal diffusivity of double epitaxial layers of n-type GaAs doped with various concentrations of Si and a p-type Be-doped GaAs layer grown on a GaAs substrate by the molecular beam epitaxial method, using the laser-induced nondestructive photothermal deflection technique. The thermal diffusivity value is evaluated from the slope of the graph of the phase of the photothermal deflection signal as a function of pump-probe offset. Analysis of the data shows that the cross-plane thermal diffusivity is less than that of the in-plane thermal diffusivity. It is also seen that the doping concentration has a great influence on the thermal diffusivity value. Measurement of p-type Be-doped samples shows that the nature of the dopant also influences the effective thermal diffusivity value. The results are interpreted in terms of a phonon-assisted heat transfer mechanism and the various scattering process involved in the propagation of phonons

Photothermal deflection measurement on heat transport in GaAs epitaxial layers

In this paper, we report the in-plane and cross-plane measurements of the thermal diffusivity of double epitaxial layers of n-type GaAs doped with various concentrations of Si and a p-type Be-doped GaAs layer grown on a GaAs substrate by the molecular beam epitaxial method, using the laser-induced nondestructive photothermal deflection technique. The thermal diffusivity value is evaluated from the slope of the graph of the phase of the photothermal deflection signal as a function of pump-probe offset. Analysis of the data shows that the cross-plane thermal diffusivity is less than that of the in-plane thermal diffusivity. It is also seen that the doping concentration has a great influence on the thermal diffusivity value. Measurement of p-type Be-doped samples shows that the nature of the dopant also influences the effective thermal diffusivity value. The results are interpreted in terms of a phonon-assisted heat transfer mechanism and the various scattering process involved in the propagation of phonons

Thermal characterization of nanocrystalline porous CePO4 ceramics

A laser-induced photoacoustic technique was employed to investigate thermal transport through nanocrystalline CePO4 samples prepared via the sol–gel route. Evaluation of thermal diffusivity was carried out using the one-dimensional model of Rosencwaig and Gersho for the reflection configuration of the photoacoustic method. Structural analyses of samples revealed that they are nanoporous in nature, possessing micron-sized grains. Analysis of results shows that thermal diffusivity value varies with sintering temperature. Results are explained in terms of the variation in porosity with sintering temperature and the effects of various scattering mechanisms on the propagation of phonons through the nanoporous ceramic matrix. Further analyses confirm that apart from porosity, grain boundary resistance and interface thermal resistance influence the effective value of thermal diffusivity of the samples under investigation.