Absorbing boundary conditions for the Westervelt equation

The focus of this work is on the construction of a family of nonlinear absorbing boundary conditions for the Westervelt equation in one and two space dimensions. The principal ingredient used in the design of such conditions is pseudo-differential calculus. This approach enables to develop high order boundary conditions in a consistent way which are typically more accurate than their low order analogs. Under the hypothesis of small initial data, we establish local well-posedness for the Westervelt equation with the absorbing boundary conditions. The performed numerical experiments illustrate the efficiency of the proposed boundary conditions for different regimes of wave propagation

A modified time-of-flight method for precise determination of high speed ratios in molecular beams

Time-of-flight (TOF) is a standard experimental technique for determining, among others, the speed ratio S (velocity spread) of a molecular beam. The speed ratio is a measure for the monochromaticity of the beam and an accurate determination of S is crucial for various applications, for example, for characterising chromatic aberrations in focussing experiments related to helium microscopy or for precise measurements of surface phonons and surface structures in molecular beam scattering experiments. For both of these applications, it is desirable to have as high a speed ratio as possible. Molecular beam TOF measurements are typically performed by chopping the beam using a rotating chopper with one or more slit openings. The TOF spectra are evaluated using a standard deconvolution method. However, for higher speed ratios, this method is very sensitive to errors related to the determination of the slit width and the beam diameter. The exact sensitivity depends on the beam diameter, the number of slits, the chopper radius, and the chopper rotation frequency. We present a modified method suitable for the evaluation of TOF measurements of high speed ratio beams. The modified method is based on a systematic variation of the chopper convolution parameters so that a set of independent measurements that can be fitted with an appropriate function are obtained. We show that with this modified method, it is possible to reduce the error by typically one order of magnitude compared to the standard method

Parameter identification in a semilinear hyperbolic system

We consider the identification of a nonlinear friction law in a one-dimensional damped wave equation from additional boundary measurements. Well-posedness of the governing semilinear hyperbolic system is established via semigroup theory and contraction arguments. We then investigte the inverse problem of recovering the unknown nonlinear damping law from additional boundary measurements of the pressure drop along the pipe. This coefficient inverse problem is shown to be ill-posed and a variational regularization method is considered for its stable solution. We prove existence of minimizers for the Tikhonov functional and discuss the convergence of the regularized solutions under an approximate source condition. The meaning of this condition and some arguments for its validity are discussed in detail and numerical results are presented for illustration of the theoretical findings

The Iteratively Regularized Gau{\ss}-Newton Method with Convex Constraints and Applications in 4Pi-Microscopy

This paper is concerned with the numerical solution of nonlinear ill-posed operator equations involving convex constraints. We study a Newton-type method which consists in applying linear Tikhonov regularization with convex constraints to the Newton equations in each iteration step. Convergence of this iterative regularization method is analyzed if both the operator and the right hand side are given with errors and all error levels tend to zero. Our study has been motivated by the joint estimation of object and phase in 4Pi microscopy, which leads to a semi-blind deconvolution problem with nonnegativity constraints. The performance of the proposed algorithm is illustrated both for simulated and for three-dimensional experimental data

Iteratively regularized Newton-type methods for general data misfit functionals and applications to Poisson data

We study Newton type methods for inverse problems described by nonlinear operator equations $F(u)=g$ in Banach spaces where the Newton equations $F'(u_n;u_{n+1}-u_n) = g-F(u_n)$ are regularized variationally using a general data misfit functional and a convex regularization term. This generalizes the well-known iteratively regularized Gauss-Newton method (IRGNM). We prove convergence and convergence rates as the noise level tends to 0 both for an a priori stopping rule and for a Lepski{\u\i}-type a posteriori stopping rule. Our analysis includes previous order optimal convergence rate results for the IRGNM as special cases. The main focus of this paper is on inverse problems with Poisson data where the natural data misfit functional is given by the Kullback-Leibler divergence. Two examples of such problems are discussed in detail: an inverse obstacle scattering problem with amplitude data of the far-field pattern and a phase retrieval problem. The performence of the proposed method for these problems is illustrated in numerical examples

Annihilation of low energy antiprotons in silicon

The goal of the AE$\mathrm{\bar{g}}$IS experiment at the Antiproton Decelerator (AD) at CERN, is to measure directly the Earth's gravitational acceleration on antimatter. To achieve this goal, the AE$\mathrm{\bar{g}}$IS collaboration will produce a pulsed, cold (100 mK) antihydrogen beam with a velocity of a few 100 m/s and measure the magnitude of the vertical deflection of the beam from a straight path. The final position of the falling antihydrogen will be detected by a position sensitive detector. This detector will consist of an active silicon part, where the annihilations take place, followed by an emulsion part. Together, they allow to achieve 1$$ precision on the measurement of $\bar{g}$ with about 600 reconstructed and time tagged annihilations. We present here, to the best of our knowledge, the first direct measurement of antiproton annihilation in a segmented silicon sensor, the first step towards designing a position sensitive silicon detector for the AE$\mathrm{\bar{g}}$IS experiment. We also present a first comparison with Monte Carlo simulations (GEANT4) for antiproton energies below 5 MeVComment: 21 pages in total, 29 figures, 3 table

Prospects for measuring the gravitational free-fall of antihydrogen with emulsion detectors

The main goal of the AEgIS experiment at CERN is to test the weak equivalence principle for antimatter. AEgIS will measure the free-fall of an antihydrogen beam traversing a moir\'e deflectometer. The goal is to determine the gravitational acceleration g for antihydrogen with an initial relative accuracy of 1% by using an emulsion detector combined with a silicon micro-strip detector to measure the time of flight. Nuclear emulsions can measure the annihilation vertex of antihydrogen atoms with a precision of about 1 - 2 microns r.m.s. We present here results for emulsion detectors operated in vacuum using low energy antiprotons from the CERN antiproton decelerator. We compare with Monte Carlo simulations, and discuss the impact on the AEgIS project.Comment: 20 pages, 16 figures, 3 table

Building the impedance model of a real machine

A reliable impedance model of a particle accelerator can be built by combining the beam coupling impedances of all the components. This is a necessary step to be able to evaluate the machine performance limitations, identify the main contributors in case an impedance reduction is required, and study the interaction with other mechanisms such as optics nonlinearities, transverse damper, noise, space charge, electron cloud, beam-beam (in a collider). The main phases to create a realistic impedance model, and verify it experimentally, will be reviewed, highlighting the main challenges. Some examples will be presented revealing the levels of precision of machine impedance models that have been achieved

A Moiré Deflectometer for Antimatter

The precise measurement of forces is one way to obtain deep insight into the fundamental interactions present in nature. In the context of neutral antimatter, the gravitational interaction is of high interest, potentially revealing new forces that violate the weak equivalence principle. Here we report on a successful extension of a tool from atom optics - the moirè deflectometer - for a measurement of the acceleration of slow antiprotons. The setup consists of two identical transmission gratings and a spatially resolving emulsion detector for antiproton annihilations. Absolute referencing of the observed antimatter pattern with a photon pattern experiencing no deflection allows the direct inference of forces present. The concept is also straightforwardly applicable to antihydrogen measurements as pursued by the AEgIS collaboration. The combination of these very different techniques from high energy and atomic physics opens a very promising route to the direct detection of the gravitational acceleration of neutral antimatter

Measuring the free fall of antihydrogen

After the first production of cold antihydrogen by the ATHENA and ATRAP experiments ten years ago, new second-generation experiments are aimed at measuring the fundamental properties of this anti-atom. The goal of AEGIS (Antimatter Experiment: Gravity, Interferometry, Spectroscopy) is to test the weak equivalence principle by studying the gravitational interaction between matter and antimatter with a pulsed, cold antihydrogen beam. The experiment is currently being assembled at CERN's Antiproton Decelerator. In AEGIS, antihydrogen will be produced by charge exchange of cold antiprotons with positronium excited to a high Rydberg state (n > 20). An antihydrogen beam will be produced by controlled acceleration in an electric-field gradient (Stark acceleration). The deflection of the horizontal beam due to its free fall in the gravitational field of the earth will be measured with a moire deflectometer. Initially, the gravitational acceleration will be determined to a precision of 1%, requiring the detection of about 105 antihydrogen atoms. In this paper, after a general description, the present status of the experiment will be reviewed