Three Dimensional Relativistic Electromagnetic Sub-cycle Solitons

Three dimensional (3D) relativistic electromagnetic sub-cycle solitons were observed in 3D Particle-in-Cell simulations of an intense short laser pulse propagation in an underdense plasma. Their structure resembles that of an oscillating electric dipole with a poloidal electric field and a toroidal magnetic field that oscillate in-phase with the electron density with frequency below the Langmuir frequency. On the ion time scale the soliton undergoes a Coulomb explosion of its core, resulting in ion acceleration, and then evolves into a slowly expanding quasi-neutral cavity.Comment: 5 pages, 6 figures; http://www.ile.osaka-u.ac.jp/research/TSI/Timur/soliton/index.htm

Radiation Spectral Synthesis of Relativistic Filamentation

Radiation from many astrophysical sources, e.g. gamma-ray bursts and active galactic nuclei, is believed to arise from relativistically shocked collisionless plasmas. Such sources often exhibit highly transient spectra evolving rapidly, compared with source lifetimes. Radiation emitted from these sources is typically associated with non-linear plasma physics, complex field topologies and non-thermal particle distributions. In such circumstances a standard synchrotron paradigm may fail to produce accurate conclusions regarding the underlying physics. Simulating spectral emission and spectral evolution numerically in various relativistic shock scenarios is then the only viable method to determine the detailed physical origin of the emitted spectra. In this Letter we present synthetic radiation spectra representing the early stage development of the filamentation (streaming) instability of an initially unmagnetized plasma, which is relevant for both collisionless shock formation and reconnection dynamics in relativistic astrophysical outflows, as well as for laboratory astrophysics experiments. Results were obtained using a highly efficient "in situ" diagnostics method, based on detailed particle-in-cell modeling of collisionless plasmas. The synthetic spectra obtained here are compared with those predicted by a semi-analytical model for jitter radiation from the filamentation instability, the latter including self-consistent generated field topologies and particle distributions obtained from the simulations reported upon here. Spectra exhibit dependence on the presence - or absence - of an inert plasma constituent, when comparing baryonic plasmas (i.e. containing protons) with pair plasmas. The results also illustrate that considerable care should be taken when using lower-dimensional models to obtain information about the astrophysical phenomena generating observed spectra.Comment: 5 pages, 5 figures, accepted in Astrophysical Journal Letter

Computationally efficient methods for modelling laser wakefield acceleration in the blowout regime

Electron self-injection and acceleration until dephasing in the blowout regime is studied for a set of initial conditions typical of recent experiments with 100 terawatt-class lasers. Two different approaches to computationally efficient, fully explicit, three-dimensional particle-in-cell modelling are examined. First, the Cartesian code VORPAL using a perfect-dispersion electromagnetic solver precisely describes the laser pulse and bubble dynamics, taking advantage of coarser resolution in the propagation direction, with a proportionally larger time step. Using third-order splines for macroparticles helps suppress the sampling noise while keeping the usage of computational resources modest. The second way to reduce the simulation load is using reduced-geometry codes. In our case, the quasi-cylindrical code CALDER-CIRC uses decomposition of fields and currents into a set of poloidal modes, while the macroparticles move in the Cartesian 3D space. Cylindrical symmetry of the interaction allows using just two modes, reducing the computational load to roughly that of a planar Cartesian simulation while preserving the 3D nature of the interaction. This significant economy of resources allows using fine resolution in the direction of propagation and a small time step, making numerical dispersion vanishingly small, together with a large number of particles per cell, enabling good particle statistics. Quantitative agreement of the two simulations indicates that they are free of numerical artefacts. Both approaches thus retrieve physically correct evolution of the plasma bubble, recovering the intrinsic connection of electron self-injection to the nonlinear optical evolution of the driver

Steady-State Ab Initio Laser Theory for N-level Lasers

We show that Steady-state Ab initio Laser Theory (SALT) can be applied to find the stationary multimode lasing properties of an N-level laser. This is achieved by mapping the N-level rate equations to an effective two-level model of the type solved by the SALT algorithm. This mapping yields excellent agreement with more computationally demanding N-level time domain solutions for the steady state

Ion dynamics and coherent structure formation following laser pulse self-channeling

The propagation of a superintense laser pulse in an underdense, inhomogeneous plasma has been studied numerically by two-dimensional particle-in-cell simulations on a time scale extending up to several picoseconds. The effects of the ion dynamics following the charge-displacement self-channeling of the laser pulse have been addressed. Radial ion acceleration leads to the ``breaking'' of the plasma channel walls, causing an inversion of the radial space-charge field and the filamentation of the laser pulse. At later times a number of long-lived, quasi-periodic field structures are observed and their dynamics is characterized with high resolution. Inside the plasma channel, a pattern of electric and magnetic fields resembling both soliton- and vortex-like structures is observed.Comment: 10 pages, 5 figures (visit http://www.df.unipi.it/~macchi to download a high-resolution version), to appear in Plasma Physics and Controlled Fusion (Dec. 2007), special issue containing invited papers from the 34th EPS Conference on Plasma Physics (Warsaw, July 2007

Adiabatic Formation of Rydberg Crystals with Chirped Laser Pulses

Ultracold atomic gases have been used extensively in recent years to realize textbook examples of condensed matter phenomena. Recently, phase transitions to ordered structures have been predicted for gases of highly excited, 'frozen' Rydberg atoms. Such Rydberg crystals are a model for dilute metallic solids with tunable lattice parameters, and provide access to a wide variety of fundamental phenomena. We investigate theoretically how such structures can be created in four distinct cold atomic systems, by using tailored laser-excitation in the presence of strong Rydberg-Rydberg interactions. We study in detail the experimental requirements and limitations for these systems, and characterize the basic properties of small crystalline Rydberg structures in one, two and three dimensions.Comment: 23 pages, 10 figures, MPIPKS-ITAMP Tandem Workshop, Cold Rydberg Gases and Ultracold Plasmas (CRYP10), Sept. 6-17, 201

Marangoni driven turbulence in high energy surface melting processes

Experimental observations of high-energy surface melting processes, such as laser welding, have revealed unsteady, often violent, motion of the free surface of the melt pool. Surprisingly, no similar observations have been reported in numerical simulation studies of such flows. Moreover, the published simulation results fail to predict the post-solidification pool shape without adapting non-physical values for input parameters, suggesting the neglect of significant physics in the models employed. The experimentally observed violent flow surface instabilities, scaling analyses for the occurrence of turbulence in Marangoni driven flows, and the fact that in simulations transport coefficients generally have to be increased by an order of magnitude to match experimentally observed pool shapes, suggest the common assumption of laminar flow in the pool may not hold, and that the flow is actually turbulent. Here, we use direct numerical simulations (DNS) to investigate the role of turbulence in laser melting of a steel alloy with surface active elements. Our results reveal the presence of two competing vortices driven by thermocapillary forces towards a local surface tension maximum. The jet away from this location at the free surface, separating the two vortices, is found to be unstable and highly oscillatory, indeed leading to turbulence-like flow in the pool. The resulting additional heat transport, however, is insufficient to account for the observed differences in pool shapes between experiment and simulations