Hamiltonian structure of propagation equations for ultrashort optical pulses

A Hamiltonian framework is developed for a sequence of ultrashort optical pulses propagating in a nonlinear dispersive medium. To this end a second-order nonlinear wave equation is first simplified using an unidirectional approximation. All non-resonant nonlinear terms are then rigorously eliminated using a suitable change of variables in the spirit of the canonical perturbation theory. The derived propagation equation operates with a properly defined complexification of the real electric field. It accounts for arbitrary dispersion, four-wave mixing processes, weak absorption, and arbitrary pulse duration. Thereafter the so called normal variables, i.e., classical fields corresponding to the quantum creation and annihilation operators, are introduced. Neglecting absorption we finally derive the Hamiltonian formulation. The latter yields the most essential integrals of motion for the pulse propagation. These integrals reflect the time-averaged fluxes of energy, momentum, and classical photon number transferred by the pulse. The conservation laws are further used to control the numerical solutions when calculating supercontinuum generation by an ultrashort optical pulse

Hamiltonian structure of propagation equations for ultrashort optical pulses

A Hamiltonian framework is developed for a sequence of ultrashort optical pulses propagating in a nonlinear dispersive medium. To this end a second-order nonlinear wave equation is first simplified using an unidirectional approximation. All non-resonant nonlinear terms are then rigorously eliminated using a suitable change of variables in the spirit of the canonical perturbation theory. The derived propagation equation operates with a properly defined complexification of the real electric field. It accounts for arbitrary dispersion, four-wave mixing processes, weak absorption, and arbitrary pulse duration. Thereafter the so called normal variables, i.e., classical fields corresponding to the quantum creation and annihilation operators, are introduced. Neglecting absorption we finally derive the Hamiltonian formulation. The latter yields the most essential integrals of motion for the pulse propagation. These integrals reflect the time-averaged fluxes of energy, momentum, and classical photon number transferred by the pulse. The conservation laws are further used to control the numerical solutions when calculating supercontinuum generation by an ultrashort optical pulse

Weyl Pseudodifferential Calculus and the Heisenberg Group in New Settings

We examine both Weyl pseudodifferential calculus and the Heisenberg group in new settings, through two published papers and a chapter representing work towards a construction of Haar bases of unimodular LC groups. We finish with some remarks about an extension of both Heisenberg group and pseudodifferential calculus techniques to non-Abelian settings, detailing possible links with representation theory and the Langlands program. The first of the two papers achieves the following "We give a simple proof of the fact that the classical Ornstein-Uhlenbeck operator L is R-sectorial of angle arcsin|1-2/p| on $L^p(\mathbb{R}^n, \exp(-|x|^2/2)dx)$ (for $1<p<\infty$). Applying the abstract holomorphic functional calculus theory of Kalton and Weis, this immediately gives a new proof of the fact that L has a bounded $H^\infty$ functional calculus with this optimal angle." In the second paper, "We construct a Weyl pseudodifferential calculus tailored to studying boundedness of operators on weighted $L^p$ spaces over $\mathbb{R}^d$ with weights of the form $exp(-\phi(x))$, for $\phi(x)$ a $C^2$ function, a setting in which the operator associated to the weighted Dirichlet form typically has only holomorphic functional calculus. A symbol class giving rise to bounded operators on $L^p$ is determined, and its properties analysed. This theory is used to calculate an upper bounded on the $H^\infty$ angle of relevant operators, and deduces known optimal results in some cases. Finally, the symbol class is enriched and studied under an algebraic viewpoint." The construction of Haar bases on unimodular LC groups proceeds via the tools of fractal tilings. After a review of these concepts we prove the key results required to obtain a "good" Haar basis, namely that the boundary of the relevant tilings is of measure $0$. We explain how the constructed Haar bases can be used to study Fourier multipliers in such settings, detailing future work

States of Low Energy on Bianchi I spacetimes

States of Low Energy are a class of exact Hadamard states for free quantum fields on cosmological spacetimes whose structure is fixed at {\it all} scales by a minimization principle. The original construction was for Friedmann-Lema\^{i}tre geometries and is here generalized to anisotropic Bianchi I geometries relevant to primordial cosmology. In addition to proving the Hadamard property systematic series expansions in the infrared and ultraviolet are developed. The infrared expansion is convergent and induces in the massless case a leading spatial long distance decay that is always Minkowski-like but anisotropy modulated. For the ultraviolet expansion a non-recursive formula for the coefficients is presented.Comment: 44 page

Institute for Computational Mechanics in Propulsion (ICOMP)

The Institute for Computational Mechanics in Propulsion (ICOMP) is operated by the Ohio Aerospace Institute (OAI) and funded under a cooperative agreement by the NASA Lewis Research Center in Cleveland, Ohio. The purpose of ICOMP is to develop techniques to improve problem-solving capabilities in all aspects of computational mechanics related to propulsion. This report describes the activities at ICOMP during 1994

Selected developments in computational electromagnetics for radio engineering

This thesis deals with the development and application of two simulation methods commonly used in radio engineering, namely the Finite-Difference Time-Domain method (FDTD) and the Finite Element Method (FEM). The main emphasis of this thesis is in FDTD. FDTD has become probably the most popular computational technique in radio engineering. It is a well established, fairly accurate and easy-to-implement method. Being a time-domain method, it can provide wide-band information in a single simulation. It simulates physical wave propagation in the computational volume, and is thus especially useful for educational purposes and for gaining engineering insight into complicated wave interaction and coupling phenomena. In this thesis, numerical dispersion taking place in the FDTD algorithm is analyzed, and a novel dispersion reduction procedure is described, based on artificial anisotropy. As a result, larger cells can be used to obtain the same accuracy in terms of dispersion error. Simulation experiments suggest that typically the dispersion reduction allows roughly doubling the cell size in each coordinate direction, without sacrificing the accuracy. The obtainable advantage is, however, dependent on the problem. In the open literature, a few other procedures are also presented to reduce the dispersion error. However, the rather dominating effect of unequal grid resolution along different coordinate directions has been neglected in previous studies. The so-called Perfectly Matched Layer (PML) has proven to be a very useful absorbing boundary condition (ABC) in FDTD simulations. It is reliable, works well in wide frequency band and is easy to implement. The most notable deficiency of PML is that it enlarges the computational volume - in open 3-D structures easily by a factor of two. However, due to its advantages, PML has become a standard ABC. In this thesis, the operation of PML in FDTD has been studied theoretically, and some interesting properties of it not known before are uncovered. For example, it is shown that, surprisingly, PML can absorb perfectly (i.e. with zero reflection) plane waves propagating towards almost arbitrary given direction at given frequency. Optimizing the conductivity profile allows reduction of the PML thickness. A typical application of the FDTD method is the design of a mobile handset antenna. An improved coaxial probe model has been developed for antenna simulations. The well-known resistive voltage source (RVS) model has also been discussed. A reference plane transformation is proposed to correct the simulated input impedance. A popular thin-wire model in 2-D FDTD is discussed, and it is shown to be based on erroneous reasoning. The error has been corrected by a simple procedure, and the corrected model has been demonstrated to simulate infinite long thin wires much better than the commonly used model. A novel way to implement singular basis functions in FEM is discussed. It is shown theoretically and demonstrated by examples that if a waveguide propagation mode contains field singularities, then explicit inclusion of singularities in finite element analysis is crucial in order to obtain accurate cut-off wavenumbers.reviewe