Resonant properties, of modified triangular plasmonic nanoparticles with higher field concentration

In this paper, we present an analysis of the resonant response of modified triangular metallic nanoparticles with polynomial sides. The particles are illuminated by an incident plane wave and the method of moments is used to solve numerically the electromagnetic scattering problem. We investigate spectral response and near field distribution in function of the length and polynomial order of the nanoparticles. Our results show that in the analyzed wavelength range (0.5-1.8) µm these particles possess smaller number of resonances and their resonant wavelengths, near field enhancement and field confinement are higher than those of the conventional triangular particle with linear sides.

Analysis of modified bowtie nanoantennas in the excitation and emission regimes

In this work, we analyze modified bowtie nanoantennas with polynomial sides in the excitation and emission regimes. In the excitation regime, the antennas are illuminated by an incident plane wave, and in the emission regime, the excitation is fulfilled by infinitesimal electric dipole positioned in the gap of the nanoantennas. Several antennas with different sizes and polynomial order were numerically analyzed by method of moments. The results show that these novel antennas possess a controllable resonance by the polynomial order and good characteristics of near field enhancement and confinement for applications in enhancement of spontaneous emission of a single molecule.

Lensometry by two-laser holography with photorefractive Bi12TiO20

Refractive and profilometric measurements of lenses were performed through holography with a photorefractive Bi12TiO20 crystal as the recording medium. Two properly aligned diode lasers emitting in the red region were employed as light sources. Both lasers were tuned in order to provide millimetric and sub-millimetric synthetic wavelengths. The surfaces of the test lens were covered by a 25-μm opaque plastic tape in order to allow the lens profilometry upon illuminating them with a collimated beam. The resulting holographic images appear covered by interference fringes corresponding to the wavefront geometry of the wave scattered by the lens. For refractive index measurement a diffusely scattering flat surface was positioned behind the uncovered lens which was also illuminated by a plane wave. The resulting contour interferogram describes the form of the wavefront after the beam traveled back and forth through the lens. The fringe quantitative evaluation was carried out through the four-stepping technique and the resulting phase map and the Branch-cut method was employed for phase unwrapping. The only non-optical procedure for lens characterization was the thickness measurement, made by a dial caliper. Exact ray tracing calculation was performed in order to establish a relation between the output wavefront geometry and the lens parameters like radii of curvature, thickness and refractive index. By quantitatively comparing the theoretical wavefront geometry with the experimental results relative uncertainties bellow 3% for refractive index and 1 % for focal length were obtained. © 2008 American Institute of Physics.

Site localization of Cd impurities in sapphire

By combining first-principles electronic structure calculations and existing time-differential gamma-gamma perturbed-angular-correlation experiments we studied the site localization, the local environment, and the electronic structure of Cd impurities in sapphire (alpha-Al2O3) single crystals in different charged states. The ab initio calculations were performed with the full-potential augmented plane wave plus local orbitals method and the projector augmented wave method. Comparing the calculated electric-field-gradient tensor at the Cd nuclei in the alpha-Al2O3 host lattice and the corresponding available experimental values, we have seen that it is equally possible for Cd to replace an Al atom (in a negative charge state) or to be placed in an interstitial site (in a neutral charge state). To finally address the issue of the Cd impurity localization, we performed formation energy calculations. These results have shown that Cd placed in the substitutional Al site, in the negatively charged state, is the most probable configuration.

Study of the vertical transport in p-doped superlattices based on group III-V semiconductors

The electrical conductivity σ has been calculated for p-doped GaAs/Al0.3Ga0.7As and cubic GaN/Al0.3Ga0.7N thin superlattices (SLs). The calculations are done within a self-consistent approach to the k → ⋅ p → theory by means of a full six-band Luttinger-Kohn Hamiltonian, together with the Poisson equation in a plane wave representation, including exchange correlation effects within the local density approximation. It was also assumed that transport in the SL occurs through extended minibands states for each carrier, and the conductivity is calculated at zero temperature and in low-field ohmic limits by the quasi-chemical Boltzmann kinetic equation. It was shown that the particular minibands structure of the p-doped SLs leads to a plateau-like behavior in the conductivity as a function of the donor concentration and/or the Fermi level energy. In addition, it is shown that the Coulomb and exchange-correlation effects play an important role in these systems, since they determine the bending potential.

First principles calculations for vacancies and antisites in PbSe and PbTe: bulk and nanowire.

The energetic stability and the electronic properties of vacancies (VX) and antisites (XY) in PbSe and PbTe are investigated. PbSe and PbTe are narrow band gap semiconductors and have the potential to be used in infrared detectors, laser, and diodes. They are also of special interest for thermoelectric devices (TE). The calculations are based in the Density Functional Theory (DFT) and the General Gradient Approximation (GGA) for the exchange-correlation term, as implemented in the VASP code. The core and valence electrons are described by the Projected Augmented Wave (PAW) and the Plane Wave (PW) methods, respectively. The defects are studied in the bulk and nanowire (NW) system. Our results show that intrinsec defects (vacancies and antisites) in PbTe have lower formation energies in the NW as compared to the bulk and present a trend in migrate to the surface of the NW. For the PbSe we obtain similar results when compare the formation energy for the bulk and NW. However, the Pb vacancy and the antisites are more stable in the core of the NW. The intrinsec defects are shallow defects for the bulk system. For both PbSe and PbTe VPb is a shallow acceptor defect and VSe and VT e are shallow donor defects for the PbSe and PbTe, respectively. Similar electronic properties are observed for the antisites. For the Pb in the anion site we obtain an n-type semiconductor for both PbSe and PbTe, SeP b is a p-type for the PbSe, and T eP b is a n-type for PbTe. Due the quantum con¯nement effects present in the NW (the band gap open), these defects have different electronic properties for the NW as compared to the bulk. Now these defects give rise to electronic levels in the band gap of the PbTe NW and the VT e present a metallic character. For the PbSe NW a p-type and a n-type semiconductor is obtained for the VP b and P bSe, respectively. On the other hand, deep electronic levels are present in the band gap for the VSe and SePb. These results show that due an enhanced in the electronic density of states (DOS) near the Fermi energy, the defective PbSe and PbTe are candidates for efficient TE devices.

Untersuchung der Kernstruktur von 3 He mittels Polarisationsobservablen

Durch die Möglichkeit, gleichzeitig mehrere Polarisationsfreiheitsgradernin der quasi-elastischen Elektronstreuung an $^3mathrm{He}$ zurnmessen, bietet sich ein neuer experimenteller Zugang zu kleinen, aber rnwichtigen Partialwellenbeiträgen ($S'$, $D$-Welle) desrn$^3mathrm{He}$-Grundzustands. Dies ermöglicht nicht nur ein tieferesrnVerständnis des Drei-Körper-Systems, sondern bietet auch diernMöglichkeit, Erkenntnisse über die $^3mathrm{He}$-Struktur undrnDynamik zu erlangen. Mit Hilfe dieser Informationen lassen sich abrninitio Rechnungen testen, sowie Korrekturen berechnen, die für anderernExperimente (z.B. Messung von $G_{en}$) benötigt werden. rnrnModerne Faddeev-Rechnungen liefern nicht nur eine quantitativernBeschreibung des $^3mathrm{He}$-Grundzustands, sondern geben auchrneinen Einblick in die sogenannten spinabhängigenrnImpulsverteilungen. Eine gründliche experimentelle Untersuchung ist in rndiesem Zusammenhang nötig, um eine solide Basis für die Üperprüfungrnder theoretische Modelle zu liefern. EinrnDreifach-Polarisationsexperiment liefert hier zum einen wichtigernDaten, zum anderen kann damit untersucht werden, ob mit der Methoderndes glqq Deuteron-Tagginggrqq polarisiertes $^3mathrm{He}$ alsrneffektives polarisiertes Protonentarget verwendet werden kann. rnrnDas hier vorgestellte Experiment kombiniert erstmals Strahl- undrnTargetpolarisation sowie die Messung der Polarisation des auslaufendenrnProtons. Das Experiment wurde im Sommer 2007 an derrnDrei-Spektrometer-Anlage der A1-Kollaboration an MAMI rndurchgeführt. Dabei wurde mit einer Strahlenergie vonrn$E=855,mathrm{MeV}$ bei $q^2=-0,14,(mathrm{GeV/c})^2$rn$(omega=0,13,mathrm{GeV}$, $q=0,4,mathrm{GeV/c})$ gemessen.rnrnDie bestimmten Wirkungsquerschnitte, sowie die Strahl-Target- und diernDreifach-Asymmetrie werden mit theoretischen Modellrechnungen vonrnJ. Golak (Plane Wave Impuls Approximation PWIA, sowie ein Modell mitrnEndzustandswechselwirkung) verglichen. Zudem wurde das Modell von dernForest verwendet, welches den Wirkungsquerschnitt über eine gemessenernSpektralfunktion berechnet. Der Vergleich mit den Modellrechnungenrnzeigt, dass sowohl der Wirkungsquerschnitt, als auch die Doppel- undrnDreifach-Asymmetrie gut mit den theoretischen Rechnungenrnübereinstimmen. rnrnDie Ergebnisse dieser Arbeit bestätigen, dass polarisiertesrn$^3mathrm{He}$ nicht nur als polarisiertes Neutronentarget, sondernrndurch Nachweis des Deuterons ebenfalls als polarisiertesrnProtonentarget verwendet werden kann.

Development of a high performance parallel computing platform and its use in the study of nanostructures : clusters, sheets and tubes .

Small clusters of gallium oxide, technologically important high temperature ceramic, together with interaction of nucleic acid bases with graphene and small-diameter carbon nanotube are focus of first principles calculations in this work. A high performance parallel computing platform is also developed to perform these calculations at Michigan Tech. First principles calculations are based on density functional theory employing either local density or gradient-corrected approximation together with plane wave and gaussian basis sets. The bulk Ga2O3 is known to be a very good candidate for fabricating electronic devices that operate at high temperatures. To explore the properties of Ga2O3 at nonoscale, we have performed a systematic theoretical study on the small polyatomic gallium oxide clusters. The calculated results find that all lowest energy isomers of GamOn clusters are dominated by the Ga-O bonds over the metal-metal or the oxygen-oxygen bonds. Analysis of atomic charges suggest the clusters to be highly ionic similar to the case of bulk Ga2O3. In the study of sequential oxidation of these slusters starting from Ga2O, it is found that the most stable isomers display up to four different backbones of constituent atoms. Furthermore, the predicted configuration of the ground state of Ga2O is recently confirmed by the experimental result of Neumark's group. Guided by the results of calculations the study of gallium oxide clusters, performance related challenge of computational simulations, of producing high performance computers/platforms, has been addressed. Several engineering aspects were thoroughly studied during the design, development and implementation of the high performance parallel computing platform, rama, at Michigan Tech. In an attempt to stay true to the principles of Beowulf revolutioni, the rama cluster was extensively customized to make it easy to understand, and use - for administrators as well as end-users. Following the results of benchmark calculations and to keep up with the complexity of systems under study, rama has been expanded to a total of sixty four processors. Interest in the non-covalent intereaction of DNA with carbon nanotubes has steadily increased during past several years. This hybrid system, at the junction of the biological regime and the nanomaterials world, possesses features which make it very attractive for a wide range of applicatioins. Using the in-house computational power available, we have studied details of the interaction between nucleic acid bases with graphene sheet as well as high-curvature small-diameter carbon nanotube. The calculated trend in the binding energies strongly suggests that the polarizability of the base molecules determines the interaction strength of the nucleic acid bases with graphene. When comparing the results obtained here for physisorption on the small diameter nanotube considered with those from the study on graphene, it is observed that the interaction strength of nucleic acid bases is smaller for the tube. Thus, these results show that the effect of introducing curvature is to reduce the binding energy. The binding energies for the two extreme cases of negligible curvature (i.e. flat graphene sheet) and of very high curvature (i.e. small diameter nanotube) may be considered as upper and lower bounds. This finding represents an important step towards a better understanding of experimentally observed sequence-dependent interaction of DNA with Carbon nanotubes.

Reduction of truncation errors in planar, cylindrical and partial spherical near-field antenna measurements

A method to reduce truncation errors in near-field antenna measurements is presented. The method is based on the Gerchberg-Papoulis iterative algorithm used to extrapolate band-limited functions and it is able to extend the valid region of the calculated far-field pattern up to the whole forward hemisphere. The extension of the valid region is achieved by the iterative application of a transformation between two different domains. After each transformation, a filtering process that is based on known information at each domain is applied. The first domain is the spectral domain in which the plane wave spectrum (PWS) is reliable only within a known region. The second domain is the field distribution over the antenna under test (AUT) plane in which the desired field is assumed to be concentrated on the antenna aperture. The method can be applied to any scanning geometry, but in this paper, only the planar, cylindrical, and partial spherical near-field measurements are considered. Several simulation and measurement examples are presented to verify the effectiveness of the method.

Performance Comparison Between Serrated Edge and Rolled Edge Reflectors Inside CATR Facilities

CATR facilities are attractive antenna measurement facilities. Main reasons which contribute to this fact lie on its inherent reduced volume, on-the-fly measurements and the extension of both to a wide range of frequencies. However, these features rely on the assumption that the field collimation scheme is able to generate a plane wave distribution (quiet zone) where the AUT is to be placed and operated in RX mode. Unfortunately, electromagnetic theory states that this field distribution is not possible to be generated by a finite size scatterer operated as the collimator of a nonzero wavelength time-harmonic propagating field. This is the background of this paper, where two well-known electromagnetic field collimators will be discussed: the serrated edge reflector and the blended rolled edge reflector. To reach this purpose, electromagnetic hybrid analysis techniques developed at Technical University of Madrid will be applied.

Monitoring maximum peak power density values as a more reliable marker for measurement procedures to elucidate RF EMF health effects: Leganes project

Considering the measurement procedures recommended by the ICNIRP, this communication is a proposal for a measurement procedure based in the maximum peak values of equivalent plane wave power density. This procedure has been included in a project being developed in Leganés, Spain. The project plans to deploy a real time monitoring system for RF to provide the city with a useful tool to adapt the environmental EM conditions to the new regulations approved. A first stage consisting of 105 measurement points has been finished and all the values are under the threshold of the regulation.

Development of a Uniform Theory of Diffraction for Scattered and Surface Electromagnetic Field Analysis on a Cylinder

Con esta tesis ”Desarrollo de una Teoría Uniforme de la Difracción para el Análisis de los Campos Electromagnéticos Dispersados y Superficiales sobre un Cilindro” hemos iniciado una nueva línea de investigación que trata de responder a la siguiente pregunta: ¿cuál es la impedancia de superficie que describe una estructura de conductor eléctrico perfecto (PEC) convexa recubierta por un material no conductor? Este tipo de estudios tienen interés hoy en día porque ayudan a predecir el campo electromagnético incidente, radiado o que se propaga sobre estructuras metálicas y localmente convexas que se encuentran recubiertas de algún material dieléctrico, o sobre estructuras metálicas con pérdidas, como por ejemplo se necesita en determinadas aplicaciones aeroespaciales, marítimas o automovilísticas. Además, desde un punto de vista teórico, la caracterización de la impedancia de superficie de una estructura PEC recubierta o no por un dieléctrico es una generalización de varias soluciones que tratan ambos tipos de problemas por separado. En esta tesis se desarrolla una teoría uniforme de la difracción (UTD) para analizar el problema canónico del campo electromagnético dispersado y superficial en un cilindro circular eléctricamente grande con una condición de contorno de impedancia (IBC) para frecuencias altas. Construir una solución basada en UTD para este problema canónico es crucial en el desarrollo de un método UTD para el caso más general de una superficie arbitrariamente convexa, mediante el uso del principio de localización de los campos electromagnéticos a altas frecuencias. Esta tesis doctoral se ha llevado a cabo a través de una serie de hitos que se enumeran a continuación, enfatizando las contribuciones a las que ha dado lugar. Inicialmente se realiza una revisión en profundidad del estado del arte de los métodos asintóticos con numerosas referencias. As í, cualquier lector novel puede llegar a conocer la historia de la óptica geométrica (GO) y la teoría geométrica de la difracción (GTD), que dieron lugar al desarrollo de la UTD. Después, se investiga ampliamente la UTD y los trabajos más importantes que pueden encontrarse en la literatura. As í, este capítulo, nos coloca en la posición de afirmar que, hasta donde nosotros conocemos, nadie ha intentado antes llevar a cabo una investigación rigurosa sobre la caracterización de la impedancia de superficie de una estructura PEC recubierta por un material dieléctrico, utilizando para ello la UTD. Primero, se desarrolla una UTD para el problema canónico de la dispersión electromagnética de un cilindro circular eléctricamente grande con una IBC uniforme, cuando es iluminado por una onda plana con incidencia oblicua a frecuencias altas. La solución a este problema canónico se construye a partir de una solución exacta mediante una expansión de autofunciones de propagación radial. Entonces, ésta se convierte en una nueva expansión de autofunciones de propagación circunferencial muy apropiada para cilindros grandes, a través de la transformación de Watson. De esta forma, la expresión del campo se reduce a una integral que se evalúa asintóticamente, para altas frecuencias, de manera uniforme. El resultado se expresa según el trazado de rayos descrito en la UTD. La solución es uniforme porque tiene la importante propiedad de mantenerse continua a lo largo de la región de transición, a ambos lados de la superficie del contorno de sombra. Fuera de la región de transición la solución se reduce al campo incidente y reflejado puramente ópticos en la región iluminada del cilindro, y al campo superficial difractado en la región de sombra. Debido a la IBC el campo dispersado contiene una componente contrapolar a causa de un acoplamiento entre las ondas TEz y TMz (donde z es el eje del cilindro). Esta componente contrapolar desaparece cuando la incidencia es normal al cilindro, y también en la región iluminada cuando la incidencia es oblicua donde el campo se reduce a la solución de GO. La solución UTD presenta una muy buena exactitud cuando se compara numéricamente con una solución de referencia exacta. A continuación, se desarrolla una IBC efectiva para el cálculo del campo electromagnético dispersado en un cilindro circular PEC recubierto por un dieléctrico e iluminado por una onda plana incidiendo oblicuamente. Para ello se derivan dos impedancias de superficie en relación directa con las ondas creeping y de superficie TM y TE que se excitan en un cilindro recubierto por un material no conductor. Las impedancias de superficie TM y TE están acopladas cuando la incidencia es oblicua, y dependen de la geometría del problema y de los números de onda. Además, se ha derivado una impedancia de superficie constante, aunque con diferente valor cuando el observador se encuentra en la zona iluminada o en la zona de sombra. Después, se presenta una solución UTD para el cálculo de la dispersión de una onda plana con incidencia oblicua sobre un cilindro eléctricamente grande y convexo, mediante la generalización del problema canónico correspondiente al cilindro circular. La solución asintótica es uniforme porque se mantiene continua a lo largo de la región de transición, en las inmediaciones del contorno de sombra, y se reduce a la solución de rayos ópticos en la zona iluminada y a la contribución de las ondas de superficie dentro de la zona de sombra, lejos de la región de transición. Cuando se usa cualquier material no conductor se excita una componente contrapolar que tiende a desaparecer cuando la incidencia es normal al cilindro y en la región iluminada. Se discuten ampliamente las limitaciones de las fórmulas para la impedancia de superficie efectiva, y se compara la solución UTD con otras soluciones de referencia, donde se observa una muy buena concordancia. Y en tercer lugar, se presenta una aproximación para una impedancia de superficie efectiva para el cálculo de los campos superficiales en un cilindro circular conductor recubierto por un dieléctrico. Se discuten las principales diferencias que existen entre un cilindro PEC recubierto por un dieléctrico desde un punto de vista riguroso y un cilindro con una IBC. Mientras para un cilindro de impedancia se considera una impedancia de superficie constante o uniforme, para un cilindro conductor recubierto por un dieléctrico se derivan dos impedancias de superficie. Estas impedancias de superficie están asociadas a los modos de ondas creeping TM y TE excitadas en un cilindro, y dependen de la posición y de la orientación del observador y de la fuente. Con esto en mente, se deriva una solución UTD con IBC para los campos superficiales teniendo en cuenta las dependencias de la impedancia de superficie. La expansión asintótica se realiza, mediante la transformación de Watson, sobre la representación en serie de las funciones de Green correspondientes, evitando as í calcular las derivadas de orden superior de las integrales de tipo Fock, y dando lugar a una solución rápida y precisa. En los ejemplos numéricos realizados se observa una muy buena precisión cuando el cilindro y la separación entre el observador y la fuente son grandes. Esta solución, junto con el método de los momentos (MoM), se puede aplicar para el cálculo eficiente del acoplamiento mutuo de grandes arrays conformados de antenas de parches. Los métodos propuestos basados en UTD para el cálculo del campo electromagnético dispersado y superficial sobre un cilindro PEC recubierto de dieléctrico con una IBC efectiva suponen un primer paso hacia la generalización de una solución UTD para superficies metálicas convexas arbitrarias cubiertas por un material no conductor e iluminadas por una fuente electromagnética arbitraria. ABSTRACT With this thesis ”Development of a Uniform Theory of Diffraction for Scattered and Surface Electromagnetic Field Analysis on a Cylinder” we have initiated a line of investigation whose goal is to answer the following question: what is the surface impedance which describes a perfect electric conductor (PEC) convex structure covered by a material coating? These studies are of current and future interest for predicting the electromagnetic (EM) fields incident, radiating or propagating on locally smooth convex parts of highly metallic structures with a material coating, or by a lossy metallic surfaces, as for example in aerospace, maritime and automotive applications. Moreover, from a theoretical point of view, the surface impedance characterization of PEC surfaces with or without a material coating represents a generalization of independent solutions for both type of problems. A uniform geometrical theory of diffraction (UTD) is developed in this thesis for analyzing the canonical problem of EM scattered and surface field by an electrically large circular cylinder with an impedance boundary condition (IBC) in the high frequency regime, by means of a surface impedance characterization. The construction of a UTD solution for this canonical problem is crucial for the development of the corresponding UTD solution for the more general case of an arbitrary smooth convex surface, via the principle of the localization of high frequency EM fields. The development of the present doctoral thesis has been carried out through a series of landmarks that are enumerated as follows, emphasizing the main contributions that this work has given rise to. Initially, a profound revision is made in the state of art of asymptotic methods where numerous references are given. Thus, any reader may know the history of geometrical optics (GO) and geometrical theory of diffraction (GTD), which led to the development of UTD. Then, the UTD is deeply investigated and the main studies which are found in the literature are shown. This chapter situates us in the position to state that, as far as we know, nobody has attempted before to perform a rigorous research about the surface impedance characterization for material-coated PEC convex structures via UTD. First, a UTD solution is developed for the canonical problem of the EM scattering by an electrically large circular cylinder with a uniform IBC, when it is illuminated by an obliquely incident high frequency plane wave. A solution to this canonical problem is first constructed in terms of an exact formulation involving a radially propagating eigenfunction expansion. The latter is converted into a circumferentially propagating eigenfunction expansion suited for large cylinders, via the Watson transformation, which is expressed as an integral that is subsequently evaluated asymptotically, for high frequencies, in a uniform manner. The resulting solution is then expressed in the desired UTD ray form. This solution is uniform in the sense that it has the important property that it remains continuous across the transition region on either side of the surface shadow boundary. Outside the shadow boundary transition region it recovers the purely ray optical incident and reflected ray fields on the deep lit side of the shadow boundary and to the modal surface diffracted ray fields on the deep shadow side. The scattered field is seen to have a cross-polarized component due to the coupling between the TEz and TMz waves (where z is the cylinder axis) resulting from the IBC. Such cross-polarization vanishes for normal incidence on the cylinder, and also in the deep lit region for oblique incidence where it properly reduces to the GO or ray optical solution. This UTD solution is shown to be very accurate by a numerical comparison with an exact reference solution. Then, an effective IBC is developed for the EM scattered field on a coated PEC circular cylinder illuminated by an obliquely incident plane wave. Two surface impedances are derived in a direct relation with the TM and TE surface and creeping wave modes excited on a coated cylinder. The TM and TE surface impedances are coupled at oblique incidence, and depend on the geometry of the problem and the wave numbers. Nevertheless, a constant surface impedance is found, although with a different value when the observation point lays in the lit or in the shadow region. Then, a UTD solution for the scattering of an obliquely incident plane wave on an electrically large smooth convex coated PEC cylinder is introduced, via a generalization of the canonical circular cylinder problem. The asymptotic solution is uniform because it remains continuous across the transition region, in the vicinity of the shadow boundary, and it recovers the ray optical solution in the deep lit region and the creeping wave formulation within the deep shadow region. When a coating is present a cross-polar field term is excited, which vanishes at normal incidence and in the deep lit region. The limitations of the effective surface impedance formulas are discussed, and the UTD solution is compared with some reference solutions where a very good agreement is met. And in third place, an effective surface impedance approach is introduced for determining surface fields on an electrically large coated metallic circular cylinder. Differences in analysis of rigorouslytreated coated metallic cylinders and cylinders with an IBC are discussed. While for the impedance cylinder case a single constant or uniform surface impedance is considered, for the coated metallic cylinder case two surface impedances are derived. These are associated with the TM and TE creeping wave modes excited on a cylinder and depend on observation and source positions and orientations. With this in mind, a UTD based method with IBC is derived for the surface fields by taking into account the surface impedance variation. The asymptotic expansion is performed, via the Watson transformation, over the appropriate series representation of the Green’s functions, thus avoiding higher-order derivatives of Fock-type integrals, and yielding a fast and an accurate solution. Numerical examples reveal a very good accuracy for large cylinders when the separation between the observation and the source point is large. Thus, this solution could be efficiently applied in mutual coupling analysis, along with the method of moments (MoM), of large conformal microstrip array antennas. The proposed UTD methods for scattered and surface EM field analysis on a coated PEC cylinder with an effective IBC are considered the first steps toward the generalization of a UTD solution for large arbitrarily convex smooth metallic surfaces covered by a material coating and illuminated by an arbitrary EM source.

Modelización de las propiedades radiativas e interacción de láseres ultraintensos con la materia

Desde el año 2004 el código ARWEN ha sido utilizado con éxito para simular y diseñar experimentos relacionados con blancos para fusión por confinamiento inercial [146], astrofísica de laboratorio [145], plasmas como amplificadores de láseres de rayos X [107] o plasmas creados por láser para la medición de espectros de transmisión. Para la realización de estas simulaciones es necesario, además de métodos de alto orden precisos y que presenten buenas propiedades conservativas, conocer ciertas propiedades de los plasmas. En el caso de la fluidodinámica y la conducción electrónica necesitaremos conocer la ecuación de estado [94, 49, 36], y para el transporte de la radiación será preciso disponer de los datos de absorción y emisión [104, 95, 40]. Hasta el año 2009 ARWEN dependía de códigos externos para la generación de estas tablas de opacidad, careciendo de control sobre los métodos empleados para su generación. Además estos códigos asumían equilibrio local termodinámico (LTE), limitando su validez a rangos de alta densidad y baja temperatura. En el marco de esta tesis se ha desarrollado el código BIGBART para la generación de tablas detalladas de opacidad y emisividad para su uso en el módulo de transporte de radiación. De esta forma el grupo dispondrá de su propia herramienta de generación de propiedades radiativas. El código desarrollado es capaz de tratar plasmas en estado fuera de equilibrio (non-LTE) mediante el modelo colisional-radiativo, extendiendo así el rango de validez de las tablas generadas. El trabajo desarrollado para implementar un código LTE/non-LTE estacionario es el siguiente Cálculo de estructura y datos atómicos. Se ha acoplado en código FAC a BIGBART, incorporando la capacidad para generar potenciales atómicos para una configuración y el cálculo de funciones de onda de electrones en orbitales ligados y libres. Aproximaciones y métodos para la obtención de tasas y secciones eficaces de procesos. Se han incluido y programado los modelos implementados en FAC para el cálculo de secciones eficaces de fotoionización, y tasas de decaimiento de emisión espontánea y autoionización. Además se ha incluido el modelo Plane-Wave Born (PWBA) para el cálculo de las secciones eficaces de ionización y excitación colisional. Modelos para la obtención de la distribución de estados iónicos dentro del plasma. Se ha programado un solver LTE basado en la ecuación de Saha-Boltzmann con efectos de ionización por presión debida a los iones adyacentes. También se ha implementado un modelo non-LTE colisionalradiativo para la resolución del sistema de ecuaciones que nos permite obtener la densidad de estados iónicos fuera de equilibrio. Modelo non-LTE RADIOM. Se ha implementado el modelo RADIOM para aproximar efectos de no-equilibrio mediante cálculos LTE a una temperatura equivalente, menor o igual que la temperatura electrónica real. Cálculo de las propiedades espectrales de absorción y emisión. Se han implementado los modelos para el cálculo de los perfiles espectrales de absorción y emisión para procesos entre niveles ligados, ligado-libre y librelibre. Aprovechando el trabajo realizado en este sentido, durante el transcurso de esta tesis se amplió el código BIGBART para tratar problemas con dependencia temporal. La extensión para tratar este tipo de problemas se orientó a la simulación numérica de la interacción de láseres ultra intensos en el rango XUV/rayos X. Para ello, además de adaptar el modelo non-LTE colisionalradiativo se incluyeron procesos adicionales asociados a la interacción de la materia con fotones altamente energéticos. También se han incluido modelos para el cálculo de las propiedades ópticas, y por ende las propiedades dieléctricas de la materia irradiada, de gran interés en algunas aplicaciones novedosas de estos láseres intensos. Debido a la naturaleza fuertemente fuera de equilibrio en la interacción de fotones de alta energía con la materia, se incluyó el tratamiento de la distribución de electrones libres fuera de equilibrio en la aproximación de Fokker-Planck, tanto para condiciones degeneradas como no degeneradas. El trabajo desarrollado en el código non-LTE con dependencia temporal es el siguiente Procesos asociados a láseres intensos XUV/rayos X. Se ha implementado el cálculo de procesos radiativos estimulados de absorción y emisión por el láser. También se han incluido procesos asociados a la creación de vacantes en capas internas electrónicas (Shake), además de doble autoionización y doble fotoionización. Cálculo de propiedades ópticas y dieléctricas en blancos sólidos. Se ha implementado un modelo para la absorción por bremsstrahlung inverso en blancos en estado sólido. Con el coeficiente de extinción debido a procesos de fotoabsorción resonante, fotoionización y bremsstrahlung inverso se obtiene el ´ındice de refracción mediante la relación de Kronig-Kramers. Electrones fuera de equilibrio. Se ha tratado la evolución de la distribución de electrones, cuando no está justificado asumir que es Maxwelliana o de Fermi-Dirac, mediante la aproximación de Fokker-Planck para la colisión entre electrones libres. En la resolución de la ecuación de Fokker-Planck se han incluido los procesos inelásticos por colisiones con iones y términos fuente por interacción con el láser y otros procesos. ABSTRACT Since 2004 the ARWEN code has been successfully used to simulate and design targets for inertial confinement fusion experiments [146], laboratory astrophysics [145], plasmas as X-ray lasers amplifiers [107] or laser created plasmas for measuring transmission spectra. To perform these simulations it is necessary, in addition to high order precise methods with good conservative properties, to know certain properties of plasmas. For fluid dynamic and electronic conduction we need to know the equation of state [94, 49, 36], and for radiation transport it will be necessary to have the data of the absorption and emission [104, 95, 40]. Until 2009 ARWEN depended on external codes to generate these opacity tables, lacking of control over the methods used for their generation. Besides, these codes assumed local thermodynamic equilibrium (LTE), limiting their validity ranges to high densities and low temperatures. As part of this thesis it has been developed the BIGBART code for generating detailed opacity and emissivity tables for use in the radiation transport module. This group will have its own tool for the generation of radiative properties. The developed code is capable of treating plasmas out of equilibrium (non-LTE) by means of a collisional-radiative model, extending the range of validity of the generated tables. The work to implement an LTE/non-LTE steady-state code is as follows Calculation of structure and atomic data. the FAC code was coupled to BIGBART, incorporating the ability to generate atomic potentials for calculating configuration wave functions for bound and free electrons. Approaches and methods for obtaining cross sections and processes rates. We have included and reprogrammed in Fortran the models implemented in FAC for calculation of photoionization cross sections and decay rates of spontaneous emission and autoionization. We also included the Plane- Wave Born (PWBA) model to calculate the cross sections of ionization and collisional excitation. Models for the obtention of the distribution of ionic states within the plasma. We programmed a LTE solver based on the Saha-Boltzmann equation with pressure ionization effects due to adjacent ions. It has also been implemented a non-LTE collisional-radiative model for solving the system of equations that allows us to obtain the density of ionic states out of equilibrium. Non-LTE RADIOM model. We have implemented the non-LTE RADIOM model to approximate non-equilibrium effects with LTE data at an equivalent temperature, lower or equal to the actual electronic temperature. Calculation of the spectral absorption and emission properties. Models have been implemented for the calculation of the spectral profiles of absorption and emission processes between bound levels, free-bound and free-free. Taking advantage of the work done in this direction throughout the course of this thesis the code BIGBART was extended to treat time-dependent problems. The extension to treat such problems is oriented to the numerical simulation of the interaction of ultra intense lasers in the XUV/X-ray range. For this range, in addition to adapting the non-LTE collisional-radiative model, additional processes associated with the interaction of matter with high energy photons. We also included models for calculation of the optical properties, and therefore the dielectric properties of the irradiated material, of great interest in some novel applications of these intense lasers. Due to the strong non-equilibrium nature of the interaction of high energy photons with matter, we included the treatment of the distribution of free electrons out of equilibrium in the Fokker-Planck approximation for both degenerate and non-degenerate conditions. The work in the non-LTE time-dependent code is as follows Processes associated with intense XUV/X-ray lasers. We have implemented the calculation of stimulated radiative processes in absorption and emission. Also we included processes associated with the creation of electronic vacancies in inner shells (Shake), double autoionization and double photoionization. Calculation of optical and dielectric properties in solid targets. We have implemented a model for inverse bremsstrahlung absorption in solid targets. With the extinction coefficient from resonant photoabsorption, photoionization and inverse bremsstrahlung the refractive index is obtained by the Kramers-Kronig relation. Electrons out of equilibrium. We treat the evolution of the electron distribution, when it is not justified to assume a Maxwellian or Fermi-Dirac distribution, by the Fokker-Planck approximation for collisions between electrons. When solving the Fokker-Planck equation we included inelastic collision processes with ions and source terms by interaction with the laser and other processes.

A source reconstruction technique for planar arrays of wide slots

A method to reconstruct the excitation coefficients of wide-slot arrays from near-field data is presented. The plane wave spectrum (PWS) is used for reconstruction, and the shape of the field distribution on a wide slot is considered in the calculation of the PWS. The proposed algorithm is validated through the reconstruction of the excitation coefficients of a wide-slot array with element failures from the simulated nearfield data. The element failures are clearly located by the proposed algorithm

Multipole expansion of strongly focussed laser beams

Multipole expansion of an incident radiation field-that is, representation of the fields as sums of vector spherical wavefunctions-is essential for theoretical light scattering methods such as the T-matrix method and generalised Lorenz-Mie theory (GLMT). In general, it is theoretically straightforward to find a vector spherical wavefunction representation of an arbitrary radiation field. For example, a simple formula results in the useful case of an incident plane wave. Laser beams present some difficulties. These problems are not a result of any deficiency in the basic process of spherical wavefunction expansion, but are due to the fact that laser beams, in their standard representations, are not radiation fields, but only approximations of radiation fields. This results from the standard laser beam representations being solutions to the paraxial scalar wave equation. We present an efficient method for determining the multipole representation of an arbitrary focussed beam. (C) 2003 Elsevier Science Ltd. All rights reserved.