Hybrid simulations of collapse of Alfvenic wave packets
Nonlinear dynamics, of large-amplitude circularly polarized Alfvenic wave packets, is investigated using a one-dimensional hybrid model.
SEARCH · Engineering Papers
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.
Nonlinear dynamics, of large-amplitude circularly polarized Alfvenic wave packets, is investigated using a one-dimensional hybrid model.
The Fano line shape, arising from the interference of pathways for the excitation of discrete and continuum states, plays a fundamental role in many branches of physics, chemistry, and materials science. Exciting the resonance with a high harmonic provides naturally a phase delay between the pathways leading to a complex asymmetry parameter. We demonstrate that its amplitude and phase can be controlled on the femtosecond and attosecond time scales, respectively. With our high-energy-resolution (10-meV) experiment, we dynamically image a resonance-enhanced electron wave packet during its temporal evolution, extracting both the amplitude and the phase. Calculations reproduce our experimental results. Our approach constitutes a method for measuring the photoionization delays of a resonance and enables the reconstruction of the electron wave packet in the time domain. This concept of an interference-controlled Fano line shape is a step toward attosecond quantum optics with potential ramifications into nanoscience and next-generation optical materials.
We study the reconstruction of a wave packet and the corresponding electron dynamics in an atom via photoelectron angular distributions (PADs) in a pump-probe scheme as a function of time delay. The method is applied to the superposition of ground and one or two excited states in helium atom representing field-free charge migrations on the attosecond timescale in form of ring currents around the core. It is based on the interference between one- and two-photon transitions from ground and excited states into the continuum. In the reconstruction predictions of first- and second-order perturbation theory are used to determine the unknown phases and amplitudes from the PADs, which we simulate via solutions of the time-dependent Schrödinger equation in single-active-electron approximation. Results of calculations show that the reconstruction technique works well for peak laser intensities less than 10 13 W/cm 2 . Knowledge of the electric field of the probe pulse is required with shot-to-shot variations of carrier-to-envelope phase and peak intensity of up to 10% and 20%, respectively. The relevance of different one- and two-photon pathways for the reconstruction as a function of peak intensity and pulse duration is analyzed—specifically their role for ultrashort probe pulses with broad bandwidths.
The relativistic motion of electrons in an intense electromagnetic wave packet propagating obliquely to a uniform magnetic field is analytically studied on the basis of the Fokker-Planck-Kolmogorov (FPK) approach. The wavepacket consists of circularly polarized electron-cyclotron waves. The dynamical system in question is shown to be reducible to one with three degrees of freedom. Within the framework of the Hamiltonian analysis the nonlinear diffusion tensor is derived, and it is shown that this tensor can be separated into zeroth-, first-, and second-order parts with respect to the relative bandwidth. The zeroth-order part describes diffusive acceleration along lines of constant unperturbed Hamiltonian. The second-order part, which corresponds to the longest time scale, describes diffusion across those lines. A possible transport theory is outlined on the basis of this separation of the time scales.
The dispersion equation of Barnes (1966) is used to study the dissipation of asymptotic wave packets generated by localized periodic sources. The solutions of the equation are linear waves, damped by Landau and transit-time processes, in a collisionless warm plasma. For the case of an ideal MHD system, most of the waves emitted from a source are shown to cancel asympotically through destructive interference. The modes transporting significant flux to asymptotic distances are found to be Alfven waves and fast waves with theta (the angle between the magnetic field and the characteristics of the far-field waves) of about 0 and about pi/2.
The creation of structured electronic wave packets (EWPs) energetically close to Fano resonances has been achieved with ultrafast extreme ultraviolet coherent light sources. However, direct real-time observations of EWP evolution and full reconstructions of the quantum properties of EWPs, including both amplitude and phase, are lacking. Here we introduce and demonstrate a comprehensive approach for the direct measurement and complete characterization of structured EWPs created within a prototypical Fano resonance. Because of its analogy with frequency-resolved optical gating (FROG), we named the method photoelectron FROG. The correlated EWP is initiated by a carefully engineered extreme UV pump pulse. A weak near-infrared laser field, serving as a probe pulse, samples the evolution of the EWPs in the time domain, as well as in the frequency domain. The amplitude and phase of the EWPs are obtained via a time-dependent reconstruction algorithm based on a short-time Fourier transformation. Given the excellent agreement between our experimental results and time-dependent reconstructions, we expect this method to be broadly applicable to the study of ultrafast processes, especially electronic ones, in complex systems, as well as the coherent control of such systems on their fundamental timescales.
Three time-dependent wave packet methods for performing elastic scattering calculations from screened Coulomb potentials are compared. The three methods are the time-dependent amplitude density method (TDADM), a Cayley-transform method (CTM), and the Chebyshev propagation method of Tal-Ezer and Kosloff. Both the TDADM and the CTM are based on a time-dependent integral equation for the wave function. In the first, the time-dependent amplitude density is propagated, while in the other two, the wave function is propagated. As a numerical example, phase shifts and cross sections are calculated using a screened Coulomb, Yukawa type potential over the range 200-1000 eV. It is found that, in most cases, all three methods yield comparable accuracy and are about equally efficient computationally. However for l = 0, where the Coulomb well is not screened by the centrifugal potential, the TDADM requires smaller grid spacings to maintain accuracy.
Particle-in-cell (PIC) simulations and Zakharov's partial differential equations (PDEs) are used to investigate the collapse of isolated Langmuir wave packets in two dimensions. Collapse thresholds are determined numerically, and the roles of enhanced Langmuir wave damping and nonlinearities not included in the standard Zakharov equations are discussed. The Langmuir wave and ion dynamics in PIC simulations are compared with the predictions of PDE simulations incorporating enhanced Langmuir damping. Electron heating and coherent acceleration in the PIC simulations are discussed and compared with predictions of the transit-time theory.
We use the Fokker-Planck-Kolmogorov equation to study the statistical aspects of stochastic dynamics of the radiation belt (RB) protons driven by nonlinear electromagnetic ion-cyclotron (EMIC) wave packets. We obtain the spectra of keV protons scattered by these waves that showsteeping near the gyroresonance, the signature of resonant wave-particle interaction that cannot be described by a simple power law. The most likely mechanism for proton precipitation events in RBs is shown to be nonlinear wave-particle interaction, namely, the scattering of RB protons into the loss cone by EMIC waves.
Photoelectron interferometry with femtosecond and attosecond light pulses is a powerful probe of the fast electron wave-packet dynamics, albeit it has practical limitations on the energy resolution. We show that one can simultaneously obtain both high temporal and spectral resolution by stimulating Raman interferences with one light pulse and monitoring the modification of the electron yield in a separate step. Applying this spectroscopic approach to the autoionizing states of argon, we experimentally resolved its electronic composition and time evolution in exquisite detail. Theoretical calculations show remarkable agreement with the observations and shed light on the light-matter interaction parameters. Using appropriate Raman probing and delayed detection steps, this technique enables highly sensitive probing and control of electron dynamics in complex systems.
Antiferromagnets are promising platforms for transduction and transmission of quantum information via magnons—the quanta of spin waves—and they offer advantages over ferromagnets in regard to dissipation, speed of response and robustness to external fields. Recently, transduction was shown in a van der Waals antiferromagnet, where strong spin-exciton coupling enables readout of the amplitude and phase of coherent magnons by photons of visible light. This discovery shifts the focus of research to transmission, specifically to exploring the non-local interactions that enable magnon wave packets to propagate. Here we demonstrate that magnon propagation is mediated by long-range dipole–dipole interaction. This coupling is an inevitable consequence of fundamental electrodynamics and, as such, will likely mediate the propagation of spin at long wavelengths in the entire class of van der Waals magnets currently under investigation. In conclusion, successfully identifying the mechanism of spin propagation provides a set of optimization rules, as well as caveats, that are essential for any future applications of these promising systems.
A numerically exact spectral method for solving the time-dependent Schroedinger equation in spherical coordinates is described. The angular dependence of the wave function is represented on a two-dimensional grid of evenly spaced points. The fast Fourier transform algorithm is used to transform between the angle space representation of the wave function and its conjugate representation in momentum space. The time propagation of the wave function is evaluated using an expansion of the time evolution operator as a series of Chebyshev polynomials. Calculations performed for a model system representing H2 scattering from a rectangular corrugated surface yield transition probabilities that are in excellent agreement with those obtained using the close-coupling wave packet (CCWP) method. However, the new method is found to require substantially more computation time than the CCWP method because of the large number of grid points needed to represent the angular dependence of the wave function and the variation in the number of terms required in the Chebyshev representation of the time evolution operator.
The method of multiple scales is used to analyze the nonlinear effects of the gas motion and the acoustic lining material on the propagation and attenuation of a wave packet in a two-dimensional duct of uniform cross section. The partial differential equations describing the space and time variation of amplitudes and phases are obtained and used to show that both the monochromatic waves and the pure amplitude modulated waves are stable. The spatial attenuation of the pure amplitude modulated waves is found to be lower than that of monochromatic waves, while the temporal attenuation of the former waves has a minimum value near the resonant frequency. The nonlinearity shifts the wavenumber and frequency to higher values without changing phase speed for the pure phase modulated waves.
We present a method to extract the phase shift of a scattering process using the real-time evolution in the early and intermediate stages of the collision in order to estimate the time delay of a wave packet. This procedure is convenient when using noisy quantum computers for which the asymptotic out-state behavior is unreachable. We demonstrate that the challenging Fourier transforms involved in the state preparation and measurements can be implemented in 1+1 dimensions with current trapped ion devices and IBM quantum computers. We compare quantum computations of the time delays obtained in the one-particle quantum mechanics limit and the scalable quantum field theory formulation with accurate numerical results. We discuss the finite volume effects in the Wigner formula connecting time delays to phase shifts. The results reported involve two- and four-qubit calculations, and we discuss the possibility of larger scale computations in the near future.
In this work, we experimentally study the angle-dependent single ionization of carbon dioxide (CO 2 ) by linearly and circularly polarized pulses. The angle dependence of the ionization probability by linearly polarized pulses extracted from time-domain measurements on an impulsively excited rotational wave packet is compared with data obtained from a direct angle-scan measurement. The results from the measurement with linear and circular polarization are consistent with the adiabatic ionization approximation. Further, we extend the time-domain method to extract the dependence of the asymptotic momentum distribution of fragment ions on the orientation of the molecular axis, and apply it to investigate dissociative double ionization of CO 2 . We show that such measurements can directly test the validity of the axial recoil approximation.
A practical difficulty of traditional Wentzel-Kramer-Brillouin (WKB) theory is that it suffers from non-physical infinities at caustics. A theoretical difficulty is that it does not possess any reasonable covariance or invariance properties under transformations in phase space. It turns out that the solution of one of these problems also solves the other, and leads to a version of WKB theory in phase space. The new theory is most easily implemented by using wave packets.
Extreme ultraviolet (XUV) transient absorption spectroscopy has emerged as a sensitive tool for mapping the real-time structural and electronic evolution of molecules. Here, attosecond XUV transient absorption is used to track dynamics in the A-band of methyl iodide (CH 3 I). Gaseous CH 3 I molecules are excited to the A-band by a UV pump (277 nm, ~20 fs) and probed by attosecond XUV pulses targeting iodine I(4d) core-to-valence transitions. Owing to the excellent temporal resolution of the technique, passage through a conical intersection is mapped through spectral signatures of nonadiabatic wave packet bifurcation observed to occur at 15 ± 4 fs following UV photoexcitation. The observed XUV signatures and time dynamics are in agreement with previous simulations [H. Wang, M. Odelius, and D. Prendergast, J. Chem. Phys. 151, 124106 (2019)]. Furthermore, due to the short duration of the UV pump pulse, coherent vibrational motion in the CH 3 I ground state along the C-I stretch mode (538 ± 7 cm -1 ) launched by resonant impulsive stimulated Raman scattering and dynamics in multiphoton excited states of CH 3 I are also detected.
In this work, a model for multimode perturbations subject to the Richtmyer–Meshkov (RM) instability is presented and compared with simulations and experiments for conditions relevant to inertial confinement fusion. The model utilizes the single mode response to the RM impulse whereby its amplitude h(k, t) first grows with an initial velocity V 0 ∝ kh(k, 0) that eventually decays in time as 1/kV 0 t. Both the growth and saturation stages are subject to nonlinearities since they depend explicitly on the initial amplitude. However, rather than using the individual mode amplitude h(k, t), nonlinearity is taken to occur when the root-mean-square amplitude h rms (k, t) of a wave-packet within wavenumbers k ± δk becomes comparable to 1/k. This is done because nearby sidebands can act in unison for an auto-correlation distance 1/δk beyond nonlinearity as observed in the beam-plasma instability. Thus, the nonlinear saturation amplitude for each mode is reduced from the usual 1/k by a phase space factor that depends on the physical dimensionality, as in the Haan model for the Rayleigh–Taylor instability. In addition, for RM, the average value of kh rms for the initial spectrum is used to calculate a nonlinear factor F NL that reduces V 0 , as observed for single modes. For broadband perturbations, the model describes self-similar growth ∝t θ as successively longer wavelength modes reach saturation. The growing and saturated modes must be discerned because only the former promote θ and are enhanced by reshock and spherical convergence. All of these flows are described here by the model in good agreement with simulations and experiments.