Search NASA⌕ Search

Engineering topics

Miller, Kyle G.

Publications and source records attributed to Miller, Kyle G..

Dephasingless two-color terahertz generation

A laser pulse composed of a fundamental and an appropriately phased second harmonic can drive a time-dependent current of photoionized electrons that generates broadband THz radiation. Over the propagation distances relevant to many experiments, dispersion causes the relative phase between the harmonics to evolve. This “dephasing” slows the accumulation of THz energy and results in a multi-cycle THz pulse with significant angular dispersion. Here, we introduce a novel optical configuration that compensates the relative phase evolution, allowing for the formation of a half-cycle THz pulse with almost no angular dispersion. The configuration uses the spherical aberration of an axilens to map a prescribed radial phase variation in the near field to a desired longitudinal phase variation in the far field. Simulations that combine this configuration with an ultrashort flying focus demonstrate the formation of a half-cycle THz pulse with a controlled emission angle and 1/4 the angular divergence of the multi-cycle pulse created by a conventional optical configuration.

47 OTHER INSTRUMENTATION↗

Generation of attosecond gigawatt soft x-ray pulses through coherent Thomson backscattering

Collision between relativistic electron sheets and counterpropagating laser pulses is recognized as a promising way to produce intense attosecond x rays through coherent Thomson backscattering (TBS). In a double-layer scheme, the electrons in an ultrathin solid foil are first pushed out by an intense laser driver and then interact with the laser reflected off a second foil to form a high-density relativistic electron sheet with vanishing transverse momentum. However, the repulsion between these concentrated electrons can increase the thickness of the layer, reducing both its density and subsequently the coherent TBS. Here, we present a systematic study on the evolution of the flying electron layer and find that its resulting thickness is determined by the interplay between the intrinsic space-charge expansion and the velocity compression induced by the drive laser. How the laser driver, the target areal density, the reflector, and the collision laser intensity affect the properties of the produced x rays is explored. Multidimensional particle-in-cell simulations indicate that employing this scheme in the nonlinear regime has the potential to stably produce soft x rays with several gigawatt peak power in hundreds of terawatt ultrafast laser facilities. The pulse duration can be tuned to tens of attoseconds. Here, this compact and intense attosecond x-ray source may have broad applications in attosecond science.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dephasingless laser wakefield acceleration in the bubble regime

Abstract Laser wakefield accelerators (LWFAs) have electric fields that are orders of magnitude larger than those of conventional accelerators, promising an attractive, small-scale alternative for next-generation light sources and lepton colliders. The maximum energy gain in a single-stage LWFA is limited by dephasing, which occurs when the trapped particles outrun the accelerating phase of the wakefield. Here, we demonstrate that a single space–time structured laser pulse can be used for ionization injection and electron acceleration over many dephasing lengths in the bubble regime. Simulations of a dephasingless laser wakefield accelerator driven by a 6.2-J laser pulse show 25 pC of injected charge accelerated over 20 dephasing lengths (1.3 cm) to a maximum energy of 2.1 GeV. The space–time structured laser pulse features an ultrashort, programmable-trajectory focus. Accelerating the focus, reducing the focused spot-size variation, and mitigating unwanted self-focusing stabilize the electron acceleration, which improves beam quality and leads to projected energy gains of 125 GeV in a single, sub-meter stage driven by a 500-J pulse.

43 PARTICLE ACCELERATORS↗

Accurate simulation of direct laser acceleration in a laser wakefield accelerator

In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The traditional particle-in-cell (PIC) algorithm, although often the tool of choice to study DLA, contains inherent errors due to numerical dispersion and the time staggering of the electric and magnetic fields. Furthermore, conventional PIC implementations cannot reliably disentangle the fields of the plasma wave and laser pulse, which obscures interpretation of the dominant acceleration mechanism. Here, a customized field solver that reduces errors from both numerical dispersion and time staggering is used in conjunction with a field decomposition into azimuthal modes to perform PIC simulations of DLA in an LWFA. Comparisons with traditional PIC methods, model equations, and experimental data show improved accuracy with the customized solver and convergence with an order-of-magnitude fewer cells. Furthermore, the azimuthal-mode decomposition reveals that the most energetic electrons receive comparable energy from DLA and LWFA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extended particle absorber for efficient modeling of intense laser-solid interactions

Here, an extended thermal particle boundary condition is devised to more efficiently and accurately model laser–plasma interactions in overdense plasmas. Particle-in-cell simulations of such interactions require many particles per cell, and a large region of background plasma is often necessary to correctly mimic a semi-infinite plasma and avoid electron refluxing from a truncated plasma. For long-pulse lasers of many picoseconds, such constraints can become prohibitively expensive. Here, an extended particle boundary condition (absorber) is designed that instantaneously stops and re-emits energetic particles streaming toward the simulation boundary over a defined region, allowing sufficient time and space for a suitably cool return current to develop in the background plasma. Tunable parameters of the absorber are explained, and simulations using the absorber with a 3-ps laser are shown to accurately reproduce those of a causally separated boundary while requiring only 20% the number of particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dynamic load balancing with enhanced shared-memory parallelism for particle-in-cell codes

Furthering our understanding of many of today’s interesting problems in plasma physics – including plasma based acceleration and magnetic reconnection with pair production due to quantum electrodynamic effects – requires large-scale kinetic simulations using particle-in-cell (PIC) codes. However, these simulations are extremely demanding, requiring that contemporary PIC codes be designed to efficiently use a new fleet of exascale computing architectures. To this end, the key issue of parallel load balance across computational nodes must be addressed. We discuss the implementation of dynamic load balancing by dividing the simulation space into many small, self-contained regions or ‘‘tiles,’’ along with shared-memory (e.g., OpenMP) parallelism both over many tiles and within single tiles. The load balancing algorithm can be used with three different topologies, including two space-filling curves. Here, we tested this implementation in the code Osiris and show low overhead and improved scalability with OpenMP thread number on simulations with both uniform load and severe load imbalance. Compared to other load-balancing techniques, our algorithm gives order-of-magnitude improvement in parallel scalability for simulations with severe load imbalance issues.

97 MATHEMATICS AND COMPUTING↗

A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm

A customized finite-difference field solver for the particle-in-cell (PIC) algorithm that provides higher fidelity for wave-particle interactions in intense electromagnetic waves is presented. In many problems of interest, particles with relativistic energies interact with intense electromagnetic fields that have phase velocities near the speed of light. Numerical errors can arise due to (1) dispersion errors in the phase velocity of the wave, (2) the staggering in time between the electric and magnetic fields and between particle velocity and position and (3) errors in the time derivative in the momentum advance. Errors of the first two kinds are analyzed in detail. It is shown that by using field solvers with different -space operators in Faraday’s and Ampere’s law, the dispersion errors and magnetic field time-staggering errors in the particle pusher can be simultaneously removed for electromagnetic waves moving primarily in a specific direction. Here, the new algorithm was implemented into Osiris by using customized higher-order finite-difference operators. Schemes using the proposed solver in combination with different particle pushers are compared through PIC simulation. It is shown that the use of the new algorithm, together with an analytic particle pusher (assuming constant fields over a time step), can lead to accurate modeling of the motion of a single electron in an intense laser field with normalized vector potentials, eA / mc 2 , exceeding for typical cell sizes and time steps.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Analyzing data from a PDV microphone

Photon Doppler Velocimetry (PDV) is used at Los Alamos National Laboratory to measure extremely fast velocities (on the order of 0.01 - 5 km/s) of projectiles during material testing procedures. This kind of testing is extremely important to understand how materials behave in extreme conditions (such as in nuclear explosions). PDV data is created by summing an incident beam with a Doppler shifted reflection, forming a beat frequency. The beat frequency is directly proportional to the velocity of the reflecting surface. We explored a new way of using PDV: measuring slow velocities (less than 1 m/s). The principle method we used was to measure vibrations in a Mylar sheet due to incident sound waves. I analyzed these files using both my own coding as well as a program designed at NSTec for this purpose. Analysis included performing Fourier Transforms, extracting velocity as a function of time, and playing the velocity function through a speaker to recreate the initial sound. We were able to create a successful microphone, though playback quality was poor. The limiting factor was the low velocity of the Mylar sheet, which was about 0.2 m/s. Other methods are now being explored to extract velocity data in a more efficient manner. During this project I gained extensive experience in programming, data analysis, problem solving, and communication. I learned that one must be creative in order to solve complex problems, approaching them from various angles. I also gained a lot of practice working on a team, which consisted of a high school student, another college intern, and our mentor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗