Search NASA⌕ Search

SEARCH · Search NASA

Results for “density wave oscillation”

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.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Role of Dislocations in the Anelasticity of the Upper Mantle

Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work.

Hein, Diede [Univ. of Minnesota, Minneapolis, MN (↗
DOE OSTI2021

Characterization of hollow cathode plasma turbulence using coherent Thomson scattering

Turbulence in hollow cathodes used for space propulsion is believed to play an important role in anomalous electron transport and ion heating. In this work, the implementation of coherent Thomson scattering to identify and characterize MHz-frequency ion acoustic turbulence and kHz-frequency oscillations in the plume of a hollow cathode is achieved. In the presence of a background magnetic field of a Hall thruster, a number of unstable modes are observed. A directive ion acoustic mode propagating predominantly within a restricted angle around the magnetic field is found, exhibiting an energy scaling with wavenumber k of the form k -5.2±0.58 , which differs from the classic Kadomtsev k -3 scaling for unmagnetized conditions. Bi-directional ion acoustic mode fluctuations propagating over a range of angles with respect to the magnetic field have been measured, possibly signifying the existence of a large-amplitude plasma wave, similar to the Buneman instability. Finally, electron density fluctuations in the kHz-frequency range, a possible consequence of drift-driven instabilities in the plane perpendicular to the magnetic field, have also been identified. Furthermore, these results not only are an indication of the diversity of wave types that exist in hollow cathode plumes but also point to the key role played by the presence of, and the configuration of, the magnetic field in their appearance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2024

Accelerated Steam Methane Reforming by Dynamically Applied Charges

Catalyst design has traditionally focused on tuning active site properties to optimally bind reaction intermediates and balance the kinetic requirements of multiple competing chemical processes, as necessitated by the Sabatier principle. It has recently been proposed that for reactions following certain potential energy landscapes, the activity limit imposed by the Sabatier principle may be overcome by using programmed oscillations of surface electron density at the timescales of surface reactions (i.e., “catalytic resonance”). Here, we use a combination of density functional theory (DFT) simulations and transient kinetic models (TKMs) to simulate the kinetics of steam methane reforming (SMR) on Ru(211) surfaces under statically and dynamically applied charges. DFT-calculated binding energies of SMR intermediates and transition states exhibit strong sensitivity to positively applied charges and follow unique scaling relationships that deviate from linear periodic trends across transition metals. Our simulations demonstrate that applying a small positive charge to Ru dramatically enhances the steady-state turnover frequency (TOF) of SMR by up to 5 orders of magnitude above the TOF observed over neutral Ru. Thus, statically charging Ru catalysts may be an effective strategy to lower the temperature requirements for SMR. Dynamic square-wave oscillations in charge resulted in SMR catalytic resonance with an onset frequency f ∼ 106 Hz and the corresponding average TOFs exceeding the statically charged Ru surface by an additional 15%. Here, based on sensitivity analyses performed for the two end points of oscillation, we propose that dynamic TOF improvement beyond the Sabatier maximum can be expected when the system oscillates between two kinetic regimes that are uniquely controlled by distinct elementary steps.

Catalysts↗
DOE OSTI2021

System-on-Chip approach microwave diagnostics development for burning plasma

System-on-chip (SoC) millimeter wave integrated circuit technology has been developed and employed in chip heterodyne radiometers in Electron Cyclotron Emission Imaging (ECEI) and Microwave Imaging Reflectometer (MIR) on the DIII-D tokamak for 2D electron temperature and density fluctuation diagnostics. The W-band (75-110 GHz) ECEI system has been employed with 20 poloidal horn receiver modules with customized monolithic integrated circuit chips and internal local oscillator coupling path. The new SoC approach demonstrates > 30 dB additional gain on the receiver and 95% electronics noise suppression compared to the previous Schottky diode mixer array . An F-band (110-140 GHz) receiver chip has been developed to measure MHD and turbulence behavior on the DIII-D high-field side, such as Alfven eigenmodes. Furthermore, it can image pedestal structure evolution on fusion facilities with stronger magnetic fields, such as ASDEX-U, EAST, HL-2M, and LHD. The F-band chip performance has been tested in the laboratory. Also, V-band (55-75 GHz) transmitter and receiver modules (55-75 GHz) have been developed for DIII-D and NSTX-U Furthermore, microwave Imaging Reflectometer (MIR) application suitability for 2D density fluctuation imaging on the pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Experiments and Simulations on the Turbulent, Rarefaction Wave Driven Rayleigh–Taylor Instability

Experiments were performed to observe the growth of the turbulent, Rayleigh–Taylor unstable mixing layer generated between air and SF 6 , with an Atwood number of A=(ρ 2 –ρ 1 )/(ρ 2 +ρ 1 ) = 0.64, where ρ 1 and ρ 2 are the densities of air and SF 6 , respectively. A nonconstant acceleration with an average value of 2300g 0 , where g 0 is the acceleration due to gravity, was generated by interaction of the interface between the two gases with a rarefaction wave. Three-dimensional, multimode perturbations were generated on the diffuse interface, with a diffusion layer thickness of δ = 3.6 mm, using a membraneless vertical oscillation technique, and 20 experiments were performed to establish a statistical ensemble. The average perturbation from this ensemble was extracted and used as input for a numerical simulation using the Lawrence Livermore National Laboratory (LLNL) Miranda code. Good qualitative agreement between the experiment and simulation was observed, while quantitative agreement was best at early to intermediate times. Several methods were used to extract the turbulent growth constant α from experiments and simulations while accounting for time varying acceleration. Experimental, average bubble and spike asymptotic self-similar growth rate values range from α = 0.022 to α = 0.032 depending on the method used, and accounting for variable acceleration. Values found from the simulations range from α = 0.024 to α = 0.041. Here, the values of α measured in the experiments are lower than what are typically measured in the literature but are more in line with those found in recent simulations.

42 ENGINEERING↗
DOE OSTI2021

Dynamical friction from scalar dark matter in the relativistic regime

Light bosonic scalars (e.g., axions) may form clouds around black holes via superradiant instabilities or via accretion if they form some component of the dark matter. It has been suggested that their presence may lead to a distinctive dephasing of the gravitational wave signal when a small compact object spirals into a larger black hole. Motivated by this, we have studied numerically the dynamical friction force on a black hole moving at relativistic velocities in a background scalar field with an asymptotically homogeneous energy density. Here we show that the relativistic scaling is analogous to that found for supersonic collisional fluids, assuming an approximate expression for the pressure correction which depends on the velocity and scalar mass. While we focus on a complex scalar field, our results confirm the expectation that real scalars would exert a force which oscillates between positive and negative values in time with a frequency set by the scalar mass. The complex field describes the time averaged value of this force, but in a real scalar, the rapid force oscillations could, in principle, leave an imprint on the trajectory. The approximation we obtain can be used to inform estimates of dephasing in the final stages of an extreme mass ratio inspiral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Lightwave terahertz quantum manipulation of nonequilibrium superconductor phases and their collective modes

In this work, we present a gauge-invariant density matrix description of nonequilibrium superconductor (SC) states with spatial and temporal correlations driven by intense terahertz (THz) lightwaves. We derive superconductor Bloch-Maxwell equations of motion that extend Anderson pseudospin models to include the Cooper pair center-of-mass motion and electromagnetic propagation effects. We thus describe quantum control of dynamical phases, collective modes, quasiparticle coherence, and high nonlinearities during cycles of carrier wave oscillations, which relates to our recent experiments. Coherent photogeneration of a nonlinear supercurrent with a dc component, achieved via condensate acceleration by an effective lightwave field, dynamically breaks the equilibrium inversion symmetry. Experimental signatures include high harmonic light emission at equilibrium-symmetry-forbidden frequencies, Rabi-Higgs collective modes and quasiparticle coherence, and nonequilibrium moving condensate states tuned by few-cycle THz fields. We use such lightwaves as an oscillating accelerating force that drives strong nonlinearities and anisotropic quasiparticle populations to control and amplify different classes of collective modes, e.g., damped oscillations, persistent oscillations, and overdamped dynamics via Rabi flopping. Recent phase-coherent nonlinear spectroscopy experiments can be modeled by solving the full nonlinear quantum dynamics including self-consistent light-matter coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Isospin equilibration in neutron star mergers

We analyze the isospin equilibration properties of neutrinoless nuclear (npe) matter in the temperature and density range that is relevant to neutron star mergers. Our analysis incorporates neutrino-transparency corrections to the isospin (“beta”) equilibrium condition which become noticeable at T ≳ 1 MeV. We find that the isospin relaxation rate rises rapidly as temperature rises, and at T ≈ 5 MeV it is comparable to the timescale of the density oscillations that occur immediately after the merger. This produces a resonant peak in the bulk viscosity at T ≈ 5 MeV, which causes density oscillations to be damped on the timescale of the merger. Furthermore, our calculations suggest that isospin relaxation dynamics may also be relevant when neutrinos are treated more accurately via neutrino transport schemes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Computing molecular excited states on a D-Wave quantum annealer

Abstract The possibility of using quantum computers for electronic structure calculations has opened up a promising avenue for computational chemistry. Towards this direction, numerous algorithmic advances have been made in the last five years. The potential of quantum annealers, which are the prototypes of adiabatic quantum computers, is yet to be fully explored. In this work, we demonstrate the use of a D-Wave quantum annealer for the calculation of excited electronic states of molecular systems. These simulations play an important role in a number of areas, such as photovoltaics, semiconductor technology and nanoscience. The excited states are treated using two methods, time-dependent Hartree–Fock (TDHF) and time-dependent density-functional theory (TDDFT), both within a commonly used Tamm–Dancoff approximation (TDA). The resulting TDA eigenvalue equations are solved on a D-Wave quantum annealer using the Quantum Annealer Eigensolver (QAE), developed previously. The method is shown to reproduce a typical basis set convergence on the example $$\hbox {H}_2$$ H 2 molecule and is also applied to several other molecular species. Characteristic properties such as transition dipole moments and oscillator strengths are computed as well. Three potential energy profiles for excited states are computed for $$\hbox {NH}_3$$ NH 3 as a function of the molecular geometry. Similar to previous studies, the accuracy of the method is dependent on the accuracy of the intermediate meta-heuristic software called qbsolv.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Measuring the properties of f - mode oscillations of a protoneutron star by third-generation gravitational-wave detectors

Core-collapse supernovae are among the astrophysical sources of gravitational waves that could be detected by third-generation gravitational-wave detectors. Here, we analyze the gravitational-wave strain signals from two- and three-dimensional simulations of core-collapse supernovae generated using the code FORNAX. A subset of the two-dimensional simulations has non-zero core rotation at the core bounce. A dominant source of time changing quadrupole moment is the l = 2 fundamental mode (f- mode) oscillation of the proto-neutron star. From the time-frequency spectrogram of the gravitational-wave strain we see that, starting ~ 400 ms after the core bounce, most of the power lies within a narrow track that represents the frequency evolution of the f-mode oscillations. The f-mode frequencies obtained from linear perturbation analysis of the angle averaged profile of the protoneutron star corroborate what we observe in the spectrograms of the gravitational wave signal. We explore the measurability of the f-mode frequency evolution of protoneutron star for a supernova signal observed in the third-generation gravitational-wave detectors. Measurement of the frequency evolution can reveal information about the masses, radii, and densities of the proto-neutron stars. We find that if the third generation detectors observe a supernova within 10 kpc, we can measure these frequencies to within ~90% accuracy. Here, we can also measure the energy emitted in the fundamental f-mode using the spectrogram data of the strain signal. We find that the energy in the f-mode can be measured to within 20% error for signals observed by Cosmic Explorer using simulations with successful explosion, assuming source distances within 10 kpc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

System-on-chip approach microwave imaging reflectometer on DIII-D tokamak

System-on-chip millimeter wave integrated circuit technology is used on the two-dimensional millimeter-wave imaging reflectometer (MIR) upgrade for density fluctuation imaging on the DIII-D tokamak fusion plasma. Customized CMOS chips have been successfully developed for the transmitter module and receiver module array, covering the 55–75 GHz working band. The transmitter module has the capability of simultaneously launching eight tunable probe frequencies (>0 dBm output power each). The receiver enclosure contains 12 receiver modules in two vertical lines. The quasi-optical local oscillator coupling of previous MIR systems has been replaced with an internal active frequency multiplier chain for improved local oscillator power delivery and flexible installation in a narrow space together with improved shielding against electromagnetic interference. The 55–75 GHz low noise amplifier, used between the receiver antenna and the first-stage mixer, significantly improves module sensitivity and suppresses electronics noise. The receiver module has a 20 dB gain improvement compared with the mini-lens approach and better than –75 dBm sensitivity, and its electronics noise temperature has been reduced from 55 000 K down to 11 200 K. Further, the V-band MIR system is developed for co-located multi-field investigation of MHD-scale fluctuations in the pedestal region with W-band electron cyclotron emission imaging on DIII-D tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory of the Tertiary Instability and the Dimits Shift from Reduced Drift-Wave Models

Tertiary modes in electrostatic drift-wave turbulence are localized near extrema of the zonal velocity U(x) with respect to the radial coordinate x. We argue that these modes can be described as quantum harmonic oscillators with complex frequencies, so their spectrum can be readily calculated. The corresponding growth rate $γ$ TI is derived within the modified Hasegawa-Wakatani model. We show that $γ$TI equals the primary-instability growth rate plus a term that depends on the local U''; hence, the instability threshold is shifted compared to that in homogeneous turbulence. This provides a generic explanation of the well-known yet elusive Dimits shift, which we find explicitly in the Terry-Horton limit. Linearly unstable tertiary modes either saturate due to the evolution of the zonal density or generate radially propagating structures when the shear |U'| is sufficiently weakened by viscosity. The Dimits regime ends when such structures are generated continuously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory of the tertiary instability and the Dimits shift from reduced drift-wave models

Tertiary modes in electrostatic drift-wave turbulence are localized near extrema of the zonal velocity $U(x)$ with respect to the radial coordinate $x$. We argue that these modes can be described as quantum harmonic oscillators with complex frequencies, so their spectrum can be readily calculated. The corresponding growth rate $\gamma_{\rm TI}$ is derived within the modified Hasegawa--Wakatani model. We show that $\gamma_{\rm TI}$ equals the primary-instability growth rate plus a term that depends on the local $U''$; hence, the instability threshold is shifted compared to that in homogeneous turbulence. This provides a generic explanation of the well-known yet elusive Dimits shift, which we find explicitly in the Terry--Horton limit. Linearly unstable tertiary modes either saturate due to the evolution of the zonal density or generate radially propagating structures when the shear $|U'|$ is sufficiently weakened by viscosity. The Dimits regime ends when such structures are generated continuously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A New Plasma Radar Concept for Simultaneous Magnetic and Density Measurements

An innovative, compact 288GHz interferometer has been fabricated, tested, installed and successfully demonstrated on the LAPD-U magnetized plasma at UCLA. The system takes advantage of frequency modulated (FM) radar techniques to deliver a compact heterodyne system. In addition, the reflected power from the source is taken advantage of to eliminate the need for additional quasi-optical components. Electron density in LAPD-U plasma has recently been increased substantially thereby requiring a higher frequency/shorter wavelength interferometer to avoid deleterious refractive effects. This system satisfies those needs. The system uses a 96GHz varactor tuned Gunn oscillator which passes to a passive tripler. This tripler has ~3% conversion efficiency. The 288GHz radiation is then coupled to free space using a so-called dual-mode or Pickett horn. The output 288GHz beam is then coupled to an aspheric lens manufactured from low-loss, high-density polyethylene. This lens is employed to collimate the emerging beam. Small axial adjustment of the lens position can also be used to create a slowly focusing beam so as to optimize the measured signal. In addition, up-down or side-to-side adjustment of the lens can be utilized to steer the beam vertically or horizontally – again to optimize alignment. The propagating beam passes through a beam splitter and then through a water-free, bubble-free fused quartz window into the LAPD-U vacuum vessel. The beam-splitter is a thin sheet of G10 which reflects a small fraction of the incident power (~5 %) towards a zero-bias detector optimized for the frequency range from 220 to 300GHz. Note that waveguides at this frequency have dimensions of ~0.9mm x 0.45mm and so have very large conductive losses. This drives the use of quasi-optical propagation. The detector requires no DC bias and is very responsive (> 1V/mW into 1MΩ). Radiation is coupled to the detector via a similar lens-horn arrangement used for the launch. This reflected beam acts as the local oscillator or reference millimeter-wave beam for the detector. The remainder of the launched source beam then enters the LAPD-U vacuum vessel and passes through the plasma at the mid-plane until reaching the opposing port which is closed off with an aluminum flange. This flange is used as a mirror to retroreflect the incident 288GHz beam back along its path. The retroreflected beam exits the input port but does NOT couple directly into the zero-bias detector. Instead, the majority of the return power continues towards the 288GHz source. As mentioned above the transmitted beam enters the source a second time. This would appear undesirable. However, at these frequencies multipliers are highly non-linear elements which results in a significant portion of the return beam (~20%) re-emerging from the multiplier and horn and then coupling via the G10 beam-splitter to the zero-bias detector. This approach eliminated the need for a second quasi-optical beam-splitter. The system is extremely compact measuring approximately 28 inches x 20 inches. The above did not explain how heterodyne operation was achieved. As mentioned above the Gunn oscillator is able to be varactor tuned. This allows a low voltage to be applied to control the operating frequency of the Gunn oscillator. During heterodyne operation a sawtooth shaped voltage is applied to the varactor at 750kHz using an 80MHz Arbitrary Waveform Generator (AWG). This voltage changes the Gunn frequency linearly during the up-sweep which is then reset abruptly at the sawtooth crash to be immediately followed by another linear sweep. Passage through the 288GHz multiplier triples the frequency change experienced by the electromagnetic wave. These frequency changes are small – tens of megahertz. This FM radar approach results in the launched electromagnetic wave frequencies at the detector for the reference and plasma wave to be different. The approximately 10 ns delay propagation delay for the plasma beam results in the local oscillator and plasma beams NOT having an identical frequency – there is in fact a fixed difference frequency. The frequency tuning level of the Gunn oscillator is then adjusted so that there is ONE cycle of this difference frequency during each linear ramp. During the sawtooth crash or downward re-sweep this one cycle replays in reverse but on a very fast timescale. The process then repeats. Low-pass filtering eliminates the fast re-sweep to leave a pure sine wave heterodyne signal. When the plasma is present it introduces a phase delay in the sine wave (caused by the extremely small Doppler shift resulting from the optical path length change). Of course, to measure this phase change we need a reference. This is simply obtained from an arbitrary waveform generator which provides a synchronized output pulse train which again is low pass filtered to obtain a 750kHz sinusoidal voltage reference for the interferometer. The interferometer was installed on LAPD-U where it has worked reliably and has established that electron densities exceeding 1x10 13 cm -3 are routinely achieved. In addition, the system sensitivity was able to easily observe density fluctuation at frequencies up to 50kHz. FM Radar techniques have enabled a full demonstration of a compact, sensitive, high frequency (288GHz/1mm) heterodyne interferometer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Nontrivial Fermi surface topology of the kagome superconductor CsV 3 Sb 5 probed by de Haas–van Alphen oscillations

Here, we have investigated the normal state Fermi-surface properties of the kagome superconductor CsV 3 Sb 5 using torque magnetometry with applied fields (H) up to 35 T. The torque signal shows clear de Haas-van Alphen (dHvA) oscillations above 15 T. The oscillations are smooth and consist of seven distinct frequencies with values from similar to ~18 T to 2135 T. The presence of higher frequencies in CsV 3 Sb 5 is further confirmed by carrying out additional measurements using the tunnel diode oscillator technique. All frequencies measured at different tilt angles (θ) of the field direction with respect to the c axis show a 1/cosθ dependence, implying that the Fermi surfaces corresponding to these frequencies are two dimensional (2D). The absence of dHvA oscillations at θ = 90° further supports the presence of 2D Fermi surfaces. The Berry phase (φ) determined from the Landau level fan diagram for all frequencies is similar to ~ 0.4. This value is close to the theoretical value of φ = 0.5 for a nontrivial system, which strongly supports the nontrivial topology of the Fermi surfaces of these frequencies. Several quantities characterizing the Fermi surface are calculated employing the Lifshitz-Kosevich theory. These findings are crucial for exploring the interplay between nontrivial band topology, charge-density wave, and unconventional superconductivity of CsV 3 Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Nonlinear evolution, propagation, electron-trapping, and damping effects of ion-acoustic solitons using fully kinetic PIC simulations

We investigate ion acoustic solitary waves (solitons) of varying amplitudes in a one-dimensional plasma using fully kinetic particle-in-cell simulations. The initial soliton conditions are based on the Korteweg–de Vries (KdV) equation, treating ions as a cold species and electrons with finite temperature. Our findings reveal that KdV solitons evolve nonlinearly to a saturated state at higher amplitude, deviating from KdV predictions for ion density and electric potential, and from the Boltzmann relation for electron density. At this saturated state, the KdV model cannot accurately describe the soliton behavior. For small amplitudes, Sagdeev's model describes the saturated state, but not the soliton width; for larger amplitudes, it models the width accurately, but not the amplitude. These discrepancies arise from assuming a Boltzmann relation for electron density, while electron trapping creates non-Boltzmann densities—a deviation that increases with soliton amplitude. Additionally, we observe that the soliton amplitude oscillates roughly at the electron bounce frequency. The soliton is better described by Schamel's electron density formulation and a modified KdV equation incorporating electron trapping. The soliton velocity matches best with predictions from Sagdeev's and Schamel's models. Moreover, the soliton speed–amplitude relationship differs from existing theoretical predictions. Finally, we find minimal ion and electron Landau damping effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Scattering of electron from a disk in 2D electron gas: Full cross section, transport cross section, and the interaction correction

It is known that the presence of the Fermi sea modifies the scattering of an electron from a point-like impurity. This is due to the Friedel oscillations of the electron density around the impurity. These oscillations create an additional scattering potential for incident electrons. The closer the energy of the incident electron to the Fermi level, the stronger the additional scattering. Here, we study this effect for the case when the impurity is not point-like but rather a hard disk, with a radius much bigger than the de Broglie wavelength. We start with a careful examination of the full and transport cross sections from an extended target. Both cross sections approach their limiting values upon increasing the wave vector of the incident electron. We establish that the transport cross section saturates much faster than the full cross section. With regard to the interaction correction, we establish that it vanishes for the full cross section, while for the transport cross section, it is enhanced compared to the case of a point-like scatterer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗