RF-sputter fabrication of magnetic garnet thin films and simulation modeling for 1-D magnetic photonic crystal waveguide devices

One dimensional magnetic photonic crystals (1D-MPC) are promising structures for integrated optical isolator applications. Rare earth substituted garnet thin films with proper Faraday rotation are required to fabricate planar 1D-MPCs. In this thesis, flat-top response 1D-MPC was proposed and spectral responses and Faraday rotation were modeled. Bismuth substituted iron garnet films were fabricated by RF magnetron sputtering and structures, compositions, birefringence and magnetooptical properties were studied. Double layer structures for single mode propagation were also fabricated by sputtering for the first time. Multilayer stacks with multiple defects (phase shift) composed of Ce-YIG and GGG quarter-wave plates were simulated by the transfer matrix method. The transmission and Faraday rotation characteristics were theoretically studied. It is found that flat-top response, with 100% transmission and near 45o rotation is achievable by adjusting the inter-defect spacing, for film structures as thin as 30 to 35 μm. This is better than 3-fold reduction in length compared to the best Ce-YIG films for comparable rotations, thus allows a considerable reduction in size in manufactured optical isolators. Transmission bands as wide as 7nm were predicted, which is considerable improvement over 2 defects structure. Effect of repetition number and ratio factor on transmission and Faraday rotation ripple factors for the case of 3 and 4 defects structure has been discussed. Diffraction across the structure corresponds to a longer optical path length. Thus the use of guided optics is required to minimize the insertion losses in integrated devices. This part is discussed in chapter 2 in this thesis. Bismuth substituted iron garnet thin films were prepared by RF magnetron sputtering. We investigated or measured the deposition parameters optimization, crystallinity, surface morphologies, composition, magnetic and magnetooptical properties. A very high crystalline quality garnet film with smooth surface has been heteroepitaxially grown on (111) GGG substrate for films less than 1μm. Dual layer structures with two distinct XRD peaks (within a single sputtered film) start to develop when films exceed this thickness. The development of dual layer structure was explained by compositional gradient across film thickness, rather than strain gradient proposed by other authors. Lower DC self bias or higher substrate temperature is found to help to delay the appearance of the 2nd layer. The deposited films show in-plane magnetization, which is advantageous for waveguide devices application. Propagation losses of fabricated waveguides can be decreased by annealing in an oxygen atmosphere from 25dB/cm to 10dB/cm. The Faraday rotation at λ=1.55μm were also measured for the waveguides. FR is small (10° for a 3mm long waveguide), due to the presence of linear birefringence. This part is covered in chapter 4. We also investigated the elimination of linear birefringence by thickness tuning method for our sputtered films. We examined the compressively and tensilely strained films and analyze the photoelastic response of the sputter deposited garnet films. It has been found that the net birefringence can be eliminated under planar compressive strain conditions by sputtering. Bi-layer GGG on garnet thin film yields a reduced birefringence. Temperature control during the sputter deposition of GGG cover layer is critical and strongly influences the magnetization and birefringence level in the waveguide. High temperature deposition lowers the magnetization and increases the linear birefringence in the garnet films. Double layer single mode structures fabricated by sputtering were also studied. The double layer, which shows an in-plane magnetization, has an increased RMS roughness upon upper layer deposition. The single mode characteristic was confirmed by prism coupler measurement. This part is discussed in chapter 5.

Yang-Baxter deformations of Minkowski spacetime

We study Yang-Baxter deformations of 4D Minkowski spacetime. The Yang-Baxter sigma model description was originally developed for principal chiral models based on a modified classical Yang-Baxter equation. It has been extended to coset curved spaces and models based on the usual classical Yang-Baxter equation. On the other hand, for flat space, there is the obvious problem that the standard bilinear form degenerates if we employ the familiar coset Poincaré group/Lorentz group. Instead we consider a slice of AdS5 by embedding the 4D Poincaré group into the 4D conformal group SO(2, 4) . With this procedure we obtain metrics and B-fields as Yang-Baxter deformations which correspond to well-known configurations such as T-duals of Melvin backgrounds, Hashimoto-Sethi and Spradlin-Takayanagi-Volovich backgrounds, the T-dual of Grant space, pp-waves, and T-duals of dS4 and AdS4. Finally we consider a deformation with a classical r-matrix of Drinfeld-Jimbo type and explicitly derive the associated metric and B-field which we conjecture to correspond to a new integrable system.

A computer program for the analysis of the dynamic bending-torsion coupling in bridges using a mini-computer

The analysis of modes and natural frequencies is of primary interest in the computation of the response of bridges. In this article the transfer matrix method is applied to this problem to provide a computer code to calculate the natural frequencies and modes of bridge-like structures. The Fortran computer code is suitable for running on small computers and results are presented for a railway bridge.

Measurements and Prediction of Pedestrian Walking Loads on an Aluminium Catwalk

During the last two decades the topic of human induced vibration has attracted a lot of attention among civil engineering practitioners and academics alike. Usually this type of problem may be encountered in pedestrian footbridges or floors of paperless offices. Slender designs are becoming increasingly popular, and as a consequence, the importance of paying attention to vibration serviceability also increases. This paper resumes the results obtained from measurements taken at different points of an aluminium catwalk which is 6 m in length by 0.6 m in width. Measurements were carried out when subjecting the structure to different actions:1)Static test: a steel cylinder of 35 kg was placed in the middle of the catwalk; 2)Dynamic test: this test consists of exciting the structure with singles impulses; 3)Dynamic test: people walking on the catwalk. Identification of the mechanical properties of the structure is an achievement of the paper. Indirect methods were used to estimate properties including the support stiffness, the beam bending stiffness, the mass of the structure (using Rayleigh method and iterative matrix method), the natural frequency (using the time domain and frequency domain analysis) and the damping ratio (by calculating the logarithmic decrement). Experimental results and numerical predictions for the response of an aluminium catwalk subjected to walking loads have been compared. The damping of this light weight structure depends on the amplitude of vibration which complicates the tuning of a structural model. In the light of the results obtained it seems that the used walking load model is not appropriate as the predicted transient vibration values (TTVs) are much higher than the measured ones.

Análisis, desarrollo y optimización de un sistema para el diseño de redes topográficas valorando el modelo digital del terreno

El presente trabajo tiene como objetivo general el análisis de las técnicas de diseño y optimización de redes topográficas, observadas mediante topografía convencional (no satelital) el desarrollo e implementación de un sistema informático capaz de ayudar a la definición de la geometría más fiable y precisa, en función de la orografía del terreno donde se tenga que ubicar. En primer lugar se realizará un estudio de la metodología del ajuste mediante mínimos cuadrados y la propagación de varianzas, para posteriormente analizar su dependencia de la geometría que adopte la red. Será imprescindible determinar la independencia de la matriz de redundancia (R) de las observaciones y su total dependencia de la geometría, así como la influencia de su diagonal principal (rii), números de redundancia, para garantizar la máxima fiabilidad interna de la misma. También se analizará el comportamiento de los números de redundancia (rii) en el diseño de una red topográfica, la variación de dichos valores en función de la geometría, analizando su independencia respecto de las observaciones así como los diferentes niveles de diseño en función de los parámetros y datos conocidos. Ha de señalarse que la optimización de la red, con arreglo a los criterios expuestos, está sujeta a los condicionantes que impone la necesidad de que los vértices sean accesibles, y además sean visibles entre sí, aquellos relacionados por observaciones, situaciones que dependen esencialmente del relieve del terreno y de los obstáculos naturales o artificiales que puedan existir. Esto implica la necesidad de incluir en el análisis y en el diseño, cuando menos de un modelo digital del terreno (MDT), aunque lo más útil sería la inclusión en el estudio del modelo digital de superficie (MDS), pero esta opción no siempre será posible. Aunque el tratamiento del diseño esté basado en un sistema bidimensional se estudiará la posibilidad de incorporar un modelo digital de superficie (MDS); esto permitirá a la hora de diseñar el emplazamiento de los vértices de la red la viabilidad de las observaciones en función de la orografía y los elementos, tanto naturales como artificiales, que sobre ella estén ubicados. Este sistema proporcionaría, en un principio, un diseño óptimo de una red constreñida, atendiendo a la fiabilidad interna y a la precisión final de sus vértices, teniendo en cuenta la orografía, lo que equivaldría a resolver un planteamiento de diseño en dos dimensiones y media1; siempre y cuando se dispusiera de un modelo digital de superficie o del terreno. Dado que la disponibilidad de obtener de manera libre el MDS de las zonas de interés del proyecto, hoy en día es costoso2, se planteará la posibilidad de conjuntar, para el estudio del diseño de la red, de un modelo digital del terreno. Las actividades a desarrollar en el trabajo de esta tesis se describen en esta memoria y se enmarcan dentro de la investigación para la que se plantean los siguientes objetivos globales: 1. Establecer un modelo matemático del proceso de observación de una red topográfica, atendiendo a todos los factores que intervienen en el mismo y a su influencia sobre las estimaciones de las incógnitas que se obtienen como resultado del ajuste de las observaciones. 2. Desarrollar un sistema que permita optimizar una red topográfica en sus resultados, aplicando técnicas de diseño y simulación sobre el modelo anterior. 3. Presentar una formulación explícita y rigurosa de los parámetros que valoran la fiabilidad de una red topográfica y de sus relaciones con el diseño de la misma. El logro de este objetivo se basa, además de en la búsqueda y revisión de las fuentes, en una intensa labor de unificación de notaciones y de construcción de pasos intermedios en los desarrollos matemáticos. 4. Elaborar una visión conjunta de la influencia del diseño de una red, en los seis siguientes factores (precisiones a posteriori, fiabilidad de las observaciones, naturaleza y viabilidad de las mismas, instrumental y metodología de estacionamiento) como criterios de optimización, con la finalidad de enmarcar el tema concreto que aquí se aborda. 5. Elaborar y programar los algoritmos necesarios para poder desarrollar una aplicación que sea capaz de contemplar las variables planteadas en el apartado anterior en el problema del diseño y simulación de redes topográficas, contemplando el modelo digital de superficie. Podrían considerarse como objetivos secundarios, los siguientes apartados: Desarrollar los algoritmos necesarios para interrelacionar el modelo digital del terreno con los propios del diseño. Implementar en la aplicación informática la posibilidad de variación, por parte del usuario, de los criterios de cobertura de los parámetros (distribución normal o t de Student), así como los grados de fiabilidad de los mismos ABSTRACT The overall purpose of this work is the analysis of the techniques of design and optimization for geodetic networks, measured with conventional survey methods (not satellite), the development and implementation of a computational system capable to help on the definition of the most liable and accurate geometry, depending on the land orography where the network has to be located. First of all, a study of the methodology by least squares adjustment and propagation of variances will be held; then, subsequently, analyze its dependency of the geometry that the network will take. It will be essential to determine the independency of redundancy matrix (R) from the observations and its absolute dependency from the network geometry, as well as the influence of the diagonal terms of the R matrix (rii), redundancy numbers, in order to ensure maximum re liability of the network. It will also be analyzed first the behavior of redundancy numbers (rii) in surveying network design, then the variation of these values depending on the geometry with the analysis of its independency from the observations, and finally the different design levels depending on parameters and known data. It should be stated that network optimization, according to exposed criteria, is subject to the accessibility of the network points. In addition, common visibility among network points, which of them are connected with observations, has to be considered. All these situations depends essentially on the terrain relief and the natural or artificial obstacles that should exist. Therefore, it is necessary to include, at least, a digital terrain model (DTM), and better a digital surface model (DSM), not always available. Although design treatment is based on a bidimensional system, the possibility of incorporating a digital surface model (DSM) will be studied; this will allow evaluating the observations feasibility based on the terrain and the elements, both natural and artificial, which are located on it, when selecting network point locations. This system would provide, at first, an optimal design of a constrained network, considering both the internal reliability and the accuracy of its points (including the relief). This approach would amount to solving a “two and a half dimensional”3 design, if a digital surface model is available. As the availability of free DSM4 of the areas of interest of the project today is expensive, the possibility of combining a digital terrain model will arise. The activities to be developed on this PhD thesis are described in this document and are part of the research for which the following overall objectives are posed: 1. To establish a mathematical model for the process of observation of a survey network, considering all the factors involved and its influence on the estimates of the unknowns that are obtained as a result of the observations adjustment. 2. To develop a system to optimize a survey network results, applying design and simulation techniques on the previous model. 3. To present an explicit and rigorous formulation of parameters which assess the reliability of a survey network and its relations with the design. The achievement of this objective is based, besides on the search and review of sources, in an intense work of unification of notation and construction of intermediate steps in the mathematical developments. 4. To develop an overview of the influence on the network design of six major factors (posterior accuracy, observations reliability, viability of observations, instruments and station methodology) as optimization criteria, in order to define the subject approached on this document. 5. To elaborate and program the algorithms needed to develop an application software capable of considering the variables proposed in the previous section, on the problem of design and simulation of surveying networks, considering the digital surface model. It could be considered as secondary objectives, the following paragraphs: To develop the necessary algorithms to interrelate the digital terrain model with the design ones. To implement in the software application the possibility of variation of the coverage criteria parameters (normal distribution or Student t test) and therefore its degree of reliability.

Terahertz Response of Excitons in Nanorod Heterostructures

Quantum-confined systems are one of the most promising ways to enable us to control a material's interactions with light. Nanorods in particular offer the right dimensions for exploring and manipulating the terahertz region of the spectrum. In this thesis, we model excitons confined inside a nanorod using the envelope function approximation. A region-matching transfer matrix method allows us to simulate excitonic states inside arbitrary heterostructures grown along the length of the rod. We apply the method to colloidal CdSe rods 70 nm in length and under 10 nm in diameter, capped with ligands of DDPA and pyridine. We extend past studies on these types of rods by taking into account their dielectric permittivity mismatch. Compared to previous calculations and experimentally measured terahertz absorption, we predict a higher energy main 1S$z$ to 2P$z$ transition peak. This indicates that the rods are likely larger in diameter than previously thought. We also investigate a nanorod with GaAs/Al$_{0.3}$Ga$_{0.7}$As coupled double dots. The excitonic transitions were found to be manipulable by varying the strength of an applied electric field. We employ quasi-static state population distributions to simulate the effects of exciton relaxation from optically active states to dim ground states. A critical value of the applied field, corresponding to the exciton binding energy of ~18 meV, was found to dramatically alter the terahertz absorption due to state mixing. Above this critical field, more nuanced shifts in transition energies were observed, and gain from radiative relaxation to the ground state is predicted.

Some twisted results

The Drinfeld twist for the opposite quasi-Hopf algebra, H-COP, is determined and is shown to be related to the (second) Drinfeld twist on a quasi-Hopf algebra. The twisted form of the Drinfeld twist is investigated. In the quasi-triangular case, it is shown that the Drinfeld u-operator arises from the equivalence of H-COP to the quasi-Hopf algebra induced by twisting H with the R-matrix. The Altschuler-Coste u-operator arises in a similar way and is shown to be closely related to the Drinfeld u-operator. The quasi-cocycle condition is introduced and is shown to play a central role in the uniqueness of twisted structures on quasi-Hopf algebras. A generalization of the dynamical quantum Yang-Baxter equation, called the quasi-dynamical quantum Yang-Baxter equation, is introduced.

Effect of anode material and cavity design on the performance of microcavity OLEDs

We report on the effect of the replacement of the conventional ITO anode with the semitransparent metallic material on the performance of microcavity OLEDs. We performed comprehensive simulations of the emission from microcavity OLEDs consisting of widely used organic materials, N,N′-di(naphthalene-1- yl)-N,N′-diphenylbenzidine (NPB) as a hole transport layer and tris (8-hydroxyquinoline) (Alq3) as emitting and electron transporting layer. Silver and LiF/Al were considered as a cathode, while metallic (Au and Ag) anode was used and simulations were performed on devices with both the metallic and conventional ITO anode. The electroluminescence emission spectra, electric field distribution inside the device, carrier density, recombination rate and exciton density were calculated as a function of the position of the emission layer. The results show that the metallic anode enhances light output and that optimum emission from a microcavity OLED is achieved when the position of the recombination region is aligned with the antinode of the standing wave inside the cavity. The microcavity OLED devices with Ag/Ag and Ag/Au mirrors were fabricated and characterized. The experimental results have been compared to the simulations and the influence of the different anode, emission region width and position on the performance of microcavity OLEDs was discussed.

Electromagnetic interference from variable speed motor drives

A methodology is presented which can be used to produce the level of electromagnetic interference, in the form of conducted and radiated emissions, from variable speed drives, the drive that was modelled being a Eurotherm 583 drive. The conducted emissions are predicted using an accurate circuit model of the drive and its associated equipment. The circuit model was constructed from a number of different areas, these being: the power electronics of the drive, the line impedance stabilising network used during the experimental work to measure the conducted emissions, a model of an induction motor assuming near zero load, an accurate model of the shielded cable which connected the drive to the motor, and finally the parasitic capacitances that were present in the drive modelled. The conducted emissions were predicted with an error of +/-6dB over the frequency range 150kHz to 16MHz, which compares well with the limits set in the standards which specify a frequency range of 150kHz to 30MHz. The conducted emissions model was also used to predict the current and voltage sources which were used to predict the radiated emissions from the drive. Two methods for the prediction of the radiated emissions from the drive were investigated, the first being two-dimensional finite element analysis and the second three-dimensional transmission line matrix modelling. The finite element model took account of the features of the drive that were considered to produce the majority of the radiation, these features being the switching of the IGBT's in the inverter, the shielded cable which connected the drive to the motor as well as some of the cables that were present in the drive.The model also took account of the structure of the test rig used to measure the radiated emissions. It was found that the majority of the radiation produced came from the shielded cable and the common mode currents that were flowing in the shield, and that it was feasible to model the radiation from the drive by only modelling the shielded cable. The radiated emissions were correctly predicted in the frequency range 30MHz to 200MHz with an error of +10dB/-6dB. The transmission line matrix method modelled the shielded cable which connected the drive to the motor and also took account of the architecture of the test rig. Only limited simulations were performed using the transmission line matrix model as it was found to be a very slow method and not an ideal solution to the problem. However the limited results obtained were comparable, to within 5%, to the results obtained using the finite element model.

Study of Periodic and Quasi-Periodic Structures in Silicon-on-Insulator

In this thesis, a numerical design approach has been proposed and developed based on the transmission matrix method in order to characterize periodic and quasi-periodic photonic structures in silicon-on-insulator. The approach and its performance have been extensively tested with specific structures in 2D and its validity has been verified in 3D.

Acousto-optic modulation of optical reflection properties in fiber Bragg gratings and the application in fiber lasers

This thesis aims to investigate the interaction of acoustic waves and fiber Bragg gratings (FBGs) in standard and suspended-core fibers (SCFs), to evaluate the influence of the fiber, grating and modulator design on the increase of the modulation efficiency, bandwidth and frequency. Initially, the frequency response and the resonant acoustic modes of a low frequency acousto-optic modulator (f < 1.2 MHz) are numerically investigated by using the finite element method. Later, the interaction of longitudinal acoustic waves and FBGs in SCFs is also numerically investigated. The fiber geometric parameters are varied and the strain and grating properties are simulated by means of the finite element method and the transfer matrix method. The study indicates that the air holes composing the SCF cause a significant reduction of the amount of silica in the fiber cross section increasing acousto-optic interaction in the core. Experimental modulation of the reflectivity of FBGs inscribed in two distinct SCFs indicates evidences of this increased interaction. Besides, a method to acoustically induce a dynamic phase-shift in a chirped FBG employing an optimized design of modulator is shown. Afterwards, a combination of this modulator and a FBG inscribed in a three air holes SCF is applied to mode-lock an ytterbium doped fiber laser. To improve the modulator design for future applications, two other distinct devices are investigated to increase the acousto-optic interaction, bandwidth and frequency (f > 10 MHz). A high reflectivity modulation has been achieved for a modulator based on a tapered fiber. Moreover, an increased modulated bandwidth (320 pm) has been obtained for a modulator based on interaction of a radial long period grating (RLPG) and a FBG inscribed in a standard fiber. In summary, the results show a considerable reduction of the grating/fiber length and the modulator size, indicating possibilities for compact and faster acousto-optic fiber devices. Additionally, the increased interaction efficiency, modulated bandwidth and frequency can be useful to shorten the pulse width of future all-fiber mode-locked fiber lasers, as well, to other photonic devices which require the control of the light in optical fibers by electrically tunable acoustic waves.

Search for the Higgs boson produced in association with Z→ℓ+ℓ- using the matrix element method at CDF II

We present a search for associated production of the standard model (SM) Higgs boson and a $Z$ boson where the $Z$ boson decays to two leptons and the Higgs decays to a pair of $b$ quarks in $p\bar{p}$ collisions at the Fermilab Tevatron. We use event probabilities based on SM matrix elements to construct a likelihood function of the Higgs content of the data sample. In a CDF data sample corresponding to an integrated luminosity of 2.7 fb$^{-1}$ we see no evidence of a Higgs boson with a mass between 100 GeV$/c^2$ and 150 GeV$/c^2$. We set 95% confidence level (C.L.) upper limits on the cross-section for $ZH$ production as a function of the Higgs boson mass $m_H$; the limit is 8.2 times the SM prediction at $m_H = 115$ GeV$/c^2$.

Top Quark Mass Measurement in the lepton+jets Channel Using a Matrix Element Method and in situ Jet Energy Calibration

A precision measurement of the top quark mass m_t is obtained using a sample of ttbar events from ppbar collisions at the Fermilab Tevatron with the CDF II detector. Selected events require an electron or muon, large missing transverse energy, and exactly four high-energy jets, at least one of which is tagged as coming from a b quark. A likelihood is calculated using a matrix element method with quasi-Monte Carlo integration taking into account finite detector resolution and jet mass effects. The event likelihood is a function of m_t and a parameter DJES to calibrate the jet energy scale /in situ/. Using a total of 1087 events, a value of m_t = 173.0 +/- 1.2 GeV/c^2 is measured.

Top quark mass measurement in the lepton plus jets channel using a modified matrix element method

We report a measurement of the top quark mass, m_t, obtained from ppbar collisions at sqrt(s) = 1.96 TeV at the Fermilab Tevatron using the CDF II detector. We analyze a sample corresponding to an integrated luminosity of 1.9 fb^-1. We select events with an electron or muon, large missing transverse energy, and exactly four high-energy jets in the central region of the detector, at least one of which is tagged as coming from a b quark. We calculate a signal likelihood using a matrix element integration method, with effective propagators to take into account assumptions on event kinematics. Our event likelihood is a function of m_t and a parameter JES that determines /in situ/ the calibration of the jet energies. We use a neural network discriminant to distinguish signal from background events. We also apply a cut on the peak value of each event likelihood curve to reduce the contribution of background and badly reconstructed events. Using the 318 events that pass all selection criteria, we find m_t = 172.7 +/- 1.8 (stat. + JES) +/- 1.2 (syst.) GeV/c^2.