Showing 1–50 of 76 results for author: Goddard, B

Efficient nonlocal linear image denoising: Bilevel optimization with Nonequispaced Fast Fourier Transform and matrix-free preconditioning

Authors: Andrés Miniguano-Trujillo, John W. Pearson, Benjamin D. Goddard

Abstract: We present a new approach for nonlocal image denoising, based around the application of an unnormalized extended Gaussian ANOVA kernel within a bilevel optimization algorithm. A critical bottleneck when solving such problems for finely-resolved images is the solution of huge-scale, dense linear systems arising from the minimization of an energy term. We tackle this using a Krylov subspace approach… ▽ More We present a new approach for nonlocal image denoising, based around the application of an unnormalized extended Gaussian ANOVA kernel within a bilevel optimization algorithm. A critical bottleneck when solving such problems for finely-resolved images is the solution of huge-scale, dense linear systems arising from the minimization of an energy term. We tackle this using a Krylov subspace approach, with a Nonequispaced Fast Fourier Transform utilized to approximate matrix-vector products in a matrix-free manner. We accelerate the algorithm using a novel change of basis approach to account for the (known) smallest eigenvalue-eigenvector pair of the matrices involved, coupled with a simple but frequently very effective diagonal preconditioning approach. We present a number of theoretical results concerning the eigenvalues and predicted convergence behavior, and a range of numerical experiments which validate our solvers and use them to tackle parameter learning problems. These demonstrate that very large problems may be effectively and rapidly denoised with very low storage requirements on a computer. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 30 pages, 8 figures, 3 tables

Stability of stationary states for mean field models with multichromatic interaction potentials

Authors: Benedetta Bertoli, Benjamin D. Goddard, Grigorios A. Pavliotis

Abstract: We consider weakly interacting diffusions on the torus, for multichromatic interaction potentials. We consider interaction potentials that are not H-stable, leading to phase transitions in the mean field limit. We show that the mean field dynamics can exhibit multipeak stationary states, where the number of peaks is related to the number of nonzero Fourier modes of the interaction. We also conside… ▽ More We consider weakly interacting diffusions on the torus, for multichromatic interaction potentials. We consider interaction potentials that are not H-stable, leading to phase transitions in the mean field limit. We show that the mean field dynamics can exhibit multipeak stationary states, where the number of peaks is related to the number of nonzero Fourier modes of the interaction. We also consider the effect of a confining potential on the structure of non-uniform steady states. We approach the problem by means of analysis, perturbation theory and numerical simulations for the interacting particle systems and the PDEs. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 25 pages, 15 figures

Experimental and Theoretical Bulk Phase Diagram and Interfacial Tension of Ouzo

Authors: Andrew J. Archer, Benjamin D. Goddard, David N. Sibley, James T. Rawlings, Ross Broadhurst, Fouzia F. Ouali, David J. Fairhurst

Abstract: Ouzo is a well-known drink in Mediterranean countries, with ingredients water, alcohol and trans-anethole oil. The oil is insoluble in water, but completely soluble in alcohol, so when water is added to the spirit, the available alcohol is depleted and the mixture exhibits spontaneous emulsification. This process is commonly known as the louche or Ouzo effect. Although the phase boundaries of this… ▽ More Ouzo is a well-known drink in Mediterranean countries, with ingredients water, alcohol and trans-anethole oil. The oil is insoluble in water, but completely soluble in alcohol, so when water is added to the spirit, the available alcohol is depleted and the mixture exhibits spontaneous emulsification. This process is commonly known as the louche or Ouzo effect. Although the phase boundaries of this archetypal ternary mixture are well known, the properties of coexisting phases have not previously been studied. Here, we present a detailed experimental investigation into the phase behaviour, including tie-lines connecting coexisting phases, determination of the critical point (also called the plait point in ternary systems) and measurements of the surface tension and density for varying alcohol concentrations. Additionally, we present a theory for the thermodynamics and phase diagram of the system. With suitable selection of the interaction parameters, the theory captures nearly all features of the experimental work. This simple model can be used to determine both bulk and non-uniform (e.g. interfacial) properties, paving the way for a wide range of future applications of the model to ternary mixtures in general. We show how our accurate equilibrium phase diagram can be used to provide improved understanding of non-equilibrium phenomena. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 15 pages, 9 figures

Dynamic Density Functional Theory with Inertia and Background Flow

Authors: Rory D. Mills-Williams, Benjamin D. Goddard, Andrew J. Archer

Abstract: We present dynamic density functional theory (DDFT) incorporating general inhomogeneous, incompressible, time dependent background flows and inertia, describing externally driven passive colloidal systems out of equilibrium. We start by considering the underlying nonequilibrium Langevin dynamics, including the effect of the local velocity of the surrounding liquid bath, to obtain the nonlinear, no… ▽ More We present dynamic density functional theory (DDFT) incorporating general inhomogeneous, incompressible, time dependent background flows and inertia, describing externally driven passive colloidal systems out of equilibrium. We start by considering the underlying nonequilibrium Langevin dynamics, including the effect of the local velocity of the surrounding liquid bath, to obtain the nonlinear, nonlocal partial differential equations governing the evolution of the (coarse--grained) density and velocity fields describing the dynamics of colloids. Additionally, we show both with heuristic arguments, and by numerical solution, that our equations and solutions agree with existing DDFTs in the overdamped (high friction) limit. We provide numerical solutions that model the flow of hard spheres, in both unbounded and confined domains, and compare to previously--derived DDFTs with and without the background flow. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 13 pages, 3 figures

MSC Class: 82-10; 82C05; 35R09

Binding potential and wetting behaviour of binary liquid mixtures on surfaces

Authors: Mounirah Areshi, Dmitri Tseluiko, Uwe Thiele, Benjamin D. Goddard, Andrew J. Archer

Abstract: We present a theory for the interfacial wetting phase behaviour of binary liquid mixtures on rigid solid substrates, applicable to both miscible and immiscible mixtures. In particular, we calculate the binding potential as a function of the adsorptions, i.e. the excess amounts of each of the two liquids at the substrate. The binding potential fully describes the corresponding interfacial thermodyn… ▽ More We present a theory for the interfacial wetting phase behaviour of binary liquid mixtures on rigid solid substrates, applicable to both miscible and immiscible mixtures. In particular, we calculate the binding potential as a function of the adsorptions, i.e. the excess amounts of each of the two liquids at the substrate. The binding potential fully describes the corresponding interfacial thermodynamics. Our approach is based on classical density functional theory. Binary liquid mixtures can exhibit complex bulk phase behaviour, including both liquid-liquid and vapour-liquid phase separation, depending on the nature of the interactions between all the particles of the two different liquids, the temperature and the chemical potentials. Here we show that the interplay between the bulk phase behaviour of the mixture and the properties of the interactions with the substrate gives rise to a wide variety of interfacial phase behaviours, including mixing and demixing situations. We find situations where the final state is a coexistence of up to three different phases. We determine how the liquid density profiles close to the substrate change as the interaction parameters are varied and how these determine the form of the binding potential, which in certain cases can be a multi-valued function of the adsorptions. We also present profiles for sessile droplets of both miscible and immiscible binary liquids. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 22 pages, 23 figures

Journal ref: Phys. Rev. E 109, 024801 (2024)

Stability of nanoparticle laden aerosol liquid droplets

Authors: A. J. Archer, B. D. Goddard, R. Roth

Abstract: We develop a model for the thermodynamics and evaporation dynamics of aerosol droplets of a liquid such as water, surrounded by the gas. When the temperature and the chemical potential (or equivalently the humidity) are such that the vapour phase is the thermodynamic equilibrium state, then of course droplets of the pure liquid evaporate over a relatively short time. However, if the droplets also… ▽ More We develop a model for the thermodynamics and evaporation dynamics of aerosol droplets of a liquid such as water, surrounded by the gas. When the temperature and the chemical potential (or equivalently the humidity) are such that the vapour phase is the thermodynamic equilibrium state, then of course droplets of the pure liquid evaporate over a relatively short time. However, if the droplets also contain nanoparticles or any other non-volatile solute, then the droplets can become thermodynamically stable. We show that the equilibrium droplet size depends strongly on the amount and solubility of the nanoparticles within, i.e. on the nature of the particle interactions with the liquid, and of course also on the vapour temperature and chemical potential. We develop a simple thermodynamic model for such droplets and compare predictions with results from a lattice density functional theory that takes as input the same particle interaction properties, finding very good agreement. We also use dynamical density functional theory to study the evaporation/condensation dynamics of liquid from/to droplets as they equilibrate with the vapour, thereby demonstrating droplet stability. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 25 pages, 9 figures

Journal ref: J. Chem. Phys. 159, 194503 (2023)

A Numerical Optimisation Framework for Parameter Identification of the SIRD Model

Authors: Andrés Miniguano-Trujillo, John W. Pearson, Benjamin D. Goddard

Abstract: We consider a numerical framework tailored to identifying optimal parameters in the context of modelling disease propagation. Our focus is on understanding the behaviour of optimisation algorithms for such problems, where the dynamics are described by a system of ordinary differential equations associated with the epidemiological SIRD model. We examine properties of the solution operator and deter… ▽ More We consider a numerical framework tailored to identifying optimal parameters in the context of modelling disease propagation. Our focus is on understanding the behaviour of optimisation algorithms for such problems, where the dynamics are described by a system of ordinary differential equations associated with the epidemiological SIRD model. We examine properties of the solution operator and determine existence of optimal parameters for the problem considered. Further, first-order optimality conditions are derived, the solution of which provides a certificate of goodness of fit, which is not always guaranteed with parameter tuning techniques. We then propose strategies for the numerical solution of such problems, based on projected gradient descent, Fast Iterative Shrinkage-Thresholding Algorithm (FISTA), and limited memory BFGS trust region approaches. We carry out a thorough computational study for a range of problems of interest, determining the relative performance of these numerical methods. Our results provide insights into the efficacy of these strategies, contributing to ongoing research into optimising parameters for accurate and reliable disease spread modelling. Moreover, our approach paves the way for calibration of more intricate compartmental models. △ Less

Submitted 23 February, 2024; v1 submitted 17 October, 2023; originally announced October 2023.

Comments: 26 pages, 7 figures, 3 tables

Design and performance of the ENUBET monitored neutrino beam

Authors: F. Acerbi, I. Angelis, L. Bomben, M. Bonesini, F. Bramati, A. Branca, C. Brizzolari, G. Brunetti, M. Calviani, S. Capelli, S. Carturan, M. G. Catanesi, S. Cecchini, N. Charitonidis, F. Cindolo, G. Cogo, G. Collazuol, F. Dal Corso, C. Delogu, G. De Rosa, A. Falcone, B. Goddard, A. Gola, D. Guffanti, L. Halić , et al. (47 additional authors not shown)

Abstract: The ENUBET project is aimed at designing and experimentally demonstrating the concept of monitored neutrino beams. These novel beams are enhanced by an instrumented decay tunnel, whose detectors reconstruct large-angle charged leptons produced in the tunnel and give a direct estimate of the neutrino flux at the source. These facilities are thus the ideal tool for high-precision neutrino cross-sect… ▽ More The ENUBET project is aimed at designing and experimentally demonstrating the concept of monitored neutrino beams. These novel beams are enhanced by an instrumented decay tunnel, whose detectors reconstruct large-angle charged leptons produced in the tunnel and give a direct estimate of the neutrino flux at the source. These facilities are thus the ideal tool for high-precision neutrino cross-section measurements at the GeV scale because they offer superior control of beam systematics with respect to existing facilities. In this paper, we present the first end-to-end design of a monitored neutrino beam capable of monitoring lepton production at the single particle level. This goal is achieved by a new focusing system without magnetic horns, a 20 m normal-conducting transfer line for charge and momentum selection, and a 40 m tunnel instrumented with cost-effective particle detectors. Employing such a design, we show that percent precision in cross-section measurements can be achieved at the CERN SPS complex with existing neutrino detectors. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: 24 pages, 33 figures

On the study of slow-fast dynamics, when the fast process has multiple invariant measures

Authors: B. D. Goddard, M. Ottobre, K. J. Painter, I. Souttar

Abstract: Motivated by applications to mathematical biology, we study the averaging problem for slow-fast systems, {\em in the case in which the fast dynamics is a stochastic process with multiple invariant measures}. We consider both the case in which the fast process is decoupled from the slow process and the case in which the two components are fully coupled. We work in the setting in which the slow proc… ▽ More Motivated by applications to mathematical biology, we study the averaging problem for slow-fast systems, {\em in the case in which the fast dynamics is a stochastic process with multiple invariant measures}. We consider both the case in which the fast process is decoupled from the slow process and the case in which the two components are fully coupled. We work in the setting in which the slow process evolves according to an Ordinary Differential Equation (ODE) and the fast process is a continuous time Markov Process with finite state space and show that, in this setting, the limiting (averaged) dynamics can be described as a random ODE (that is, an ODE with random coefficients.) Keywords. Multiscale methods, Processes with multiple equilibria, Averaging, Collective Navigation, Interacting Piecewise Deterministic Markov Processes. △ Less

Submitted 15 August, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: 24 pages

MSC Class: 34C29; 34D05; 34E10; 34F05; 37H10; 60K35

Automatic setup of 18 MeV electron beamline using machine learning

Authors: Francesco Maria Velotti, Brennan Goddard, Verena Kain, Rebecca Ramjiawan, Giovanni Zevi Della Porta, Simon Hirlaender

Abstract: To improve the performance-critical stability and brightness of the electron bunch at injection into the proton-driven plasma wakefield at the AWAKE CERN experiment, automation approaches based on unsupervised Machine Learning (ML) were developed and deployed. Numerical optimisers were tested together with different model-free reinforcement learning agents. In order to avoid any bias, reinforcemen… ▽ More To improve the performance-critical stability and brightness of the electron bunch at injection into the proton-driven plasma wakefield at the AWAKE CERN experiment, automation approaches based on unsupervised Machine Learning (ML) were developed and deployed. Numerical optimisers were tested together with different model-free reinforcement learning agents. In order to avoid any bias, reinforcement learning agents have been trained also using a completely unsupervised state encoding using auto-encoders. To aid hyper-parameter selection, a full synthetic model of the beamline was constructed using a variational auto-encoder trained to generate surrogate data from equipment settings. This paper describes the novel approaches based on deep learning and reinforcement learning to aid the automatic setup of a low energy line, as the one used to deliver beam to the AWAKE facility. The results obtained with the different ML approaches, including automatic unsupervised feature extraction from images using computer vision are presented. The prospects for operational deployment and wider applicability are discussed. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Improving tobacco social contagion models using agent-based simulations on networks

Authors: Adarsh Prabhakaran, Valerio Restocchi, Benjamin D. Goddard

Abstract: Over the years, population-level tobacco control policies have considerably reduced smoking prevalence worldwide. However, the rate of decline of smoking prevalence is slowing down. Therefore, there is a need for models that capture the full complexity of the smoking epidemic. These models can then be used as test-beds to develop new policies to limit the spread of smoking. Current models of smoki… ▽ More Over the years, population-level tobacco control policies have considerably reduced smoking prevalence worldwide. However, the rate of decline of smoking prevalence is slowing down. Therefore, there is a need for models that capture the full complexity of the smoking epidemic. These models can then be used as test-beds to develop new policies to limit the spread of smoking. Current models of smoking dynamics mainly use ordinary differential equation (ODE) models, where studying the effect of an individual's contact network is challenging. They also do not consider all the interactions between individuals that can lead to changes in smoking behaviour, implying that they do not consider valuable information on the spread of smoking behaviour. In this context, we develop an agent-based model (ABM), calibrate and then validate it on historical trends observed in the US and UK. Our ABM considers spontaneous terms, interactions between agents, and the agent's contact network. To explore the effect of the underlying network on smoking dynamics, we test the ABM on six different networks, both synthetic and real-world. In addition, we also compare the ABM with an ODE model. Our results suggest that the dynamics from the ODE model are similar to the ABM only when the network structure is fully connected (FC). The FC network performs poorly in replicating the empirical trends in the data, while the real-world network best replicates it amongst the six networks. Further, when information on the real-world network is unavailable, our ABM on Lancichinetti-Fortunato-Radicchi benchmark networks (or networks with a similar average degree as the real-world network) can be used to model smoking behaviour. These results suggest that networks are essential for modelling smoking behaviour and that our ABM can be used to develop network-based intervention strategies and policies for tobacco control. △ Less

Submitted 17 July, 2022; originally announced July 2022.

Comments: Main paper: 17 pages, 5 Figures

MultiShape: A Spectral Element Method, with Applications to Dynamic Density Functional Theory and PDE-Constrained Optimization

Authors: Jonna C. Roden, Rory D. Mills-Williams, John W. Pearson, Benjamin D. Goddard

Abstract: A numerical framework is developed to solve various types of PDEs on complicated domains, including steady and time-dependent, non-linear and non-local PDEs, with different boundary conditions that can also include non-linear and non-local terms. This numerical framework, called MultiShape, is a class in Matlab, and the software is open source. We demonstrate that MultiShape is compatible with oth… ▽ More A numerical framework is developed to solve various types of PDEs on complicated domains, including steady and time-dependent, non-linear and non-local PDEs, with different boundary conditions that can also include non-linear and non-local terms. This numerical framework, called MultiShape, is a class in Matlab, and the software is open source. We demonstrate that MultiShape is compatible with other numerical methods, such as differential--algebraic equation solvers and optimization algorithms. The numerical implementation is designed to be user-friendly, with most of the set-up and computations done automatically by MultiShape and with intuitive operator definition, notation, and user-interface. Validation tests are presented, before we introduce three examples motivated by applications in Dynamic Density Functional Theory and PDE-constrained optimization, illustrating the versatility of the method. △ Less

Submitted 12 July, 2022; originally announced July 2022.

MSC Class: 35Q70; 35Q93; 65M70; 82C22; 82M22

Non-mean-field Vicsek-type models for collective behaviour

Authors: P. Buttà, B. Goddard, T. M. Hodgson, M. Ottobre, K. J. Painter

Abstract: We consider interacting particle dynamics with Vicsek type interactions, and their macroscopic PDE limit, in the non-mean-field regime; that is, we consider the case in which each particle/agent in the system interacts only with a prescribed subset of the particles in the system (for example, those within a certain distance). In this non-mean-field regime the influence between agents (i.e. the int… ▽ More We consider interacting particle dynamics with Vicsek type interactions, and their macroscopic PDE limit, in the non-mean-field regime; that is, we consider the case in which each particle/agent in the system interacts only with a prescribed subset of the particles in the system (for example, those within a certain distance). In this non-mean-field regime the influence between agents (i.e. the interaction term) can be normalised either by the total number of agents in the system (\textit{global scaling}) or by the number of agents with which the particle is effectively interacting (\textit{local scaling}). We compare the behaviour of the globally scaled and the locally scaled systems in many respects, considering for each scaling both the PDE and the corresponding particle model. In particular we observe that both the locally and globally scaled particle system exhibit pattern formation (i.e. formation of travelling-wave-like solutions) within certain parameter regimes, and generally display similar dynamics. The same is not true of the corresponding PDE models. Indeed, while both PDE models have multiple stationary states, for the globally scaled PDE such (space-homogeneous) equilibria are unstable for certain parameter regimes, with the instability leading to travelling wave solutions, while they are always stable for the locally scaled one, which never produces travelling waves. This observation is based on a careful numerical study of the model, supported by further analysis. △ Less

Submitted 17 May, 2022; originally announced June 2022.

Comments: 44 pages

Poisson Equations with locally-Lipschitz coefficients and Uniform in Time Averaging for Stochastic Differential Equations via Strong Exponential Stability

Authors: Dan Crisan, Paul Dobson, Ben Goddard, Michela Ottobre, Iain Souttar

Abstract: We study averaging for Stochastic Differential Equations (SDEs) and Poisson equations. We succeed in obtaining a uniform in time (UiT) averaging result, with a rate, for fully coupled SDE models with super-linearly growing coefficients. This is the main result of this paper and is, to the best of our knowledge, the first UiT multiscale result with a rate. Very few UiT averaging results exist in th… ▽ More We study averaging for Stochastic Differential Equations (SDEs) and Poisson equations. We succeed in obtaining a uniform in time (UiT) averaging result, with a rate, for fully coupled SDE models with super-linearly growing coefficients. This is the main result of this paper and is, to the best of our knowledge, the first UiT multiscale result with a rate. Very few UiT averaging results exist in the literature, and they almost exclusively apply to multiscale systems of Ordinary Differential Equations. Among these few, none of those we are aware of comes with a rate of convergence. The UiT nature of this result and the rate of convergence given by the main theorem, make it important as theoretical underpinning for a range of applications, such as applications to statistical methodology, molecular dynamics etc. Key to obtaining both our UiT averaging result and to enable dealing with the super-linear growth of the coefficients is conquering exponential decay in time of the space-derivatives of appropriate Markov semigroups. We refer to this property as being Strongly Exponentially Stable (SES). The analytic approach to proving averaging results we take requires studying a family of Poisson problems associated with the generator of the (fast component of the) SDE dynamics. The study of Poisson equations in non-compact state space is notoriously difficult, with current literature mostly covering the case when the coefficients of the Partial Differential Equation (PDE) are either bounded or satisfy linear growth assumptions. In this paper we treat Poisson equations on non-compact state spaces for coefficients that can grow super-linearly. We demonstrate how SES can be employed not only to prove the UiT result for the slow-fast system but also to overcome some of the technical hurdles in the analysis of Poisson problems, which is of independent interest as well. △ Less

Submitted 5 April, 2024; v1 submitted 6 April, 2022; originally announced April 2022.

Comments: 70 pages, 3 figures

MSC Class: 60J60; 60H10; 35B30; 34K33; 34D20; 47D07; 65M75

The Mean Field Fokker-Planck Equation with Nonlinear No-flux Boundary Conditions

Authors: R. D. Mills-Williams, B. D. Goddard, G. A. Pavliotis

Abstract: We consider the mean field Fokker-Planck equation subject to nonlinear no-flux boundary conditions, which necessarily arise when subjecting a system of Brownian particles interacting via a pair potential in a bounded domain. With the additional presence of an external potential $V_1$, we show, by analysing the linearised Fokker-Planck operator, that the spectral properties of the equilibrium densi… ▽ More We consider the mean field Fokker-Planck equation subject to nonlinear no-flux boundary conditions, which necessarily arise when subjecting a system of Brownian particles interacting via a pair potential in a bounded domain. With the additional presence of an external potential $V_1$, we show, by analysing the linearised Fokker-Planck operator, that the spectral properties of the equilibrium densities can differ considerably when compared with previous studies, e.g., with periodic boundary conditions. Amongst other mean field models of complex many-body particle systems, we present numerical experiments encompassing in a wide range of physical applications, including: generalised exponential models (Gaussian, Morse); a Kuramoto model, for noisy coupled oscillators; and an Onsager model for liquid crystals. We showcase our results by using the numerical methods developed in the pseudospectral collocation scheme 2DChebClass. △ Less

Submitted 29 March, 2022; v1 submitted 26 March, 2022; originally announced March 2022.

Comments: The paper was withdrawn by the authors due to a critical error in the definition of the system of Langevin equations

Changing the flow profile and resulting drying pattern of dispersion droplets via contact angle modification

Authors: Carmen Morcillo Perez, Marcel Rey, Benjamin D. Goddard, Job H. J. Thijssen

Abstract: Spilling tea or coffee leads to a tell-tale circular stain after the drying of the droplet. This phenomenon was termed after the latter example as the "coffee ring effect". The evaporation of suspension droplets is a complex physical process, and prediction and control over particle deposit patterns obtained from sessile droplet evaporation are essential for many industrial processes such as ink-j… ▽ More Spilling tea or coffee leads to a tell-tale circular stain after the drying of the droplet. This phenomenon was termed after the latter example as the "coffee ring effect". The evaporation of suspension droplets is a complex physical process, and prediction and control over particle deposit patterns obtained from sessile droplet evaporation are essential for many industrial processes such as ink-jet printing or crop-care applications. In this article, we present a systematic investigation of the effect of surface wettability on the evaporation dynamics of a particle-laden droplet, including the effect on the contact line stick-slip, the hydrodynamic flow of the suspended particles and the resulting particle deposit after evaporation. We tune the wettability of glass slides using silanisation and quantify the internal flow during the evaporation by tracking fluorescent tracer particles. We find that the internal flow shifts from a predominantly outward flow towards the contact line for low contact angles to an inward flow for large contact angles. Additionally, the corresponding deposit gradually changes from the typical coffee-ring to a central stain upon increasing the hydrophobicity of the substrate. Last, we corroborate these experimental findings with dynamic density functional theory, modelling the droplet evaporation process and stick-slip behaviour of the contact line. Our investigation suggests that the wettability of the substrate can substantially alter hydrodynamic flow within drying droplets and therefore the resulting particle deposit. △ Less

Submitted 27 December, 2021; v1 submitted 31 October, 2021; originally announced November 2021.

Comments: 9 pages, 6 figures

Sensitivity of the SHiP experiment to dark photons decaying to a pair of charged particles

Authors: SHiP Collaboration, C. Ahdida, A. Akmete, R. Albanese, A. Alexandrov, A. Anokhina, S. Aoki, G. Arduini, E. Atkin, N. Azorskiy, J. J. Back, A. Bagulya, F. Baaltasar Dos Santos, A. Baranov, F. Bardou, G. J. Barker, M. Battistin, J. Bauche, A. Bay, V. Bayliss, G. Bencivenni, A. Y. Berdnikov, Y. A. Berdnikov, M. Bertani, C. Betancourt , et al. (309 additional authors not shown)

Abstract: Dark photons are hypothetical massive vector particles that could mix with ordinary photons. The simplest theoretical model is fully characterised by only two parameters: the mass of the dark photon m$_{γ^{\mathrm{D}}}$ and its mixing parameter with the photon, $\varepsilon$. The sensitivity of the SHiP detector is reviewed for dark photons in the mass range between 0.002 and 10 GeV. Different pro… ▽ More Dark photons are hypothetical massive vector particles that could mix with ordinary photons. The simplest theoretical model is fully characterised by only two parameters: the mass of the dark photon m$_{γ^{\mathrm{D}}}$ and its mixing parameter with the photon, $\varepsilon$. The sensitivity of the SHiP detector is reviewed for dark photons in the mass range between 0.002 and 10 GeV. Different production mechanisms are simulated, with the dark photons decaying to pairs of visible fermions, including both leptons and quarks. Exclusion contours are presented and compared with those of past experiments. The SHiP detector is expected to have a unique sensitivity for m$_{γ^{\mathrm{D}}}$ ranging between 0.8 and 3.3$^{+0.2}_{-0.5}$ GeV, and $\varepsilon^2$ ranging between $10^{-11}$ and $10^{-17}$. △ Less

Submitted 1 March, 2021; v1 submitted 10 November, 2020; originally announced November 2020.

Experimental study of extended timescale dynamics of a plasma wakefield driven by a self-modulated proton bunch

Authors: J. Chappell, E. Adli, R. Agnello, M. Aladi, Y. Andrebe, O. Apsimon, R. Apsimon, A. -M. Bachmann, M. A. Baistrukov, F. Batsch, M. Bergamaschi, P. Blanchard, P. N. Burrows, B. Buttenschön, A. Caldwell, E. Chevallay, M. Chung, D. A. Cooke, H. Damerau, C. Davut, G. Demeter, L. H. Deubner, A. Dexter, G. P. Djotyan, S. Doebert , et al. (74 additional authors not shown)

Abstract: Plasma wakefield dynamics over timescales up to 800 ps, approximately 100 plasma periods, are studied experimentally at the Advanced Wakefield Experiment (AWAKE). The development of the longitudinal wakefield amplitude driven by a self-modulated proton bunch is measured using the external injection of witness electrons that sample the fields. In simulation, resonant excitation of the wakefield cau… ▽ More Plasma wakefield dynamics over timescales up to 800 ps, approximately 100 plasma periods, are studied experimentally at the Advanced Wakefield Experiment (AWAKE). The development of the longitudinal wakefield amplitude driven by a self-modulated proton bunch is measured using the external injection of witness electrons that sample the fields. In simulation, resonant excitation of the wakefield causes plasma electron trajectory crossing, resulting in the development of a potential outside the plasma boundary as electrons are transversely ejected. Trends consistent with the presence of this potential are experimentally measured and their dependence on wakefield amplitude are studied via seed laser timing scans and electron injection delay scans. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 11 pages, 8 figures

PDE-Constrained Optimization Models and Pseudospectral Methods for Multiscale Particle Dynamics

Authors: Mildred Aduamoah, Benjamin D. Goddard, John W. Pearson, Jonna C. Roden

Abstract: We derive novel algorithms for optimization problems constrained by partial differential equations describing multiscale particle dynamics, including non-local integral terms representing interactions between particles. In particular, we investigate problems where the control acts as an advection 'flow' vector or a source term of the partial differential equation, and the constraint is equipped wi… ▽ More We derive novel algorithms for optimization problems constrained by partial differential equations describing multiscale particle dynamics, including non-local integral terms representing interactions between particles. In particular, we investigate problems where the control acts as an advection 'flow' vector or a source term of the partial differential equation, and the constraint is equipped with boundary conditions of Dirichlet or no-flux type. After deriving continuous first-order optimality conditions for such problems, we solve the resulting systems by developing a link with computational methods for statistical mechanics, deriving pseudospectral methods in both space and time variables, and utilizing variants of existing fixed point methods. Numerical experiments indicate the effectiveness of our approach for a range of problem set-ups, boundary conditions, as well as regularization and model parameters. △ Less

Submitted 8 September, 2021; v1 submitted 21 September, 2020; originally announced September 2020.

MSC Class: 49M41

A primary electron beam facility at CERN -- eSPS Conceptual design report

Authors: M. Aicheler, T. Akesson, F. Antoniou, A. Arnalich, P. A. Arrutia Sota, P. Bettencourt Moniz Cabral, D. Bozzini, M. Brugger, O. Brunner, P. N. Burrows, R. Calaga, M. J. Capstick, R. Corsini, S. Doebert, L. A. Dougherty, Y. Dutheil, L. A. Dyks, O. Etisken, L. Evans, A. Farricker, R. Fernandez Ortega, M. A. Fraser, J. Gall, S. J. Gessner, B. Goddard , et al. (30 additional authors not shown)

Abstract: The design of a primary electron beam facility at CERN is described. The study has been carried out within the framework of the wider Physics Beyond Colliders study. It re-enables the Super Proton Synchrotron (SPS) as an electron accelerator, and leverages the development invested in Compact Linear Collider (CLIC) technology for its injector and as an accelerator research and development infrastru… ▽ More The design of a primary electron beam facility at CERN is described. The study has been carried out within the framework of the wider Physics Beyond Colliders study. It re-enables the Super Proton Synchrotron (SPS) as an electron accelerator, and leverages the development invested in Compact Linear Collider (CLIC) technology for its injector and as an accelerator research and development infrastructure. The facility would be relevant for several of the key priorities in the 2020 update of the European Strategy for Particle Physics, such as an electron-positron Higgs factory, accelerator R\&D, dark sector physics, and neutrino physics. In addition, it could serve experiments in nuclear physics. The electron beam delivered by this facility would provide access to light dark matter production significantly beyond the targets predicted by a thermal dark matter origin, and for natures of dark matter particles that are not accessible by direct detection experiments. It would also enable electro-nuclear measurements crucial for precise modelling the energy dependence of neutrino-nucleus interactions, which is needed to precisely measure neutrino oscillations as a function of energy. The implementation of the facility is the natural next step in the development of X-band high-gradient acceleration technology, a key technology for compact and cost-effective electron/positron linacs. It would also become the only facility with multi-GeV drive bunches and truly independent electron witness bunches for plasma wakefield acceleration. A second phase capable to deliver positron witness bunches would make it a complete facility for plasma wakefield collider studies. [...] △ Less

Submitted 21 December, 2020; v1 submitted 15 September, 2020; originally announced September 2020.

Noisy bounded confidence models for opinion dynamics: the effect of boundary conditions on phase transitions

Authors: Benjamin D. Goddard, Beth Gooding, Grigoris A. Pavliotis, Hannah Short

Abstract: We study SDE and PDE models for opinion dynamics under bounded confidence, for a range of different boundary conditions, with and without the inclusion of a radical population. We perform exhaustive numerical studies with pseudospectral methods to determine the effects of the boundary conditions, suggesting that the no-flux case most faithfully reproduce the underlying mechanisms in the associated… ▽ More We study SDE and PDE models for opinion dynamics under bounded confidence, for a range of different boundary conditions, with and without the inclusion of a radical population. We perform exhaustive numerical studies with pseudospectral methods to determine the effects of the boundary conditions, suggesting that the no-flux case most faithfully reproduce the underlying mechanisms in the associated deterministic models of Hegselmann and Krause. We also compare the SDE and PDE models, and use tools from analysis to study phase transitions, including a systematic description of an order parameter. △ Less

Submitted 23 September, 2021; v1 submitted 7 September, 2020; originally announced September 2020.

MSC Class: 82C22

Proton beam defocusing in AWAKE: comparison of simulations and measurements

Authors: A. A. Gorn, M. Turner, E. Adli, R. Agnello, M. Aladi, Y. Andrebe, O. Apsimon, R. Apsimon, A. -M. Bachmann, M. A. Baistrukov, F. Batsch, M. Bergamaschi, P. Blanchard, P. N. Burrows, B. Buttenschon, A. Caldwell, J. Chappell, E. Chevallay, M. Chung, D. A. Cooke, H. Damerau, C. Davut, G. Demeter, L. H. Deubner, A. Dexter , et al. (74 additional authors not shown)

Abstract: In 2017, AWAKE demonstrated the seeded self-modulation (SSM) of a 400 GeV proton beam from the Super Proton Synchrotron (SPS) at CERN. The angular distribution of the protons deflected due to SSM is a quantitative measure of the process, which agrees with simulations by the two-dimensional (axisymmetric) particle-in-cell code LCODE. Agreement is achieved for beam populations between $10^{11}$ and… ▽ More In 2017, AWAKE demonstrated the seeded self-modulation (SSM) of a 400 GeV proton beam from the Super Proton Synchrotron (SPS) at CERN. The angular distribution of the protons deflected due to SSM is a quantitative measure of the process, which agrees with simulations by the two-dimensional (axisymmetric) particle-in-cell code LCODE. Agreement is achieved for beam populations between $10^{11}$ and $3 \times 10^{11}$ particles, various plasma density gradients ($-20 ÷20\%$) and two plasma densities ($2\times 10^{14} \text{cm}^{-3}$ and $7 \times 10^{14} \text{cm}^{-3}$). The agreement is reached only in the case of a wide enough simulation box (at least five plasma wavelengths). △ Less

Submitted 26 August, 2020; originally announced August 2020.

Comments: 8 pages, 9 figures, 1 table

The Large Hadron-Electron Collider at the HL-LHC

Authors: P. Agostini, H. Aksakal, S. Alekhin, P. P. Allport, N. Andari, K. D. J. Andre, D. Angal-Kalinin, S. Antusch, L. Aperio Bella, L. Apolinario, R. Apsimon, A. Apyan, G. Arduini, V. Ari, A. Armbruster, N. Armesto, B. Auchmann, K. Aulenbacher, G. Azuelos, S. Backovic, I. Bailey, S. Bailey, F. Balli, S. Behera, O. Behnke , et al. (312 additional authors not shown)

Abstract: The Large Hadron electron Collider (LHeC) is designed to move the field of deep inelastic scattering (DIS) to the energy and intensity frontier of particle physics. Exploiting energy recovery technology, it collides a novel, intense electron beam with a proton or ion beam from the High Luminosity--Large Hadron Collider (HL-LHC). The accelerator and interaction region are designed for concurrent el… ▽ More The Large Hadron electron Collider (LHeC) is designed to move the field of deep inelastic scattering (DIS) to the energy and intensity frontier of particle physics. Exploiting energy recovery technology, it collides a novel, intense electron beam with a proton or ion beam from the High Luminosity--Large Hadron Collider (HL-LHC). The accelerator and interaction region are designed for concurrent electron-proton and proton-proton operation. This report represents an update of the Conceptual Design Report (CDR) of the LHeC, published in 2012. It comprises new results on parton structure of the proton and heavier nuclei, QCD dynamics, electroweak and top-quark physics. It is shown how the LHeC will open a new chapter of nuclear particle physics in extending the accessible kinematic range in lepton-nucleus scattering by several orders of magnitude. Due to enhanced luminosity, large energy and the cleanliness of the hadronic final states, the LHeC has a strong Higgs physics programme and its own discovery potential for new physics. Building on the 2012 CDR, the report represents a detailed updated design of the energy recovery electron linac (ERL) including new lattice, magnet, superconducting radio frequency technology and further components. Challenges of energy recovery are described and the lower energy, high current, 3-turn ERL facility, PERLE at Orsay, is presented which uses the LHeC characteristics serving as a development facility for the design and operation of the LHeC. An updated detector design is presented corresponding to the acceptance, resolution and calibration goals which arise from the Higgs and parton density function physics programmes. The paper also presents novel results on the Future Circular Collider in electron-hadron mode, FCC-eh, which utilises the same ERL technology to further extend the reach of DIS to even higher centre-of-mass energies. △ Less

Submitted 12 April, 2021; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: 373 pages, many figures, to be published by J. Phys. G

Report number: CERN-ACC-Note-2020-0002

Journal ref: J.Phys.G 48 (2021) 11, 110501

Improved Short-Baseline Neutrino Oscillation Search and Energy Spectrum Measurement with the PROSPECT Experiment at HFIR

Authors: M. Andriamirado, A. B. Balantekin, H. R. Band, C. D. Bass, D. E. Bergeron, D. Berish, N. S. Bowden, J. P. Brodsky, C. D. Bryan, T. Classen, A. J. Conant, G. Deichert, M. V. Diwan, M. J. Dolinski, A. Erickson, B. T. Foust, J. K. Gaison, A. Galindo-Uribarri, C. E. Gilbert, B. W. Goddard, B. T. Hackett, S. Hans, A. B. Hansell, K. M. Heeger, B. Heffron , et al. (39 additional authors not shown)

Abstract: We present a detailed report on sterile neutrino oscillation and U-235 antineutrino energy spectrum measurement results from the PROSPECT experiment at the highly enriched High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory. In 96 calendar days of data taken at an average baseline distance of 7.9 m from the center of the 85 MW HFIR core, the PROSPECT detector has observed more than 5… ▽ More We present a detailed report on sterile neutrino oscillation and U-235 antineutrino energy spectrum measurement results from the PROSPECT experiment at the highly enriched High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory. In 96 calendar days of data taken at an average baseline distance of 7.9 m from the center of the 85 MW HFIR core, the PROSPECT detector has observed more than 50,000 interactions of antineutrinos produced in beta decays of U-235 fission products. New limits on the oscillation of antineutrinos to light sterile neutrinos have been set by comparing the detected energy spectra of ten reactor-detector baselines between 6.7 and 9.2 meters. Measured differences in energy spectra between baselines show no statistically significant indication of antineutrinos to sterile neutrino oscillation and disfavor the Reactor Antineutrino Anomaly best-fit point at the 2.5$σ$ confidence level. The reported U-235 antineutrino energy spectrum measurement shows excellent agreement with energy spectrum models generated via conversion of the measured U-235 beta spectrum, with a $χ^2$/DOF of 31/31. PROSPECT is able to disfavor at 2.4$σ$ confidence level the hypothesis that U-235 antineutrinos are solely responsible for spectrum discrepancies between model and data obtained at commercial reactor cores. A data-model deviation in PROSPECT similar to that observed by commercial core experiments is preferred with respect to no observed deviation, at a 2.2$σ$ confidence level. △ Less

Submitted 1 July, 2020; v1 submitted 19 June, 2020; originally announced June 2020.

Comments: 42 pages, 52 Figures. Submitted to Phys. Rev. D. Supplementary Material Included

Journal ref: Phys. Rev. D 103, 032001 (2021)

Modelling inelastic granular media using Dynamical Density Functional Theory

Authors: B. D. Goddard, T. D. Hurst, R. Ocone

Abstract: We construct a new mesoscopic model for granular media using Dynamical Density Functional Theory (DDFT). The model includes both a collision operator to incorporate inelasticity and the Helmholtz free energy functional to account for external potentials, interparticle interactions and volume exclusion. We use statistical data from event-driven microscopic simulations to determine the parameters no… ▽ More We construct a new mesoscopic model for granular media using Dynamical Density Functional Theory (DDFT). The model includes both a collision operator to incorporate inelasticity and the Helmholtz free energy functional to account for external potentials, interparticle interactions and volume exclusion. We use statistical data from event-driven microscopic simulations to determine the parameters not given analytically by the closure relations used to derive the DDFT. We numerically demonstrate the crucial effects of each term in the DDFT, and the importance of including an accurately parametrised pair correlation function. △ Less

Submitted 16 March, 2020; originally announced March 2020.

Well-Posedness and Equilibrium Behaviour of Overdamped Dynamic Density Functional Theory

Authors: B. D. Goddard, R. D. Mills-Williams, M. Ottobre, G. Pavliotis

Abstract: We establish the global well-posedness of overdamped dynamic density functional theory (DDFT): a nonlinear, nonlocal integro-partial differential equation used in statistical mechanical models of colloidal fluids, and other applications including nonlinear reaction-diffusion systems and opinion dynamics. With nonlinear no-flux boundary conditions, we determine the existence and uniqueness of the w… ▽ More We establish the global well-posedness of overdamped dynamic density functional theory (DDFT): a nonlinear, nonlocal integro-partial differential equation used in statistical mechanical models of colloidal fluids, and other applications including nonlinear reaction-diffusion systems and opinion dynamics. With nonlinear no-flux boundary conditions, we determine the existence and uniqueness of the weak density and flux, subject to two-body hydrodynamic interactions (HI). We also show that the density is Lyapunov stable with respect to the usual (Helmholtz) free energy functional. Principally, this is done by rewriting the dynamics for the density in an implicit gradient flow form, resembling the classical Smoluchowski equation but with spatially inhomogeneous diffusion and advection tensors. We also rigorously show that the stationary density is independent of the HI tensors, and prove exponentially fast convergence to equilibrium. △ Less

Submitted 14 September, 2021; v1 submitted 26 February, 2020; originally announced February 2020.

Comments: 33 pages

SND@LHC

Authors: SHiP Collaboration, C. Ahdida, A. Akmete, R. Albanese, A. Alexandrov, M. Andreini, A. Anokhina, S. Aoki, G. Arduini, E. Atkin, N. Azorskiy, J. J. Back, A. Bagulya, F. Baaltasar Dos Santos, A. Baranov, F. Bardou, G. J. Barker, M. Battistin, J. Bauche, A. Bay, V. Bayliss, G. Bencivenni, A. Y. Berdnikov, Y. A. Berdnikov, M. Bertani , et al. (319 additional authors not shown)

Abstract: We propose to build and operate a detector that, for the first time, will measure the process $pp\toνX$ at the LHC and search for feebly interacting particles (FIPs) in an unexplored domain. The TI18 tunnel has been identified as a suitable site to perform these measurements due to very low machine-induced background. The detector will be off-axis with respect to the ATLAS interaction point (IP1)… ▽ More We propose to build and operate a detector that, for the first time, will measure the process $pp\toνX$ at the LHC and search for feebly interacting particles (FIPs) in an unexplored domain. The TI18 tunnel has been identified as a suitable site to perform these measurements due to very low machine-induced background. The detector will be off-axis with respect to the ATLAS interaction point (IP1) and, given the pseudo-rapidity range accessible, the corresponding neutrinos will mostly come from charm decays: the proposed experiment will thus make the first test of the heavy flavour production in a pseudo-rapidity range that is not accessible by the current LHC detectors. In order to efficiently reconstruct neutrino interactions and identify their flavour, the detector will combine in the target region nuclear emulsion technology with scintillating fibre tracking layers and it will adopt a muon identification system based on scintillating bars that will also play the role of a hadronic calorimeter. The time of flight measurement will be achieved thanks to a dedicated timing detector. The detector will be a small-scale prototype of the scattering and neutrino detector (SND) of the SHiP experiment: the operation of this detector will provide an important test of the neutrino reconstruction in a high occupancy environment. △ Less

Submitted 20 February, 2020; originally announced February 2020.

Comments: Letter of Intent

Report number: CERN-LHCC-2020-002, LHCC-I-035

Predicting the Trajectory of a Relativistic Electron Beam for External Injection in Plasma Wakefields

Authors: Felipe Peña Asmus, Francesco Maria Velotti, Marlene Turner, Spencer Gessner, Mikhail Martyanov, Chiara Bracco, Brennan Goddard, Patric Muggli

Abstract: We use beam position measurements over the first part of the AWAKE electron beamline, together with beamline modeling, to deduce the beam average momentum and to predict the beam position in the second part of the beamline. Results show that using only the first five beam position monitors leads to much larger differences between predicted and measured positions at the last two monitors than when… ▽ More We use beam position measurements over the first part of the AWAKE electron beamline, together with beamline modeling, to deduce the beam average momentum and to predict the beam position in the second part of the beamline. Results show that using only the first five beam position monitors leads to much larger differences between predicted and measured positions at the last two monitors than when using the first eight beam position monitors. These last two positions can in principle be used with ballistic calculations to predict the parameters of closest approach of the electron bunch with the proton beam. In external injection experiments of the electron bunch into plasma wakefields driven by the proton bunch, only the first five beam position monitors measurements remain un-affected by the presence of the much higher charge proton bunch. Results with eight beam position monitors show the prediction method works in principle to determine electron and proton beams closest approach within the wakefields width ($<$1\,mm), corresponding to injection of electrons into the wakefields. Using five beam position monitors is not sufficient. △ Less

Submitted 2 February, 2020; originally announced February 2020.

Comments: seven pages, five figures, submitted for EAAC 2019 Proceedings

SPS Beam Dump Facility -- Comprehensive Design Study

Authors: C. Ahdida, R. G. Alia, G. Arduini, A. Arnalich, P. Avigni, F. Bardou, M. Battistin, J. Bauche, M. Brugger, J. Busom, M. Calviani, M. Casolino, N. Colonna, L. Dougherty, Y. Dutheil, E. Fornasiere, M. A. Fraser, L. Gatignon, J. Gall, S. Gilardoni, B. Goddard, J-L. Grenard, D. Grenier, C. Hessler, R. Jacobsson , et al. (23 additional authors not shown)

Abstract: The proposed Beam Dump Facility (BDF) is foreseen to be located at the North Area of the SPS. It is designed to be able to serve both beam dump like and fixed target experiments. The SPS and the new facility would offer unique possibilities to enter a new era of exploration at the intensity frontier. Possible options include searches for very weakly interacting particles predicted by Hidden Sector… ▽ More The proposed Beam Dump Facility (BDF) is foreseen to be located at the North Area of the SPS. It is designed to be able to serve both beam dump like and fixed target experiments. The SPS and the new facility would offer unique possibilities to enter a new era of exploration at the intensity frontier. Possible options include searches for very weakly interacting particles predicted by Hidden Sector models, and flavour physics measurements. In the first instance, exploitation of the facility, in beam dump mode, is envisaged to be for the Search for Hidden Particle (SHiP) experiment. Following the first evaluation of the BDF in 2014-2016, CERN management launched a Comprehensive Design Study over three years for the facility. The BDF study team has since executed an in-depth feasibility study of proton delivery to target, the target complex, and the underground experimental area, including prototyping of key sub-systems and evaluations of the radiological aspects and safety. A first iteration of detailed integration and civil engineering studies have been performed in order to produce a realistic schedule and cost. This document gives a detailed overview of the proposed facility together with the results of the studies, and draws up a possible road map for a three-year Technical Design Report phase, followed by a 5 to 6 year construction phase. △ Less

Submitted 13 December, 2019; originally announced December 2019.

Comments: 589 pages

Report number: CERN-PBC-REPORT-2019-005; CERN-2020-002

A derivation of the Liouville equation for hard particle dynamics with non-conservative interactions

Authors: T. D. Hurst, B. D. Goddard, M. Wilkinson

Abstract: The Liouville equation is of fundamental importance in the derivation of continuum models for physical systems which are approximated by interacting particles. However, when particles undergo instantaneous interactions such as collisions, the derivation of the Liouville equation must be adapted to exclude non-physical particle positions, and include the effect of instantaneous interactions. We pre… ▽ More The Liouville equation is of fundamental importance in the derivation of continuum models for physical systems which are approximated by interacting particles. However, when particles undergo instantaneous interactions such as collisions, the derivation of the Liouville equation must be adapted to exclude non-physical particle positions, and include the effect of instantaneous interactions. We present the weak formulation of the Liouville equation for interacting particles with general particle dynamics and interactions, and discuss the results using an example. △ Less

Submitted 15 October, 2019; originally announced October 2019.

A primary electron beam facility at CERN

Authors: T. Åkesson, R. Corsini, Y. Dutheil, L. Evans, B. Goddard, A. Grudiev, A. Latina, Y. Papaphilippou, S. Stapnes

Abstract: This paper describes the concept of a primary electron beam facility at CERN, to be used for dark gauge force and light dark matter searches. The electron beam is produced in three stages: A Linac accelerates electrons from a photo-cathode up to 3.5 GeV. This beam is injected into the Super Proton Synchrotron, SPS, and accelerated up to a maximum energy of 16 GeV. Finally, the accelerated beam is… ▽ More This paper describes the concept of a primary electron beam facility at CERN, to be used for dark gauge force and light dark matter searches. The electron beam is produced in three stages: A Linac accelerates electrons from a photo-cathode up to 3.5 GeV. This beam is injected into the Super Proton Synchrotron, SPS, and accelerated up to a maximum energy of 16 GeV. Finally, the accelerated beam is slowly extracted to an experiment, possibly followed by a fast dump of the remaining electrons to another beamline. The beam parameters are optimized using the requirements of the Light Dark Matter eXperiment (LDMX) as benchmark. △ Less

Submitted 18 May, 2019; originally announced May 2019.

Comments: 3 pages, 3 figures

Design of a high power production target for the Beam Dump Facility at CERN

Authors: E. Lopez Sola, M. Calviani, P. Avigni, M. Battistin, J. Busom Descarrega, J. Canhoto Espadanal, M. A. Fraser, S. Gilardoni, B. Goddard, D. Grenier, R. Jacobsson, K. Kershaw, M. Lamont, A. Perillo-Marcone, M. Pandey, B. Riffaud, S. Sgobba, V. Vlachoudis, L. Zuccalli

Abstract: The Beam Dump Facility (BDF) project is a proposed general-purpose facility at CERN, dedicated to beam dump and fixed target experiments. In its initial phase, the facility is foreseen to be exploited by the Search for Hidden Particles (SHiP) experiment. Physics requirements call for a pulsed 400 GeV/c proton beam as well as the highest possible number of protons on target (POT) each year of opera… ▽ More The Beam Dump Facility (BDF) project is a proposed general-purpose facility at CERN, dedicated to beam dump and fixed target experiments. In its initial phase, the facility is foreseen to be exploited by the Search for Hidden Particles (SHiP) experiment. Physics requirements call for a pulsed 400 GeV/c proton beam as well as the highest possible number of protons on target (POT) each year of operation, in order to search for feebly interacting particles. The target/dump assembly lies at the heart of the facility, with the aim of safely absorbing the full high intensity Super Proton Synchrotron (SPS) beam, while maximizing the production of charmed and beauty mesons. High-Z materials are required for the target/dump, in order to have the shortest possible absorber and reduce muon background for the downstream experiment. The high average power deposited on target (305 kW) creates a challenge for heat removal. During the BDF facility Comprehensive Design Study (CDS), launched by CERN in 2016, extensive studies have been carried out in order to define and assess the target assembly design. These studies are described in the present contribution, which details the proposed design of the BDF production target, as well as the material selection process and the optimization of the target configuration and beam dilution. One of the specific challenges and novelty of this work is the need to consider new target materials, such as a molybdenum alloy (TZM) as core absorbing material and Ta2.5W as cladding. Thermo-structural and fluid dynamics calculations have been performed to evaluate the reliability of the target and its cooling system under beam operation. In the framework of the target comprehensive design, a preliminary mechanical design of the full target assembly has also been carried out, assessing the feasibility of the whole target system. △ Less

Submitted 30 October, 2019; v1 submitted 5 April, 2019; originally announced April 2019.

Comments: 17 pages, 18 figures

Journal ref: Phys. Rev. Accel. Beams 22, 113001 (2019)

The ENUBET narrow band neutrino beam

Authors: ENUBET Collaboration, M. Tenti, F. Acerbi, G. Ballerini, M. Bonesini, C. Brizzolari, G. Brunetti M. Calviani, S. Carturan, M. G. Catanesi, S. Cecchini, F. Cindolo, G. Collazuol, E. Conti F. Dal Corso, G. De Rosa, C. Delogu, A. Falcone, B. Goddard, A. Gola, R. A. Intonti, C. Jollet, V. Kain, B. Klicek, Y. Kudenko, M. Laveder, A. Longhin , et al. (32 additional authors not shown)

Abstract: The narrow band beam of ENUBET is the first implementation of the "monitored neutrino beam" technique proposed in 2015. ENUBET has been designed to monitor lepton production in the decay tunnel of neutrino beams and to provide a 1% measurement of the neutrino flux at source. In particular, the three body semi-leptonic decay of kaons monitored by large angle positron production offers a fully contr… ▽ More The narrow band beam of ENUBET is the first implementation of the "monitored neutrino beam" technique proposed in 2015. ENUBET has been designed to monitor lepton production in the decay tunnel of neutrino beams and to provide a 1% measurement of the neutrino flux at source. In particular, the three body semi-leptonic decay of kaons monitored by large angle positron production offers a fully controlled $ν_{e}$ source at the GeV scale for a new generation of short baseline experiments. In this contribution the performances of the positron tagger prototypes tested at CERN beamlines in 2016-2018 are presented. △ Less

Submitted 27 March, 2019; originally announced March 2019.

Comments: Poster presented at NuPhys2018 (London 19-21 December 2018). 5 pages, 3 figures

Report number: NuPhys2018-Tenti

The ENUBET Beamline

Authors: ENUBET Collaboration, G. Brunetti, F. Acerbi, G. Ballerini, M. Bonesini, A. Branca, C. Brizzolari, M. Calviani, S. Carturan, M. G. Catanesi, S. Cecchini, F. Cindolo, G. Collazuol, E. Conti, F. Dal Corso, G. De Rosa, C. Delogu, A. Falcone, B. Goddard, A. Gola, R. A. Intonti, C. Jollet, V. Kain, B. Klicek, Y. Kudenko , et al. (34 additional authors not shown)

Abstract: The ENUBET ERC project (2016-2021) is studying a narrow band neutrino beam where lepton production can be monitored at single particle level in an instrumented decay tunnel. This would allow to measure $ν_μ$ and $ν_{e}$ cross sections with a precision improved by about one order of magnitude compared to present results. In this proceeding we describe a first realistic design of the hadron beamline… ▽ More The ENUBET ERC project (2016-2021) is studying a narrow band neutrino beam where lepton production can be monitored at single particle level in an instrumented decay tunnel. This would allow to measure $ν_μ$ and $ν_{e}$ cross sections with a precision improved by about one order of magnitude compared to present results. In this proceeding we describe a first realistic design of the hadron beamline based on a dipole coupled to a pair of quadrupole triplets along with the optimisation guidelines and the results of a simulation based on G4beamline. A static focusing design, though less efficient than a horn-based solution, results several times more efficient than originally expected. It works with slow proton extractions reducing drastically pile-up effects in the decay tunnel and it paves the way towards a time-tagged neutrino beam. On the other hand a horn-based transferline would ensure higher yields at the tunnel entrance. The first studies conducted at CERN to implement the synchronization between a few ms proton extraction and a horn pulse of 2-10 ms are also described. △ Less

Submitted 26 November, 2020; v1 submitted 21 March, 2019; originally announced March 2019.

Comments: Poster presented at NuPhys2018 (London 19-21 December 2018). 4 pages, 3 figures

Report number: NuPhys2018-Brunetti

Gamma Factory at CERN -- novel research tools made of light

Authors: W. Placzek, A. Abramov, S. E. Alden, R. Alemany Fernandez, P. S. Antsiferov, A. Apyan, H. Bartosik, E. G. Bessonov, N. Biancacci, J. Bieron, A. Bogacz, A. Bosco, R. Bruce, D. Budker, K. Cassou, F. Castelli, I. Chaikovska, C. Curatolo, P. Czodrowski, A. Derevianko, K. Dupraz, Y. Dutheil, K. Dzierzega, V. Fedosseev, N. Fuster Martinez , et al. (37 additional authors not shown)

Abstract: We discuss the possibility of creating novel research tools by producing and storing highly relativistic beams of highly ionised atoms in the CERN accelerator complex, and by exciting their atomic degrees of freedom with lasers to produce high-energy photon beams. Intensity of such photon beams would be by several orders of magnitude higher than offered by the presently operating light sources, in… ▽ More We discuss the possibility of creating novel research tools by producing and storing highly relativistic beams of highly ionised atoms in the CERN accelerator complex, and by exciting their atomic degrees of freedom with lasers to produce high-energy photon beams. Intensity of such photon beams would be by several orders of magnitude higher than offered by the presently operating light sources, in the particularly interesting gamma-ray energy domain of 0.1-400 MeV. In this energy range, the high-intensity photon beams can be used to produce secondary beams of polarised electrons, polarised positrons, polarised muons, neutrinos, neutrons and radioactive ions. New research opportunities in a wide domain of fundamental and applied physics can be opened by the Gamma Factory scientific programme based on the above primary and secondary beams. △ Less

Submitted 14 June, 2019; v1 submitted 21 March, 2019; originally announced March 2019.

Comments: 12 pages; presented by W. Placzek at the XXV Cracow Epiphany Conference on Advances in Heavy Ion Physics, 8-11 January 2019, Cracow, Poland

Summary Report of Physics Beyond Colliders at CERN

Authors: R. Alemany, C. Burrage, H. Bartosik, J. Bernhard, J. Boyd, M. Brugger, M. Calviani, C. Carli, N. Charitonidis, D. Curtin, A. Dainese, A. de Roeck, M. Diehl, B. Döbrich, L. Evans, J. L. Feng, M. Ferro-Luzzi, L. Gatignon, S. Gilardoni, S. Gninenko, G. Graziani, E. Gschwendtner, B. Goddard, A. Hartin, I. Irastorza , et al. (39 additional authors not shown)

Abstract: Physics Beyond Colliders is an exploratory study aimed at exploiting the full scientific potential of CERN's accelerator complex and its scientific infrastructure in the next two decades through projects complementary to the LHC, HL-LHC and other possible future colliders. These projects should target fundamental physics questions that are similar in spirit to those addressed by high-energy collid… ▽ More Physics Beyond Colliders is an exploratory study aimed at exploiting the full scientific potential of CERN's accelerator complex and its scientific infrastructure in the next two decades through projects complementary to the LHC, HL-LHC and other possible future colliders. These projects should target fundamental physics questions that are similar in spirit to those addressed by high-energy colliders, but that require different types of beams and experiments. A kick-off workshop held in September 2016 identified a number of areas of interest and working groups have been set-up to study and develop these directions. All projects currently under consideration are presented including physics motivation, a brief outline of the experimental set-up and the status of the corresponding beam and detector technological studies. The proposals are also put in context of the worldwide landscape and their implementation issues are discussed. △ Less

Submitted 1 February, 2019; originally announced February 2019.

Comments: This document (66 pages, 19 figures) is the summary document of the Physics Beyond Colliders (PBC) study. It follows the PBC mandate and draws on a whole set of documents produced in the context of PBC (https://pbc.web.cern.ch), in particular the reports of the QCD and BSM working groups also available at arXiv:1901.04482 and arXiv:1901.09966

Report number: CERN-PBC-REPORT-2018-003

External Electron Injection for the AWAKE Experiment

Authors: Marlene Turner, Chiara Bracco, Spencer Gessner, Brennan Goddard, Edda Gschwendtner, Patric Muggli, Felipe Pena Asmus, Francesco Velotti, Livio Verra

Abstract: We summarize and explain the realization of witness particle injection into wakefields for the AWAKE experiment. In AWAKE, the plasma wakefields are driven by a self-modulating relativistic proton bunch. To demonstrate that these wakefields can accelerate charged particles, we inject a \unit[10-20]{MeV} electron bunch produced by a photo-injector. We summarize the experimental challenges of this i… ▽ More We summarize and explain the realization of witness particle injection into wakefields for the AWAKE experiment. In AWAKE, the plasma wakefields are driven by a self-modulating relativistic proton bunch. To demonstrate that these wakefields can accelerate charged particles, we inject a \unit[10-20]{MeV} electron bunch produced by a photo-injector. We summarize the experimental challenges of this injection process and present our plans for the near future. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: 4 pages, 3 figures

The Singular Hydrodynamic Interactions Between Two Spheres In Stokes Flow

Authors: B. D. Goddard, R. D. Mills-Williams, J. Sun

Abstract: We study exact solutions for the slow viscous flow of an infinite liquid caused by two rigid spheres approaching each either along or parallel to their line of centres, valid at all separations. This goes beyond the applicable range of existing solutions for singular hydrodynamic interactions (HIs) which, for practical applications, are limited to the near-contact or far field region of the flow.… ▽ More We study exact solutions for the slow viscous flow of an infinite liquid caused by two rigid spheres approaching each either along or parallel to their line of centres, valid at all separations. This goes beyond the applicable range of existing solutions for singular hydrodynamic interactions (HIs) which, for practical applications, are limited to the near-contact or far field region of the flow. For the normal component of the HI, by use of a bipolar coordinate system, we derive the stream function for the flow as $Re\to 0$ and a formula for the singular (squeeze) force between the spheres as an infinite series. We also obtain the asymptotic behaviour of the forces as the nondimensional separation between the spheres goes to zero and infinity, rigorously confirming and improving upon known results relevant to a widely accepted lubrication theory. Additionally, we recover the force on a sphere moving perpendicularly to a plane as a special case. For the tangential component, again by using a bipolar coordinate system, we obtain the corresponding infinite series expression of the (shear) singular force between the spheres. All results hold for retreating spheres, consistent with the reversibility of Stokes flow. We demonstrate substantial differences in numerical simulations of colloidal fluids when using the present theory compared with existing multipole methods. Furthermore, we show that the present theory preserves positive definiteness of the resistance matrix $\boldsymbol{R}$ in a number of situations in which positivity is destroyed for multipole/perturbative methods. △ Less

Submitted 27 March, 2020; v1 submitted 14 September, 2018; originally announced September 2018.

Comments: 28 pages, 12 Figures

The PROSPECT Reactor Antineutrino Experiment

Authors: PROSPECT Collaboration, J. Ashenfelter, A. B. Balantekin, C. Baldenegro, H. R. Band, C. D. Bass, D. E. Bergeron, D. Berish, L. J. Bignell, N. S. Bowden, J. Boyle, J. Bricco, J. P. Brodsky, C. D. Bryan, A. Bykadorova Telles, J. J. Cherwinka, T. Classen, K. Commeford, A. Conant, A. A. Cox, D. Davee, D. Dean, G. Deichert, M. V. Diwan, M. J. Dolinski , et al. (64 additional authors not shown)

Abstract: The Precision Reactor Oscillation and Spectrum Experiment, PROSPECT, is designed to make both a precise measurement of the antineutrino spectrum from a highly-enriched uranium reactor and to probe eV-scale sterile neutrinos by searching for neutrino oscillations over meter-long baselines. PROSPECT utilizes a segmented $^6$Li-doped liquid scintillator detector for both efficient detection of reacto… ▽ More The Precision Reactor Oscillation and Spectrum Experiment, PROSPECT, is designed to make both a precise measurement of the antineutrino spectrum from a highly-enriched uranium reactor and to probe eV-scale sterile neutrinos by searching for neutrino oscillations over meter-long baselines. PROSPECT utilizes a segmented $^6$Li-doped liquid scintillator detector for both efficient detection of reactor antineutrinos through the inverse beta decay reaction and excellent background discrimination. PROSPECT is a movable 4-ton antineutrino detector covering distances of 7m to 13m from the High Flux Isotope Reactor core. It will probe the best-fit point of the $\barν_e$ disappearance experiments at 4$σ$ in 1 year and the favored regions of the sterile neutrino parameter space at more than 3$σ$ in 3 years. PROSPECT will test the origin of spectral deviations observed in recent $θ_{13}$ experiments, search for sterile neutrinos, and address the hypothesis of sterile neutrinos as an explanation of the reactor anomaly. This paper describes the design, construction, and commissioning of PROSPECT and reports first data characterizing the performance of the PROSPECT antineutrino detector. △ Less

Submitted 21 August, 2019; v1 submitted 31 July, 2018; originally announced August 2018.

Comments: 30 pages, 33 figures; updated with journal revision and reference

Journal ref: Nucl.Instrum.Meth. A922 (2019) 287-309

First search for short-baseline neutrino oscillations at HFIR with PROSPECT

Authors: J. Ashenfelter, A. B. Balantekin, C. Baldenegro, H. R. Band, C. D. Bass, D. E. Bergeron, D. Berish, L. J. Bignell, N. S. Bowden, J. Bricco, J. P. Brodsky, C. D. Bryan, A. Bykadorova Telles, J. J. Cherwinka, T. Classen, K. Commeford, A. J. Conant, A. A. Cox, D. Davee, D. Dean, G. Deichert, M. V. Diwan, M. J. Dolinski, A. Erickson, M. Febbraro , et al. (63 additional authors not shown)

Abstract: This Letter reports the first scientific results from the observation of antineutrinos emitted by fission products of $^{235}$U at the High Flux Isotope Reactor. PROSPECT, the Precision Reactor Oscillation and Spectrum Experiment, consists of a segmented 4 ton $^6$Li-doped liquid scintillator detector covering a baseline range of 7-9 m from the reactor and operating under less than 1 meter water e… ▽ More This Letter reports the first scientific results from the observation of antineutrinos emitted by fission products of $^{235}$U at the High Flux Isotope Reactor. PROSPECT, the Precision Reactor Oscillation and Spectrum Experiment, consists of a segmented 4 ton $^6$Li-doped liquid scintillator detector covering a baseline range of 7-9 m from the reactor and operating under less than 1 meter water equivalent overburden. Data collected during 33 live-days of reactor operation at a nominal power of 85 MW yields a detection of 25461 $\pm$ 283 (stat.) inverse beta decays. Observation of reactor antineutrinos can be achieved in PROSPECT at 5$σ$ statistical significance within two hours of on-surface reactor-on data-taking. A reactor-model independent analysis of the inverse beta decay prompt energy spectrum as a function of baseline constrains significant portions of the previously allowed sterile neutrino oscillation parameter space at 95% confidence level and disfavors the best fit of the Reactor Antineutrino Anomaly at 2.2$σ$ confidence level. △ Less

Submitted 27 September, 2018; v1 submitted 7 June, 2018; originally announced June 2018.

Comments: 7 pages, 5 figures; v3: Added additional supplemental files

Journal ref: Phys. Rev. Lett. 121, 251802 (2018)

Non-adiabatic transitions in multiple dimensions

Authors: V. Betz, B. D. Goddard, T. Hurst

Abstract: We consider non-adiabatic transitions in multiple dimensions, which occur when the Born-Oppenheimer approximation breaks down. We present a general, multi-dimensional algorithm which can be used to accurately and efficiently compute the transmitted wavepacket at an avoided crossing. The algorithm requires only one-level Born-Oppenheimer dynamics and local knowledge of the potential surfaces. Cruci… ▽ More We consider non-adiabatic transitions in multiple dimensions, which occur when the Born-Oppenheimer approximation breaks down. We present a general, multi-dimensional algorithm which can be used to accurately and efficiently compute the transmitted wavepacket at an avoided crossing. The algorithm requires only one-level Born-Oppenheimer dynamics and local knowledge of the potential surfaces. Crucially, in contrast to standard methods in the literature, we compute the whole wavepacket, including its phase, rather than simply the transition probability. We demonstrate the excellent agreement with full quantum dynamics for a range of examples in two dimensions. We also demonstrate surprisingly good agreement for a system with a full conical intersection. △ Less

Submitted 12 April, 2018; originally announced April 2018.

Comments: 20 pages, 9 figures

Multiple Superadiabatic Transitions and Landau-Zener Formulas

Authors: Benjamin D. Goddard, Tim Hurst

Abstract: We consider nonadiabatic systems in which the classical Born-Oppenheimer approximation breaks down. We present a general theory that accurately captures the full transmitted wavepacket after multiple transitions through either a single or distinct avoided crossings, including phase information and associated interference effects. Under suitable approximations, we recover both the celebrated Landau… ▽ More We consider nonadiabatic systems in which the classical Born-Oppenheimer approximation breaks down. We present a general theory that accurately captures the full transmitted wavepacket after multiple transitions through either a single or distinct avoided crossings, including phase information and associated interference effects. Under suitable approximations, we recover both the celebrated Landau-Zener formula and standard surface-hopping algorithms. Our algorithm shows excellent agreement with the full quantum dynamics for a range of avoided crossing systems, and can also be applied to single full crossings with similar accuracy. △ Less

Submitted 12 April, 2018; originally announced April 2018.

Comments: 31 pages, 10 figures

Experimental results of crystal-assisted slow extraction at the SPS

Authors: M. A. Fraser, F. Addesa, G. Cavoto, F. Galluccio, S. Gilardoni, B. Goddard, F. Iacoangeli, V. Kain, D. Mirarchi, S. Montesano, F. Murtas, S. Petrucci, S. Redaelli, R. Rossi, W. Scandale, L. Stoel

Abstract: The possibility of extracting highly energetic particles from the Super Proton Synchrotron (SPS) by means of silicon bent crystals has been explored since the 1990's. The channelling effect of a bent crystal can be used to strongly deflect primary protons and eject them from the synchrotron. Many studies and experiments have been carried out to investigate crystal channelling effects. The extracti… ▽ More The possibility of extracting highly energetic particles from the Super Proton Synchrotron (SPS) by means of silicon bent crystals has been explored since the 1990's. The channelling effect of a bent crystal can be used to strongly deflect primary protons and eject them from the synchrotron. Many studies and experiments have been carried out to investigate crystal channelling effects. The extraction of 120 and 270 GeV proton beams has already been demonstrated in the SPS with dedicated experiments located in the ring. Presently in the SPS, the UA9 experiment is performing studies to evaluate the possibility to use bent silicon crystals to steer particle beams in high energy accelerators. Recent studies on the feasibility of extraction from the SPS have been made using the UA9 infrastructure with a longer-term view of using crystals to help mitigate slow extraction induced activation of the SPS. In this paper, the possibility to eject particles into the extraction channel in LSS2 using the bent crystals already installed in the SPS is presented. Details of the concept, simulations and measurements carried out with beam are presented, before the outlook for the future is discussed. △ Less

Submitted 17 July, 2017; originally announced July 2017.

Comments: 4 pages, 7 figures, submitted to to International Particle Accelerator Conference (IPAC) 2017 in Copenhagen, Denmark

Status of the HIE-ISOLDE project at CERN

Authors: M. A. Fraser, Y. Kadi, A. P. Bernardes, Y. Blumenfeld, E. Bravin, S. Calatroni, R. Catherall, B. Goddard, D. Parchet, E. Siesling, W. Venturini Delsolaro, G. Vandoni, D. Voulot, L. R. Williams

Abstract: The HIE-ISOLDE project represents a major upgrade of the ISOLDE nuclear facility with a mandate to significantly improve the quality and increase the intensity and energy of radioactive nuclear beams produced at CERN. The project will expand the experimental nuclear physics programme at ISOLDE by focusing on an upgrade of the existing Radioactive ion beam EXperiment (REX) linac with a 40 MV superc… ▽ More The HIE-ISOLDE project represents a major upgrade of the ISOLDE nuclear facility with a mandate to significantly improve the quality and increase the intensity and energy of radioactive nuclear beams produced at CERN. The project will expand the experimental nuclear physics programme at ISOLDE by focusing on an upgrade of the existing Radioactive ion beam EXperiment (REX) linac with a 40 MV superconducting linac comprising thirty-two niobium-on-copper sputter-coated quarter-wave resonators housed in six cryomodules. The new linac will raise the energy of post-accelerated beams from 3 MeV/u to over 10 MeV/u. The upgrade will be staged to first deliver beam energies of 5.5 MeV/u using two high-$β$ cryomodules placed downstream of REX, before the energy variable section of the existing linac is replaced with two low-$β$ cryomodules and two additional high-$β$ cryomodules are installed to attain over 10 MeV/u with full energy variability above 0.45 MeV/u. An overview of the project including a status summary of the different R&D activities and the schedule will outlined. △ Less

Submitted 17 July, 2017; originally announced July 2017.

Comments: 7 pages, 12 figures, submitted to the Heavy Ion Accelerator Technology conference (HIAT) 2012, in Chicago

Physics Opportunities with the FCC-hh Injectors

Authors: B. Goddard, G. Isidori, F. Teubert, M. Bai, A. Ball, B. Batell, T. Bowcock, G. Cavoto, A. Ceccucci, M. Chrzaszcz, A. Golutvin, W. Herr, J. Jowett, M. Moulson, T. Nakada, J. Rojo, Y. Semertzidis

Abstract: In this chapter we explore a few examples of physics opportunities using the existing chain of accelerators at CERN, including potential upgrades. In this context the LHC ring is also considered as a part of the injector system. The objective is to find examples that constitute sensitive probes of New Physics that ideally cannot be done elsewhere or can be done significantly better at theCERN acce… ▽ More In this chapter we explore a few examples of physics opportunities using the existing chain of accelerators at CERN, including potential upgrades. In this context the LHC ring is also considered as a part of the injector system. The objective is to find examples that constitute sensitive probes of New Physics that ideally cannot be done elsewhere or can be done significantly better at theCERN accelerator complex. Some of these physics opportunities may require a more flexible injector complex with additional functionality than that just needed to inject protons into the FCC-hh at the right energy, intensity and bunch structure. Therefore it is timely to discuss these options concurrently with the conceptual design of the FCC-hh injector system. △ Less

Submitted 22 June, 2017; originally announced June 2017.

Comments: 13 pages, chapter 5 in Physics at the FCC-hh, a 100 TeV pp collider

Journal ref: CERN Yellow Report CERN 2017-003-M, pp. 693-705

Injection and Dumping Systems

Authors: B. Goddard, A. Lechner, J. Uythoven

Abstract: Chapter 14 in High-Luminosity Large Hadron Collider (HL-LHC) : Preliminary Design Report. The Large Hadron Collider (LHC) is one of the largest scientific instruments ever built. Since opening up a new energy frontier for exploration in 2010, it has gathered a global user community of about 7,000 scientists working in fundamental particle physics and the physics of hadronic matter at extreme tempe… ▽ More Chapter 14 in High-Luminosity Large Hadron Collider (HL-LHC) : Preliminary Design Report. The Large Hadron Collider (LHC) is one of the largest scientific instruments ever built. Since opening up a new energy frontier for exploration in 2010, it has gathered a global user community of about 7,000 scientists working in fundamental particle physics and the physics of hadronic matter at extreme temperature and density. To sustain and extend its discovery potential, the LHC will need a major upgrade in the 2020s. This will increase its luminosity (rate of collisions) by a factor of five beyond the original design value and the integrated luminosity (total collisions created) by a factor ten. The LHC is already a highly complex and exquisitely optimised machine so this upgrade must be carefully conceived and will require about ten years to implement. The new configuration, known as High Luminosity LHC (HL-LHC), will rely on a number of key innovations that push accelerator technology beyond its present limits. Among these are cutting-edge 11-12 tesla superconducting magnets, compact superconducting cavities for beam rotation with ultra-precise phase control, new technology and physical processes for beam collimation and 300 metre-long high-power superconducting links with negligible energy dissipation. The present document describes the technologies and components that will be used to realise the project and is intended to serve as the basis for the detailed engineering design of HL-LHC. △ Less

Submitted 26 May, 2017; originally announced May 2017.

Comments: 8 pages, chapter 14 in High-Luminosity Large Hadron Collider (HL-LHC) : Preliminary Design Report

Journal ref: CERN Yellow Report CERN 2015-005, pp.221-228

PERLE: Powerful Energy Recovery Linac for Experiments - Conceptual Design Report

Authors: D. Angal-Kalinin, G. Arduini, B. Auchmann, J. Bernauer, A. Bogacz, F. Bordry, S. Bousson, C. Bracco, O. Brüning, R. Calaga, K. Cassou, V. Chetvertkova, E. Cormier, E. Daly, D. Douglas, K. Dupraz, B. Goddard, J. Henry, A. Hutton, E. Jensen, W. Kaabi, M. Klein, P. Kostka, F. Marhauser, A. Martens , et al. (17 additional authors not shown)

Abstract: A conceptual design is presented of a novel ERL facility for the development and application of the energy recovery technique to linear electron accelerators in the multi-turn, large current and large energy regime. The main characteristics of the powerful energy recovery linac experiment facility (PERLE) are derived from the design of the Large Hadron electron Collider, an electron beam upgrade u… ▽ More A conceptual design is presented of a novel ERL facility for the development and application of the energy recovery technique to linear electron accelerators in the multi-turn, large current and large energy regime. The main characteristics of the powerful energy recovery linac experiment facility (PERLE) are derived from the design of the Large Hadron electron Collider, an electron beam upgrade under study for the LHC, for which it would be the key demonstrator. PERLE is thus projected as a facility to investigate efficient, high current (> 10 mA) ERL operation with three re-circulation passages through newly designed SCRF cavities, at 801.58 MHz frequency, and following deceleration over another three re-circulations. In its fully equipped configuration, PERLE provides an electron beam of approximately 1 GeV energy. A physics programme possibly associated with PERLE is sketched, consisting of high precision elastic electron-proton scattering experiments, as well as photo-nuclear reactions of unprecedented intensities with up to 30 MeV photon beam energy as may be obtained using Fabry-Perot cavities. The facility has further applications as a general technology test bed that can investigate and validate novel superconducting magnets (beam induced quench tests) and superconducting RF structures (structure tests with high current beams, beam loading and transients). Besides a chapter on operation aspects, the report contains detailed considerations on the choices for the SCRF structure, optics and lattice design, solutions for arc magnets, source and injector and on further essential components. A suitable configuration derived from the here presented design concept may next be moved forward to a technical design and possibly be built by an international collaboration which is being established. △ Less

Submitted 24 May, 2017; originally announced May 2017.

Comments: 112 pages, 63 figures

Pseudospectral methods for density functional theory in bounded and unbounded domains

Authors: Andreas Nold, Benjamin D. Goddard, Peter Yatsyshin, Nikos Savva, Serafim Kalliadasis

Abstract: Classical Density Functional Theory (DFT) is a statistical-mechanical framework to analyze fluids, which accounts for nanoscale fluid inhomogeneities and non-local intermolecular interactions. DFT can be applied to a wide range of interfacial phenomena, as well as problems in adsorption, colloidal science and phase transitions in fluids. Typical DFT equations are highly non-linear, stiff and conta… ▽ More Classical Density Functional Theory (DFT) is a statistical-mechanical framework to analyze fluids, which accounts for nanoscale fluid inhomogeneities and non-local intermolecular interactions. DFT can be applied to a wide range of interfacial phenomena, as well as problems in adsorption, colloidal science and phase transitions in fluids. Typical DFT equations are highly non-linear, stiff and contain several convolution terms. We propose a novel, efficient pseudo-spectral collocation scheme for computing the non-local terms in real space with the help of a specialized Gauss quadrature. Due to the exponential accuracy of the quadrature and a convenient choice of collocation points near interfaces, we can use grids with a significantly lower number of nodes than most other reported methods. We demonstrate the capabilities of our numerical methodology by studying equilibrium and dynamic two-dimensional test cases with single- and multispecies hard-sphere and hard-disc particles modelled with fundamental measure theory, with and without van der Waals attractive forces, in bounded and unbounded physical domains. We show that our results satisfy statistical mechanical sum rules. △ Less

Submitted 22 January, 2017; originally announced January 2017.

Sum of two repdigits a square

Authors: Bart Goddard, Jeremy Rouse

Abstract: A \emph{repdigit} is a natural number greater than 10 which has all of its base-10 digits the same. In this paper we find all examples of two repdigits adding to a square. The proofs lead to interesting questions about consecutive quadratic residues and non-residues, and provide an elementary application of elliptic curves. A \emph{repdigit} is a natural number greater than 10 which has all of its base-10 digits the same. In this paper we find all examples of two repdigits adding to a square. The proofs lead to interesting questions about consecutive quadratic residues and non-residues, and provide an elementary application of elliptic curves. △ Less

Submitted 22 July, 2016; originally announced July 2016.

Comments: 9 pages

MSC Class: 11D61 (Primary); 11G05; 11A63 (Secondary)

Journal ref: Integers 17 (2017), Paper No. A24, 8 pp

Dynamical density functional theory for orientable colloids including inertia and hydrodynamic interactions

Authors: Miguel A. Durán-Olivencia, Benjamin D. Goddard, Serafim Kalliadasis

Abstract: Over the last few decades, classical density-functional theory (DFT) and its dynamic extensions (DDFTs) have become powerful tools in the study of colloidal fluids. Recently, previous DDFTs for spherically-symmetric particles have been generalised to take into account both inertia and hydrodynamic interactions, two effects which strongly influence non-equilibrium properties. The present work furth… ▽ More Over the last few decades, classical density-functional theory (DFT) and its dynamic extensions (DDFTs) have become powerful tools in the study of colloidal fluids. Recently, previous DDFTs for spherically-symmetric particles have been generalised to take into account both inertia and hydrodynamic interactions, two effects which strongly influence non-equilibrium properties. The present work further generalises this framework to systems of anisotropic particles. Starting from the Liouville equation and utilising Zwanzig's projection-operator techniques, we derive the kinetic equation for the Brownian particle distribution function, and by averaging over all but one particle, a DDFT equation is obtained. Whilst this equation has some similarities with DDFTs for spherically-symmetric colloids, it involves a translational-rotational coupling which affects the diffusivity of the (asymmetric) particles. We further show that, in the overdamped (high friction) limit, the DDFT is considerably simplified and is in agreement with a previous DDFT for colloids with arbitrary shape particles. △ Less

Submitted 8 July, 2016; originally announced July 2016.

Comments: dynamical density functional theory ; colloidal fluids ; arbitrary-shape particles ; orientable colloids

Journal ref: J Stat Phys (2016) 164: 785