Transport of laser accelerated proton beams and isochoric heating of matter

The acceleration of intense proton and ion beams by ultra-intense lasers has matured to a point where applications in basic research and technology are being developed. Crucial for harvesting the unmatched beam parameters driven by the relativistic electron sheath is the precise control of the beam. We report on recent experiments using the PHELIX laser at GSI, the VULCAN laser at RAL and the TRIDENT laser at LANL to control and use laser accelerated proton beams for applications in high energy density research. We demonstrate efficient collimation of the proton beam using high field pulsed solenoid magnets, a prerequisite to capture and transport the beam for applications. Furthermore we report on two campaigns to use intense, short proton bunches to isochorically heat solid targets up to the warm dense matter state. The temporal profile of the proton beam allows for rapid heating of the target, much faster than the hydrodynamic response time thereby creating a strongly coupled plasma at solid density. The target parameters are then probed by X-ray Thomson scattering (XRTS) to reveal the density and temperature of the heated volume. This combination of two powerful techniques developed during the past few years allows for the generation and investigation of macroscopic samples of matter in states present in giant planets or the interior of the earth.

Observation of the inhomogeneous spatial distribution of MeV ions accelerated by the hydrodynamic ambipolar expansion of clusters

An inhomogeneous spatial distribution of laser accelerated carbon/oxygen ions produced via the hydrodynamic ambipolar expansion of CO₂ clusters has been measured by using CR-39 detectors. An inhomogeneous etch pits spatial distribution has appeared on the etched CR-39 detector installed on the laser propagation direction, while homogeneous ones are appeared on those installed at 45°and 90°from the laser propagation direction. From the range of ions in CR-39 obtained by using the multi-step etching technique, the averaged energies of carbon/oxygen ions for all directions are determined as 0.78 ± 0.09 MeV/n. The number of ions in the laser propagation direction is about 1.5 times larger than those in other directions. The inhomogeneous etch pits spatial distribution in the laser propagation direction could originate from an ion beam collimation and modulation by the effect of electromagnetic structures created in the laser plasma.

Picosecond metrology of laser-driven proton bursts

Tracking primary radiation-induced processes in matter requires ultrafast sources and high precision timing. While compact laser-driven ion accelerators are seeding the development of novel high instantaneous flux applications, combining the ultrashort ion and laser pulse durations with their inherent synchronicity to trace the real-time evolution of initial damage events has yet to be realized. Here we report on the absolute measurement of proton bursts as short as 3.5±0.7 ps from laser solid target interactions for this purpose. Our results verify that laser-driven ion acceleration can deliver interaction times over a factor of hundred shorter than those of state-of-the-art accelerators optimized for high instantaneous flux. Furthermore, these observations draw ion interaction physics into the field of ultrafast science, opening the opportunity for quantitative comparison with both numerical modelling and the adjacent fields of ultrafast electron and photon interactions in matter.

Buffered high charge spectrally-peaked proton beams in the relativistic-transparency regime

Spectrally-peaked proton beams of high charge (E_p ≈ 8 MeV, ΔE ≈ 4 MeV, N ≈ 50 nC ) have been observed from the interaction of an intense laser (>10¹⁹Wcm^-2) with ultrathinCHfoils, as measured by spectrally-resolved full beam profiles. These beams are reproducibly generated for foil thicknesses 5-100 nm, and exhibit narrowing divergence with decreasing target thickness down to ≈8° for 5 nm. Simulations demonstrate that the narrow energy spread feature is a result of buffered acceleration of protons. The radiation pressure at the front of the target results in asymmetric sheath fields which permeate throughout the target, causing preferential forward acceleration. Due to their higher chargeto-mass ratio, the protons outrun a carbon plasma driven in the relativistic transparency regime.

Location and characteristics of the reconnection X line deduced from low-altitude satellite and ground-based observations: 1. Theory

A method for estimating both the Alfvén speed and the field-aligned flow of the magnetosheath at the magnetopause reconnection site is presented. The method employs low-altitude cusp ion observations and requires the identification of a feature in the cusp ion spectra near the low-energy cutoff which will often be present for a low-latitude dayside reconnection site. The appearance of these features in data of limited temporal, energy, and pitch angle resolution is illustrated by using model calculations of cusp ion distribution functions. These are based on the theory of ion acceleration at the dayside magnetopause and allow for the effects on the spectrum of flight times of ions precipitating down newly opened field lines. In addition, the variation of the reconnection rate can be evaluated, and comparison with ground-based observations of the corresponding sequence of transient events allows the field-aligned distance from the ionosphere to the reconnection site to be estimated.

Studio teorico e modellizzazione numerica degli effetti di preplasma in regime di accelerazione di ioni tramite laser

Argomento di questo lavoro di tesi è l’accelerazione di ioni tramite interazione laser-plasma. Nel particolare, uno studio teorico sui principali meccanismi di accelerazione e una campagna di simulazioni numeriche volte ad analizzare gli effetti sullo spettro energetico dei protoni accelerati dovuti al preplasma sono stati svolti. Quando laser ad impulsi brevi, ad alta intensità e a contrasto finito, interagiscono con un campione solido è presente un preimpulso che causa la formazione di una regione di preplasma antistante il bersaglio che può rendere più efficace l’accelerazione dei protoni. Vengono dunque introdotti nel primo capitolo i concetti e le tecniche base per studiare a livello teorico e numericamente le interazioni laser-plasma. Nel secondo capitolo sono trattati analiticamente alcuni meccanismi di accelerazione di ioni. Nel terzo ed ultimo capitolo vengono descritti gli effetti di preplasma ed esposti risultati di simulazioni PIC volte a trovare i parametri ottimali per ottenere ioni più energetici.

Código híbrido avanzado de motores de plasma de efecto Hall

La propulsión eléctrica constituye hoy una tecnología muy competitiva y de gran proyección de futuro. Dentro de los diversos motores de plasma existentes, el motor de efecto Hall ha adquirido una gran madurez y constituye un medio de propulsión idóneo para un rango amplio de misiones. En la presente Tesis se estudian los motores Hall con geometría convencional y paredes dieléctricas. La compleja interacción entre los múltiples fenómenos físicos presentes hace que sea difícil la simulación del plasma en estos motores. Los modelos híbridos son los que representan un mejor compromiso entre precisión y tiempo de cálculo. Se basan en utilizar un modelo fluido para los electrones y algoritmos de dinámica de partículas PIC (Particle-In- Cell) para los iones y los neutros. Permiten hacer uso de la hipótesis de cuasineutralidad del plasma, a cambio de resolver separadamente las capas límite (o vainas) que se forman en torno a las paredes de la cámara. Partiendo de un código híbrido existente, llamado HPHall-2, el objetivo de la Tesis doctoral ha sido el desarrollo de un código híbrido avanzado que mejorara la simulación de la descarga de plasma en un motor de efecto Hall. Las actualizaciones y mejoras realizadas en las diferentes partes que componen el código comprenden tanto aspectos teóricos como numéricos. Fruto de la extensa revisión de la algoritmia del código HPHall-2 se han conseguido reducir los errores de precisión un orden de magnitud, y se ha incrementado notablemente su consistencia y robustez, permitiendo la simulación del motor en un amplio rango de condiciones. Algunos aspectos relevantes a destacar en el subcódigo de partículas son: la implementación de un nuevo algoritmo de pesado que permite determinar de forma más precisa el flujo de las magnitudes del plasma; la implementación de un nuevo algoritmo de control de población, que permite tener suficiente número de partículas cerca de las paredes de la cámara, donde los gradientes son mayores y las condiciones de cálculo son más críticas; las mejoras en los balances de masa y energía; y un mejor cálculo del campo eléctrico en una malla no uniforme. Merece especial atención el cumplimiento de la condición de Bohm en el borde de vaina, que en los códigos híbridos representa una condición de contorno necesaria para obtener una solución consistente con el modelo de interacción plasma-pared, y que en HPHall-2 aún no se había resuelto satisfactoriamente. En esta Tesis se ha implementado el criterio cinético de Bohm para una población de iones con diferentes cargas eléctricas y una gran dispersión de velocidades. En el código, el cumplimiento de la condición cinética de Bohm se consigue por medio de un algoritmo que introduce una fina capa de aceleración nocolisional adyacente a la vaina y mide adecuadamente el flujo de partículas en el espacio y en el tiempo. Las mejoras realizadas en el subcódigo de electrones incrementan la capacidad de simulación del código, especialmente en la región aguas abajo del motor, donde se simula la neutralización del chorro del plasma por medio de un modelo de cátodo volumétrico. Sin abordar el estudio detallado de la turbulencia del plasma, se implementan modelos sencillos de ajuste de la difusión anómala de Bohm, que permiten reproducir los valores experimentales del potencial y la temperatura del plasma, así como la corriente de descarga del motor. En cuanto a los aspectos teóricos, se hace especial énfasis en la interacción plasma-pared y en la dinámica de los electrones secundarios libres en el interior del plasma, cuestiones que representan hoy en día problemas abiertos en la simulación de los motores Hall. Los nuevos modelos desarrollados buscan una imagen más fiel a la realidad. Así, se implementa el modelo de vaina de termalización parcial, que considera una función de distribución no-Maxwelliana para los electrones primarios y contabiliza unas pérdidas energéticas más cercanas a la realidad. Respecto a los electrones secundarios, se realiza un estudio cinético simplificado para evaluar su grado de confinamiento en el plasma, y mediante un modelo fluido en el límite no-colisional, se determinan las densidades y energías de los electrones secundarios libres, así como su posible efecto en la ionización. El resultado obtenido muestra que los electrones secundarios se pierden en las paredes rápidamente, por lo que su efecto en el plasma es despreciable, no así en las vainas, donde determinan el salto de potencial. Por último, el trabajo teórico y de simulación numérica se complementa con el trabajo experimental realizado en el Pnnceton Plasma Physics Laboratory, en el que se analiza el interesante transitorio inicial que experimenta el motor en el proceso de arranque. Del estudio se extrae que la presencia de gases residuales adheridos a las paredes juegan un papel relevante, y se recomienda, en general, la purga completa del motor antes del modo normal de operación. El resultado final de la investigación muestra que el código híbrido desarrollado representa una buena herramienta de simulación de un motor Hall. Reproduce adecuadamente la física del motor, proporcionando resultados similares a los experimentales, y demuestra ser un buen laboratorio numérico para estudiar el plasma en el interior del motor. Abstract Electric propulsion is today a very competitive technology and has a great projection into the future. Among the various existing plasma thrusters, the Hall effect thruster has acquired a considerable maturity and constitutes an ideal means of propulsion for a wide range of missions. In the present Thesis only Hall thrusters with conventional geometry and dielectric walls are studied. The complex interaction between multiple physical phenomena makes difficult the plasma simulation in these engines. Hybrid models are those representing a better compromise between precision and computational cost. They use a fluid model for electrons and Particle-In-Cell (PIC) algorithms for ions and neutrals. The hypothesis of plasma quasineutrality is invoked, which requires to solve separately the sheaths formed around the chamber walls. On the basis of an existing hybrid code, called HPHall-2, the aim of this doctoral Thesis is to develop an advanced hybrid code that better simulates the plasma discharge in a Hall effect thruster. Updates and improvements of the code include both theoretical and numerical issues. The extensive revision of the algorithms has succeeded in reducing the accuracy errors in one order of magnitude, and the consistency and robustness of the code have been notably increased, allowing the simulation of the thruster in a wide range of conditions. The most relevant achievements related to the particle subcode are: the implementation of a new weighing algorithm that determines more accurately the plasma flux magnitudes; the implementation of a new algorithm to control the particle population, assuring enough number of particles near the chamber walls, where there are strong gradients and the conditions to perform good computations are more critical; improvements in the mass and energy balances; and a new algorithm to compute the electric field in a non-uniform mesh. It deserves special attention the fulfilment of the Bohm condition at the edge of the sheath, which represents a boundary condition necessary to match consistently the hybrid code solution with the plasma-wall interaction, and remained as a question unsatisfactory solved in the HPHall-2 code. In this Thesis, the kinetic Bohm criterion has been implemented for an ion particle population with different electric charges and a large dispersion in their velocities. In the code, the fulfilment of the kinetic Bohm condition is accomplished by an algorithm that introduces a thin non-collisional layer next to the sheaths, producing the ion acceleration, and measures properly the flux of particles in time and space. The improvements made in the electron subcode increase the code simulation capabilities, specially in the region downstream of the thruster, where the neutralization of the plasma jet is simulated using a volumetric cathode model. Without addressing the detailed study of the plasma turbulence, simple models for a parametric adjustment of the anomalous Bohm difussion are implemented in the code. They allow to reproduce the experimental values of the plasma potential and the electron temperature, as well as the discharge current of the thruster. Regarding the theoretical issues, special emphasis has been made in the plasma-wall interaction of the thruster and in the dynamics of free secondary electrons within the plasma, questions that still remain unsolved in the simulation of Hall thrusters. The new developed models look for results closer to reality, such as the partial thermalization sheath model, that assumes a non-Maxwellian distribution functions for primary electrons, and better computes the energy losses at the walls. The evaluation of secondary electrons confinement within the chamber is addressed by a simplified kinetic study; and using a collisionless fluid model, the densities and energies of free secondary electrons are computed, as well as their effect on the plasma ionization. Simulations show that secondary electrons are quickly lost at walls, with a negligible effect in the bulk of the plasma, but they determine the potential fall at sheaths. Finally, numerical simulation and theoretical work is complemented by the experimental work carried out at the Princeton Plasma Physics Laboratory, devoted to analyze the interesting transitional regime experienced by the thruster in the startup process. It is concluded that the gas impurities adhered to the thruster walls play a relevant role in the transitional regime and, as a general recomendation, a complete purge of the thruster before starting its normal mode of operation it is suggested. The final result of the research conducted in this Thesis shows that the developed code represents a good tool for the simulation of Hall thrusters. The code reproduces properly the physics of the thruster, with results similar to the experimental ones, and represents a good numerical laboratory to study the plasma inside the thruster.

Ion cascade acceleration from the interaction of a relativistic femtosecond laser pulse with a narrow thin target

Particle-in-cell simulations are performed to study the acceleration of ions due to the interaction of a relativistic femtosecond laser pulse with a narrow thin target. The numerical results show that ions can be accelerated in a cascade by two electrostatic fields if the width of the target is smaller than the laser beam waist. The first field is formed in front of the target by the central part of the laser beam, which pushes the electron layer inward. The major part of the abaxial laser energy propagates along the edges to the rear side of the target and pulls out some hot electrons from the edges of the target, which form another electrostatic field at the rear side of the target. The ions from the front surface are accelerated stepwise by these two electrostatic fields to high energies at the rear side of the target. The simulations show that the largest ion energy gain for a narrow target is about four times higher than in the case of a wide target. (c) 2006 American Institute of Physics.

Radiation-Pressure Acceleration of Ion Beams from Nanofoil Targets: The Leaky Light-Sail Regime

A new ion radiation-pressure acceleration regime, the "leaky light sail," is proposed which uses sub-skin-depth nanometer foils irradiated by circularly polarized laser pulses. In the regime, the foil is partially transparent, continuously leaking electrons out along with the transmitted laser field. This feature can be exploited by a multispecies nanofoil configuration to stabilize the acceleration of the light ion component, supplementing the latter with an excess of electrons leaked from those associated with the heavy ions to avoid Coulomb explosion. It is shown by 2D particle-in-cell simulations that a monoenergetic proton beam with energy 18 MeV is produced by circularly polarized lasers at intensities of just 10(19) W/cm(2). 100 MeV proton beams are obtained by increasing the intensities to 2 x 10(20) W/cm(2).

Superthermal ion signatures of auroral acceleration processes

The retarding ion mass spectrometer on the Dynamics Explorer 1 spacecraft has generated a unique data set which documents, among other things, the occurrence of non-Maxwellian superthermal features in the auroral topside ionosphere distribution functions. In this paper, we provide a representative sampling of the observed features and their spatial morphology as observed at altitudes in the range from a few thousand kilometers to a few earth radii. At lower altitudes, these features appear at auroral latitudes separating regions of polar cap and subauroral light ion polar wind. The most common signature is the appearance of an upgoing energetic tail having conical lobes representing significant ion heat and number flux in all species, including O+. Transverse ion heating below the observation point at several thousand kilometers is clearly associated with O+ outflows. In some events observed, transverse acceleration apparently involves nearly the entire thermal plasma, the distribution function becomes highly anisotropic with T⊥ > T∥, and may actually develop a minimum at zero velocity, i.e., become a torus having as its axis the local magnetic field direction. At higher altitudes, the localized dayside source region appears as a field aligned flow which is dispersed tailward across the polar cap according to parallel velocity by antisunward convective flow, so that upflowing low energy O+ ions appear well within the polar cap region. While this flow can appear beamlike in a given location, the energy dispersion observed implies a very broad energy distribution at the source, extending from a few tenths of an eV to in excess of 50 eV. On the nightside, upgoing ion beams are found to be latitudinally bounded by regions of ion conics whose half angles increase with increasing separation from the beam region, indicating low altitude transverse acceleration in immediate proximity to, and below, the parallel acceleration region. These observations reveal a clear distinction between classical polar wind ion outflow and O+ enhanced superthermal flows, and confirm the importance of low altitude transverse acceleration in ionospheric plasma transport, as suggested by previous observations.

Thermal ion flows in the topside auroral ionosphere and the effects of low-altitude, transverse acceleration

Topside ionospheric profiles are used to study the upward field-aligned flow of thermal O+ at high latitudes. On the majority of the field lines outside the plasmasphere, the mean flux is approximately equal to the mean polar wind measured by spacecraft at greater altitudes. This is consistent with the theory of thermal light ion escape supported, via charge exchange, by upward O+ flow at lower heights. Events of larger O+ flow are detected at auroral latitudes and their occurrence is found to agree with that of transversely accelerated ions within the topside ionosphere and the magnetosphere. The effects of low altitude heating of O+ by oxygen cyclotron waves, driven by downward field-aligned currents, are considered as a possible common cause of these two types of event.

Shock acceleration in the solar corona

In this thesis acceleration of energetic particles at collisionless shock waves in space plasmas is studied using numerical simulations, with an emphasis on physical conditions applicable to the solar corona. The thesis consists of four research articles and an introductory part that summarises the main findings reached in the articles and discusses them with respect to theory of diffusive shock acceleration and observations. This thesis gives a brief review of observational properties of solar energetic particles and discusses a few open questions that are currently under active research. For example, in a few large gradual solar energetic particle events the heavy ion abundance ratios and average charge states show characteristics at high energies that are typically associated with flare-accelerated particles, i.e. impulsive events. The role of flare-accelerated particles in these and other gradual events has been discussed a lot in the scientific community, and it has been questioned if and how the observed features can be explained in terms of diffusive shock acceleration at shock waves driven by coronal mass ejections. The most extreme solar energetic particle events are the so-called ground level enhancements where particle receive so high energies that they can penetrate all the way through Earth's atmosphere and increase radiation levels at the surface. It is not known what conditions are required for acceleration into GeV/nuc energies, and the presence of both very fast coronal mass ejections and X-class solar flares makes it difficult to determine what is the role of these two accelerators in ground level enhancements. The theory of diffusive shock acceleration is reviewed and its predictions discussed with respect to the observed particle characteristics. We discuss how shock waves can be modeled and describe in detail the numerical model developed by the author. The main part of this thesis consists of the four scientific articles that are based on results of the numerical shock acceleration model developed by the author. The novel feature of this model is that it can handle complex magnetic geometries which are found, for example, near active regions in the solar corona. We show that, according to our simulations, diffusive shock acceleration can explain the observed variations in abundance ratios and average charge states, provided that suitable seed particles and magnetic geometry are available for the acceleration process in the solar corona. We also derive an injection threshold for diffusive shock acceleration that agrees with our simulation results very well, and which is valid under weakly turbulent conditions. Finally, we show that diffusive shock acceleration can produce GeV/nuc energies under suitable coronal conditions, which include the presence of energetic seed particles, a favourable magnetic geometry, and an enhanced level of ambient turbulence.