Cette thèse traite des équations de Maxwell en domaine temporel. Le principal objectif est de proposer des méthodes de type éléments finis d'ordre élevé pour les équations de Maxwell et des schémas d'intégration en temps efficaces sur des maillages localement raffinés. Nous considérons des méthodes GDDT (Galerkine Discontinues en Domaine Temporel) s'appuyant sur une interpolation polynomiale d'ordre arbitrairement élevé des composantes du champ électromagnétique. Les méthodes GDDT pour les équations de Maxwell s'appuient le plus souvent sur des schémas d'intégration en temps explicites dont la condition de stabilité peut être très restrictive pour des maillages raffinés. Pour surmonter cette limitation, nous considérons des schémas en temps qui consistent à appliquer un schéma implicite localement, dans les régions raffinées, tout en préservant un schéma explicite sur le reste du maillage. Nous présentons une étude théorique complète et une comparaison de deux méthodes GDDT localement implicites. Des expériences numériques en 2D et 3D illustrent l'utilité des schémas proposés. Le traitement numérique de milieux de propagation complexes est également l'un des objectifs. Nous considérons l'interaction des ondes électromagnétiques avec les tissus biologiques qui est au cœur de nombreuses applications dans le domaine biomédical. La modélisation numérique nécessite alors de résoudre le système de Maxwell avec des modèles appropriés de dispersion. Nous formulons une méthode GDDT localement implicite pour le modèle de Debye et proposons une analyse théorique et numérique complète du schéma.

Le présent projet de maîtrise a pour but d’étudier les interactions optiques des films de nanotubes de carbone (FNTCs) avec les ondes THz. Des expériences d’absorption térahertz faites par spectroscopie THz dans le domaine temporel ont été entreprises sur les films dont l’épaisseur varie. Les films d’épaisseurs allant de 14 à 145 nm, sont des couches minces de nanotubes de carbone (NTCs) empilés les uns sur les autres et sont déposés sur substrats (GaAs et silicium). Une caractérisation comparative des épaisseurs des films est entreprise dans un premier temps par AFM et par ellipsométrie spectroscopique. À cause de la rugosité de la surface et de porosité des films qui compliquent les interactions de la lumière avec les films, les épaisseurs déterminées par AFM sont gardées au détriment de celles d’ellipsométrie. La relation entre les épaisseurs mesurées par AFM en fonction des épaisseurs nominales s’est révélée linéaire. Les couleurs des FNTC sont aussi caractérisées en fonction de leurs épaisseurs. L’expérience d’absorption THz sur les films consiste à enregistrer la transmission d’une impulsion THz à large bande à travers les échantillons. Sur les spectres, on détecte aussi l’impulsion de réflexion, l’écho de réflexion de l’impulsion principale THz à l’intérieur du substrat séparé par un délai temporel. La diminution du pic de l’impulsion principale THz en fonction de l’épaisseur est non linéaire et atteint une saturation pour les films les plus épais. Ce résultat est en lien direct avec les mesures quatre pointes de conductivité dc des films où l’inverse de la résistivité de feuille sature à partir des mêmes épaisseurs de film. L’écho de réflexion de l’impulsion principale à l’intérieur du substrat perd de l’amplitude plus rapidement en fonction de l’épaisseur à cause de près de deux passages supplémentaires de l’impulsion dans le film au moment de la réflexion. Finalement, une disparition de l’impulsion de réflexion à une épaisseur particulière de film (100 nm pour le GaAs et 60 nm pour le Si) démontre les propriétés antiréfléchissantes des FNTCs. / In the present masters project, the goal is to study the optical interactions of carbon nanotube films (CNTFs) with terahertz (THz) waves. The THz absorption experiments made by time domain THz spectroscopy have been undertaken on thickness-variable films. CNTFs, which have their thicknesses range from 14 to 145 nm, are thin CNT layers that are piled one on another are deposited on a substrate (GaAs or silicon). First, a comparative characterization of film thicknesses is undertaken with AFM and with spectroscopic ellipsometry. Because of surface rugosity and film porosity which has the effect of complexifying the interaction of light with the films, AFM thicknesses are held for the rest of the analysis instead of those determined with ellipsometry. AFM measured thicknesses scale linearly with respect to nominal thicknesses that are proportional to the CNT density. CNTFs’ colors reveal to be correlated with their thicknesses. THz absorption experiments consist of taking the transmission spectrum of a broad band THz pulse through the samples. On the spectra, we also detect the reflection pulse, which is the echo of the main THz pulse inside the substrate separated by a time delay. The decrease of the main THz pulse with respect to the film thickness is non linear and reaches a saturation plateau for the thickest films. This finding is in direct relationship with four-point probe sheet conductivity measurements made on the films where a saturation is also observed from the same thicknesses. The reflection pulse loses amplitude more rapidly as the film thickness increases because of two additional wave passages in the film during reflection. Lastly, a quenching of the reflection pulse which is observed at a particular film thickness (100 nm for GaAs and 60 nm for silicon) demonstrates antireflection properties for the CNTFs. / Les expériences de spectroscopie ont été réalisées en collaboration avec Jean-François Allard du groupe de Denis Morris de l'Université de Sherbrooke.

Cette thèse traite des équations de Maxwell en domaine temporel. Le principal objectif est de proposer des méthodes de type éléments finis d'ordre élevé pour les équations de Maxwell et des schémas d'intégration en temps efficaces sur des maillages localement raffinés. Nous considérons des méthodes GDDT (Galerkine Discontinues en Domaine Temporel) s'appuyant sur une interpolation polynomiale d'ordre arbitrairement élevé des composantes du champ électromagnétique. Les méthodes GDDT pour les équations de Maxwell s'appuient le plus souvent sur des schémas d'intégration en temps explicites dont la condition de stabilité peut être très restrictive pour des maillages raffinés. Pour surmonter cette limitation, nous considérons des schémas en temps qui consistent à appliquer un schéma implicite localement, dans les régions raffinées, tout en préservant un schéma explicite sur le reste du maillage. Nous présentons une étude théorique complète et une comparaison de deux méthodes GDDT localement implicites. Des expériences numériques en 2D et 3D illustrent l'utilité des schémas proposés. Le traitement numérique de milieux de propagation complexes est également l'un des objectifs. Nous considérons l'interaction des ondes électromagnétiques avec les tissus biologiques qui est au cœur de nombreuses applications dans le domaine biomédical. La modélisation numérique nécessite alors de résoudre le système de Maxwell avec des modèles appropriés de dispersion. Nous formulons une méthode GDDT localement implicite pour le modèle de Debye et proposons une analyse théorique et numérique complète du schéma.

Les lasers à cascade quantique (LCQ) sont des nanostructures de semiconducteurs se basant sur des transitions intersousbandes entre états confinés de la bande de conduction. Ils ont permis de combler un manque de sources puissantes et compactes d'abord dans le moyen infrarouge (MIR) puis dans le térahertz (THz). Cette thèse présente deux études en rapport avec ces lasers. La première partie présente les propriétés d'optique non-linéaire résonantes des LCQ. Il s'agit de démontrer la génération de différence de fréquences entre un faisceau proche infrarouge (IR) et le champ THz du LCQ. L'excitation proche IR est résonante avec les transitions interbandes des puits quantiques qui composent le LCQ. Ceci exalte la susceptibilité non-linéaire du milieu. Le champ THz intense dans la cavité combiné à cette excitation résonante permet d'obtenir de bonnes efficacités (jusqu'à 0.13%) et de générer des harmoniques supérieures jusqu'à l'ordre 3. Ces interactions non-linéaires ont également été étudiées dans les LCQ MIR ce qui a permis d'augmenter la température de fonctionnement jusqu'à 275 K. Une deuxième partie traite du contrôle de la phase du champ THz d'un LCQ au moyen d'un montage de spectroscopie THz dans le domaine temporel. L'originalité résulte dans l'usage d'un LCQ ayant un guide double métal. Ces guides permettent d'avoir de meilleures performances en température mais ont des dimensions largement sous longueur d'onde. Ces dimensions compliquent le couplage d'une onde THz externe nécessaire à l'amorçage du champ THz du LCQ sur une phase fixe. Des antennes en forme de V sont déposées à la surface du LCQ pour faire une adaptation d'impédance et favoriser ainsi le couplage.

Ce travail de thèse a été réalisé dans le cadre d'un contrat avec l'ANR (ANR-09-VERS-013 « THID ») et porte sur le développement d'une nouvelle génération de tags Chipless à bas coûts fonctionnant dans le domaine THz, pour des applications d'identification et/ou authentification unitaire des articles commerciaux, des papiers d'identités, des personnes pour le contrôle d'accès... Les structures proposées, constituées d'un empilement périodique de couches diélectriques d'indices de réfraction différents, utilisent les propriétés particulières des cristaux photoniques 1D de présenter une réponse électromagnétique entrecoupée de bandes interdites photoniques (BIP). Toute perturbation de la périodicité de la structure engendre des pics dans les bandes interdites qui sont utilisés pour coder une information binaire. Cette structuration particulière des matériaux permet donc de manipuler précisément une signature électromagnétique. Pour des raisons liées à l'industrialisation (facilité de fabrication en masse) et aussi de coût, nous avons retenu des matériaux de base déjà couramment utilisés dans l'industrie papetière : le papier et le polyéthylène. Le choix de ces matériaux, qui doivent allier contraste d'indice élevé et faible absorption, représente une étape cruciale dans ce travail. Ainsi, à partir des résultats expérimentaux obtenus par spectroscopie THz dans le domaine temporel (THz-TDS) sur un grand nombre de matériaux, nous avons pu concevoir deux familles de tags sur la base de ces différents matériaux. Par ailleurs, nous avons développé deux méthodes de codage d'une information binaire, toutes deux basées sur l'absence ou la présence de pics dans une BIP, pics dont la position et le nombre dépendent bien évidemment des défauts de périodicité introduits. Pour des applications liées à l'identification, des capacités de codage de près de 20 bits ont été démontrées. Nous avons aussi montré que la richesse d'information contenue dans la réponse électromagnétique de ces Tags THz peut être utilisée pour les applications liées à l'authentification unitaire, en utilisant comme critère de discrimination le coefficient d'autocorrélation. Nous avons ainsi pu évaluer les performances d'un test d'authentification basé sur ce critère dans différents domaines d'analyse : temporel, fréquentiel et temps-fréquence. Nous avons montré qu'une étude du spectrogramme (combinant temps et fréquence) est ainsi bien plus pertinente qu'une étude dans les seuls domaines temporel ou fréquentiel. / This thesis work deals with the development of a new generation of low-cost Chipless tags operating in the THz frequency domain, it has been supported by the french national agency for research (ANR-09-VERS-013 « THID » ). It covers a wide area of applications such as the identification and/or unitary authentication of commercial items, identity papers, access control…To manufacture these tags, we proposed to use a periodic stack of dielectric material layers with different refractive index and whose thickness is of the order of the wavelength, commonly known as a one dimensional photonic crystal. The electromagnetic signature of such a structure exhibits photonic bandgaps (PBG), i.e. frequency windows in which light propagation is prohibited. We suggested modifying the periodicity of the crystal to create defect levels (peaks) for example in the 1st PBG to encode binary information. This particular structure allows to precisely tuned an electromagnetic signature. To ensure a mass and cost effective industrialization, we retained basic materials which are widely used in the pulp and paper industry: paper and polyethylene. The choice of these materials, which must combine high index contrast and low absorption, represents the first and a crucial step in this work. We characterize a wide range of materials using classical THz time domain spectroscopy (THz-TDS) and we propose two families of tags based on paper and polyethylene. Furthermore, we developed two methods to encode binary information, both based on the absence or presence of peaks in a PBG, peaks whose number and position depend on the introduced defects of periodicity. In a real identification test, a coding capacity of nearly 20-bit has been demonstrated. We also showed that the information contained in the electromagnetic response of these THz tags can be used for other applications related to the unitary authentication and by using the correlation coefficient as criterion for discrimination of the different signatures. Therefore, we evaluate the performance of an authentication test based on this criterion in various analysis domains: time, frequency and time-frequency. We showed that a study of the spectrogram (combining time and frequency representation) is much more relevant than a study in the only time or frequency domain.

La tomographie optique diffuse (TOD) est une modalité d’imagerie biomédicale 3D peu dispendieuse et non-invasive qui permet de reconstruire les propriétés optiques d’un tissu biologique. Le processus de reconstruction d’images en TOD est difficile à réaliser puisqu’il nécessite de résoudre un problème non-linéaire et mal posé. Les propriétés optiques sont calculées à partir des mesures de surface du milieu à l’étude. Dans ce projet, deux méthodes de reconstruction non-linéaire pour la TOD ont été développées. La première méthode utilise un modèle itératif, une approche encore en développement qu’on retrouve dans la littérature. L’approximation de la diffusion est le modèle utilisé pour résoudre le problème direct. Par ailleurs, la reconstruction d’image à été réalisée dans différents régimes, continu et temporel, avec des mesures intrinsèques et de fluorescence. Dans un premier temps, un algorithme de reconstruction en régime continu et utilisant des mesures multispectrales est développé pour reconstruire la concentration des chromophores qui se trouve dans différents types de tissus. Dans un second temps, un algorithme de reconstruction est développé pour calculer le temps de vie de différents marqueurs fluorescents à partir de mesures optiques dans le domaine temporel. Une approche innovatrice a été d’utiliser la totalité de l’information du signal temporel dans le but d’améliorer la reconstruction d’image. Par ailleurs, cet algorithme permettrait de distinguer plus de trois temps de vie, ce qui n’a pas encore été démontré en imagerie de fluorescence. La deuxième méthode qui a été développée utilise l’apprentissage machine et plus spécifiquement l’apprentissage profond. Un modèle d’apprentissage profond génératif est mis en place pour reconstruire la distribution de sources d’émissions de fluorescence à partir de mesures en régime continu. Il s’agit de la première utilisation d’un algorithme d’apprentissage profond appliqué à la reconstruction d’images en TOD de fluorescence. La validation de la méthode est réalisée avec une mire aux propriétés optiques connues dans laquelle sont inséres des marqueurs fluorescents. La robustesse de cette méthode est démontrée même dans les situations où le nombre de mesures est limité et en présence de bruit. / Abstract : Diffuse optical tomography (DOT) is a low cost and noninvasive 3D biomedical imaging technique to reconstruct the optical properties of biological tissues. Image reconstruction in DOT is inherently a difficult problem, because the inversion process is nonlinear and ill-posed. During DOT image reconstruction, the optical properties of the medium are recovered from the boundary measurements at the surface of the medium. In this work, two approaches are proposed for non-linear DOT image reconstruction. The first approach relies on the use of iterative model-based image reconstruction, which is still under development for DOT and that can be found in the literature. A 3D forward model is developed based on the diffusion equation, which is an approximation of the radiative transfer equation. The forward model developed can simulate light propagation in complex geometries. Additionally, the forward model is developed to deal with different types of optical data such as continuous-wave (CW) and time-domain (TD) data for both intrinsic and fluorescence signals. First, a multispectral image reconstruction algorithm is developed to reconstruct the concentration of different tissue chromophores simultaneously from a set of CW measurements at different wavelengths. A second image reconstruction algorithm is developed to reconstruct the fluorescence lifetime (FLT) of different fluorescent markers from time-domain fluorescence measurements. In this algorithm, all the information contained in full temporal curves is used along with an acceleration technique to render the algorithm of practical use. Moreover, the proposed algorithm has the potential of being able to distinguish more than 3 FLTs, which is a first in fluorescence imaging. The second approach is based on machine learning techniques, in particular deep learning models. A deep generative model is proposed to reconstruct the fluorescence distribution map from CW fluorescence measurements. It is the first time that such a model is applied for fluorescence DOT image reconstruction. The performance of the proposed algorithm is validated with an optical phantom and a fluorescent marker. The proposed algorithm recovers the fluorescence distribution even from very noisy and sparse measurements, which is a big limitation in fluorescence DOT imaging.

Dans le contexte de l’amélioration des modèles de prévision en acoustique extérieure, ces travaux de thèse se focalisent sur la modélisation des effets des irrégularités de la surface du sol sur la propagation acoustique. Pour ce faire, des méthodes numériques temporelles sont utilisées : d’une part, la méthode FDTD basée sur la résolution des équations de l’acoustique par des schémas aux différences finies, et d’autre part, la méthode des lignes de transmission (TLM). La modélisation des effets de la rugosité de surface est abordée en considérant le formalisme de l’impédance effective. Deux modèles d’impédance effective sont étudiés : le premier caractérise les effets d’une rugosité déterministe constituée de diffuseurs de géométrie constante, le second caractérise les effets moyens d’une rugosité aléatoire définie par un spectre de rugosité. Ce second modèle est validé expérimentalement par une campagne de mesures en salle semi-anéchoïque, audessus de surfaces rugueuses dont la rugosité a été définie très précisément. Les deux modèles d’impédance effective sont également validés par des simulations numériques FDTD et TLM. La possibilité d’implémenter ces conditions d’impédance effective dans les deux codes temporels est ainsi montrée, ce qui permet de modéliser les effets de la rugosité sans avoir à réaliser un maillage précis du profil des hauteurs de la surface du sol. Une campagne de mesures de l’impédance de différents terrains est réalisée afin d’étudier les effets de la variabilité spatiale et saisonnière de l’impédance sur la prévision des niveaux sonores. Les impédances mesurées lors de cette campagne sont également utilisées comme des données d’entrée réalistes pour le code TLM, afin de simuler et d’étudier les effets de la propagation acoustique au-dessus d’un sol hétérogène présentant une impédance spatialement variable. / In the context of prediction models improvement for outdoor sound propagation, this work focuses on the modelling of the effects of ground irregularities on sound propagation. Time-domain numerical methods are used: on one hand, the solving of the governing equations by finite difference schemes (FDTD method), and on the other hand, the transmission line matrix (TLM) method. Effective impedance is considered to model the effects of surface roughness. Two effective impedance models are studied: the first one takes into account the effects of a deterministic roughness formed by scatterers of constant geometry, the second one takes into account the mean effects of a random roughness defined by a roughness spectrum. This second model is validated experimentally by a measurements campaign carried out in a semi-anechoic chamber, above rough surfaces whose roughness profiles were precisely designed. The two effective impedance models are also validated by FDTD and TLM simulations. The possibility to use the effective impedances directly into the numerical methods is then shown, allowing the modelling of roughness effects without meshing the exact height profile of the ground surface. A measurements campaign of the impedance of different grounds is performed in order to assess the effects of space and seasonal variability of the impedance on the sound levels predictions. The results of this campaign are also used as realistic entry data for the TLM code, and the propagation above a heterogeneous ground showing spatially variable impedance is simulated.

Les travaux présentés dans ce mémoire s’inscrivent dans la thématique de la transmission d’énergie sans fil, appliquée à l’alimentation à distance de capteurs, de réseaux de capteurs et d’actionneurs à faible consommation. Cette étude porte sur la conception,l’optimisation, la réalisation et la mesure de circuits Rectennas (Rectifying antennas)compacts, à faible coût et à haut rendement de conversion RF-DC.Un outil d’analyse globale, basé sur la méthode des Différences Finies dans le Domaine Temporel (FDTD), a été développé et utilisé pour prédire avec précision la sortie DC des rectennas étudiées. Les résultats numériques obtenus se sont avérés plus précis et plus complets que ceux de simulations à base d’outils commerciaux. La diode Schottky a été rigoureusement modélisée, en tenant compte de ses éléments parasites et de son boîtier SOT23, et introduite dans le calcul itératif FDTD.Trois rectennas innovantes, en technologie micro-ruban, ont été développées,optimisées et caractérisées expérimentalement. Elles fonctionnent à 2.45 GHz et elles ne contiennent ni filtre d’entrée HF ni vias de retour à la masse. Des rendements supérieurs à 80% ont pu être mesurés avec une densité surfacique de puissance de l’ordre de 0.21 mW/cm²(E = 28 V/m). Une tension DC de 3.1 V a été mesurée aux bornes d’une charge optimale de1.05 k_, lorsque le niveau du champ électrique est égal à 34 V/m (0.31 mW/cm²).Des réseaux de rectennas connectées en série et en parallèle ont été développés. Les tensions et les puissances DC ont été doublées et quadruplées à l’aide de deux et de quatre éléments, respectivement. / The work presented in this thesis is included within the theme of wireless power transmission, applied to wireless powering of sensors, sensor nodes and actuators with low consumption. This study deals with the design, optimization, fabrication and experimental characterization of compact, low cost and efficient Rectennas (Rectifying antennas).A global analysis tool, based on the Finite Difference Time Domain method (FDTD),has been developed and used to predict with a good precision the DC output of studied rectennas. The packaged Schottky diode has been rigorously modeled, taking into account the parasitic elements, and included in the iterative FDTD calculation.Three new rectennas, with microstrip technology, have been developed and measured.They operate at 2.45 GHz and they don’t need neither input HF filter nor via hole connections. Efficiencies more than 80 % have been measured when the power density is 0.21mW/cm² (E = 28 V/m). An output DC voltage of about 3.1 V has been measured with anoptimal load of 1.05 k_, when the power density is equal to 0.31mW/cm² (34 V/m).Rectenna arrays, with series and parallel interconnections, have been developed and measured. Output DC voltages and powers have been doubled and quadrupled using two andfour rectenna elements, respectively.

L’objectif de cette thèse est de développer une méthode Galerkine discontinue d’ordre élevé capable de prendre en considération des simulations réalistes liées à la nanophotonique. Au cours des dernières décennies, l’évolution des techniques de lithographie a permis la création de structure géométriques de tailles nanométriques, révélant ainsi une large gamme de phénomènes nouveaux nés de l’interaction lumière-matière à ces échelles. Ces effets apparaissent généralement pour des objets de taille égale ou (très) inférieure à la longueur d’onde du champ incident. Ce travail repose sur le développement et l’implémentation de modèles de dispersion appropriés (principalement pour les métaux), ainsi que sur un large éventail de méthodes computationnelles classiques. Deux développements méthodologiques majeurs sont présentés et étudiés en détails: (i) les éléments courbes, et (ii) l’ordre d’approximation local. Ces études sont accompagnées de plusieurs cas-tests réalistes tirés de la nanophotonique. / The goal of this thesis is to develop a discontinuous Galerkin time-domain method to be able to handle realistic nanophotonics computations. During the last decades, the evolution of lithography techniques allowed the creation of geometrical structures at the nanometer scale, thus unveiling a variety of new phenomena arising from light-matter interactions at such levels. These effects usually occur when the device is of comparable size or (much) smaller than the wavelength of the incident field. This work relies on the development and implementation of appropriate models for dispersive materials (mostly metals), as well as on a large panel of classical computational techniques. Two major methodological developments are presented and studied in details: (i) curvilinear elements, and (ii) local order of approximation. This work is complemented with several physical studies of real-life nanophotonics applications.

La conduction thermique devient moins efficace à mesure que la taille des struc-tures diminuent en desous du micron, car la diffusion de surface des phononsdevient prédominante et limite plus efficacement les phonons que la diffusionphonon-phonon Umklapp. Des études récentes ont indiqué que les phonon po-laritons de surface (SPhPs), qui sont les ondes électromagnétiques évanescentesgénérées par l’hybridation des phonons optiques et des photons et se propageantà la surface d’une surface diélectrique polaire, pourraient servir de nouveauxvecteurs de chaleur pour améliorer les performances thermiques dans des dis-positifs micro- et nano-métriques. Nous étudions l’état des SPhPs existantdans un film submicronique diélectrique dans une large gamme de fréquences.Le calcul de la conductivité thermique des SPhPs basé sur l’équation de trans-port de Boltzmann (BTE) montre que le flux de chaleur transporté par lesSPhPs est supérieur à celui des phonons. Nous effectuons également une mesurede réflectance thermique dans le domaine temporel (TDTR) de films submi-croniques deSiNet démontrons que la conductivité thermique due aux SPhPsà haute température augmente lorsque l’épaisseur du film dimine. Les résultatsprésentés dans cette thèse ont des applications potentielles dans le domaine dutransfert de chaleur, de la gestion thermique, du rayonnement en champ proche et de la polaritoniques. / Thermal conduction becomes less efficient as structures scale down into submicron sizes since phonon-boundary scattering becomes predominant and impede phonons more efficiently than Umklapp scattering. Recent studies indicated that the surface phonon polaritons (SPhPs), which are the evanescent electromagnetic waves generated by the hybridation of the optical phonons and the photons and propagating at the surface of a polar dielectric material surface, potentially serve as novel heat carriers to enhance the thermal performance in micro- and nanoscale devices. We study the condition of SPhPs existing in a dielectric submicron film with a broad frequency range. The calculaton of SPhPs thermal conductivity based on Boltzmann transport equation (BTE) demonstrates that the heat flux carried by SPhPs exceeds the one carried by phonons. We also conduct a time-domain-thermal-reflectance (TDTR) measurement of $SiN$ submicron films and demonstrate that the thermal conductivity due to the SPhPs at high temperatures increases by decreasing the film thickness. The results presented in this thesis have potential applications in the field of heat transfer, thermal management, near-field radiation and polaritonics.