Cavity dumping of an injection-locked free-electron laser

This letter reports cavity dumping of an electrostatic-accelerator-driven free-electron laser (FEL) while it is injection-locked to a frequency-stabilized 240 GHz solid-state source. Cavity dumping enhances the FEL output power by a factor of $\sim$8, and abruptly cuts off the end of the FEL pulse. The cavity-dumped, injection-locked FEL output is used in a 240 GHz pulsed electron spin resonance (ESR) experiment.Comment: 8 pages including 3 figure

Optimising Sparse Matrix Vector multiplication for large scale FEM problems on FPGA

Sparse Matrix Vector multiplication (SpMV) is an important kernel in many scientific applications. In this work we propose an architecture and an automated customisation method to detect and optimise the architecture for block diagonal sparse matrices. We evaluate the proposed approach in the context of the spectral/hp Finite Element Method, using the local matrix assembly approach. This problem leads to a large sparse system of linear equations with block diagonal matrix which is typically solved using an iterative method such as the Preconditioned Conjugate Gradient. The efficiency of the proposed architecture combined with the effectiveness of the proposed customisation method reduces BRAM resource utilisation by as much as 10 times, while achieving identical throughput with existing state of the art designs and requiring minimal development effort from the end user. In the context of the Finite Element Method, our approach enables the solution of larger problems than previously possible, enabling the applicability of FPGAs to more interesting HPC problems

Cosmological Constraints from Moments of the Thermal Sunyaev-Zel'dovich Effect

In this paper, we explain how moments of the thermal Sunyaev-Zel'dovich (tSZ) effect can constrain both cosmological parameters and the astrophysics of the intracluster medium (ICM). As the tSZ signal is strongly non-Gaussian, higher moments of tSZ maps contain useful information. We first calculate the dependence of the tSZ moments on cosmological parameters, finding that higher moments scale more steeply with sigma_8 and are sourced by more massive galaxy clusters. Taking advantage of the different dependence of the variance and skewness on cosmological and astrophysical parameters, we construct a statistic, ||/^1.4, which cancels much of the dependence on cosmology (i.e., sigma_8) yet remains sensitive to the astrophysics of intracluster gas (in particular, to the gas fraction in low-mass clusters). Constraining the ICM astrophysics using this statistic could break the well-known degeneracy between cosmology and gas physics in tSZ measurements, allowing for tight constraints on cosmological parameters. Although detailed simulations will be needed to fully characterize the accuracy of this technique, we provide a first application to data from the Atacama Cosmology Telescope and the South Pole Telescope. We estimate that a Planck-like full-sky tSZ map could achieve a <1% constraint on sigma_8 and a 1-sigma error on the sum of the neutrino masses that is comparable to the existing lower bound from oscillation measurements.Comment: 11 pages, 12 figures, to be submitted to Phys. Rev. D; v2: 14 pages, 16 figures, matches PRD accepted version (changes from v1 include additional calculations with primordial non-Gaussianity and a new appendix discussing the tSZ kurtosis

Dealiasing techniques for high-order spectral element methods on regular and irregular grids

High-order methods are becoming increasingly attractive in both academia and industry, especially in the context of computational fluid dynamics. However, before they can be more widely adopted, issues such as lack of robustness in terms of numerical stability need to be addressed, particularly when treating industrial-type problems where challenging geometries and a wide range of physical scales, typically due to high Reynolds numbers, need to be taken into account. One source of instability is aliasing effects which arise from the nonlinearity of the underlying problem. In this work we detail two dealiasing strategies based on the concept of consistent integration. The first uses a localised approach, which is useful when the nonlinearities only arise in parts of the problem. The second is based on the more traditional approach of using a higher quadrature. The main goal of both dealiasing techniques is to improve the robustness of high order spectral element methods, thereby reducing aliasing-driven instabilities. We demonstrate how these two strategies can be effectively applied to both continuous and discontinuous discretisations, where, in the latter, both volumetric and interface approximations must be considered. We show the key features of each dealiasing technique applied to the scalar conservation law with numerical examples and we highlight the main differences in terms of implementation between continuous and discontinuous spatial discretisations