Beam Dynamics in High Intensity Cyclotrons Including Neighboring Bunch Effects: Model, Implementation and Application

Space charge effects, being one of the most significant collective effects, play an important role in high intensity cyclotrons. However, for cyclotrons with small turn separation, other existing effects are of equal importance. Interactions of radially neighboring bunches are also present, but their combined effects has not yet been investigated in any great detail. In this paper, a new particle in cell based self-consistent numerical simulation model is presented for the first time. The model covers neighboring bunch effects and is implemented in the three-dimensional object-oriented parallel code OPAL-cycl, a flavor of the OPAL framework. We discuss this model together with its implementation and validation. Simulation results are presented from the PSI 590 MeV Ring Cyclotron in the context of the ongoing high intensity upgrade program, which aims to provide a beam power of 1.8 MW (CW) at the target destination

Two decades of numerical modelling to understand long term fluvial archives: Advances and future perspectives

The development and application of numerical models to investigate fluvial sedimentary archives has increased during the last decades resulting in a sustained growth in the number of scientific publications with keywords, 'fluvial models', 'fluvial process models' and 'fluvial numerical models'. In this context we compile and review the current contributions of numerical modelling to the understanding of fluvial archives. In particular, recent advances, current limitations, previous unexpected results and future perspectives are all discussed. Numerical modelling efforts have demonstrated that fluvial systems can display non-linear behaviour with often unexpected dynamics causing significant delay, amplification, attenuation or blurring of externally controlled signals in their simulated record. Numerical simulations have also demonstrated that fluvial records can be generated by intrinsic dynamics without any change in external controls. Many other model applications demonstrate that fluvial archives, specifically of large fluvial systems, can be convincingly simulated as a function of the interplay of (palaeo) landscape properties and extrinsic climate, base level and crustal controls. All discussed models can, after some calibration, produce believable matches with real world systems suggesting that equifinality - where a given end state can be reached through many different pathways starting from different initial conditions and physical assumptions - plays an important role in fluvial records and their modelling. The overall future challenge lies in the development of new methodologies for a more independent validation of system dynamics and research strategies that allow the separation of intrinsic and extrinsic record signals using combined fieldwork and modelling

Two-neutron transfer reaction mechanisms in 12^{12}C(6^6He,4^{4}He)14^{14}C using a realistic three-body 6^{6}He model

The reaction mechanisms of the two-neutron transfer reaction $^{12}$C($^6$He,$^4$He) have been studied at 30 MeV at the TRIUMF ISAC-II facility using the SHARC charged-particle detector array. Optical potential parameters have been extracted from the analysis of the elastic scattering angular distribution. The new potential has been applied to the study of the transfer angular distribution to the 2$^+_2$ 8.32 MeV state in $^{14}$C, using a realistic 3-body $^6$He model and advanced shell model calculations for the carbon structure, allowing to calculate the relative contributions of the simultaneous and sequential two-neutron transfer. The reaction model provides a good description of the 30 MeV data set and shows that the simultaneous process is the dominant transfer mechanism. Sensitivity tests of optical potential parameters show that the final results can be considerably affected by the choice of optical potentials. A reanalysis of data measured previously at 18 MeV however, is not as well described by the same reaction model, suggesting that one needs to include higher order effects in the reaction mechanism.Comment: 9 pages, 9 figure

The Acceleration and Storage of Radioactive Ions for a Beta-Beam Facility

The term beta-beam has been coined for the production of a pure beam of electron neutrinos or their antiparticles through the decay of radioactive ions circulating in a storage ring. This concept requires radioactive ions to be accelerated to as high Lorentz gamma as 150. The neutrino source itself consists of a storage ring for this energy range, with long straight sections in line with the experiment(s). Such a decay ring does not exist at CERN today, nor does a high-intensity proton source for the production of the radioactive ions. Nevertheless, the existing CERN accelerator infrastructure could be used as this would still represent an important saving for a beta-beam facility.Comment: beta-beam working group website at http://cern.ch/beta-bea

Two-Neutron Transfer Reaction Mechanisms in \u3csup\u3e12\u3c/sup\u3eC(\u3csup\u3e6\u3c/sup\u3eHe, \u3csup\u3e4\u3c/sup\u3eHe) \u3csup\u3e14\u3c/sup\u3eC using a Realistic Three-Body \u3csup\u3e6\u3c/sup\u3eHe Model

The reaction mechanisms of the two-neutron transfer reaction 12C(6He,4He) have been studied at Elab=30 MeV at the TRIUMF ISAC-II facility using the Silicon Highly-segmented Array for Reactions and Coulex (SHARC) charged-particle detector array. Optical potential parameters have been extracted from the analysis of the elastic scattering angular distribution. The new potential has been applied to the study of the transfer angular distribution to the 2+2 8.32 MeV state in 14C, using a realistic three-body 6He model and advanced shell-model calculations for the carbon structure, allowing to calculate the relative contributions of the simultaneous and sequential two-neutron transfer. The reaction model provides a good description of the 30-MeV data set and shows that the simultaneous process is the dominant transfer mechanism. Sensitivity tests of optical potential parameters show that the final results can be considerably affected by the choice of optical potentials. A reanalysis of data measured previously at Elab=18 MeV, however, is not as well described by the same reaction model, suggesting that one needs to include higher-order effects in the reaction mechanism

Impact of instrument and data characteristics in the interferometric reconstruction of the 21 cm power spectrum

Combining the visibilities measured by an interferometer to form a cosmological power spectrum is a complicated process. In a delay-based analysis, the mapping between instrumental and cosmological space is not a one-to-one relation. Instead, neighbouring modes contribute to the power measured at one point, with their respective contributions encoded in the window functions. To better understand the power measured by an interferometer, we assess the impact of instrument characteristics and analysis choices on these window functions. Focusing on the Hydrogen Epoch of Reionization Array (HERA) as a case study, we find that long-baseline observations correspond to enhanced low-k tails of the window functions, which facilitate foreground leakage, whilst an informed choice of bandwidth and frequency taper can reduce said tails. With simple test cases and realistic simulations, we show that, apart from tracing mode mixing, the window functions help accurately reconstruct the power spectrum estimator of simulated visibilities. The window functions depend strongly on the beam chromaticity and less on its spatial structure – a Gaussian approximation, ignoring side lobes, is sufficient. Finally, we investigate the potential of asymmetric window functions, down-weighting the contribution of low-k power to avoid foreground leakage. The window functions presented here correspond to the latest HERA upper limits for the full Phase I data. They allow an accurate reconstruction of the power spectrum measured by the instrument and will be used in future analyses to confront theoretical models and data directly in cylindrical space

Characterization of inpaint residuals in interferometric measurements of the epoch of reionization

To mitigate the effects of Radio Frequency Interference (RFI) on the data analysis pipelines of 21 cm interferometric instruments, numerous inpaint techniques have been developed. In this paper, we examine the qualitative and quantitative errors introduced into the visibilities and power spectrum due to inpainting. We perform our analysis on simulated data as well as real data from the Hydrogen Epoch of Reionization Array (HERA) Phase 1 upper limits. We also introduce a convolutional neural network that is capable of inpainting RFI corrupted data. We train our network on simulated data and show that our network is capable of inpainting real data without requiring to be retrained. We find that techniques that incorporate high wavenumbers in delay space in their modelling are best suited for inpainting over narrowband RFI. We show that with our fiducial parameters discrete prolate spheroidal sequences (DPSS) and CLEAN provide the best performance for intermittent RFI while Gaussian progress regression (GPR) and least squares spectral analysis (LSSA) provide the best performance for larger RFI gaps. However, we caution that these qualitative conclusions are sensitive to the chosen hyperparameters of each inpainting technique. We show that all inpainting techniques reliably reproduce foreground dominated modes in the power spectrum. Since the inpainting techniques should not be capable of reproducing noise realizations, we find that the largest errors occur in the noise dominated delay modes. We show that as the noise level of the data comes down, CLEAN and DPSS are most capable of reproducing the fine frequency structure in the visibilities

Improved Constraints on the 21 cm EoR Power Spectrum and the X-Ray Heating of the IGM with HERA Phase I Observations

We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits, we find at 95% confidence that Δ2(k = 0.34 h Mpc−1) ≤ 457 mK2 at z = 7.9 and that Δ2(k = 0.36 h Mpc−1) ≤ 3496 mK2 at z = 10.4, an improvement by a factor of 2.1 and 2.6, respectively. These limits are mostly consistent with thermal noise over a wide range of k after our data quality cuts, despite performing a relatively conservative analysis designed to minimize signal loss. Our results are validated with both statistical tests on the data and end-to-end pipeline simulations. We also report updated constraints on the astrophysics of reionization and the cosmic dawn. Using multiple independent modeling and inference techniques previously employed by HERA Collaboration, we find that the intergalactic medium must have been heated above the adiabatic cooling limit at least as early as z = 10.4, ruling out a broad set of so-called “cold reionization” scenarios. If this heating is due to high-mass X-ray binaries during the cosmic dawn, as is generally believed, our result’s 99% credible interval excludes the local relationship between soft X-ray luminosity and star formation and thus requires heating driven by evolved low-metallicity stars