Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave function”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Perturbative model for the saturation of energetic-particle-driven modes limited by self-generated zonal modes

We present a simplified energy-conserving approach to incorporate wave–wave nonlinear effects within the framework commonly used to describe wave–particle nonlinearities. In particular, the effects of zonal mode (ZM) generation on the determination of the saturation amplitude of energetic particle (EP)-driven Alfvénic instabilities is studied. The model assumes that the zonal perturbations grow at a rate twice that of the original (pump) wave, consistent with a beat-driven (or force-driven) generation mechanism. The evolution and saturation of the mode amplitude are investigated both analytically and numerically within our reduced model assumptions, in both the collisionless and scattering-dominated regimes. These studies underscore the crucial role of sources and sinks in capturing the impact and the role of beat-driven zonal perturbations on mode evolution. In the realistic case of saturation set by sources and sinks, we discuss the role of a finite amplitude ZM in reducing microturbulent particle scattering, thus limiting the energy source for the resonant mode. We then discuss comparisons between the model’s predictions and simulation results. The model reproduces key features observed in gyrokinetic simulations as the reduction in saturated mode amplitude and the onset of wave–wave nonlinear effects as functions of mode growth rate and amplitude. Thanks to its simplicity, it can be readily implemented into codes based on reduced models, thereby improving their predictive capability for strongly driven instabilities.

Alfvén eigenmodes↗
DOE OSTI2021

Identification of a low-energy metastable 1 T -type phase for monolayer VSe 2

Elucidating the polymorphism of transition metal dichalcogenide layers and the interplay between structure and properties is a key challenge for the application of these materials. We identify a low-energy metastable phase of monolayer VSe 2 and elucidate its magnetic and electronic properties. This structure is distinct from the previously identified charge density wave (CDW) phase. However, while having rather distinct properties it is very close in energy to the CDW phase and is likely to be realized in experiments. Importantly, local bonding instabilities, as characterized by reconstruction of the electronic structure over a wide energy range, are important for this distortion, which includes both V off-centering in the octahedral coordination cages and a partial disproportionation into two distinct types of V. The phase does not have a ferromagnetic ground state. Furthermore, the results show that the physics of 1T–VSe 2 are richer than previously known with an interplay of Fermi surface instabilities and local bonding effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Intrinsic and doping-enhanced superconductivity in monolayer 1 H − TaS 2 : Critical role of charge ordering and spin-orbit coupling

The interplay of superconductivity with charge density wave (CDW) in metallic transition-metal dichalcogenides has been widely debated, and viable strategies manipulating these quantum states in the two-dimensional (2D) limit remain unclear. Using the ab initio anisotropic Migdal-Eliashberg theory, we successfully explain the superconductivity observed in monolayer 1H-TaS 2 by simultaneously determining its precise CDW structure and treating the marked modification of electron-phonon interaction and critical temperature T c by spin-orbit coupling effects. With this paradigm, we further show that electron doping weakens the CDW order leading to increased T c up to 11 K, along with a single-gap to two-gap superconductivity transition due to the suppression of the CDW gap. By contrast, a low hole doping barely affects the CDW but still yields a significantly enhanced superconducting order, implying their good coexistence. Combined with the synergistic behavior of CDW and superconductivity, which cooperate upon TaS 2 thickness reduction causing an unusual rise of T c , our results unravel diversified interactions between the two collective orders in ultrathin TaS 2 , being competition, coexistence or cooperation depending on external stimuli, which provide key clues for controlling correlated states in devices based on 2D CDW superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multiband Fermi surface in $1T–$$\mathrm{VSe_2}$ and its implication for the charge density wave phase

Angle-resolved photoemission spectroscopy experiments reveal a surprisingly richer surface electronic structure in $1T–$$\mathrm{VSe_2}$ than previously predicted or probed. Earlier claims supporting a charge density wave phase in this material are reexamined in terms of these findings and are found to be untenable. Here, the Fermi surface is found to be gapless, while band warping effects, currently attributed to three-dimensional lattice distortion, result from the simultaneous dispersion of the closely lying multiple bands. Based on these findings, a charge density wave scenario in $1T–$$\mathrm{VSe_2}$ is unlikely. On the other side, the presence of multiple states crossing the Fermi level should constitute relevant constraints for any viable microscopic model of the structural phase transition of $\mathrm{VSe_2}$.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Elastic scattering of 3 He + 4 He with the SONIK scattering chamber

Measurements of the elastic scattering cross section of 3 He and 4 He are important in order to improve constraints on theoretical models of 4 He ⁢( 3 He ,𝛾)⁢ 7 Be, a key reaction in Big Bang nucleosynthesis and solar neutrino production. The astrophysical 𝑆 factor for this reaction is a significant source of uncertainty in the standard-solar-model prediction of the 7 Be and 8 B solar neutrino fluxes. The elastic scattering measurements reported in the literature do not extend to low energies and lack proper uncertainty quantification. A new measurement of the 4 He ⁢( 3 He, 3 He)⁢ 4 He reaction has been made at center-of-mass energies 𝐸c.m.=0.38−3.13 MeV using the Scattering of Nuclei in Inverse Kinematics (SONIK) scattering chamber: a windowless, extended gas target surrounded by an array of 30 collimated silicon charged particle detectors situated at TRIUMF. This is the first elastic scattering measurement of 3 He + 4 He made below 500 keV and it has greater angular range and better precision than previous measurements. The elastic scattering data were analyzed using both 𝑅-matrix and halo effective field theory frameworks, and values of the 𝑠-wave scattering length and effective range were extracted. Here, the resulting improvement in knowledge of the 𝑠-wave effective-range function at low energies reduces the overall uncertainty in 𝑆 34 at solar energies.

Elastic scattering reactions↗
DOE OSTI2022

Fermi level tuning and double-dome superconductivity in the kagome metal CsV 3 Sb 5 – x Sn x

The recently reported AV 3 Sb 5 (A = K , Rb, Cs) family of kagome metals are candidates for unconventional superconductivity and chiral charge density wave (CDW) order; both potentially arise from nested saddle points in their band structures close to the Fermi energy. Here, we use chemical substitution to introduce holes into CsV 3 Sb 5 and unveil an unconventional coupling of the CDW and superconducting states. Specifically, we generate a phase diagram for CsV 3 Sb 5 – x Sn x that illustrates the impact of hole doping the system and lifting the nearest van Hove singularity toward and above E F . Superconductivity exhibits a nonmonotonic evolution with the introduction of holes, resulting in two “domes” peaked at 3.6 and 4.1 K and the rapid suppression of three-dimensional CDW order. Further, the evolution of CDW and superconducting order is compared with the evolution of the electronic band structure of CsV 3 Sb 5 – x Sn x , where the complete suppression of superconductivity seemingly coincides with an electronlike band comprised of Sb p z orbitals pushed above E F .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Rupture Model of the 5 April 2024 Tewksbury, New Jersey, Earthquake Based on Regional Lg -Wave Data

On 5 April 2024, an earthquake of magnitude 4.8 occurred in Tewksbury, New Jersey. It was the largest instrumentally recorded event since 1900 in New Jersey and southern New York. Millions of people around New York City, ~65 km east–northeast of Tewksbury, felt the shaking from the mainshock, but the epicentral area experienced no known significant property damages. We determine the focal mechanism, which is oblique faulting, and retrieve the Lg–wave relative source time functions (RSTFs) from the stations at regional distances to understand rupture processes and ground motions. Our fault–slip models well explain azimuthal variations of the RSTFs. The models show the rupture propagating toward the east–northeast (~50° to 60°), not along the fault strike. The slip distribution on the nodal plane striking north and dipping to the east shows a slip area of 1.1 km radius with the rupture propagating down–dip. The down–dip rupture may account for the observed lack of strong shaking in the epicentral area.

58 GEOSCIENCES↗
DOE OSTI2024

Spectral Analysis of Regular Material Point Method and its Application to Study High Pressure Reverse Osmosis Membrane Compaction and Embossing

Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.

material point method↗
DOE OSTI2024

Comparison of interband related optical transitions and excitons in ZnGeN 2 and GaN

The optical dielectric function of ZnGeN 2 is calculated from the interband transitions using the energy bands calculated in the quasiparticle self-consistent (QS)GW^ method using two different levels of approximation: the independent particle approximation (IPA) and the Bethe-Salpeter Equation (BSE) approach. The first allows us to relate peaks in ε 2 to specific bands and k-points but does not include electron-hole interaction effects. The second includes electron-hole interaction or excitonic effects. The corresponding changes in the shape of ε 2 (ω) are found to be similar to those in GaN. The screened Coulomb interaction W^ is here calculated including electron-hole interactions in the polarization function and gives a band structure already going beyond the random phase approximation. The static dielectric constants including only electronic screening, commonly referred to as ε∞ , were calculated separately by extrapolating the wave vector dependent macroscopic dielectric function in the long-wavelength limit. Below the quasiparticle gap, we find three bound excitons optically active for different polarization. The convergence of these bound excitons with respect to the density of the k-mesh used in the BSE is studied and found to require a fine mesh. It is also found that these bound excitons originate from only the lowest conduction band and the top three valence bands. To incorporate the lattice screening, we include a scaling factor (ε∞/ε 0 ) 2 , which allows us to obtain exciton binding energies of the correct order of magnitude similar to those in GaN. Here, the excitons are related to each of the three fold split valence bands and the splittings of the latter are also studied as function of strain. Finally, a relation between the anisotropic effective masses and the valence band splitting is pointed out and explained.

36 MATERIALS SCIENCE↗
DOE OSTI2022

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. So, we describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Revisiting Competing Paths in Electrochemical CO 2 Reduction on Copper via Embedded Correlated Wavefunction Theory

We re-evaluate two key steps in the mechanism of CO 2 reduction on copper at a higher level of theory capable of correcting inherent errors in density functional theory (DFT) approximations, namely, embedded correlated wavefunction (ECW) theory. Here, we consider the CO reduction step on Cu(111), which is critical to understanding reaction selectivity. Here, we optimize embedding potentials at the periodic plane-wave DFT level using density functional embedding theory (DFET). All possible adsorption sites (adsites) for each adsorbate then are screened with ECW theory at the catalytically active site to refine the local electronic structure. Unsurprisingly, DFT and ECW theory predict different adsite preferences, largely because of DFT’s inability to properly situate the CO 2π* level. Differing preferred adsites suggest that different reaction pathways could emerge from DFT versus ECW theory. Starting from these preferred ECW theory adsites, we then obtain reaction pathways at the plane-wave DFT level using the climbing-image nudged elastic band method to determine minimum energy paths. Thereafter, we perform ECW calculations at the catalytically active site to correct the energetics at each interpolated structure (image) along the reaction pathways. Via this approach, we confirm that the first step in CO reduction via hydrogen transfer on Cu(111) is to form hydroxymethylidyne (*COH) instead of formyl (*CHO). Although the prediction to preferentially form *COH is consistent with that of DFT, the two theories predict quite different structural and mechanistic behaviors, suggesting that verification is needed for other parts of the mechanism of CO 2 reduction, which is the subject of ongoing work.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

$S$ Hmax orientation in the Alpine region from observations of stress-induced anisotropy of nonlinear elasticity

The orientation of $S$ Hmax is commonly estimated from in situ borehole breakouts and earthquake focal mechanisms. Borehole measurements are expensive, and therefore sparse, and earthquake measurements can only be made in regions with many well-characterized earthquakes. Here, we derive the stress-field orientation using stress-induced anisotropy in nonlinear elasticity. In this method, we measure the strain derivative of velocity as a function of azimuth. We use a natural pump-probe (NPP) approach which consists of measuring elastic wave speed using empirical Green’s functions (probe) at different points of the earth tidal strain cycle (pump). The approach is validated using a larger data set in the Northern Alpine Foreland region where the orientation of maximum horizontal compressive stress is known from borehole breakouts and drilling-induced fractures. The technique resolves NNW-SSW to N-S directed $S$ Hmax which is in good agreement with conventional methods and the recent crustal stress model. We confirm that the NPP method can be applied to dense large-scale seismic arrays. The technique is then applied to the Southern Alps to understand the contemporary stress pattern associated with the ongoing deformation due to counterclockwise rotation of the Adriatic plate with respect to the European plate. Our results explain why the two major faults in Northeastern Italy, the Giudicarie Fault and the Periadriatic Line (Pustertal–Gailtal Fault) are currently inactive, while the currently acting stress field allows faults in Slovenia to deform actively. We have demonstrated that the pump-probe method has the potential to fill in the measurement gap left by conventional approaches, both in terms of regional coverage and in depth.

58 GEOSCIENCES↗
DOE OSTI2025

​​Hydrogen Transportation and Distributed Energy Systems Seismic Risk Assessment for Cascadia Subduction Zone Airport Facilities​

Portland International Airport in Oregon is exploring operating a fleet of 28 fuel cell electric buses to support airport operations and possibly provide backup power during outages. In this report, we evaluate the risk of hydrogen deployment at an airport considering the potential for seismic activity. We present a seismic risk assessment for pressurized piping for a generic hydrogen refueling station at the airport. Only pressurized piping was considered in the risk analysis because of its vulnerability and availability of fragility information unlike the other components. Seismic capacity of other components should be included as more information becomes available. To characterize the seismic hazard, we considered the initiating event frequency of a Cascadia Subduction Zone earthquake event based on data from the literature. Seismic stressors that translate to the site from an earthquake event are expressed as peak ground acceleration. Site disturbance is characterized as a function of soil conditions, represented by different shear wave velocity values, which affect how seismic waves propagate and impact structures. Pipe fragility curves as a function of seismic capacity were used to correlate ground acceleration and failure probability. These fragility characteristics were combined with site-specific soil conditions to calculate the probability of pipe rupture during a seismic event. It should be noted that multiple simultaneous failures due to an earthquake event as a common cause are not considered. Some epistemic uncertainties such as self-exciting shocks, aftershocks, and the time variant nature of the ground motion are not considered either.

08 HYDROGEN↗
DOE OSTI2025

Data for "Plasma shape effects on the Alfvén eigenmode spectrum through Alfvén slow-magnetosonic wave coupling"

electron density profile as a function of normalized poloidal flux: psinorm ne(10^20/m^3) dne/dpsiN Pressure profile as a function of normalized poloidal flux: psinorm ptot(KPa) dptot/dpsiN scaled to oftain beta=5% at magnetic axis Magnetic safety factor profile as a function of normalized poloidal flux: psinorm q_profile All needed to perform the NOVA calculations.

Alfven waves↗
DOE OSTI2021

First-principles characterization of the magnetic properties of Cu 2 (OH) 3 Br

Low-dimensional spin-1/2 transition metal oxides and oxyhalides continue to be at the forefront of research investigating nonclassical phases such as quantum spin liquids. In this work, we examine the magnetic properties of the oxyhalide Cu 2 (OH) 3 Br in the botallackite structure using first-principles density functional theory, linear spin-wave theory, and exact diagonalization calculations. This quasi-two-dimensional system consists of Cu 2+ S = 1/2 moments arranged on a distorted triangular lattice. Our exact diagonalization calculations, which rely on a first-principles-based magnetic model, generate spectral functions consistent with inelastic neutron scattering data. By performing computational experiments to disentangle the chemical and steric effects of the halide ions, we find that the dominant effect of the halogen ions is steric in the Cu 2 (OH) 3 X series of compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Space-Time Structured Plasma Waves

Electrostatic waves play a critical role in nearly every branch of plasma physics from fusion to advanced accelerators, to astro, solar, and ionospheric physics. The properties of planar electrostatic waves are fully determined by the plasma conditions, such as density, temperature, ionization state, or details of the distribution functions. Furthermore we demonstrate that electrostatic wave packets structured with space-time correlations can have properties that are independent of the plasma conditions. For instance, an appropriately structured electrostatic wave packet can travel at any group velocity, even backward with respect to its phase fronts, while maintaining a localized energy density. These linear, propagation-invariant wave packets can be constructed with or without orbital angular momentum by superposing natural modes of the plasma and can be ponderomotively excited by space-time structured laser pulses like the flying focus.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Amplitude analysis of the ${D}_s^{+}$ → π - π + π + decay

A Dalitz plot analysis of the ${D}_s^{+}$ → π - π + π + decay is presented. The analysis is based on proton-proton collision data recorded by the LHCb experiment at a centre of-mass energy of 8 TeV and corresponding to an integrated luminosity of 1.5 fb -1 . The resonant structure of the decay is obtained using a quasi-model-independent partial-wave analysis, in which the π + π - S-wave amplitude is parameterised as a generic complex function determined by a fit to the data. The S-wave component is found to be dominant, followed by the contribution from spin-2 resonances and a small contribution from spin-1 resonances. The latter includes the first observation of the ${D}_s^{+}$ → ω(782)π + channel in the ${D}_s^{+}$ → π - π + π + decay. The resonant structures of the ${D}_s^{+}$ → π - π + π + and ${D}^{+}$ → π - π + π + decays are compared, providing information about the mechanisms for the hadron formation in these decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Unveiling the mechanism of surface corrugation formation on a quasi free-standing bi-layer graphene via experimental and modeling investigations

Extensive studies have been conducted to accomplish the controls over the formation of surface corrugations on chemical vapor deposited (CVD) graphene. However, the underlying mechanism of the surface corrugation formation on graphene remains elusive, due to the difficulty delineating the growth-induced surface corrugations from wrinkles formed during graphene transfer with the use of a supporting layer casted on top of the sample. Here, in this study, “quasi free-standing” graphene was realized using the weakly-interacting h-BN layer during a wet transfer process to study the surface corrugation formation on graphene. The weak interaction of graphene with the supporting h-BN layer allows for the growth-induced surface corrugations to be released, enabling the investigation of the primary mechanism for the graphene corrugation behavior during the transfer process. The experimental measurements of surface corrugations were conducted using atomic force microscopy (AFM) to reveal that surface corrugations are formed following the graphene’s crystallographic directions. Density functional theory calculations using the projector augmented-wave (PAW) method and Perdew-Burke-Ernzerhof (PBE) exchange correlation functional were also conducted for two models of “atomic-scale ripples” and “nanoscale wrinkles” and showed that all of them are formed along the crystallographic axis of graphene, which is in good agreement with the experimental observations.

36 MATERIALS SCIENCE↗