A Rigorous Finite-Element Domain Decomposition Method for Electromagnetic Near Field Simulations

Rigorous computer simulations of propagating electromagnetic fields have become an important tool for optical metrology and design of nanostructured optical components. A vectorial finite element method (FEM) is a good choice for an accurate modeling of complicated geometrical features. However, from a numerical point of view solving the arising system of linear equations is very demanding even for medium sized 3D domains. In numerics, a domain decomposition method is a commonly used strategy to overcome this problem. Within this approach the overall computational domain is split up into smaller domains and interface conditions are used to assure continuity of the electromagnetic field. Unfortunately, standard implementations of the domain decomposition method as developed for electrostatic problems are not appropriate for wave propagation problems. In an earlier paper we therefore proposed a domain decomposition method adapted to electromagnetic field wave propagation problems. In this paper we apply this method to 3D mask simulation.Comment: 9 pages, 7 figures, SPIE conference Advanced Lithography / Optical Microlithography XXI (2008

JCMmode: An Adaptive Finite Element Solver for the Computation of Leaky Modes

We present our simulation tool JCMmode for calculating propagating modes of an optical waveguide. As ansatz functions we use higher order, vectorial elements (Nedelec elements, edge elements). Further we construct transparent boundary conditions to deal with leaky modes even for problems with inhomogeneous exterior domains as for integrated hollow core Arrow waveguides. We have implemented an error estimator which steers the adaptive mesh refinement. This allows the precise computation of singularities near the metal's corner of a Plasmon-Polariton waveguide even for irregular shaped metal films on a standard personal computer.Comment: 11 page

The Two-Screen Measurement Setup to Indirectly Measure Proton Beam Self-Modulation in AWAKE

The goal of the first phase of the AWAKE \cite{AWAKE1,AWAKE2} experiment at CERN is to measure the self-modulation \cite{SMI} of the $\sigma_z = 12\,\rm{cm}$ long SPS proton bunch into microbunches after traversing $10\,\rm{m}$ of plasma with a plasma density of $n_{pe}=7\times10^{14}\,\rm{electrons/cm}^3$. The two screen measurement setup \cite{Turner2016} is a proton beam diagnostic that can indirectly prove the successful development of the self-modulation of the proton beam by imaging protons that got defocused by the transverse plasma wakefields after passing through the plasma, at two locations downstream the end of the plasma. This article describes the design and realization of the two screen measurement setup integrated in the AWAKE experiment. We discuss the performance and background response of the system based on measurements performed with an unmodulated Gaussian SPS proton bunch during the AWAKE beam commissioning in September and October 2016. We show that the system is fully commissioned and adapted to eventually image the full profile of a self-modulated SPS proton bunch in a single shot measurement during the first phase of the AWAKE experiment.Comment: 5 pages 8 figure

Indirect Self-Modulation Instability Measurement Concept for the AWAKE Proton Beam

AWAKE, the Advanced Proton-Driven Plasma Wakefield Acceleration Experiment, is a proof-of-principle R&D experiment at CERN using a 400 GeV/c proton beam from the CERN SPS (longitudinal beam size sigma_z = 12 cm) which will be sent into a 10 m long plasma section with a nominal density of approx. 7x10^14 atoms/cm3 (plasma wavelength lambda_p = 1.2mm). In this paper we show that by measuring the time integrated transverse profile of the proton bunch at two locations downstream of the AWAKE plasma, information about the occurrence of the self-modulation instability (SMI) can be inferred. In particular we show that measuring defocused protons with an angle of 1 mrad corresponds to having electric fields in the order of GV/m and fully developed self-modulation of the proton bunch. Additionally, by measuring the defocused beam edge of the self-modulated bunch, information about the growth rate of the instability can be extracted. If hosing instability occurs, it could be detected by measuring a non-uniform defocused beam shape with changing radius. Using a 1 mm thick Chromox scintillation screen for imaging of the self-modulated proton bunch, an edge resolution of 0.6 mm and hence a SMI saturation point resolution of 1.2 m can be achieved.Comment: 4 pages, 4 figures, EAAC conference proceeding

Numerical Investigation of Light Scattering off Split-Ring Resonators

Recently, split ring-resonators (SRR's) have been realized experimentally in the near infrared (NIR) and optical regime. In this contribution we numerically investigate light propagation through an array of metallic SRR's in the NIR and optical regime and compare our results to experimental results. We find numerical solutions to the time-harmonic Maxwell's equations by using advanced finite-element-methods (FEM). The geometry of the problem is discretized with unstructured tetrahedral meshes. Higher order, vectorial elements (edge elements) are used as ansatz functions. Transparent boundary conditions and periodic boundary conditions are implemented, which allow to treat light scattering problems off periodic structures. This simulation tool enables us to obtain transmission and reflection spectra of plane waves which are incident onto the SRR array under arbitrary angles of incidence, with arbitrary polarization, and with arbitrary wavelength-dependencies of the permittivity tensor. We compare the computed spectra to experimental results and investigate resonances of the system.Comment: 9 pages, 8 figures (see original publication for images with a better resolution

Counting Hyperbolic Manifolds

No Abstrac

Bose-Condensed Gases in a 1D Optical Lattice at Finite Temperatures

We study equilibrium properties of Bose-Condensed gases in a one-dimensional (1D) optical lattice at finite temperatures. We assume that an additional harmonic confinement is highly anisotropic, in which the confinement in the radial directions is much tighter than in the axial direction. We derive a quasi-1D model of the Gross-Pitaeavkill equation and the Bogoliubov equations, and numerically solve these equations to obtain the condensate fraction as a function of the temperature.Comment: Comments: 6 pages, 3 figures, submitted to Quantum Fluids and Solids Conference (QFS 2006