SEARCH · Search NASA

Results for “wave functions”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Advanced Laboratory and Field Arrays (ALFA)/Lab Collaboration Project (LCP) for Marine Energy (Final Scientific/Technical Report)

The objective of the Advanced Laboratory and Field Arrays (ALFA) project was to reduce the Levelized Cost of Energy (LCOE) of Marine and Hydrokinetic (MHK) energy by leveraging research, development, and testing capabilities at Oregon State University, University of Washington, and the University of Alaska, Fairbanks. ALFA is a project within the Pacific Marine Energy Center (PMEC; formerly NNMREC), a multi-institution entity with a diverse funding base that focuses on research and development for marine renewables. The ALFA project aimed to accelerate the development of next-generation arrays of wave energy conversion (WEC) and tidal energy conversion (TEC) devices through a suite of field-focused R&D activities spanning a broad range of strategic opportunity areas identified in the Funding Opportunity Announcement: • Device and/or array operation and maintenance (O&M) logistics development; • High-fidelity resource characterization and/or modeling technique development and validation; • Array-specific component technology development (e.g. moorings and foundations, transmission, and other offshore grid components); • Array performance testing and evaluation; and • Novel cost-effective environmental monitoring techniques and instrumentation testing and evaluation. The objective of the Lab Collaboration Project (LCP) was to accelerate the development of next-generation marine energy conversion systems. The LCP aimed to achieve these project objectives in collaboration with the national laboratories by: • Developing concept generation and assessment tools; • Improving access to existing testing resources; • Validating collision risk models between fish and turbines; and • Advancing analysis and simulation capabilities for wave-WEC interactions and PTO analysis in nonlinear ocean waves. The ALFA portion of the project was comprised of six overarching technical tasks: • Task 1: Debris Modeling, Detection and Mitigation; • Task 2: Autonomous Monitoring & Intervention; • Task 3: Resource Characterization for Extreme Conditions; • Task 4: Robust Models for Design of Offshore Anchoring and Mooring Systems; • Task 5: Performance Enhancement for Marine Energy Converter (MEC) Arrays; and • Task 6: Evaluating Sampling Techniques for MHK Biological Monitoring. The LCP was divided into four overarching technical tasks: • Task 7: Project Management and Reporting • Task 8: Novel Design and Assessment Methodologies for Wave Energy Converter Design (Wave- SPARC) • Task 9: Testing Access for Commercial Marine Renewable Energy Technology Developers • Task 10: Quantifying Collision Risk for Fish and Turbines • Task 11: Nonlinear Ocean Waves and PTO Control Strategy Each ALFA/LCP task listed above functioned as a separate and discreet project. A final Technical Report was written for each individual task and these reports were uploaded to OSTI, after receiving DOE approval. The following document is a compilation of each of these final, approved reports arranged as individual chapters.

13 HYDRO ENERGY↗

DOE OSTI2024

Effects of low-frequency voltage on nonlinear standing wave excitation, plasma uniformity, and ion dynamics in dual-frequency asymmetric capacitive discharges

Abstract It is known that in very-high-frequency (VHF) capacitively coupled plasmas, the higher harmonics generated by nonlinear sheath motion can enhance the standing wave effect (SWE), which can lead to center-peaked plasma density profiles. In this work, an improved nonlinear electromagnetic model incorporating a transmission line model, an electron momentum balance model, a bulk plasma model, a collisionless nonlinear numerical sheath model, and an ion Monte-Carlo collision (MCC) model is developed to study the effects of low-frequency (LF) voltage V L on the nonlinear standing wave excitation, plasma uniformity, and ion energy and angular distribution functions (IEDFs and IADFs) in dual-frequency (DF) asymmetric capacitive argon discharges at relatively low pressure of 3 Pa. The plasma diffusion in the radial direction and ion dynamics within the LF oscillating sheath are self-consistently considered. The LF voltage V L at 2 MHz varies from 0 to 700 V while the HF voltage V H at 60 MHz is fixed at 100 V. Simulation results indicate that without the addition of an LF source (i.e. V L = 0 V), there are a considerable number of high-order harmonics with short wavelengths, leading to significant SWE and central peak in the radial plasma density profile. Nevertheless, the high-order harmonic excitations tend to be weakened and merely occur around the phase of the full LF sheath collapse due to a shorter characteristic damping time of the surface waves as V L increases. This, combined with increased surface wavelengths of both the driving frequency and the higher harmonics at a higher V L , leads to suppressed standing waves and improved plasma uniformity. Meanwhile, the simulations show that both the low and the high energy peaks of IEDF move towards higher energies, and the energy peak separation width ΔEbecomes wider with the increase of V L . The IEDF at the radial center of the powered electrode exhibits a broader ΔEthan that at the edge. For the IADF, an increased V L results in more ions incident on the electrode with a smaller deflection angle. Because of a thinner sheath and a higher sheath voltage at the electrode center, the peak value of IADF at the electrode center is greater than that at the edge.

Physics↗

DOE OSTI2025

A Cost-Effective Wave Energy Harvesting System with Maximum Power Point Tracking Capability

In this paper, we propose a low-cost bidirectional four-switch inverter that transfers three-phase energy between a permanent magnet synchronous generator (PMSG) driven by a wave energy collector and the power grid or a local ac load. To achieve high control performance, the machine-side inverter uses a high-gain observer for rotor position estimation and a Lyapunov-based approach for torque control. For the gridside inverter control, we estimate the filter capacitor voltage with a model reference adaptive system and control the dc-link voltage with a Lyapunov-based energy function. We also track the maximum power point from the sea waves with an integral control law that actively matches the generator impedance with that of the wave collection device. The proposed system supports both grid-following/forming modes. We validate our control design with simulation results on a 3-kW hardware setup.

wave energy, power systems integration, control de↗

DOE OSTI2024

The slow wave resonance cone in the collisional regime

In the low-density edge plasma of tokamaks, ion cyclotron range of frequencies actuators may parasitically emit slow waves. If the density is sufficiently low, which may be common in large future devices such as international thermonuclear experimental reactor (ITER), these slow waves take the form of so-called resonance cones. The traditional theoretical description of this wave mode relies on formally relating an electrostatic approximation of the frequency-domain wave equation to a time-domain wave equation and relating the cone angle to the wave speed in the time-domain wave equation. In the cold plasma collisional regime, that wave speed is complex. We investigate that scenario in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Nanoscale strain wave generation by a piezoelectric grating from polar vortices

Nanostructures formed by spontaneously broken symmetry have provided new ways to manipulate quantum states. Specifically, topological structures with periodic spatial ordering, such as polar vortices and skyrmions, can be ideal hosts for creating engineered responses in both spatial and frequency domains. So far, however, only a few examples of such hierarchical engineering have been reported in the literature. Here we demonstrate that the spatially modulated piezoelectric response of a polar vortex structure can create strain waves with a characteristic nanoscale wavefront. Using time-resolved pump–probe resonant X-ray scattering and diffraction measurements, coupled with dynamical phase-field simulations, we show that the piezoelectric modulation of the spontaneously formed polar vortex crystal functions as an acoustic diffraction grating. This system converts incoming laterally uniform strain waves into outgoing waves with a characteristic sub-terahertz frequency, driven by an intrinsic excitation of the polar vortex crystal. Moreover, our phase-field simulations suggest that the dynamic mechanical displacements exhibiting vortex textures are generated from both space- and time-varying piezoelectric responses. Our findings illustrate a new method for generating nanoscale strain waves with unique spatial textures by tuning the hierarchical order of polar topologies to engineer new collective modes, allowing for a wide range of control through the topological lattice.

ferroelectrics↗

DOE OSTI2024

Pressure-induced polyamorphic transition in CaA l 2 O 4 glass

In situ high-pressure ultrasonic velocity measurements of CaAl 2 O 4 glass reveal abrupt irreversible discontinuities in the elastic wave velocities at ~8–10 GPa. Total structure factor and pair distribution functions measured by synchrotron x-ray diffraction show a rapid change in the intermediate range structure attributed to a rearrangement of calcium ions over this narrow pressure region. Atomistic models obtained from molecular dynamics simulations reveal that this intermediate range structure is explained by a transition of Ca–O void radius distribution from a bimodal distribution with peaks at ~2.1 and ~2.4 Å to a single distribution centered at ~2.1 Å. Further, these abrupt structural changes involving the rapid increase in elastic wave velocity are markedly different to the continuous transformations observed in conventional network-forming glasses, such as SiO 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

How Does Metal Spin State Affect Electronic Communication in Mixed-Valence Dimers? Insights from Ultrafast Near-Infrared and Soft X-ray Transient Absorption Spectroscopy

Recent advancements in photocatalysis, photovoltaics, and quantum information science take advantage of electron spin, and determining how spin multiplicity affects electron transfer is key to understanding these phenomena. Here, in this study, we examine how metal spin state affects electronic communication in an organometallic mixed-valence dimer, ferrocenyl cobaltocenium ([Fe II Cp 2 Co III Cp 2 ] + ). This complex can be photoexcited from its low-spin singlet Fe II ground state to form intermediate-spin triplet Fe II and high-spin quintet Fe II excited states. Using femtosecond optical transient absorption (OTA) spectroscopy with visible (vis), near-infrared (NIR), and short-wave IR (SWIR) probes, supported by time-dependent density functional theory (TD-DFT) calculations, we measure Fe II Co III →Fe III Co II intervalence charge transfer (IVCT) bands in each of the Fe II spin states. Mulliken–Hush analysis of the excited-state IVCT bands was used to compute the electronic coupling between the metal centers in all three spin states, which increased as quintet < triplet < singlet. Meanwhile, the peak energy of the bands, and thus the ΔG of the IVCT transition, increased as triplet < quintet < singlet. Then, to directly probe the electronic structure at each metal center, we employed picosecond soft X-ray transient absorption (XTA) spectroscopy at the Fe and Co L 3 edges. Our results show that the low-spin and high-spin states of [Fe II Cp 2 Co III Cp 2 ] + are valence-localized, while the intermediate-spin state is partially delocalized. The differences in charge delocalization are attributed to differences in orbital occupation and geometry that affect the free energy and superexchange coupling.

Burke, John H. [Univ. of Illinois at Urbana-Champa↗

DOE OSTI2025

Anisotropic light-tailored RKKY interaction in two-dimensional 𝑑-wave altermagnets

Altermagnets are known in spintronics for their intrinsic spin-splitting and unconventional magnetic responses, particularly to magnetic impurities. However, effectively controlling the magnetic exchange interactions in altermagnets is challenging for practical applications. Here, in this work, we propose using circularly polarized light to tune the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction in two-dimensional 𝑑-wave altermagnets. Using the real-space retarded Green's functions approach, our results show that while the Heisenberg and Ising exchanges dominate, a notable Dzyaloshinskii–Moriya (DM) interaction also plays a key role. Furthermore, the inherent strength of altermagnetism imprints chirp-like signatures into the magnetic responses, which can be dynamically tuned via light. We mainly demonstrate that gate-induced Rashba spin-orbit coupling is essential in response to light—light selectively and anisotropically adjusts the DM interaction without affecting the other exchanges. Our findings further indicate that rotating the altermagnet by 45° relative to the light's polarization direction generates a Dirac-like dispersion and different DM interactions. We finally extract critical thresholds where light reverses DM interactions along one axis or balances both in-plane components. The anisotropic light-driven control of RKKY interactions in altermagnets not only highlights their unique properties but also opens new avenues for engineering tailored magnetic characteristics in spintronic applications.

altermagnetism↗

DOE OSTI2024

code for "Improving the QBO forcing by resolved waves with vertical grid refinement in E3SMv2"

This is all code that used in the journal article "Improving the QBO forcing by resolved waves with vertical grid refinement in E3SMv2", including: Running and post processing E3SM: Run E3SM: run_E3SM.2023.scidac.MMF.amip.py Regrid: regrid.E3SM.py post processing, diagnostics, archive: run_post.2023.scidac.MMF.py TEM calculation, including filter TEM, E3SM: calculate.TEM.v2.py calculate.TEM.monthly_mean.v1.py TEM, ERA5: ERA5_calculate_TEM_90.py merge TEM data: merge_TEM_E3SM_ERA5.py TEM, filter, E3SM: filter_TEM_E3SM.py TEM, filter, ERA5: filter_TEM_ERA5.py WK wave analysis code from (https://github.com/brianpm/wavenumber_frequency): wavenumber_frequency_functions.py calculate the wave spectra so easy to plot: wk_calculation.py Make figures some functions for plot: plot.py code for drawing paper plots: figures_for_paper.ipynb

Hannah, WalterM↗

DOE OSTI2024

Tracking Seismic Velocity Perturbations at Ridgecrest Using Ballistic Correlation Functions

Abstract We present results based on data of a dense nodal array composed of 147 stations, deployed in 2022 near the epicenter of the 2019 Mw 7.1 Ridgecrest earthquake to investigate characteristics of the seismic wavefields. Through array analyses, we identified two primary components. First, we observed far-field P waves dominating the 0.5–1.2 Hz frequency range, which are likely primarily generated by wind-driven oceanic swell activity. Second, we detected near-field body waves resulting from anthropogenic activities in the frequency range 2–8 Hz. We examined noise correlation functions derived from data of the dense deployment and regional stations to explore fault-zone seismic velocity changes using ballistic arrivals, with a focus on velocity perturbation shortly before and after the Ridgecrest earthquake sequence. Our findings exhibit distinct behavior compared to results obtained through standard coda-wave interferometry. Particularly, we observed a decrease in P-wave travel time on certain station pairs prior to the 2019 earthquake sequence. Supported by detailed investigation of the local seismic wavefields, we interpret the decreasing P-wave travel time as likely caused by a velocity increase away from the fault, possibly related to fluid migration. However, additional information is necessary to verify this hypothesis.

Geochemistry & Geophysics↗

DOE OSTI2025

Tunable magnons of an antiferromagnetic Mott insulator via interfacial metal-insulator transitions

Antiferromagnetic insulators present a promising alternative to ferromagnets due to their ultrafast spin dynamics essential for low-energy terahertz spintronic device applications. Magnons, i.e., quantized spin waves capable of transmitting information through excitations, serve as a key functional element in this paradigm. However, identifying external mechanisms to effectively tune magnon properties has remained a major challenge. Here we demonstrate that interfacial metal-insulator transitions offer an effective method for controlling the magnons of Sr 2 IrO 4 , a strongly spin-orbit coupled antiferromagnetic Mott insulator. Resonant inelastic x-ray scattering experiments reveal a significant softening of zone-boundary magnon energies in Sr 2 IrO 4 films epitaxially interfaced with metallic 4 d transition-metal oxides. Therefore, the magnon dispersion of Sr 2 IrO 4 can be tuned by metal-insulator transitions of the 4 d transition-metal oxides. We tentatively attribute this non-trivial behavior to a long-range phenomenon mediated by magnon-acoustic phonon interactions. Our experimental findings introduce a strategy for controlling magnons and underscore the need for further theoretical studies to better understand the underlying microscopic interactions between magnons and phonons.

magnetic properties and materials↗

DOE OSTI2025

Unraveling trace anomaly of supradense matter via neutron star compactness scaling

The trace anomaly Δ ≡ 1/3 −𝑃/𝜖 =1/3 −𝜙 quantifies the possibly broken conformal symmetry in supradense matter under pressure 𝑃 at energy density 𝜖. Perturbative QCD (pQCD) predicts a vanishing Δ at extremely high energy or baryon densities when the conformal symmetry is realized but its behavior at intermediate densities reachable in neutron stars (NSs) is still very uncertain. The extraction of Δ from NS observations strongly depends on the employed model for nuclear equation of state (EOS). Using the IPAD-TOV method based on an intrinsic and perturbative analysis of the dimensionless (IPAD) Tolman-Oppenheimer-Volkoff (TOV) equations that are further verified numerically by using 10 5 EOSs generated randomly with a metamodel in a very broad EOS parameter space constrained by terrestrial nuclear experiments and astrophysical observations, here we first show that the compactness 𝜉 ≡ 𝐺⁡𝑀 NS /𝑅⁢𝑐 2 ≡ 𝑀 NS /𝑅 of a NS with mass 𝑀 NS and radius 𝑅 scales very accurately with $\bar{Π}$ c ≡ $Π$ c · (1 +18⁢X/25) ≡ X/(1 +3⁢X 2 +4⁢X) · (1 +18⁢X/25) where X ≡ 𝜙 c = 𝑃 c /𝜖 c is the ratio of pressure over energy density at NS centers. The scaling of NS compactness thus enables one to readily read off the central trace anomaly Δ c = 1/3 −X directly from the observational data of either the mass-radius or red-shift measurements. Finally, we then demonstrate indeed that the available NS data themselves from recent X-ray and gravitational wave observations can determine model insensitively the trace anomaly as a function of energy density in NS cores, providing a stringent test of existing NS models and a clear guidance in a new direction for further understanding the nature and EOS of supradense matter.

nuclear astrophysics↗

DOE OSTI2024

First-order structural phase transition at low temperature in GaPt 5 P and its rapid enhancement with pressure

Single crystals of X Pt 5 ⁢P (X = Al, Ga, and In), belonging to the 1-5-1 family of compounds, were grown from a Pt-P solution at high temperatures, and measurements of the ambient pressure, temperature-dependent magnetization, resistivity, and x-ray diffraction were made. Additionally, the ambient-pressure Hall resistivity and temperature-dependent resistance under pressure were measured on GaPt 5 ⁢P. All three compounds have a tetragonal P4/mmm crystal structure at room temperature with metallic transport and weak diamagnetism over the 2–300 K temperature range. Surprisingly, at ambient pressure, both the transport and magnetization measurements on GaPt 5 ⁢P show a steplike feature in the 70–90 K region, suggesting a possible structural phase transition. Neither AlPt 5 ⁢P nor InPt 5 ⁢P have any signatures of a phase transition in their temperature-dependent electrical resistance and magnetization data. Both the hysteretic nature and sharpness of the features in the GaPt 5 ⁢P data suggest that the transition is first-order. Further, single-crystal x-ray diffraction measurements provided further details of the structural transition with a possibility of a crystal symmetry different from P⁢4/mmm below the transition temperature. The transition is characterized by anisotropic changes in the lattice parameters and a volume collapse with respect to the high-temperature tetragonal crystal structure. Furthermore, satellite peaks are observed at two distinct and nonequivalent wave vectors (0, 0, 0.5) and (0.5, 0.5, 0.5), and density functional theory calculations present phonon softening, especially at (0.5, 0.5, 0.5), as a possible driving mechanism. Additionally, we find that the structural transition temperature increases rapidly with increasing pressure, reaching room temperature by ~2.2 GPa, highlighting the high degree of pressure sensitivity of GaPt 5 ⁢P and fragile nature of its room-temperature structure. Even though the volume collapse and extreme pressure sensitivity suggest chemical pressure should drive a similar structural change in AlPt 5⁢ P, where both unit-cell dimensions and volume are smaller, its structure is found to be the same as that of the room-temperature GaPt 5 ⁢P. Overall, GaPt 5 ⁢P stands out as a sole member of the 1-5-1 family of compounds for which a temperature-driven structural change has been observed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Lorentzian OPE inversion formula: A geometric perspective

We give a new perspective on the Lorentzian operator product expansion inversion formula [S. Caron-Huot, Analyticity in spin in conformal theories, .; D. Simmons-Duffin, D. Stanford, and E. Witten, A spacetime derivation of the Lorentzian OPE inversion formula, .], building on [P. Agarwal , companion paper, Embedding space approach to Lorentzian CFT amplitudes and causal spherical functions, .]. We introduce an “auxiliary” fourpoint function that can be related to the traditionally defined ones via a Radon transform. The Mellin amplitudes associated with this auxiliary function can be shown to be equivalent to the conventional partial wave amplitudes. This has the intuitive geometrical meaning of a generalization of the projection-slice theorem. Published by the American Physical Society 2025

Agarwal, Pulkit (ORCID:0000000346581691)↗

DOE OSTI2024

A Nonlinear Least Squares Phasor Estimation Algorithm with a Trust Metric

The paper presents a separable nonlinear least squares (NLLS)-based approach for estimation of fundamental-frequency phasors from sampled point-on-wave measurements. An analytical connection is established between the NLLS cost function and the discrete Fourier transform (DFT)-based periodogram of the input signal. This periodogram-based interpretation of the cost function offers an intuitive and easy-to-implement solution for the frequency estimate. Using the residual error of the NLLS-fit, the paper also presents an insightful measure for ascertaining the quality of the phasor estimates and their validity, especially for data windows containing signal transients.

Chatterjee, Kaustav↗

DOE OSTI2026

Kramers nodal line in the charge density wave state of YTe 3 and the influence of twin domains

Recent studies have focused on the relationship between charge density wave (CDW) collective electronic ground states and nontrivial topological states. YTe 3 , a nonmagnetic quasi-two-dimensional chalcogenide, has been reported to exhibit a CDW state below 334 K. Using angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT), we establish that YTe 3 is a CDW-induced Kramers nodal line (KNL) metal, a recently proposed topological state of matter. Scanning tunneling microscopy and low energy electron diffraction reveal two orthogonal domains, each with a unidirectional CDW and a similar wave vector (𝐪 CDW ). When the influence of twin domains is considered, the effective band structure (EBS) computations that utilize DFT-calculated bands using a noncentrosymmetric structure determined by x-ray crystallography, show excellent agreement with ARPES. The noncentrosymmetry of YTe 3 is established by Raman spectroscopy. The Fermi surface and ARPES intensity plots show weak shadow bands displaced by 𝐪 CDW from the main bands. Furthermore, these are linked to CDW modulation, as the EBS calculation confirms. Bilayer split main and shadow bands suggest the existence of crossings, according to theory and experiment. DFT bands, including spin-orbit coupling, indicate existence of a KNL along the Σ direction from multiple crossings of bands dispersing perpendicular to it. Additionally, doubly degenerate bands are only found along the KNL at all energies, with some bands dispersing through the Fermi level.

Angle-resolved photoemission spectroscopy↗

DOE OSTI2025

Characterization of Gramicidin A in Triblock and Diblock Polymersomes and Hybrid Vesicles via Continuous Wave Electron Paramagnetic Resonance Spectroscopy

Studying membrane proteins in a native environment is crucial to understanding their structural and/or functional studies. Often, widely accepted mimetic systems have limitations that prevent the study of some membrane proteins. Micelles, bicelles, and liposomes are common biomimetic systems but have problems with membrane compatibility, limited lipid composition, and heterogeneity. To overcome these limitations, polymersomes and hybrid vesicles have become popular alternatives. Polymersomes form from amphiphilic triblock or diblock copolymers and are considered more robust than liposomes. Hybrid vesicles are a combination of lipids and block copolymers that form vesicles composed of a mixture of the two. These hybrid vesicles are appealing because they have the native lipid environment of bilayers but also the stability and customizability of polymersomes. Gramicidin A was incorporated into these polymersomes and characterized using continuous wave electron paramagnetic resonance (CW-EPR) and transmission electron microscopy (TEM). EPR spectroscopy is a powerful biophysical technique used to study the structure and dynamic properties of membrane proteins in their native environment. Spectroscopic studies of gramicidin A have been limited to liposomes; in this study, the membrane peptide is studied in both polymersomes and hybrid vesicles using CW-EPR spectroscopy. Lineshape analysis of spin-labeled gramicidin A revealed linewidth broadening, suggesting that the thicker polymersome membranes restrict the motion of the spin label more when compared to liposome membranes. Statement of Significance: Understanding membrane proteins’ structures and functions is critical in the study of many diseases. In order to study them in a native environment, membrane mimetics must be developed that can be suitable for obtaining superior biophysical data quality to characterize structural dynamics while maintaining their native functions and structures. Many currently widely accepted methods have limitations, such as a loss of native structure and function, heterogeneous vesicle formation, restricted lipid types for the vesicle formation for many proteins, and experimental artifacts, which leaves rooms for the development of new biomembrane mimetics. The triblock and diblock polymersomes and hybrid versicles utilized in this study may overcome these limitations and provide the stability and customizability of polymersomes, keeping the biocompatibility and functionality of liposomes for EPR studies of membrane proteins.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗