An efficient and portable SIMD algorithm for charge/current deposition in Particle-In-Cell codes

In current computer architectures, data movement (from die to network) is by far the most energy consuming part of an algorithm (10pJ/word on-die to 10,000pJ/word on the network). To increase memory locality at the hardware level and reduce energy consumption related to data movement, future exascale computers tend to use more and more cores on each compute nodes ("fat nodes") that will have a reduced clock speed to allow for efficient cooling. To compensate for frequency decrease, machine vendors are making use of long SIMD instruction registers that are able to process multiple data with one arithmetic operator in one clock cycle. SIMD register length is expected to double every four years. As a consequence, Particle-In-Cell (PIC) codes will have to achieve good vectorization to fully take advantage of these upcoming architectures. In this paper, we present a new algorithm that allows for efficient and portable SIMD vectorization of current/charge deposition routines that are, along with the field gathering routines, among the most time consuming parts of the PIC algorithm. Our new algorithm uses a particular data structure that takes into account memory alignement constraints and avoids gather/scatter instructions that can significantly affect vectorization performances on current CPUs. The new algorithm was successfully implemented in the 3D skeleton PIC code PICSAR and tested on Haswell Xeon processors (AVX2-256 bits wide data registers). Results show a factor of $\times 2$ to $\times 2.5$ speed-up in double precision for particle shape factor of order $1$ to $3$. The new algorithm can be applied as is on future KNL (Knights Landing) architectures that will include AVX-512 instruction sets with 512 bits register lengths (8 doubles/16 singles).Comment: 36 pages, 5 figure

Neutron Scattering to Characterize Cu/Mg(Li) Destabilized Hydrogen Storage Materials

Cu-Li-Mg-(H,D) was studied as an example of destabilizer of the Ti-(H,D) system. A Cu-Li-Mg alloy was prepared resulting in the formation of a system with 60.5 at% of CuLi0.08Mg1.92, 23.9 at% of CuMg 2 and 15.6 at% of Cu2Mg. Titanium was added to a fraction of this mixture so that 68.2 at% (47.3 wt%) of the final mixture was Ti. The mixture was ground and kept at 200 °C/473 K for 7h under H2 or 9h under D2 at P = 34 bar. Under those conditions, neutron powder diffraction shows the formation of TiD2, as well as of the deuteride of CuLi0.08Mg1.92. Similarly inelastic neutron scattering shows that at 10 K TiH2 is present in the sample, together with the hydride of CuLi0.08Mg1.92. Interestingly, at 10 K TiH 2 is very clearly detected and at 300 K TiH2 is still clearly present as indicated by the neutron vibrational spectrum, but CuLi 0.08Mg1.92-H is not detected anymore. These results indicate that Ti(H,D)2 is possibly formed by diffusion of hydrogen from the Cu-Li-Mg-(H,D) alloys. This is an intriguing result since TiH 2 is normally synthesized from the metal at T > 400°C/673 K (and most commonly at T ∼ 700 °C/973 K). In the presence of CuLi 0.08Mg1.92, TiH2 forms at a temperature that is 300 - 400 K lower than that needed to synthesize it just from the elements

γ-ray generation enhancement by the charge separation field in laser-target interaction in the radiation dominated regime

A new source of radiation can be created with a laser pulse of intensity 1023W/cm2 interacting with a slightly overdense plasma. Collective effects driven by the electrostatic field significantly enhance the synchrotron radiation. They impact on the laser energy repartition leading to a specific emission but also constitute a crucial element for the intense radiation production. They allow electrons to be accelerated over a length up to 10 laser wavelengths favoring emission of an intense radiation. It is shown that charge separation field depends on the ion mass and target thickness but also on laser polarization. These phenomena are studied with an one dimensional relativistic particle-in-cell code accounting for the classical radiation reaction force

Comprehensive care on paper only? The challenge for physiotherapy provision in day to day haemophilia practice

Comprehensive care on paper only? The challenge for physiotherapy provision in day to day haemophilia practic

Neutron Scattering to Characterize Cu/Mg(Li) Destabilized Hydrogen Storage Materials

Cu-Li-Mg-(H,D) was studied as an example of destabilizer of the Ti-(H,D) system. A Cu-Li-Mg alloy was prepared resulting in the formation of a system with 60.5 at% of CuLi0.08Mg1.92, 23.9 at% of CuMg 2 and 15.6 at% of Cu2Mg. Titanium was added to a fraction of this mixture so that 68.2 at% (47.3 wt%) of the final mixture was Ti. The mixture was ground and kept at 200 °C/473 K for 7h under H2 or 9h under D2 at P = 34 bar. Under those conditions, neutron powder diffraction shows the formation of TiD2, as well as of the deuteride of CuLi0.08Mg1.92. Similarly inelastic neutron scattering shows that at 10 K TiH2 is present in the sample, together with the hydride of CuLi0.08Mg1.92. Interestingly, at 10 K TiH 2 is very clearly detected and at 300 K TiH2 is still clearly present as indicated by the neutron vibrational spectrum, but CuLi 0.08Mg1.92-H is not detected anymore. These results indicate that Ti(H,D)2 is possibly formed by diffusion of hydrogen from the Cu-Li-Mg-(H,D) alloys. This is an intriguing result since TiH 2 is normally synthesized from the metal at T > 400°C/673 K (and most commonly at T ∼ 700 °C/973 K). In the presence of CuLi 0.08Mg1.92, TiH2 forms at a temperature that is 300 - 400 K lower than that needed to synthesize it just from the elements