Cosmic Sculpture: A new way to visualise the Cosmic Microwave Background

3D printing presents an attractive alternative to visual representation of physical datasets such as astronomical images that can be used for research, outreach or teaching purposes, and is especially relevant to people with a visual disability. We here report the use of 3D printing technology to produce a representation of the all-sky Cosmic Microwave Background (CMB) intensity anisotropy maps produced by the Planck mission. The success of this work in representing key features of the CMB is discussed as is the potential of this approach for representing other astrophysical data sets. 3D printing such datasets represents a highly complementary approach to the usual 2D projections used in teaching and outreach work, and can also form the basis of undergraduate projects. The CAD files used to produce the models discussed in this paper are made available.Comment: Accepted for publication in the European Journal of Physic

Controlled Shock Shells and Intracluster Fusion Reactions in the Explosion of Large Clusters

The ion phase-space dynamics in the Coulomb explosion of very large ($\sim 10^6 - 10^7$ atoms) deuterium clusters can be tailored using two consecutive laser pulses with different intensities and an appropriate time delay. For suitable sets of laser parameters (intensities and delay), large-scale shock shells form during the explosion, thus highly increasing the probability of fusion reactions within the single exploding clusters. In order to analyze the ion dynamics and evaluate the intracluster reaction rate, a one-dimensional theory is used, which approximately accounts for the electron expulsion from the clusters. It is found that, for very large clusters (initial radius $\sim$ 100 nm), and optimal laser parameters, the intracluster fusion yield becomes comparable to the intercluster fusion yield. The validity of the results is confirmed with three-dimensional particle-in-cell simulations.Comment: 25 pages, 11 figures, to appear in Physical Review

DC magnetic field generation in unmagnetized shear flows

The generation of DC magnetic fields in unmagnetized plasmas with velocity shear is predicted for non relativistic and relativistic scenarios either due to thermal effects or due to the onset of the Kelvin-Helmholtz instability (KHI). A kinetic model describes the growth and the saturation of the DC field. The predictions of the theory are confirmed by multidimensional particle-in-cell simulations, demonstrating the formation of long lived magnetic fields ($t \sim 100s \omega_{pi}^{-1}$) along the full longitudinal extent of the shear layer, with transverse width on the electron length scale ($\sqrt{\gamma_0}c/\omega_{pe}$), reaching magnitudes $eB_{\mathrm{DC}}/m_ec\omega_{pe}\sim \beta_0\sqrt{\gamma_0}$

The impact of kinetic effects on the properties of relativistic electron-positron shocks

We assess the impact of non-thermally shock-accelerated particles on the magnetohydrodynamic (MHD) jump conditions of relativistic shocks. The adiabatic constant is calculated directly from first principle particle-in-cell simulation data, enabling a semi-kinetic approach to improve the standard fluid model and allowing for an identification of the key parameters that define the shock structure. We find that the evolving upstream parameters have a stronger impact than the corrections due to non-thermal particles. We find that the decrease of the upstream bulk speed yields deviations from the standard MHD model up to 10%. Furthermore, we obtain a quantitative definition of the shock transition region from our analysis. For Weibel-mediated shocks the inclusion of a magnetic field in the MHD conservation equations is addressed for the first time

Electromagnetic field generation in the downstream of electrostatic shocks due to electron trapping

A new magnetic field generation mechanism in electrostatic shocks is found, which can produce fields with magnetic energy density as high as 0.01 of the kinetic energy density of the flows on time scales $\tilde \, 10^4 \, {\omega}_{pe}^{-1}$. Electron trapping during the shock formation process creates a strong temperature anisotropy in the distribution function, giving rise to the pure Weibel instability. The generated magnetic field is well-confined to the downstream region of the electrostatic shock. The shock formation process is not modified and the features of the shock front responsible for ion acceleration, which are currently probed in laser-plasma laboratory experiments, are maintained. However, such a strong magnetic field determines the particle trajectories downstream and has the potential to modify the signatures of the collisionless shock

Electron-scale shear instabilities: magnetic field generation and particle acceleration in astrophysical jets

Strong shear flow regions found in astrophysical jets are shown to be important dissipation regions, where the shear flow kinetic energy is converted into electric and magnetic field energy via shear instabilities. The emergence of these self-consistent fields make shear flows significant sites for radiation emission and particle acceleration. We focus on electron-scale instabilities, namely the collisionless, unmagnetized Kelvin-Helmholtz instability (KHI) and a large-scale dc magnetic field generation mechanism on the electron scales. We show that these processes are important candidates to generate magnetic fields in the presence of strong velocity shears, which may naturally originate in energetic matter outburst of active galactic nuclei and gamma-ray bursters. We show that the KHI is robust to density jumps between shearing flows, thus operating in various scenarios with different density contrasts. Multidimensional particle-in-cell (PIC) simulations of the KHI, performed with OSIRIS, reveal the emergence of a strong and large-scale dc magnetic field component, which is not captured by the standard linear fluid theory. This dc component arises from kinetic effects associated with the thermal expansion of electrons of one flow into the other across the shear layer, whilst ions remain unperturbed due to their inertia. The electron expansion forms dc current sheets, which induce a dc magnetic field. Our results indicate that most of the electromagnetic energy developed in the KHI is stored in the dc component, reaching values of equipartition on the order of $10^{-3}$ in the electron time-scale, and persists longer than the proton time-scale. Particle scattering/acceleration in the self generated fields of these shear flow instabilities is also analyzed

Classical Radiation Reaction in Particle-In-Cell Simulations

Under the presence of ultra high intensity lasers or other intense electromagnetic fields the motion of particles in the ultrarelativistic regime can be severely affected by radiation reaction. The standard particle-in-cell (PIC) algorithms do not include radiation reaction effects. Even though this is a well known mechanism, there is not yet a definite algorithm nor a standard technique to include radiation reaction in PIC codes. We have compared several models for the calculation of the radiation reaction force, with the goal of implementing an algorithm for classical radiation reaction in the Osiris framework, a state-of-the-art PIC code. The results of the different models are compared with standard analytical results, and the relevance/advantages of each model are discussed. Numerical issues relevant to PIC codes such as resolution requirements, application of radiation reaction to macro particles and computational cost are also addressed. The Landau and Lifshitz reduced model is chosen for implementation.Comment: 12 pages, 8 figure

Three-body Faddeev-Alt-Grassberger-Sandhas approach to direct nuclear reactions

Momentum space three-body Faddeev-like equations are used to calculate elastic, transfer and charge exchange reactions resulting from the scattering of deuterons on 12C and 16O or protons on 13C and 17O; 12C and 16O are treated as inert cores. All possible reactions are calculated in the framework of the same model space. Comparison with previous calculations based on approximate methods used in nuclear reaction theory is discussed.Comment: 10 pages, 13 figures, to be published in Phys. Rev.