On the motion and radiation of charged particles in strong electromagnetic waves. 1 - Motion in plane and spherical waves
Motion and radiation of charged particles in strong electromagnetic waves in plane and spherical waves
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Motion and radiation of charged particles in strong electromagnetic waves in plane and spherical waves
An exact solution of the Dirac equation in the presence of an arbitrary electromagnetic plane wave is found, which corresponds to a focused electron wave packet, with the focus of the wave packet moving at the speed of light in the opposite direction of the average momentum of the electron wave packet (unless the plane wave is so intense to reflect the electron). The photon spectrum emitted by such an electron wave packet in the presence of a linearly polarized plane wave is studied both analytically and numerically. The spectrum is also compared with the one emitted by a single-momentum, plane-wave electron in the case of the electron being initially counterpropagating (on average for the flying-focus case) with the plane wave and within the locally constant field approximation. It is found that if the electron flying-focus wave packet is focused beyond a Compton wavelength, the angular distribution of the emitted radiation along the magnetic field of the electromagnetic plane wave is broader than for an electron with definite momentum. Corresponding the maximum value of the photon yield on the transverse plane is smaller in the flying-focus electron case. This could represent an experimental signature of a laser-driven flying-focus electron wave packet.
A first principles theory is developed to study the nonlinear spatial evolution of a near-resonance triad of instability waves in boundary layer transition. This triad consists of a plane wave at fundamental frequency and a pair of symmetrical, oblique waves at the subharmonic frequency. A low frequency, high Reynolds number asymptotic scaling leads to a distinct critical layer where nonlinearity first becomes important; the development of the triad's waves is determined by the critical layer's nonlinear, viscous dynamics. The resulting theory is fully nonlinear in that all nonlinearly generated oscillatory and nonoscillatory components are accounted for. The presence of the plane wave initially causes exponential of exponential growth of the oblique waves. However, the plane wave continues to follow the linear theory, even when the oblique waves' amplitude attains the same order of magnitude as that of the plane wave. A fully interactive stage then comes into effect when the oblique waves exceed a certain level compared to that of the plane wave. The oblique waves react back on the fundamental, slowing its growth rate. The oblique waves' saturation results from their self-interaction - a mechanism that does not require the presence of the plane wave. The oblique waves' saturation level is independent of their initial level, but decreases as the obliqueness angle increases.
Flows and rotation, particularly E×B rotation, are critical to improving plasma performance, and waves are a primary tool of plasma control. Thus, it is paramount to understand under what conditions waves can drive E×B flows in plasmas. In this didactic review, an invited paper accompanying the 2023 Marshall N. Rosenbluth Doctoral Thesis Award, this question is answered in the context of momentum-conserving quasilinear theory. There are two primary frameworks for momentum-conserving quasilinear theories that can handle both resonant and nonresonant particles: Eulerian averaging theories and oscillation-center Hamiltonian theories. There are also two different paradigmatic wave problems: plane-wave initial value problems, and steady-state boundary value problems. Here, it is shown that each of these frameworks “naturally” works better with a different problem type. By using these theories, one finds a great difference in the behavior of time- vs space-dependent waves. A time-evolving plane wave can only drive flow if the electromagnetic momentum of the wave, given by the Poynting flux, changes. This result precludes flow drive by any planar electrostatic wave. In contrast, a steady-state spatially evolving wave can drive flow whenever there is divergence in the flux of Minkowski momentum, a completely different physical quantity. This review aims to provide a high-level, intuitive understanding of the very different behaviors observed for these two types of problem.
Plane elastic waves of finite amplitude in Hadamard materials and harmonic materials
Charge particle motion and radiation in strong plane and spherical electromagnetic waves with nonthermal astrophysical applications
Differential equation for passage of plane electromagnetic waves through plane-parallel inhomogeneous plasma slab
The motion of a charged particle in electromagnetic fields of various geometric configurations and arising from a variety of sources is of intrinsic interest in electromagnetic theory. The particular configuration consisting of a plane wave propagating in the presence of a static uniform magnetic field whose direction is parallel to the wave normal is examined. The analysis presented here is treated within the context of classical electromagnetic theory. A numerical solution - at least to the approximate Lorentz-Dirac equation - is obtained.
The interaction of a large‐amplitude disturbing wave with a thermally generated, highly ionized cesium plasma has been investigated using a small‐amplitude sensing wave as a probe. The subsequent interaction between the disturbing and sensing waves is related to the plasma electron temperature. The plasma was contained in a magnetic field of 500–1000 G and had an electron plasma frequency near 10 GHz and an electron gyromagnetic frequency around 2 GHz. The sensing wave was a 10 GHz plane wave of approximately 1.5 mW input power and the disturbing wave was 7.5 GHz plane wave of about 1.2 W peak power. Measurements of the sensing‐wave transmitted power were made using standard microwave techniques. A notch filter in the sensing‐wave receiving circuit reduced the scattered disturbing‐wave interference to an immeasurably low level. Langmuir probes were used to detect disturbing‐wave modulation and to measure plasma temperatures and densities. Measurements of the interaction were correlated with a solution of the Boltzmann transport equation for a perturbed Maxwellian distribution function with binary Coulomb collisions.
A modelling technique is described which is used to explore three dimensional (3D) image irradiance distributions formed by high numerical aperture (NA is greater than 0.5) lenses in homogeneous, linear films. This work uses a 3D modelling approach that is based on a plane-wave decomposition in the exit pupil. Each plane wave component is weighted by factors due to polarization, aberration, and input amplitude and phase terms. This is combined with a modified thin-film matrix technique to derive the total field amplitude at each point in a film by a coherent vector sum over all plane waves. Then the total irradiance is calculated. The model is used to show how asymmetries present in the polarized image change with the influence of a thin film through varying degrees of focus.
Solutions for diffraction of plane shock wave by two dimensional disturbances
Experimental measurements of the attenuation of plane shock waves moving over rough walls have been made in a shock tube. Measurements of the boundary-layer characteristics, including thickness and velocity distribution behind the shock, have also been made with the aid of new cal techniques which provide direct information on the local boundary-layer conditions at the rough walls. Measurements of shock speed and shock pressure ratio are presented for both smooth-wall and rough-wall flow over lengths of machined-smooth and rough strips which lined all four walls of the shock tube. A simplified theory based on Von Karman's expression for skin-friction coefficient for flow over rough walls, along with a wave-model concept and extensions to include time effects, is presented. In this theory, the shock-tube flow is assumed to be one-dimensional at all times and the wave-model concept is used to relate the local layer growth to decreases in shock strength. This concept assumes that local boundary-layer growths act as local mass-flow sinks, which give rise to expansion waves which, in turn, overtake the shock and lower its mass flow accordingly.
Obliquely incident p-polarized plane electromagnetic wave interaction with hot plasma half space, using linearized relativistic Vlasov equation and Laplace transform technique
Array deconvolution is commonly used in aeroacoustic analysis to remove the influence of a microphone array's point spread function from a conventional beamforming map. Unfortunately, the majority of deconvolution algorithms assume that the acoustic sources in a measurement are incoherent, which can be problematic for some aeroacoustic phenomena with coherent, spatially-distributed characteristics. While several algorithms have been proposed to handle coherent sources, some are computationally intractable for many problems while others require restrictive assumptions about the source field. Newer generalized inverse techniques hold promise, but are still under investigation for general use. An alternate coherent deconvolution method is proposed based on a wavespace transformation of the array data. Wavespace analysis offers advantages over curved-wave array processing, such as providing an explicit shift-invariance in the convolution of the array sampling function with the acoustic wave field. However, usage of the wavespace transformation assumes the acoustic wave field is accurately approximated as a superposition of plane wave fields, regardless of true wavefront curvature. The wavespace technique leverages Fourier transforms to quickly evaluate a shift-invariant convolution. The method is derived for and applied to ideal incoherent and coherent plane wave fields to demonstrate its ability to determine magnitude and relative phase of multiple coherent sources. Multi-scale processing is explored as a means of accelerating solution convergence. A case with a spherical wave front is evaluated. Finally, a trailing edge noise experiment case is considered. Results show the method successfully deconvolves incoherent, partially-coherent, and coherent plane wave fields to a degree necessary for quantitative evaluation. Curved wave front cases warrant further investigation. A potential extension to nearfield beamforming is proposed.
Emission of resonance fluorescence by an atom near the surface of a four-wave mixing phase conjugator is considered. The dipole radiation field, regarded as a Heisenberg-operator field, is decomposed into plane waves with the aid of Weyl's representation of the Green's function for the wave equation. Each plane-wave component which is incident on the surface of the nonlinear medium, is reflected as its phase-conjugate image. Summation of all reflected plane waves then yields the phase conjugate replica of the incident dipole radiation. This field adds to the radiation which is emitted by the atom into the direction away from the medium. The condition under which squeezing occurs in the emitted resonance fluorescence is investigated.
A calibrator, referred to as the spider design, can be used to calibrate probes incorporating multiple acoustic sensing elements. The application is an acoustic energy density probe, although the calibrator can be used for other types of acoustic probes. The calibrator relies on the use of acoustic waveguide technology to produce the same acoustic field at each of the sensing elements. As a result, the sensing elements can be separated from each other, but still calibrated through use of the acoustic waveguides. Standard calibration techniques involve placement of an individual microphone into a small cavity with a known, uniform pressure to perform the calibration. If a cavity is manufactured with sufficient size to insert the energy density probe, it has been found that a uniform pressure field can only be created at very low frequencies, due to the size of the probe. The size of the energy density probe prevents one from having the same pressure at each microphone in a cavity, due to the wave effects. The "spider" design probe is effective in calibrating multiple microphones separated from each other. The spider design ensures that the same wave effects exist for each microphone, each with an indivdual sound path. The calibrator s speaker is mounted at one end of a 14-cm-long and 4.1-cm diameter small plane-wave tube. This length was chosen so that the first evanescent cross mode of the plane-wave tube would be attenuated by about 90 dB, thus leaving just the plane wave at the termination plane of the tube. The tube terminates with a small, acrylic plate with five holes placed symmetrically about the axis of the speaker. Four ports are included for the four microphones on the probe. The fifth port is included for the pre-calibrated reference microphone. The ports in the acrylic plate are in turn connected to the probe sensing elements via flexible PVC tubes. These five tubes are the same length, so the acoustic wave effects are the same in each tube. The flexible nature of the tubes allows them to be positioned so that each tube terminates at one of the microphones of the energy density probe, which is mounted in the acrylic structure, or the calibrated reference microphone. Tests performed verify that the pressure did not vary due to bends in the tubes. The results of these tests indicate that the average sound pressure level in the tubes varied by only 0.03 dB as the tubes were bent to various angles. The current calibrator design is effective up to a frequency of approximately 4.5 kHz. This upper design frequency is largely due to the diameter of the plane-wave tubes.
A representation was described whose basis levels combine the important physical aspects of a finite set of plane waves with those of a set of Bloch tight-binding levels. The chosen combination has a particularly simple dependence on the wave vector within the Brillouin Zone, and its use in reducing the standard one-electron band structure problem to the usual secular equation has the advantage that the lattice sums involved in the calculation of the matrix elements are actually independent of the wave vector. For systems with complicated crystal structures, for which the Korringa-Kohn-Rostoker (KKR), Augmented-Plane Wave (APW) and Orthogonalized-Plane Wave (OPW) methods are difficult to apply, the present method leads to results with satisfactory accuracy and convergence.
We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.