Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RANDOM PHASE APPROXIMATION”

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.

207 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Anomalous elastic softening in ferroelectric hafnia under pressure

his study employs first-principles density-functional theory (DFT) calculations to explore the elastic and mechanical properties of ferroelectric hafnia (HfO 2 ) in its polar orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase under varying hydrostatic pressure conditions up to 30 GPa. Utilizing a plane-wave basis set and Perdew-Burke-Ernzerhof generalized-gradient approximation for solids in our DFT calculations, we investigate both pure and yttrium-substituted HfO 2 . Our findings reveal an anomalous reduction in the 𝐶 33 component of the elastic tensor with increasing pressure, which becomes significant above 15 GPa and signals a potential pressure-driven structural phase transition at higher pressure. The analysis of atomic displacements under pressure sheds light on the unusual mechanical behavior and phase stability of this material. Additionally, we observe a transition from an indirect band gap to a direct band gap with increasing pressure, which could have significant implications for optical applications. Here, the effects of yttrium substitution on the mechanical and electronic properties are further examined, revealing that yttrium substitution softens the elastic response of this material and reduces the electronic band gap. These results enhance our understanding of elastic and mechanical responses of ferroelectric hafnia and its potential for applications in microelectronics, piezoelectric devices, and nonvolatile ferroelectric random-access memories. Further experimental validation is recommended to confirm our predictions and explore the practical implications of the observed phase transitions and electronic behavior of the ferroelectric hafnia under high-pressure conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

N-particle irreducible actions for stochastic fluids

We construct one- and two-particle irreducible (1PI and 2PI) effective actions for the stochastic fluid dynamics of a conserved density undergoing diffusive motion. We compute the 1PI action in one-loop order and the 2PI action in two-loop approximation. We derive a set of Schwinger-Dyson equations and regularize the resulting equations using Pauli-Villars fields. We numerically solve the Schwinger-Dyson equations for a non-critical fluid. We find that higher-loop effects summed by the Schwinger-Dyson renormalize the non-linear coupling. We also find indications of a diffuson-cascade, the appearance of n-loop correction with smaller and smaller exponential suppression.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Fog Intermittency and Critical Behavior

The intermittency of fog occurrence (the switching between fog and no-fog) is a key stochastic feature that plays a role in its duration and the amount of moisture available. Here, fog intermittency is studied by using the visibility time series collected during the month of July 2022 on Sable Island, Canada. In addition to the visibility, time series of air relative humidity and turbulent kinetic energy, putative variables akin to the formation and breakup conditions of fog, respectively, are also analyzed in the same framework to establish links between fog intermittency and the underlying atmospheric variables. Intermittency in the time series is quantified with their binary telegraph approximations to isolate clustering behavior from amplitude variations. It is shown that relative humidity and turbulent kinetic energy bound many stochastic features of visibility, including its spectral exponent, clustering exponent, and the growth of its block entropy slope. Although not diagnostic, the visibility time series displays features consistent with Pomeau–Manneville Type-III intermittency in its quiescent phase duration PDF scaling (−3/2), power spectrum scaling (−1/2), and signal amplitude PDF scaling (−2). The binary fog time series exhibits properties of self-organized criticality in the relation between its power spectrum scaling and quiescent phase duration distribution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Sonic horizons and causality in phase transition dynamics

A system gradually driven through a symmetry-breaking phase transition is subject to the Kibble-Zurek mechanism (KZM). As a consequence of the critical slowing down, its state cannot follow local equilibrium, and its evolution becomes nonadiabatic near the critical point. In the simplest approximation, that stage can be regarded as an “impulse” where the state of the system remains unchanged. It leads to the correct KZM scaling laws. However, such a “freeze-out” might suggest that the coherence length of the nascent order parameter remains unchanged as the critical region is traversed. By contrast, the original causality-based discussion emphasized the role of the sonic horizon: Domains of the broken symmetry phase can expand with a velocity limited by the speed of the relevant sound. This effect was demonstrated in the quantum Ising chain where the dynamical exponent z = 1 and quasiparticles excited by the transition have a fixed speed of sound. To elucidate the role of the sonic horizon, in this paper we study two systems with z > 1 where the speed of sound is no longer fixed, and the fastest excited quasiparticles set the size of the sonic horizon. Their effective speed decays with increasing transition time. In the extreme case, the dynamical exponent z can diverge such as in the Griffiths region of the random Ising chain where the localization of excited quasiparticles freezes the growth of the correlation range when the critical region is traversed. Of particular interest is an example with z < 1 , the long-range extended Ising chain, where there is no upper limit to the velocity of excited quasiparticles with small momenta. Initially, the power-law tail of the correlation function grows adiabatically, but in the nonadiabatic stage it lags behind the adiabatic evolution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Social network structure and the spread of complex contagions from a population genetics perspective

Ideas, behaviors, and opinions spread through social networks. If the probability of spreading to a new individual is a non-linear function of the fraction of the individuals’ affected neighbors, such a spreading process becomes a “complex contagion”. This non-linearity does not typically appear with physically spreading infections, but instead can emerge when the concept that is spreading is subject to game theoretical considerations (e.g. for choices of strategy or behavior) or psychological effects such as social reinforcement and other forms of peer influence (e.g. for ideas, preferences, or opinions). Here we study how the stochastic dynamics of such complex contagions are affected by the underlying network structure. Motivated by simulations of complex contagions on real social networks, we present a framework for analyzing the statistics of contagions with arbitrary non-linear adoption probabilities based on the mathematical tools of population genetics. The central idea is to use an effective lower-dimensional diffusion process to approximate the statistics of the contagion. This leads to a tradeoff between the effects of ”selection” (microscopic tendencies for an idea to spread or die out), random drift, and network structure. Our framework illustrates intuitively several key properties of complex contagions: stronger community structure and network sparsity can significantly enhance the spread, while broad degree distributions dampen the effect of selection compared to random drift. Finally, we show that some structural features can exhibit critical values that demarcate regimes where global contagions become possible for networks of arbitrary size. Our results draw parallels between the competition of genes in a population and memes in a world of minds and ideas. Our tools provide insight into the spread of information, behaviors, and ideas via social influence, and highlight the role of macroscopic network structure in determining their fate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Understanding vortex dynamics in CaK(Fe,Ni) 4 As 4 and Ba(Fe,Co) 2 As 2 single crystals under the influence of random point disorder

We report on the influence of doping on vortex dynamics in 3 MeV proton-irradiated single crystals of CaK(Fe 1–x Ni x ) 4 As 4 (1144, x = 0.015, 0.025, and 0.03) and Ba(Fe 1–x Co x ) 2 As 2 (x = 0.04, 0.062, 0.066 and 0.074). Non-irradiated crystals of the 1144 system display superconducting critical temperatures ranging from 31 K for x = 0.015–20.5 K, as doping increases to 0.03. On the other hand, pristine crystals of the 122 system show T c values between 14.6 and 23.6 K, with the maximum T c occurring at intermediate doping levels. The fluence was set at 3 × 10 16 p cm –2 , resulting in a decrease in the T c by around 1.5 K for all samples and significantly affecting the vortex dynamics by reducing the flux creep relaxation compared to previously reported values for unirradiated crystals. Parameters such as vortex pinning energy U 0 and the glassy exponent μ dependencies on doping and magnetic field strength are identified. For the 1144 system, U 0 reaches values approaching 500 K for small fields in samples with T c = 29.3 K (x = 0.015), systematically decreasing to around 200 K as T c falls below 20 K. Furthermore, U 0 decreases as the field increases to 3 T for the same sample, varying from approximately 250 K to 100 K as T c decreases. These changes are typically accompanied by modifications in μ, gradually increasing from values around 1 towards 1.5, corresponding to small bundle relaxation in the collective creep theory. Despite differences in the substitutional disorder and magnetic phase diagram with respect to the 1144 system, the results for 122 single crystals follow a similar tendency in which U 0 usually reduces and μ increase rise as the applied magnetic field is increased. Due to moderate U 0 in these systems (few hundreds of kelvins), the resulting decay of persistent current at liquid helium temperatures is primarily determined by a balance between U 0 and bundle size contribution. Furthermore, these findings provide valuable insights for potential applications of these systems, particularly in the context of intrinsic superconducting parameters and the resulting pinning landscape.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Correcting for filter-based aerosol light absorption biases at the Atmospheric Radiation Measurement program's Southern Great Plains site using photoacoustic measurements and machine learning

Abstract. Measurement of light absorption of solar radiation by aerosols is vital for assessing direct aerosol radiative forcing, which affects local and global climate. Low-cost and easy-to-operate filter-based instruments, such as the Particle Soot Absorption Photometer (PSAP), that collect aerosols on a filter and measure light attenuation through the filter are widely used to infer aerosol light absorption. However, filter-based absorption measurements are subject to artifacts that are difficult to quantify. These artifacts are associated with the presence of the filter medium and the complex interactions between the filter fibers and accumulated aerosols. Various correction algorithms have been introduced to correct for the filter-based absorption coefficient measurements toward predicting the particle-phase absorption coefficient (Babs). However, the inability of these algorithms to incorporate into their formulations the complex matrix of influencing parameters such as particle asymmetry parameter, particle size, and particle penetration depth results in prediction of particle-phase absorption coefficients with relatively low accuracy. The analytical forms of corrections also suffer from a lack of universal applicability: different corrections are required for rural and urban sites across the world. In this study, we analyzed and compared 3 months of high-time-resolution ambient aerosol absorption data collected synchronously using a three-wavelength photoacoustic absorption spectrometer (PASS) and PSAP. Both instruments were operated on the same sampling inlet at the Department of Energy's Atmospheric Radiation Measurement program's Southern Great Plains (SGP) user facility in Oklahoma. We implemented the two most commonly used analytical correction algorithms, namely, Virkkula (2010) and the average of Virkkula (2010) and Ogren (2010)–Bond et al. (1999) as well as a random forest regression (RFR) machine learning algorithm to predict Babs values from the PSAP's filter-based measurements. The predicted Babs was compared against the reference Babs measured by the PASS. The RFR algorithm performed the best by yielding the lowest root mean square error of prediction. The algorithm was trained using input datasets from the PSAP (transmission and uncorrected absorption coefficient), a co-located nephelometer (scattering coefficients), and the Aerosol Chemical Speciation Monitor (mass concentration of non-refractory aerosol particles). A revised form of the Virkkula (2010) algorithm suitable for the SGP site has been proposed; however, its performance yields approximately 2-fold errors when compared to the RFR algorithm. To generalize the accuracy and applicability of our proposed RFR algorithm, we trained and tested it on a dataset of laboratory measurements of combustion aerosols. Input variables to the algorithm included the aerosol number size distribution from the Scanning Mobility Particle Sizer, absorption coefficients from the filter-based Tricolor Absorption Photometer, and scattering coefficients from a multiwavelength nephelometer. The RFR algorithm predicted Babs values within 5 % of the reference Babs measured by the multiwavelength PASS during the laboratory experiments. Thus, we show that machine learning approaches offer a promising path to correct for biases in long-term filter-based absorption datasets and accurately quantify their variability and trends needed for robust radiative forcing determination.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Experimental and computational studies on high-entropy carbide MoNbTaVWC 5 under high pressures

High-entropy carbide, MoNbTaVWC 5 , was synthesized from oxide precursors of the constituent metals, mixed with graphite powder in a microwave-generated hydrogen plasma at 26.66 kPa and 2100 °C. Ambient x-ray diffraction analysis confirms the full conversion of oxide precursors into a single-phase, face-centered cubic structure with a lattice parameter a = 4.3309 Å. Nanoindentation measured a hardness of 24.5 ± 1.3 GPa and an elastic modulus of 386 ± 22 GPa. The synthesized sample, mixed with a copper pressure marker, was studied by the radial x-ray diffraction technique with beryllium gasketing in a diamond anvil cell up to 70 GPa. The experimentally measured pressure–volume curve and shear strength were compared with theoretical predictions using the special quasi-random structure technique and density functional theory. MoNbTaVWC 5 achieved a 12% volume compression at 70 GPa and exhibited a high shear strength of 6.6 GPa. The present study demonstrates that the high-entropy carbide MoNbTaVWC 5 exhibits exceptional incompressibility and high strength under extreme conditions.

36 MATERIALS SCIENCE↗