Source author record

Andreas Adelmann

Andreas Adelmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

muEDM: Towards a search for the muon electric dipole moment at PSI using the frozen-spin technique

The search for a permanent electric dipole moment (EDM) of the muon is an excellent probe for physics beyond the Standard Model of particle physics. We propose the first dedicated muon EDM search employing the frozen-spin technique at the Paul Scherrer Institute (PSI), Switzerland, with a sensitivity of $6 \times 10^{-23}~e\!\cdot\!\mathrm{cm}$, improving the current best limit set by the E821 experiment at Brookhaven National Laboratory by more than three orders of magnitude. In preparation for a high precision experiment to measure the muon EDM, several R&D studies have been performed at PSI: the characterisation of a possible beamline to host the experiment for the muon beam injection study and the measurement of the multiple Coulomb scattering of positrons in potential detector materials at low momenta for the positron tracking scheme development. This paper discusses experimental concepts and the current status of the muEDM experiment at PSI.

preprint2022arXiv

New directions for surrogate models and differentiable programming for High Energy Physics detector simulation

The computational cost for high energy physics detector simulation in future experimental facilities is going to exceed the current available resources. To overcome this challenge, new ideas on surrogate models using machine learning methods are being explored to replace computationally expensive components. Additionally, differentiable programming has been proposed as a complementary approach, providing controllable and scalable simulation routines. In this document, new and ongoing efforts for surrogate models and differential programming applied to detector simulation are discussed in the context of the 2021 Particle Physics Community Planning Exercise (`Snowmass').

preprint2022arXiv

Report of the Snowmass'21 Workshop on High-Power Cyclotrons and FFAs

This whitepaper summarizes and the state of the field of high-power cyclotrons and FFAs as discussed by international experts during a three-day workshop of the same name. The workshop was held online from Sep 7 to Sep 9, 2021 as part of the US Snowmass'21 Community Exercise, specifically the Accelerator Frontier (AF) and the subpanel Accelerators for Neutrinos (AF02). Thus, we put emphasis on the application of high-power cyclotrons in particle physics, specifically neutrino physics, and as drivers for muon production. In the introduction, we discuss the role of cyclotrons for particle physics, and later we highlight existing and planned experiments in the corresponding sections. However, as these same accelerators have important applications in the fields of isotope production - both for research and medicine - and possibly even in energy research, by providing beam to demonstrator experiments in the areas of Accelerator Driven Systems (ADS), we include these far-reaching topics to provide a full picture of the status and applications of high-power cyclotrons. Furthermore, Fixed Field Alternating Gradient accelerators (FFAs) have recently seen renewed interest. They are in many respects (basic operating principles) similar to cyclotrons and have thus been included in this workshop and whitepaper as well. We are discussing current projects and whether FFAs have the prospect of becoming high-intensity machines.

preprint2021arXiv

Benchmarking Collective Effects of Electron Interactions in a Wiggler with OPAL-FEL

OPAL-FEL is a recently developed tool for the modeling of particle accelerators containing wigglers or undulators. It extends the well established 3D electrostatic particle-tracking code OPAL, by merging it with the finite-difference time-domain electromagnetic solver MITHRA. We present results of two benchmark cases where OPAL-FEL simulations are compared to experimental results. Both experiments concern electron beamlines where the longitudinal phase space is modulated with a short magnetic wiggler. Good agreement was found in both the space charge and radiation dominated regimes.

preprint2021arXiv

Order-of-Magnitude Beam Current Improvement in Compact Cyclotrons

There is great need for high intensity proton beams from compact particle accelerators in particle physics, medical isotope production, and materials- and energy-research. To address this need, we present, for the first time, a design for a compact isochronous cyclotron that will be able to deliver 10 mA of 60 MeV protons - an order of magnitude higher than on-market compact cyclotrons and a factor four higher than research machines. A key breakthrough is that vortex motion is incorporated in the design of a cyclotron, leading to clean extraction. Beam losses on the septa of the electrostatic extraction channels stay below 50 W (a factor four below the required safety limit), while maintaining good beam quality. We present a set of highly accurate particle-in-cell simulations, and an uncertainty quantification of select beam input parameters using machine learning, showing the robustness of the design. This design can be utilized for beams for experiments in particle and nuclear physics, materials science and medical physics as well as for industrial applications.

preprint2020arXiv

Global Sensitivity Analysis on Numerical Solver Parameters of Particle-In-Cell Models in Particle Accelerator Systems

Every computer model depends on numerical input parameters that are chosen according to mostly conservative but rigorous numerical or empirical estimates. These parameters could for example be the step size for time integrators, a seed for pseudo-random number generators, a threshold or the number of grid points to discretize a computational domain. In case a numerical model is enhanced with new algorithms and modelling techniques the numerical influence on the quantities of interest, the running time as well as the accuracy is often initially unknown. Usually parameters are chosen on a trial-and-error basis neglecting the computational cost versus accuracy aspects. As a consequence the cost per simulation might be unnecessarily high which wastes computing resources. Hence, it is essential to identify the most critical numerical parameters and to analyze systematically their effect on the result in order to minimize the time-to-solution without losing significantly on accuracy. Relevant parameters are identified by global sensitivity studies where Sobol' indices are common measures. These sensitivities are obtained by surrogate models based on polynomial chaos expansion. In this paper, we first introduce the general methods for uncertainty quantification. We then demonstrate their use on numerical solver parameters to reduce the computational costs and discuss further model improvements based on the sensitivity analysis. The sensitivities are evaluated for neighbouring bunch simulations of the existing PSI Injector II and PSI Ring as well as the proposed Daedalus Injector cyclotron and simulations of the rf electron gun of the Argonne Wakefield Accelerator.

preprint2020arXiv

Machine Learning for Orders of Magnitude Speedup in Multi-Objective Optimization of Particle Accelerator Systems

High-fidelity physics simulations are powerful tools in the design and optimization of charged particle accelerators. However, the computational burden of these simulations often limits their use in practice for design optimization and experiment planning. It also precludes their use as online models tied directly to accelerator operation. We introduce an approach based on machine learning to create nonlinear, fast-executing surrogate models that are informed by a sparse sampling of the physics simulation. The models are 10^6 to 10^7 times more efficient to execute.We also demonstrate that these models can be reliably used with multi-objective optimization to obtain orders-of-magnitude speedup in initial design studies and experiment planning. For example, we required 132 times fewer simulation evaluations to obtain an equivalent solution for our main test case, and initial studies suggest that between 330 to 550 times fewer simulation evaluations are needed when using an iterative retraining process. Our approach enables new ways for high-fidelity particle accelerator simulations to be used, at comparatively little computational cost.

preprint2020arXiv

Multiobjective optimization of the dynamic aperture for SLS 2.0 using surrogate models based on artificial neural networks

Modern synchrotron light source storage rings, such as the Swiss Light Source upgrade (SLS 2.0), use multi-bend achromats in their arc segments to achieve unprecedented brilliance. This performance comes at the cost of increased focusing requirements, which in turn require stronger sextupole and higher-order multipole fields for compensation and lead to a considerable decrease in the dynamic aperture and/or energy acceptance. In this paper, to increase these two quantities, a multi-objective genetic algorithm (MOGA) is combined with a modified version of the well-known tracking code tracy. As a first approach, a massively parallel implementation of a MOGA is used. Compared to a manually obtained solution this approach yields very good results. However, it requires a long computation time. As a second approach, a surrogate model based on artificial neural networks is used in the optimization. This improves the computation time, but the results quality deteriorates. As a third approach, the surrogate model is re-trained during the optimization. This ensures a solution quality comparable to the one obtained with the first approach while also providing an order of magnitude speedup. Finally, good candidate solutions for SLS 2.0 are shown and further analyzed.

preprint2020arXiv

Scientific opportunies for bERLinPro 2020+, report with ideas and conclusions from bERLinProCamp 2019

The Energy Recovery Linac (ERL) paradigm offers the promise to generate intense electron beams of superior quality with extremely small six-dimensional phase space for many applications in the physical sciences, materials science, chemistry, health, information technology and security. Helmholtz-Zentrum Berlin started in 2010 an intensive R\&D programme to address the challenges related to the ERL as driver for future light sources by setting up the bERLinPro (Berlin ERL Project) ERL with 50 MeV beam energy and high average current. The project is close to reach its major milestone in 2020, acceleration and recovery of a high brightness electron beam. The goal of bERLinProCamp 2019 was to discuss scientific opportunities for bERLinPro 2020+. bERLinProCamp 2019 was held on Tue, 17.09.2019 at Helmholtz-Zentrum Berlin, Berlin, Germany. This paper summarizes the main themes and output of the workshop.

preprint2020arXiv

Sparse Grids based Adaptive Noise Reduction strategy for Particle-In-Cell schemes

We propose a sparse grids based adaptive noise reduction strategy for electrostatic particle-in-cell (PIC) simulations. Our approach is based on the key idea of relying on sparse grids instead of a regular grid in order to increase the number of particles per cell for the same total number of particles, as first introduced in Ricketson and Cerfon (Plasma Phys. and Control. Fusion, 59(2), 024002). Adopting a new filtering perspective for this idea, we construct the algorithm so that it can be easily integrated into high performance large-scale PIC code bases. Unlike the physical and Fourier domain filters typically used in PIC codes, our approach automatically adapts to mesh size, number of particles per cell, smoothness of the density profile and the initial sampling technique. Thanks to the truncated combination technique, we can reduce the larger grid-based error of the standard sparse grids approach for non-aligned and non-smooth functions. We propose a heuristic based on formal error analysis for selecting the optimal truncation parameter at each time step, and develop a natural framework to minimize the total error in sparse PIC simulations. We demonstrate its efficiency and performance by means of two test cases: the diocotron instability in two dimensions, and the three-dimensional electron dynamics in a Penning trap. Our run time performance studies indicate that our new scheme can provide significant speedup and memory reduction as compared to regular PIC for achieving comparable accuracy in the charge density deposition.

preprint2019arXiv

Constrained multi-objective shape optimization of superconducting RF cavities considering robustness against geometric perturbations

High current storage rings, such as the Z-pole operating mode of the FCC-ee, require accelerating cavities that are optimized with respect to both the fundamental mode and the higher order modes. Furthermore, the cavity shape needs to be robust against geometric perturbations which could, for example, arise from manufacturing inaccuracies or harsh operating conditions at cryogenic temperatures. This leads to a constrained multi-objective shape optimization problem which is computationally expensive even for axisymmetric cavity shapes. In order to decrease the computation cost, a global sensitivity analysis is performed and its results are used to reduce the search space and redefine the objective functions. A massively parallel implementation of an evolutionary algorithm, combined with a fast axisymmetric Maxwell eigensolver and a frequency-tuning method is used to find an approximation of the Pareto front. The computed Pareto front approximation and a cavity shape with desired properties are shown. Further, the approach is generalized and applied to another type of cavity.

preprint2019arXiv

On Architecture and Performance of Adaptive Mesh Refinement in an Electrostatics Particle-In-Cell Code

This article presents a hardware architecture independent implementation of an adaptive mesh refinement Poisson solver that is integrated into the electrostatic Particle-In-Cell beam dynamics code OPAL. The Poisson solver is solely based on second generation Trilinos packages to ensure the desired hardware portability. Based on the massively parallel framework AMReX, formerly known as BoxLib, the new adaptive mesh refinement interface provides several refinement policies in order to enable precise large-scale neighbouring bunch simulations in high intensity cyclotrons. The solver is validated with a built-in multigrid solver of AMReX and a test problem with analytical solution. The parallel scalability is presented as well as an example of a neighbouring bunch simulation that covers the scale of the later anticipated physics simulation.

preprint2016arXiv

Examination of the Plasma located in PSI Ring Cyclotron

A plasma has been observed inside the vacuum chamber of the PSI Ring Cyclotron. This ionized gas cloud maybe a substantial contributor to several interior components having reduced lifetimes. The plasma's generation has been directly linked to the voltage that is applied to the Flat Top Cavity through visual confirmation using CCD cameras. A spectrometer was used to correlate the plasma's intensity and ignition to the Flat Top Cavity voltage as well as to determine the composition of the plasma. This paper reports on the analysis of the plasma using spectroscopy. The spectrometer data was analyzed to determine the composition of the plasma and that the plasma intensity (luminosity) directly corresponds to the Flat Top voltage. The results showed that the plasma was comprised of elements consistent with the cyclotrons vacuum interior

preprint2015arXiv

A Homotopy Method for Large-Scale Multi-Objective Optimization

A homotopy method for multi-objective optimization that produces uniformly sampled Pareto fronts by construction is presented. While the algorithm is general, of particular interest is application to simulation-based engineering optimization problems where economy of function evaluations, smoothness of result, and time-to-solution are critical. The presented algorithm achieves an order of magnitude improvement over other geometrically motivated methods, like Normal Boundary Intersection and Normal Constraint, with respect to solution evenness for similar computational expense. Furthermore, the resulting uniformity of solutions extends even to more difficult problems, such as those appearing in common Evolutionary Algorithm test cases.

preprint2015arXiv

The Dynamical Kernel Scheduler - Part 1

Emerging processor architectures such as GPUs and Intel MICs provide a huge performance potential for high performance computing. However developing software using these hardware accelerators introduces additional challenges for the developer such as exposing additional parallelism, dealing with different hardware designs and using multiple development frameworks in order to use devices from different vendors. The Dynamic Kernel Scheduler (DKS) is being developed in order to provide a software layer between host application and different hardware accelerators. DKS handles the communication between the host and device, schedules task execution, and provides a library of built-in algorithms. Algorithms available in the DKS library will be written in CUDA, OpenCL and OpenMP. Depending on the available hardware, the DKS can select the appropriate implementation of the algorithm. The first DKS version was created using CUDA for the Nvidia GPUs and OpenMP for Intel MIC. DKS was further integrated in OPAL (Object-oriented Parallel Accelerator Library) to speed up a parallel FFT based Poisson solver and Monte Carlo simulations for particle matter interaction used for proton therapy degrader modeling. DKS was also used together with Minuit2 for parameter fitting, where $χ^2$ and $max-log-likelihood$ functions were offloaded to the hardwareaccelerator. The concepts of the DKS, first results together with plans for the future will be shown in this paper.

preprint2014arXiv

Multilevel Monte Carlo for the Feynman-Kac Formula for the Laplace Equation

Since its formulation in the late 1940s, the Feynman-Kac formula has proven to be an effective tool for both theoretical reformulations and practical simulations of differential equations. The link it establishes between such equations and stochastic processes can be exploited to develop Monte Carlo sampling methods that are effective, especially in high dimensions. There exist many techniques of improving standard Monte Carlo sampling methods, a relatively new development being the so-called Multilevel Monte Carlo method. This paper investigates the applicability of multilevel ideas to the stochastic representation of partial differential equations by the Feynman- Kac formula, using the Walk on Sphere algorithm to generate the required random paths. We focus on the Laplace equation, the simplest elliptic PDE, while mentioning some extension possibilities.

preprint2012arXiv

A novel adaptive time stepping variant of the Boris-Buneman integrator for the simulation of particle accelerators with space charge

We show that adaptive time stepping in particle accelerator simulation is an enhancement for certain problems. The new algorithm has been implemented in the OPAL (Object Oriented Parallel Accelerator Library) framework, and is compared to the existing code. The idea is to adjust the frequency of costly self field calculations, which are needed to model Coulomb interaction (space charge) effects. In analogy to a Kepler orbit simulation that requires a higher time step resolution at the close encounter, we propose to choose the time step based on the magnitude of the space charge forces. Inspired by geometric integration techniques, our algorithm chooses the time step proportional to a function of the current phase space state instead of calculating a local error estimate like a conventional adaptive procedure. In this paper we build up on first observations made in recent work. A more profound argument is given on how exactly the time step should be chosen. An intermediate algorithm, initially built to allow a clearer analysis by introducing separate time steps for external field and self field integration, turned out to be useful in itself already for a large class of problems.

preprint2011arXiv

Perfectly Matched Layers in a Divergence Preserving ADI Scheme for Electromagnetics

For numerical simulations of highly relativistic and transversely accelerated charged particles including radiation fast algorithms are needed. While the radiation in particle accelerators has wavelengths in the order of 100 um the computational domain has dimensions roughly 5 orders of magnitude larger resulting in very large mesh sizes. The particles are confined to a small area of this domain only. To resolve the smallest scales close to the particles subgrids are envisioned. For reasons of stability the alternating direction implicit (ADI) scheme by D. N. Smithe et al. (J. Comput. Phys. 228 (2009) pp.7289-7299) for Maxwell equations has been adopted. At the boundary of the domain absorbing boundary conditions have to be employed to prevent reflection of the radiation. In this paper we show how the divergence preserving ADI scheme has to be formulated in perfectly matched layers (PML) and compare the performance in several scenarios.