High Accelerator Research Articles

Design and development of the machine protection system for HIAF

ABSTRACT The High Intensity Heavy-ion Accelerator Facility (HIAF), currently under construction, is a complex machine that couples a Continuous Wave (CW) superconducting ion Linear accelerator (iLinac) with a high-energy synchrotron to produce various stable and radioactive intense beams with energies from MeV/u to GeV/u. The machine has a versatile operation mode which requires a high flexibility and reliability to the Machine Protection System (MPS). A customized and robust MPS is designed and developed to give the readiness of the machine for operation, to mitigate and analyze faults related to the relative damage potential. To get a high speed and have a high level of reliability, all interlock signal processing is processed on radiation-tolerant Field-Programmable Gate Arrays (FPGA) with triple or dual redundancy, as well as with a fail-safe design. By implementing a multiprocessing platform system-on-chip FPGA, the HIAF MPS can be tightly integrated with other systems to maximize availability pinpoint failures for operations, and give the postmortem analysis. This paper will describe the architecture of the interlocks linking the protection systems, the strategies to manage the complexity, the detailed components, and the interlock logic of the customized HIAF MPS, as well as the test and verification of the prototype.

GPU accelerated Staggered Update Procedure (SUP)

The advancement in programmable capability of graphics hardware has paved new opportunities in the domain of high performance computing (HPC). The computational fluid dynamics (CFD) community, being a significant user of HPC, has started exploiting the inherent data parallelism in the numerical solvers to be able to make efficient use of these many-core, high throughput accelerator based processors. In the present work, we examine the process of accelerating our CPU based Staggered Update Procedure (SUP) solver, i.e., a higher order accurate cell-centred finite volume solver by off-loading the computationally most expensive region of the code pertaining to the explicit residual computation. We have adopted OpenACC, a directive based programming model to expose parallelism in the code. The framework evolved for GPU porting in the context of SUP is also of value to those intending to port their CFD solvers based on classical finite volume methodology. The performance analysis is conducted using scalar convection–diffusion equations in both two- and three-dimensions. The findings demonstrate a speedup factor of 9 (in case of 2D) and 28 (in case of 3D) when considering the explicit residual alone, achieved with a single NVIDIA Tesla V100 GPU card. In addition, we could establish superior algorithmic scalability by the way of recovering near perfect serial performance, on the heterogeneous CPU+GPU architecture. Further, overall code acceleration can be achieved by porting other parts of the solver on GPU.

Study of the effect of detector to front-end electronics distance on the spectrometric performances of solid-state detectors

Solid-state detectors are used for charge particles detection and spectrometry with front-end electronics located, usually, next to the detectors. Therefore, in some applications, due to high-level radiation, high temperature and/or electro-magnetic noises (e.g. in nuclear fusion energy experiments, high energy accelerators, etc.) there could be the need to locate the front-end electronics far away from the detector. To solve the problem a matching sensitive charge amplifier (MSCA) with AC coupled input was designed by Roberto Cardarelli (former INFN scientist) for measurements at the fusion JET tokamak (U.K.). The prototype, composed of two amplification stages, was “ad hoc” designed for use with Diamond detectors and resulted very reliable. A commercial version of the MSCA was then developed and is now available on the market. Its novelty and main advantage, with respect to the former version, is to be used also with other fast detectors e.g. Resistive Plate Chamber (RPC) with high counting speed and Silicon detectors. The commercial MSCA features input impedance very close to 50 Ohm, which can be matched to a 50-Ohm transmission line. This allows locating the amplifier at variable distances from the detector and up to 100 m. These claimed performances result very interesting for solid state detectors, therefore, up to now a few applications of MSCA are reported in the literature and a systematic study of its performances is missing. In this work, an analysis of the spectrometric performances of Diamond and Silicon detectors coupled to MSCA located at various distances (up to 48 m) from the detectors is reported. Pulse height spectra of a multi-peaks alpha source were measured using MSCA amplifier and compared with those obtained using a standard charge sensitive preamplifier (CA). Then, the effect of the distance between detectors and amplifiers was investigated using standard RG58 and double shielded RG214 coaxial cables. The results are presented and discussed.

Phase-contrast visualization of human tissues using superimposed wavefront imaging of diffraction-enhanced x-rays.

Phase-contrast computed tomography (CT) using high-brilliance, synchrotron-generated x-rays enable three-dimensional (3D) visualization of microanatomical structures within biological specimens, offering exceptionally high-contrast images of soft tissues. Traditional methods for phase-contrast CT; however, necessitate a gap between the subject and the x-ray camera, compromising spatial resolution due to penumbral blurring. Our newly developed technique, Superimposed Wavefront Imaging of Diffraction-enhanced x-rays (SWIDeX), leverages a Laue-case Si angle analyzer affixed to a scintillator to convert x-rays to visible light, capturing second-order differential phase contrast images and effectively eliminating the distance to the x-ray camera. This innovation achieves superior spatial resolution over conventional methods. In this paper, the imaging principle and CT reconstruction algorithm based on SWIDeX are presented in detail and compared with conventional analyzer-based imaging (ABI). It also shows the physical setup of SWIDeX that provides the resolution preserving second-order differential images for reconstruction. We compare the spatial resolution and the sensitivity of SWIDeX to conventional ABI. To demonstrate high-spatial resolution achievable by SWIDeX, the internal structures of four human tissues-ductal carcinoma in situ, normal stomach, normal pancreas, and intraductal papillary mucinous neoplasm of the pancreas-were visualized using an imaging system configured at the Photon Factory's BL14B beamline under the High Energy Accelerator Research Organization (KEK). Each tissue was thinly sliced after imaging, stained with hematoxylin and eosin (H&E) for conventional microscope-based pathology. A comparison of SWIDeX-CT and pathological images visually demonstrates the effectiveness of SWIDeX-CT for biological tissue imaging. SWIDeX could generate clearer 3D images than existing analyzer-based phase-contrast methods and accurately delineate tissue structures, as validated against histopathological images. SWIDeX can visualize important 3D structures in biological soft tissue with high spatial resolution and can be an important tool for providing information between the disparate scales of clinical and pathological imaging.

Phase space reconstruction technique for beam optimization in the Low Energy High Intensity Proton Accelerator (LEHIPA)

The Low Energy High Intensity Proton Accelerator (LEHIPA) has been commissioned to the design energy of 20 MeV at BARC, India. The low energy beam transport (LEBT) channel of LEHIPA consists of two solenoids and drift spaces for matching the 50 keV proton beam from the ECR ion source (ECR-IS) to the RFQ. The ion beam extracted from the ECR-IS also contains molecular species like H2 + and H3 +. Proton fraction in the beam is found to degrade slowly with time due to surface contamination of the plasma chamber and this reduction in proton beam current has implications for longitudinal and transverse beam dynamics. Hence, it becomes important to characterise the beam at different proton fraction levels to understand the end-to-end beam dynamics of LEHIPA. Computed Tomography (CT) technique has been used for the beam phase-space reconstruction in LEHIPA, using a Python program, incorporating the feature of filtering secondary species from the beam profiles measured using slit scanners in the LEBT. Simultaneous Algebraic Reconstruction Technique (SART) is used in the reconstruction since it is identified to be a better technique for a limited number of projections. The reconstruction program is benchmarked with the TraceWin beam dynamics code and implemented on the measured beam profiles to recreate the phase space distribution at the beginning of LEBT. Further, the tomographic reconstruction method is compared with the solenoid scan method and the rms emittance values are found to be in good agreement. The measured tomographic phase space distribution has then been used as TraceWin input for LEHIPA end-end beam dynamics simulations and the LEHIPA beam line parameters are re-optimized for minimum beam emittance growth. This paper presents the simulations of CT technique, benchmarking simulations with TraceWin, phase space reconstruction with measured beam profiles and beam dynamics studies of LEHIPA using the reconstructed beam distributions.

Beam dynamics simulations at the high energy accelerator section in the RAON accelerator

Beam dynamics simulations at the high energy accelerator section in the RAON accelerator

An inverse free electron laser-bunching driven shallow-angle superradiant Compton source

A coherent soft x-ray Compton source driven by a density modulated electron beam is investigated. The density modulation is attained by an inverse free electron laser prebuncher, which offers strong energy modulation at the relevant beam energies for Compton scattering in an shallow-angle oblique scattering geometry using relatively low laser beam power. In particular, we study the case of a 101.7 MeV electron beam produced by a state-of-the-art high brightness accelerator beamline. Calculations show that 1011 coherent soft x-ray photons per second can be generated in an ultra-narrow relative bandwidth of 0.002% from a 250 pC microbunched electron bunch colliding with a 23 mJ laser pulse at 10 kHz.

Beam physics studies for a high charge and high beam quality laser-plasma accelerator

Electron acceleration by laser-plasma techniques is approaching maturity and is getting ready for the construction of particle accelerators with dedicated applications. We present a general methodology showing how beam physics studies can be used to achieve a specific parameter set in a laser-plasma accelerator. Laser systems, plasma targets, and magnetic component properties are designed to optimize the electron beam so as to achieve the required performances. Beam physics in its full 6D phase space is studied from electron injection to beam delivery to the end user, through the plasma acceleration stage and transport line. As each beam parameter can only be modified by specific electric/magnetic field configurations, it is crucial to assign from the beginning specific roles to given accelerator sections in obtaining given beam parameters. These beam physics considerations were successfully applied to the design of a plasma-based electron injector for the AWAKE Run2 experiment. Electron beam parameters were calculated using a global simulation, achieving simultaneously unprecedented high charge (100 pC) and high quality (micrometric beam emittance and size). Published by the American Physical Society 2024

New developments in the MCUNED-Plus code for radiation transport and coupled transport-activation computational simulations in accelerator-based facilities

The radiation hazard is one of the major aspects of concern in the design of a facility where it exists a risk of exposure to ionizing radiation. The nuclear analysis for radioprotection purposes requires an accurate description of the radiation fields. This is especially true in accelerator-based facilities where several types of particles are producing the radiation fields. Radiation fields exist during the accelerator operation, but delayed radiation may also be significant for a long time after operation ends. The radiation fields of concern present during the operation are due to the secondary neutrons and photons emitted by the interaction of accelerated particles with the accelerator components (beam leakage) or with the target. The radiation present after the operation is the residual photon field produced by activated materials, these materials being activated by both accelerated particles and secondary neutrons. This residual radiation field is relevant in high intensity accelerators like IFMIF-DONES.In order to address these complex radiation transport simulations, the D1SUNED code has been updated and the new code release renamed as MCUNED-Plus. The new developments include improvements in the light-ion transport like the implementation of a variance reduction for the production of secondary particles, and a new kinematics to reproduce the angular distribution of secondary particles emitted after deuteron breakup reaction. The calculation of the residual photon field in accelerator facility has also been improved by allowing to evaluate both light-ions and secondary neutrons induced shutdown dose rate in a single coupled simulation.

High Performance CNN Accelerators based on Hardware and Algorithm Co-Optimization

The efficacy of convolutional neural networks (CNNs) has led to their extensive application in picture recognition and categorization. Nevertheless, embedded systems face challenges when it comes to storing the massive volumes of weight data required by CNNs in their on-chip memory. Pruning a convolutional neural network (CNN) model can reduce its size with no impact on accuracy; however, when used with a parallel architecture, the resulting model will be slower to execute. A CNN compression method that is hardware-centric is introduced in this paper. A model of a deep neural network (DNN) will include "pruning layers (P-layers)" and "no-pruning layers (NP-layers)" attached to it. For effective and parallel processing, an NP-layer employs a regular distribution for its weights. Pruning causes a P-layer to have an irregular form, yet the high compression ratio it yields makes the shape undesirable. Using uniform and gradual quantization methods, we can achieve a processing efficiency/compression ratio trade-off with a small loss of precision. Integrating multiple parallel finite impulse response (FIR) filters into an NP-layer distributed convolutional architecture, we further improve the regular model. The P-layer irregular sparse model is suggested to use an ADF-based processing element. The suggested compression strategy and hardware architecture can be utilised by a hardware/algorithm co-optimization (HACO) method to run an NP→P hybrid compressed CNN model on FPGAs. The VGG-16 network achieves a compression ratio of 27.5× with a top-5 accuracy loss of 0.44% by using a hardware accelerator on a single FPGA chip without off-chip memory. A result of 83.0 FPS, or 1.8 times quicker, was achieved when the compressed VGG-16 model was implemented on a Xilinx VCU118 evaluation board for image applications. Index Terms - CNN, Accelerators, Co-Optimization.

Activity estimates of localised gamma emitting radionuclides using a GR-135 survey meter

During the operation of high energy accelerators activated materials are commonly created. The activity and isotopes present in these materials must be characterised for their clearance and release from the facility, or to ascertain their duration of stay in a radiological storage area. An activity estimate method using a gamma detecting GR-135 survey meter, which has the ability to collect an energy spectrum, is presented. Using several reference radioactive sources the detection efficiency and dead time of the survey meter were characterised. This information combined with the physical properties of the survey meter, the counting time and the properties of the measured photon energy emissions can be used to calculate an accurate activity estimate for localised activation on accelerator components, or loose contamination on isolated waste materials.

Training-free demonstration of a 5.4 T Nb3Sn Canted–Cosine–Theta accelerator dipole impregnated with paraffin wax

Lawrence Berkeley National Laboratory is pursuing stress-managed Nb3Sn Canted–Cosine–Theta (CCT) magnet technology for high field accelerator magnets. Although promising results have been reported, improvements in the training performance are desired. This work describes the fabrication and testing campaigns of two subscale Nb3Sn CCT magnets; a baseline impregnated with National High Magnetic Field Laboratory Mix-61 and a magnet impregnated with paraffin wax. The paraffin magnet reached the short sample limit of the conductor, to within the measurement uncertainty, with no training quenches inside the magnet. In contrast, the baseline magnet reached 80% of the short sample limit after approximately 20 quenches and exhibits some loss of memory after thermal cycles. Inter-layer flexible quench antennas combined with voltage tap data show that all quenches appear identical and originate from a region corresponding to the location of peak field in the cable. Although this success should be replicated at higher fields, these first test results demonstrate the potential for training free Nb3Sn accelerator magnets operating near the short-sample limit.

Ultrafast laser triggered electron emission from ultrananocrystalline diamond pyramid tip cathode

Nitrogen-incorporated ultrananocrystalline diamond [(N)UNCD] pyramid tip cathode has been considered as a next-generation high peak current electron source for dielectric laser accelerators as well as other high peak current particle accelerator applications. In this work, we study non-linear photoemission from an (N)UNCD pyramid tip cathode using an ultrafast laser with the pulse length of 150 fs with the central wavelength of 800 nm in the peak intensity range of 109–1010W/cm2. We demonstrated that as the incident laser intensity increases, the current emitted from the nano-tip first increases as a power function with an exponent of about 5 and then starts to roll over to an exponent of 3. This roll over is attributed to the Coulomb interaction between electrons emitted from the tip also known as the space charge. We also measured the photoemission electron energy spectra that show electrons with energies as high as ∼10 eV. Based on the shape of the electron energy spectra, we conclude that the high-energy electrons are thermally emitted electrons due to ultrafast laser heating at the tip of the (N)UNCD pyramid tip cathode.

Analysis of single particle and collective beam effects in high intensity beams in a periodic quadrupole channel

In high intensity proton accelerators, there are two main mechanisms that can cause beam degradation: incoherent and coherent effects due to nonlinear space charge forces. The incoherent effects represent the single particle dynamics while the coherent effects represent the collective response of the beam. Of particular interest is the region above a zero current phase advance (σ0) of 90° where the (coherent) second-order envelope instability and the (incoherent) fourth-order particle resonance are seen to lead to emittance growth. Large emittance growth is also seen below the envelope instability region as the full current beam phase advance (σ) decreases. In the present study, we have studied the nonlinear effects in a high intensity beam propagating through a focusing-defocusing (FD) quadrupole channel for σ0 greater than 90°, both analytically, by studying the solutions of Kapchinsky-Vladimirsky (KV) equation and the particle core model, and through detailed particle-in-cell (PIC) simulations using the code. The KV envelope equation gives the collective response of the beam while the particle core model gives the contribution of the single particle effects. With pic simulations, which resemble the behavior of the real beams more closely as compared to envelope calculations or the particle core model, it is possible to study the evolution of the beam in a self-consistent manner. Our studies show, that it is possible to identify the specific process responsible for beam degradation in high intensity beams. It is also possible to identify which process dominates under different conditions. We further show that the width of the emittance increase stop band calculated from PIC simulations is wider than that calculated by the envelope equations and that the width depends on the length of the channel being studied. Published by the American Physical Society 2024

Multi-bit MRAM based high performance neuromorphic accelerator for image classification

Binary neural networks (BNNs) are the most efficient solution to bridge the design gap of the hardware implementation of neural networks in a resource-constrained environment. Spintronics is a prominent technology among emerging fields for next-generation on-chip non-volatile memory. Spin transfer torque (STT) and spin-orbit torque (SOT) based magnetic random-access memory (MRAM) offer non-volatility and negligible static power. Over the last few years, STT and SOT-based multilevel spintronic memories have emerged as a promising solution to attain high storage density. This paper presents the operation principle and performance evaluation of spintronics-based single-bit STT and SOT MRAM, dual-level cells, three-level cells (TLCs), and four-level cells. Further, multi-layer perceptron architectures have been utilized to perform MNIST image classification with these multilevel devices. The performance of the complete system level consisting of crossbar arrays with various MRAM bit cells in terms of area, energy, and latency is evaluated. The throughput efficiency of the BNN accelerator using TLCs is 26.6X, and 3.61X higher than conventional single-bit STT-MRAM, and SOT-MRAM respectively.

Low-energy ion beam diagnostics: an integrated solution

High gradient accelerator injectors have been widely studied throughout the world-leading accelerator facilities. The demand for high frequency cavities has led the Detector, Accelerator and instrumentation laboratory (KAHVELab) in Istanbul to deploy a four-vane Radio Frequency Quadrupole (RFQ) operating at 800 MHz to accelerate 20 keV protons to 2 MeV. The protons from the microwave discharge ion source are transversely matched to the RFQ via an optimized Low Energy Beam Transport (LEBT) line which also contains an integrated measurement station, called measurement station (MBOX). The MBOX is used to measure the proton beam's current along with pulse length duration, profile as well as the beam emittance upstream of the RFQ. It contains a number of home-built diagnostic tools: a Faraday cup, a scintillator screen and a pepper pot plate (PP). The analysis software is also locally developed and tested for the PP photo analysis. In this note, the design, construction and tests of the integrated measurement station are discussed. The results from various measurements, especially on beam profile and charge, are compared to the simulation predictions.

Effect of the second curing cycle on performance of superconducting magnet insulating system

In large-scale, high-field superconducting magnets used for magnetic confinement fusion, high energy accelerators, and magnetic resonance imaging, the insulating system made from glass fiber reinforced resin-based composites is the key component, which mainly plays the role of mechanical support, fixing and protecting superconducting conductors, as well as electrical insulation. Vacuum Pressure Impregnation (VPI) approach is widely used in the manufacturing of the insulation system. The second curing cycle is generally required after the first VPI and curing process. For example, after the superconducting coil is cured in the mold, the de-molding process requires the superconducting coil to be reheated according the curing temperature. Moreover, for large-scale superconducting magnets, the superconducting coil needs to undergo a second VPI process after the first VPI process to fix the coil in the coil case. In this work, the tensile and shear properties of pure epoxy resin and the glass fiber reinforced resin-based composite, were investigated at both room and cryogenic temperatures and the effect of the second curing cycle on the mechanical properties was analyzed. Additionally, the strain evolution of the Nb-Ti superconducting coil during the second curing cycle was measured using the Fiber Bragg Grating (FBG) sensors embedded in the composite. The results indicate that the second curing cycle will not introduce additional strain to the previously cured resin matrix, but the defective or weak parts of the resin matrix may be affected by the new added epoxy resin and a little extra strain has been observed.

High power electrostatic beam splitter for a proton beamline

The High Intensity Proton Accelerator facility (HIPA) delivers a 590 MeV cw (50.6 MHz) proton beam with up to 1.4 MW beam power (2.4 mA) to spallation and meson production targets serving particle physics experiments and material research. The main accelerator is the ring cyclotron, an isochronous proton machine accelerating an injected 72 MeV beam to a final 590 MeV. A few meters downstream of the ring cyclotron, an electrostatic beam splitter was installed in the 1980s and originally designed to peel off from a 200 μA beam up to 20 μA (12 kW beam power). Future initiatives will also make use of the splitter. Specifically, as part of the Isotope and Muon Production using Advanced Cyclotron and Target technologies (IMPACT) upgrade project, Targeted Alpha Tumour Therapy and Other Oncological Solutions (TATTOOS), an online isotope separation facility will allow to produce promising radionuclides for diagnosis and therapy of cancer in quantities sufficient for clinical studies. The TATTOOS facility includes a dedicated beamline intended to operate at a beam intensity of 100 μA (60 kW beam power), requiring continuous splitting of the high-power main beam via the splitter. As a step forward toward reaching the desired beam intensity, a beam study was carried out to test the viability of the existing splitter for TATTOOS. The results of this study show that a record of 90 μA (53 kW beam power) was peeled off a horizontally and vertically enlarged beam by the splitter. The successful beam strategy employed during the study as well as the results of several key measurements are presented in this paper, with particular emphasis on diagnostic measurements. Additionally, to support the measurements, a computational model of the splitter has been implemented using Monte Carlo simulation tools, including realistic geometry, electrostatic fields, beam optics, and power deposition calculations. Overall, the results of this paper show that through the combination of beam measurements and simulations, the existing splitter can be used to reach the 100−μA beam intensity requirement for TATTOOS. Published by the American Physical Society 2024

Fabrication and rejuvenation of high quantum efficiency caesium telluride photocathodes for high brightness and high average current photoinjectors

With their high quantum efficiency (QE) and excellent photoemissive properties, cesium telluride (Cs-Te) photocathodes are the current workhorse in the high average current electron accelerators around the globe. Their ability to generate high brightness and high charge electron beams has opened the doorway to numerous applications, including fundamental particle physics research, radiation therapy, high-energy physics experiments, and high repetition rate free-electron lasers (FELs). Their long-term performance is critical. In this work, we analyze the systematic production and rejuvenation of photocathodes via thermal evaporation and deposition techniques (both sequential and co-deposition), showing excellent values of QE exceeding 20% using deep ultraviolet illumination (λ=262–266 nm). To extend the photocathode lifetime, the rejuvenation process is performed via multilayer thin-film co-deposition, demonstrating the feasibility of reliably recovering the initial QE even after air exposure. We evaluated their final performance in both high gradient dc and rf gun setups and obtained consistent results of QE and lifetimes. The technique can be used at the gun level, which significantly extends the sustained accelerator operation without major interventions. These approaches significantly enhance the robustness and performance of Cs-Te photoinjectors and represent a significant advancement for reliable high average current electron accelerators and FELs. Published by the American Physical Society 2024

Beam optics distortion caused by orbit deviation at strong sextupole magnets in SuperKEKB

The SuperKEKB is an electron-positron double-ring asymmetric-energy collider at the High Energy Accelerator Research Organization (KEK) in Japan. It uses a novel collision method referred to as the nano-beam scheme to avoid the so-called hourglass effect. The highest peak luminosity achieved thus far is 4.65 × 1034 cm-2s-1, which was reached on June 8th, 2022. Although the SuperKEKB holds the record for the world's highest instantaneous luminosity, several unresolved issues prevent a further increase in luminosity. One of the major issues is inadequate stability and reproducibility of the beam orbit and optics. This paper discusses optics distortion caused by orbit deviation near strong sextupole magnets. Numerical estimations indicate that a horizontal orbit deviation of a few tens of micrometers induces a sizable distortion of the optical functions at the interaction point. Analyses of the amount of betatron tune feedback and orbit at sextupole magnets used for crab waist and vertical local chromaticity correction indicate that the beam-current dependence of the vertical tune shift is attributable to the change of the beam orbit at the sextupoles. The orbit fluctuation at the sextupoles causes beta-beating and makes stable operation at high beam currents more difficult. The results show that orbit tuning at the sextupoles improves both the injection efficiency and the detector background level. Optics deterioration within a few days is another urgent issue in high-beam-current operation. A vertical orbit change of a few tens of microns at sextupole magnets installed in arc cells is a matter of concern according to numerical estimations.