Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “covariance and correlation”

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 OSTI2026

Parity-odd four-point correlation function from the DESI data release 1 luminous red galaxy sample

The parity-odd four-point function provides a unique probe of fundamental symmetries and potential new physics in the large-scale structure of the Universe. We present measurements of the parity-odd four-point function using the Dark Energy Spectroscopic Instrument (DESI) DR1 luminous red galaxy (LRG) sample and assess its detection significance. Our analysis considers both auto- and cross-correlations, using two complementary approaches to the covariance: (i) the full analytic covariance matrix applied to the uncompressed data vector, and (ii) a compressed data vector combined with a hybrid covariance matrix constructed from simulations and analytic estimates. When using the full analytic covariance matrix without corrections, we observe apparent auto-correlation signals with significance up to 4⁢𝜎. However, this excess is also consistent with a mismatch between the statistical fluctuations estimated from the simulations and those present in the real data. Our findings therefore suggest that the parity-odd signal in the current DESI DR1 LRG sample is consistent with zero. We note, however, that the low completeness of this sample may have a non-negligible impact on the detection sensitivity. Future data releases with improved completeness will be crucial for further investigation.

Hou, Jiamin [Ludwig-Maximilians-Universität; Unive↗
DOE OSTI2024

Fermionic mean-field theory as a tool for studying spin Hamiltonians

The Jordan–Wigner transformation permits one to convert spin 1/2 operators into spinless fermion ones, or vice versa. In some cases, it transforms an interacting spin Hamiltonian into a noninteracting fermionic one, which is exactly solved at the mean-field level. Even when the resulting fermionic Hamiltonian is interacting, its mean-field solution can provide surprisingly accurate energies and correlation functions. Furthermore, Jordan–Wigner is, however, only one possible means of interconverting spin and fermionic degrees of freedom. Here, we apply several such techniques to the XXZ and J 1 –J 2 Heisenberg models, as well as to the pairing or reduced Bardeen–Cooper–Schrieffer Hamiltonian, with the aim of discovering which of these mappings is most useful in applying fermionic mean-field theory to the study of spin Hamiltonians.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Evolution of intermittent filaments in the scrape-off layer of NSTX

Filamentary structures naturally arise from background turbulence in the scrape-off layer (SOL) of plasmas, leading to significant particle and heat transport that can degrade overall plasma confinement. This enhanced transport can contribute to unacceptably high heat loads on plasma-facing components. As such, understanding the physics of SOL plasma filaments is critical for predicting and mitigating their effects in future fusion devices. On the National Spherical Torus Experiment (NSTX), plasma filaments—commonly referred to as “blobs”—were investigated using the gas-puff imaging (GPI) diagnostic in the edge and SOL regions. The analysis involved identifying, segmenting, and tracking the characteristic contours of the blobs in each frame of the GPI video sequences. Their evolution was characterized through shape descriptors, velocity, and angular velocity derived from their contour coordinates. The results indicate that as the blob area increases, their shapes become more concave and less circular, suggesting reduced structural stability in larger blobs. This result aligns with previous theoretical results where it was shown that larger blobs are more susceptible to instabilities [Krasheninnikov et al., J. Plasma Phys. 74, 679–717 (2008) and D'Ippolito et al., Phys. Plasmas 18, 060501 (2011)]. A positive correlation was observed between radial velocity and radial position, suggesting radially outward acceleration of the filaments, potentially driven by decreasing viscous drag toward the far SOL. Interestingly, blobs in background SOL turbulence exhibited minimal spinning in contrast to filaments originating from edge localized modes, which show substantial rotation during their paths [Lampert et al., Phys. Plasmas 29, 102502 (2022)]. Statistical analysis of the solidity and total curvature shape descriptors, along with their temporal evolution, revealed relatively broad, near-Gaussian distributions. This suggests that blob morphology is strongly influenced by stochastic turbulent processes in the surrounding plasma environment. Blob parameters were also compared with bulk plasma and radial profile measurements. Notable trends were found between blob rotation and poloidal velocity with collisionality and line-integrated density. These findings contribute to a deeper understanding of blob dynamics and provide valuable insights for refining SOL turbulence models.

Covariance and correlation↗
DOE OSTI2023

Surface coverage dynamics for reversible dissociative adsorption on finite linear lattices

Dissociative adsorption onto a surface introduces dynamic correlations between neighboring sites not found in non-dissociative absorption. We study surface coverage dynamics where reversible dissociative adsorption of dimers occurs on a finite linear lattice. We derive analytic expressions for the equilibrium surface coverage as a function of the number of reactive sites, N, and the ratio of the adsorption and desorption rates. Using these results, we characterize the finite size effect on the equilibrium surface coverage. For comparable N’s, the finite size effect is significantly larger when N is even than when N is odd. Moreover, as N increases, the size effect decays more slowly in the even case than in the odd case. The finite-size effect becomes significant when adsorption and desorption rates are considerably different. These finite-size effects are related to the number of accessible configurations in a finite system where the odd-even dependence arises from the limited number of accessible configurations in the even case. We confirm our analytical results with kinetic Monte Carlo simulations. We also analyze the surface-diffusion case where adsorbed atoms can hop into neighboring sites. As expected, the odd-even dependence disappears because more configurations are accessible in the even case due to surface diffusion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2006

Using CO2:CO Correlations to Improve Inverse Analyses of Carbon Fluxes

Observed correlations between atmospheric concentrations of CO2 and CO represent potentially powerful information for improving CO2 surface flux estimates through coupled CO2-CO inverse analyses. We explore the value of these correlations in improving estimates of regional CO2 fluxes in east Asia by using aircraft observations of CO2 and CO from the TRACE-P campaign over the NW Pacific in March 2001. Our inverse model uses regional CO2 and CO surface fluxes as the state vector, separating biospheric and combustion contributions to CO2. CO2-CO error correlation coefficients are included in the inversion as off-diagonal entries in the a priori and observation error covariance matrices. We derive error correlations in a priori combustion source estimates of CO2 and CO by propagating error estimates of fuel consumption rates and emission factors. However, we find that these correlations are weak because CO source uncertainties are mostly determined by emission factors. Observed correlations between atmospheric CO2 and CO concentrations imply corresponding error correlations in the chemical transport model used as the forward model for the inversion. These error correlations in excess of 0.7, as derived from the TRACE-P data, enable a coupled CO2-CO inversion to achieve significant improvement over a CO2-only inversion for quantifying regional fluxes of CO2.

Palmer, Paul I.↗
DOE OSTI2026

Correlation function metrology for warm dense matter: Recent developments and practical guidelines

X-ray Thomson scattering (XRTS) has emerged as a valuable diagnostic for matter under extreme conditions, as it captures the intricate many-body physics of the probed sample. Recent advances, such as the model-free temperature diagnostic of Dornheim et al. [Nat. Commun. 13 , 7911 (2022)], have demonstrated how much information can be extracted directly within the imaginary-time formalism. However, since the imaginary-time formalism is a concept often difficult to grasp, we provide here a systematic overview of its theoretical foundations and explicitly demonstrate its practical applications to temperature inference, including relevant subtleties. Furthermore, we present recent developments that enable the determination of the absolute normalization, Rayleigh weight, and density from XRTS measurements without reliance on uncontrolled model assumptions. Finally, we outline a unified workflow that guides the extraction of these key observables, offering a practical framework for applying the method to interpret experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Allosteric prediction via convolutional neural networks and protein structural and dynamical features

Allostery is the phenomenon whereby a binding event or covalent modification at one site in a protein modulates function at a distal site, thus changing a protein’s functional state. As such, it is a ubiquitous aspect of protein functional regulation. Computationally predicting allosteric states is important as part of the broader challenge of functional annotation, but it also has practical implications for drug development, as targeting an allosteric site often affords greater specificity compared with targeting an orthosteric site. This study introduces a machine learning approach to predict the allosteric functional state using the small G-protein KRas as the model system, due to its implication in many types of cancer and being well studied as a result with many x-ray crystallographic structures of KRas available with different mutations and ligands bound. Using structural and dynamical features that can be cast as images, namely interatomic distances, contact maps, covariance, and mutual information, supervised learning was performed using convolutional neural networks. Two pretrained convolutional neural network architectures, GoogLeNet and ResNet18, were fine-tuned to classify KRas into active or inactive states based on these features. Across training regimes, atomic contact maps emerged as the most effective structural feature, whereas linearized mutual information outperformed covariance in capturing dynamical correlations relevant to allostery. Models achieved significant validation accuracy, with atomic contact maps yielding up to 90% accuracy. In conclusion, the findings suggest that integrating global structural rearrangements and correlated motion patterns with deep learning can reliably predict protein allosteric states, offering a promising framework for understanding allosteric regulation and developing targeted therapeutics.

Rajeshwar T., Rajitha [Oak Ridge National Laborato↗
DOE OSTI2020

Efficient Langevin dynamics for “noisy” forces

Efficient Boltzmann-sampling using first-principles methods is challenging for extended systems due to the steep scaling of electronic structure methods with the system size. Additionally, stochastic approaches provide a gentler system-size dependency at the cost of introducing “noisy” forces, which could limit the efficiency of the sampling. When the forces are deterministic, the first-order Langevin dynamics (FOLD) offers efficient sampling by combining a well-chosen preconditioning matrix S with a time-step-bias-mitigating propagator [G. Mazzola and S. Sorella, Phys. Rev. Lett. 118, 015703 (2017)]. However, when forces are noisy, S is set equal to the force-covariance matrix, a procedure that severely limits the efficiency and the stability of the sampling. Here, we develop a new, general, optimal, and stable sampling approach for FOLD under noisy forces. We apply it for silicon nanocrystals treated with stochastic density functional theory and show efficiency improvements by an order-of-magnitude.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. In this study, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Dynamics of filaments during the edge-localized mode crash on NSTX

Edge localized modes (ELMs) are routinely observed in H-mode plasma regimes of the National Spherical Torus Experiment (NSTX). Due to the explosive nature of the instability, only diagnostics with high temporal and spatial resolution could provide a detailed insight into the dynamics associated with the ELMs. Gas-puff imaging at NSTX provides 2D measurements of the magnetic field aligned fluctuations (e.g., ELM filaments) in the scrape-off layer and at the plasma edge with 2.5 μs temporal and 10 mm optical resolution. A novel analysis technique was developed to estimate the frame-by-frame velocities and the spatial parameters of the dominant structures associated with the ELMs. The analysis was applied to single ELM events to characterize the ELM crash dynamics and then extended to a database of 159 ELM events. Statistical analysis was performed in order to find the characterizing dynamics of the ELM crash. The results show that on average, an ELM crash consists of a filament with a circular cross section, which is propelled outward with a characterizing peak radial velocity of ~3.3 km/s. The radial velocity was found to be linearly dependent on the distance of the filament from the separatrix, which has never been seen before. The ELM filament is characterized by propagation in the ion-diamagnetic direction poloidally with a peak velocity of 11.4 km/s. The ELM crash lasts for approximately 100 μs until the radial propulsion settles back to the pre-ELM level. Furthermore, the experimental findings were compared with analytical theory. Two possible mechanisms were identified for explaining the observations: the curvature interchange model and the current–filament interaction model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Equation-of-state, sound speed, and reshock of shock-compressed fluid carbon dioxide

Mechanical equation-of-state data of initially liquid and solid CO2 shock-compressed to terapascal conditions are reported. Diamond-sapphire anvil cells were used to vary the initial density and state of CO2 samples that were then further compressed with laser-driven shock waves, resulting in a data set from which precise derivative quantities, including Grüneisen parameter and sound speed, are determined. Reshock states are measured to 800 GPa and map the same pressure-density conditions as the single shock using different thermodynamic paths. The compressibility data reported here do not support current density-functional-theory calculations, but are better represented by tabular equation-of-state models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Uncertainty Quantification for a Multi-phase Carbon Equation of State Model

Many physics models have tunable parameters that are calibrated by matching the model output to experimental or calculated data. However, given that calibration data often contain uncertainty and that different model parameter sets might result in a very similar simulated output for a finite calibration data set, it is advantageous to provide an ensemble of parameter sets that are consistent with the calibration data. Uncertainty quantification (UQ) provides a means to generate such an ensemble in a statistically rigorous fashion. In this work, we perform UQ for a multi-phase equation of state (EOS) model for carbon containing the diamond, graphite, and liquid phases. We use a Bayesian framework for the UQ and introduce a novel strategy for including phase diagram information in the calibration. The method is highly general and accurately reproduces the calibration data without any material-specific prior knowledge of the EOS model parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Imaging dynamic exciton interactions and coupling in transition metal dichalcogenides

Transition metal dichalcogenides (TMDs) are regarded as a possible material platform for quantum information science and related device applications. In TMD monolayers, the dephasing time and inhomogeneity are crucial parameters for any quantum information application. In TMD heterostructures, coupling strength and interlayer exciton lifetimes are also parameters of interest. However, many demonstrations in TMDs can only be realized at specific spots on the sample, presenting a challenge to the scalability of these applications. Here, using multi-dimensional coherent imaging spectroscopy, we shed light on the underlying physics—including dephasing, inhomogeneity, and strain—for a MoSe 2 monolayer and identify both promising and unfavorable areas for quantum information applications. We, furthermore, apply the same technique to a MoSe 2 /WSe 2 heterostructure. Despite the notable presence of strain and dielectric environment changes, coherent and incoherent coupling and interlayer exciton lifetimes are mostly robust across the sample. This uniformity is despite a significantly inhomogeneous interlayer exciton photoluminescence distribution that suggests a bad sample for device applications. This robustness strengthens the case for TMDs as a next-generation material platform in quantum information science and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Temporal structure of blobs in NSTX

The time dependence of the blob pulse shape and the waiting time between blobs was found using data from the gas puff imaging diagnostic in NSTX. The database used was of 103 shots from 2010 as described in a previous paper (Zweben et al., Phys. Plasmas 29, 012505 (2022).]. The blob pulse shape was well fit by an exponential rise and fall where the average rise time was τ r = 9.0 ± 2.7 μs and the average fall time was τ f = 16.6 ± 5.8 μs. The waiting times between blob pulses above a threshold of three times the mean had a broad distribution with an average of τ w = 1.2 ± 0.85 ms over the database. The blob intermittency parameter γ b = τ d /τ w , where the blob pulse duration was τ d = τ r + τ f , ranged from γ b ~1% to 5% for shots in this database and increased almost linearly with the blob fraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Bayesian inferences of electrical current delivered to shocked transmission lines

Small radius Bdot measurements of electrical current delivered by pulsed power devices are routinely compromised by electrode/convolute plasmas endemic to multiterawatt transmission lines. Inferences of delivered current unaffected by these issues have recently been obtained by numerically optimizing consistency between model predictions and local experimental velocimetry data, but these are only unique for shockless velocity profiles. In the paper we describe a more general Bayesian method capable of inferring current despite the presence of shocks. Additionally, we describe uncertainty estimates and use of the technique on experimental data. This technique is the first to provide uncertainty estimates on the full current trace delivered to an inertial confinement fusion target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Virtual tower measurements during the American WAKE ExperimeNt (AWAKEN)

Dual-Doppler lidar measurements were made during the American WAKE ExperimeNt to provide height-resolved measurements of wind speed and direction at multiple locations immediately south of the leading row turbines in the King Plains wind farm in Oklahoma. These so-called virtual tower measurements were performed to characterize the inflow into the wind farm and to assess possible upwind blockage effects due to the collective action of the wind farm. The campaign was conducted from 12 November 2022 to 17 October 2023, during which time 14 unique virtual tower locations were sampled with heights ranging from 240 to 490 m AGL. The wind retrieval algorithm provided estimates of the horizontal winds and their uncertainties with a vertical resolution of about 10 m, while also accounting for the tilt of the lidar platform. The virtual tower results are compared to collocated lidar wind profiling data at the A1 site, which was located roughly 2.4 rotor diameters south of the nearest turbine. The wind speed difference between the wind profiler and the virtual tower was found to be quite sensitive to atmospheric stability and wind direction below 250 m AGL. The largest differences were observed for inflow under stable conditions, where the profiler wind speeds were observed to be about 22% lower than the virtual tower near hub height. These results suggest that there are persistent horizontal gradients in the flow upwind of the wind farm which result in biased estimates using standard ground-based lidar wind profiling methods.

17 WIND ENERGY↗
DOE OSTI2024

Data-driven prediction of scaling and ignition of inertial confinement fusion experiments

Recent advances in inertial confinement fusion (ICF) at the National Ignition Facility (NIF), including ignition and energy gain, are enabled by a close coupling between experiments and high-fidelity simulations. Neither simulations nor experiments can fully constrain the behavior of ICF implosions on their own, meaning pre- and postshot simulation studies must incorporate experimental data to be reliable. Linking past data with simulations to make predictions for upcoming designs and quantifying the uncertainty in those predictions has been an ongoing challenge in ICF research. We have developed a data-driven approach to prediction and uncertainty quantification that combines large ensembles of simulations with Bayesian inference and deep learning. The approach builds a predictive model for the statistical distribution of key performance parameters, which is jointly informed by past experiments and physics simulations. The prediction distribution captures the impact of experimental uncertainty, expert priors, design changes, and shot-to-shot variations. We have used this new capability to predict a 10× increase in ignition probability between Hybrid-E shots driven with 2.05 MJ compared to 1.9 MJ, and validated our predictions against subsequent experiments. We describe our new Bayesian postshot and prediction capabilities, discuss their application to NIF ignition and validate the results, and finally investigate the impact of data sparsity on our prediction results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Active deep kernel learning of molecular properties from structural embeddings

As vast databases of chemical identities become increasingly available, the challenge shifts to how we effectively explore and leverage these resources to study molecular properties. This paper presents an active learning approach for molecular discovery using deep kernel learning (DKL), demonstrated on the QM9 dataset. DKL links structural embeddings directly to properties, creating organized latent spaces that prioritize relevant property information. By iteratively recalculating embedding vectors in alignment with target properties, DKL uncovers concentrated maxima representing key molecular properties and reveals unexplored regions with potential for innovation. This approach underscores DKL’s potential in advancing molecular research and discovery.

Artificial neural networks↗