Search NASA⌕ Search

SEARCH · Search NASA

Results for “random fields”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multilevel Hierarchical Decomposition of Finite Element White Noise with Application to Multilevel Markov Chain Monte Carlo

In this work we develop a new hierarchical multilevel approach to generate Gaussian random field realizations in an algorithmically scalable manner that is well suited to incorporating into multilevel Markov chain Monte Carlo (MCMC) algorithms. This approach builds off of other partial differential equation (PDE) approaches for generating Gaussian random field realizations; in particular, a single field realization may be formed by solving a reaction-diffusion PDE with a spatial white noise source function as the right-hand side. While these approaches have been explored to accelerate forward uncertainty quantification tasks, e.g., multilevel Monte Carlo, the previous constructions are not directly applicable to multilevel MCMC frameworks which build fine-scale random fields in a hierarchical fashion from coarse-scale random fields. Our new hierarchical multilevel method relies on a hierarchical decomposition of the white noise source function in $L^2$ which allows us to form Gaussian random field realizations across multiple levels of discretization in a way that fits into multilevel MCMC algorithmic frameworks. After presenting our main theoretical results and numerical scaling results to showcase the utility of this new hierarchical PDE method for generating Gaussian random field realizations, this method is tested on a four-level MCMC algorithm to explore its feasibility.

algebraic multigrid↗
DOE OSTI2017

Random electric field instabilities of relaxor ferroelectrics

Relaxor ferroelectrics are complex oxide materials which are rather unique to study the effects of compositional disorder on phase transitions. Here, we study the effects of quenched cubic random electric fields on the lattice instabilities that lead to a ferroelectric transition and show that, within a microscopic model and a statistical mechanical solution, even weak compositional disorder can prohibit the development of long-range order and that a random field state with anisotropic and power-law correlations of polarization emerges from the combined effect of their characteristic dipole forces and their inherent charge disorder. As a result, we compare and reproduce several key experimental observations in the well-studied relaxor PbMg 1/3 Nb 2/3 O 3 –PbTiO 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2017

Charged Particle Diffusion in Isotropic Random Magnetic Fields

The investigation of the diffusive transport of charged particles in a turbulent magnetic field remains a subject of considerable interest. Research has most frequently concentrated on determining the diffusion coefficient in the presence of a mean magnetic field. Here we consider the diffusion of charged particles in fully three-dimensional isotropic turbulent magnetic fields with no mean field, which may be pertinent to many astrophysical situations. We identify different ranges of particle energy depending upon the ratio of Larmor radius to the characteristic outer length scale of turbulence. Two different theoretical models are proposed to calculate the diffusion coefficient, each applicable to a distinct range of particle energies. The theoretical results are compared to those from computer simulations, showing good agreement.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Effects of the Random Crystal Field on the Phase Diagrams and the Magnetic Properties of a Spin-1 Blume-Capel Nanoisland

The phase diagrams and the magnetic properties of a spin-1 Blume-Capel nanoisland are investigated by the use of the effective-field theory based on the probability distribution technique. The influence of the random crystal field and the system parameters on the critical behaviors of the nanoisland are examined. The results show that the system can exhibit a first-order phase transition, a tricritical point, and a re-entrant or even a double re-entrant phenomenon.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Random magnetic field and the Dirac Fermi surface

In this paper, we study a single two-dimensional Dirac fermion at finite density, subject to a quenched random magnetic field. At low energies and sufficiently weak disorder, the theory maps onto an infinite collection of 1D chiral fermions (associated to each point on the Fermi surface) coupled by a random vector potential. This low-energy theory exhibits an exactly solvable random fixed line, along which we directly compute various disorder-averaged observables without the need for the usual replica, supersymmetry, or Keldysh techniques. We find the longitudinal dc conductivity in the collisionless $\hbar$ω/k B T→∞ limit to be nonuniversal and to vary continuously along the fixed line.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Random Crystal Field in a Mixed Spin S = 1/2 and S = 3/2 Ising Model by Renormalization Group Theory

In this work, we study the Ising model with mixed spins S = 1/2 and S = 3/2 on the hypercubic lattice and the random crystal field at the sublattice with S = 3/2 described by a two peaks law. To achieve this, we use an approximation of position space renormalization group (PSRG) namely Migdal-Kadanoff in which we use both the free energy derivative and the flow in the parameter space of the Hamiltonian. For all values of the random probability, the critical behavior is determined via the critical exponents at the second-order fixed points while at low temperatures; the discontinuities of the of the free energy derivative provide the positions of the first-order transitions. The introduction of a minimal amount of disorder causes a change in the phase diagram showing the relevance of disorder for d = 2 and d = 3. The second-order transition remains always with the same critical exponents as those of the pure model, and a new first-order transition appears at very low temperature. Also, a comparison with other similar works is given.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A modular methodology for time-domain stochastic seismic wave propagation

Presented here is a modular methodology for time-domain stochastic seismic wave propagation analysis. Presented methodology is designed to analyse uncertain seismic motions as an input, propagating through uncertain material. Traditional approach for uncertain wave propagation relies on models that include deep bedrock, local soil site, and their random process and random field information. Such models can become quite large and computationally intractable. The modular approach proposed herein features two step approach that allows separate consideration of the deep bedrock and local site along with corresponding random field information. In this work, the first step considers an auxiliary stochastic motions problem in the bedrock. Stochastic local site response can then be simulated in a reduced domain within certain depth from the surface. Application of uncertain seismic motions at depth, for local uncertain site response is done using stochastic effective forces developed through the Domain Reduction Method. By using Hermite polynomial chaos expansion to represent the non-Gaussian random field of material parameters and non-stationary random process of seismic motion, the proposed modular methodology is formulated using intrusive stochastic Galerkin approach, as seen in the Stochastic Elastic–Plastic Finite Element Method (SEPFEM). Developed modular methodology is illustrated using a 1-D stochastic seismic wave propagation analysis with three cases, and simulation results are also verified with results from conventional approach.

58 GEOSCIENCES↗
DOE OSTI2022

Map-based cosmology inference with lognormal cosmic shear maps

ABSTRACT Most cosmic shear analyses to date have relied on summary statistics (e.g. ξ+ and ξ−). These types of analyses are necessarily suboptimal, as the use of summary statistics is lossy. In this paper, we forward-model the convergence field of the Universe as a lognormal random field conditioned on the observed shear data. This new map-based inference framework enables us to recover the joint posterior of the cosmological parameters and the convergence field of the Universe. Our analysis properly accounts for the covariance in the mass maps across tomographic bins, which significantly improves the fidelity of the maps relative to single-bin reconstructions. We verify that applying our inference pipeline to Gaussian random fields recovers posteriors that are in excellent agreement with their analytical counterparts. At the resolution of our maps – and to the extent that the convergence field can be described by the lognormal model – our map posteriors allow us to reconstruct all summary statistics (including non-Gaussian statistics). We forecast that a map-based inference analysis of LSST-Y10 data can improve cosmological constraints in the σ8–Ωm plane by $\approx\!{30}{{\ \rm per\ cent}}$ relative to the currently standard cosmic shear analysis. This improvement happens almost entirely along the $S_8=\sigma _8\Omega _{\rm m}^{1/2}$ directions, meaning map-based inference fails to significantly improve constraints on S8.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Study of the Vertical Magnetic Field in Face-on Galaxies Using Faraday Tomography

Faraday tomography allows astronomers to probe the distribution of the magnetic field along the line of sight (LOS), but that can be achieved only after the Faraday spectrum is interpreted. However, the interpretation is not straightforward, mainly because the Faraday spectrum is complicated due to a turbulent magnetic field; it ruins the one-to-one relation between the Faraday depth and the physical depth, and appears as many small-scale features in the Faraday spectrum. In this paper, by employing “simple toy models” for the magnetic field, we describe numerically as well as analytically the characteristic properties of the Faraday spectrum. We show that the Faraday spectrum along “multiple LOSs” can be used to extract the global properties of the magnetic field. Specifically, considering face-on spiral galaxies and modeling turbulent magnetic field as a random field with a single coherence length, we numerically calculate the Faraday spectrum along a number of LOSs and its shape-characterizing parameters, that is, the moments. When multiple LOSs cover a region of ≳(10 coherence length){sup 2}, the shape of the Faraday spectrum becomes smooth and the shape-characterizing parameters are well specified. With the Faraday spectrum constructed as a sum of Gaussian functions with different means and variances, we analytically show that the parameters are expressed in terms of the regular and turbulent components of the LOS magnetic field and the coherence length. We also consider the turbulent magnetic field modeled with a power-law spectrum, and study how the magnetic field is revealed in the Faraday spectrum. Our work suggests a way to obtain information on the magnetic field from a Faraday tomography study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Do Epistemic Uncertainties Allow for Replacing Microstructural Experiments with Reconstruction Algorithms?

A novel problem in computational materials modeling is addressed: “Are the computational microstructure reconstruction techniques reliable enough to replace experiments?” Here, “reliable” computations are associated with producing “expected” reconstructions that are adequately close to the experimental data. The output of computational techniques can deviate from the experimental measurements because of the epistemic uncertainties in the algorithms. In this work, an analytical formulation for quantification of epistemic uncertainties in a microstructure reconstruction algorithm based on Markov random field is presented. The method is used to predict the large-scale spatial distribution of a microstructure given an experimental input measured over a small spatial domain. However, small variations are observed on the microstructural features of the synthesized samples due to the Markov random field algorithm. The proposed analytical technique aims to quantify these uncertainties and estimate their propagation to the macroscale material properties to provide a significant understanding on how reliable it is to replace the experiments with the Markov random field model to predict microstructural maps over large spatial domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Random strain induced correlations in materials with intertwined nematic and magnetic orders

Electronic nematicity is rarely observed as an isolated instability of a correlated electron system. Instead, in iron pnictides and in certain cuprates and heavy-fermion materials, nematicity is intertwined with an underlying spin-stripe or charge-stripe state. As a result, random strain, ubiquitous in any real crystal, creates both random-field disorder for the nematic degrees of freedom and random-bond disorder for the spin or charge ones. Here, we put forward an Ashkin-Teller model with random Baxter fields to capture the dual role of random strain in nematic systems for which nematicity is a composite order arising from a stripe state. Using Monte Carlo to simulate this random Baxter-field model, we find not only the expected break-up of the system into nematic domains, but also the emergence of nontrivial disorder-promoted magnetic correlations. Such correlations enhance and tie up the fluctuations associated with the two degenerate magnetic stripe states from which nematicity arises, leaving characteristic signatures in the spatial profile of the magnetic domains, in the configurational space of the spin variables, and in the magnetic noise spectrum. We discuss possible experimental manifestations of these effects in iron-pnictide superconductors. Furthermore, our work establishes the random Baxter-field model as a more complete alternative to the random-field Ising model to describe complex electronic nematic phenomena in the presence of disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Optimal design of acoustic metamaterial cloaks under uncertainty

In this work, we consider the problem of optimal design of an acoustic cloak under uncertainty and develop scalable approximation and optimization methods to solve this problem. The design variable is taken as an infinite-dimensional spatially-varying field that represents the material property, while an additive infinite-dimensional random field represents the variability of the material property or the manufacturing error. Discretization of this optimal design problem results in high-dimensional design variables and uncertain parameters. To solve this problem, we develop a computational approach based on a Taylor approximation and an approximate Newton method for optimization, which is based on a Hessian derived at the mean of the random field. We show our approach is scalable with respect to the dimension of both the design variables and uncertain parameters, in the sense that the necessary number of acoustic wave propagations is essentially independent of these dimensions, for numerical experiments with up to one million design variables and half a million uncertain parameters. Additionally, we demonstrate that, using our computational approach, an optimal design of the acoustic cloak that is robust to material uncertainty is achieved in a tractable manner. The optimal design under uncertainty problem is posed and solved for the classical circular obstacle surrounded by a ring-shaped cloaking region, subjected to both a single-direction single-frequency incident wave and multiple-direction multiple-frequency incident waves. Finally, we apply the method to a deterministic large-scale optimal cloaking problem with complex geometry, to demonstrate that the approximate Newton method’s Hessian computation is viable for large, complex problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Artificial generation of representative single Li-ion electrode particle architectures from microscopy data

Accurately capturing the architecture of single lithium-ion electrode particles is necessary for understanding their performance limitations and degradation mechanisms through multi-physics modeling. Information is drawn from multimodal microscopy techniques to artificially generate LiNi 0.5 Mn 0.3 Co 0.2 O 2 particles with full sub-particle grain detail. Statistical representations of particle architectures are derived from X-ray nano-computed tomography data supporting an ‘outer shell’ model, and sub-particle grain representations are derived from focused-ion beam electron backscatter diffraction data supporting a ‘grain’ model. A random field model used to characterize and generate the outer shells, and a random tessellation model used to characterize and generate grain architectures, are combined to form a multi-scale model for the generation of virtual electrode particles with full-grain detail. This work demonstrates the possibility of generating representative single electrode particle architectures for modeling and characterization that can guide synthesis approaches of particle architectures with enhanced performance.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Optimization, Inference and Learning for Advanced Networks

This piece of software aims at providing an implementation of methods to perform efficient optimization, probabilistic inference and model selection and learning for various type of networked systems. The software contains a suite of methods and functions that: 1) perform probabilistic inference and random sampling on networked models 2) select or learn networked models from data 3) solve optimization tasks with a networked structure Network models that are supported by these features include, but are not limited to, factor graphs and Markov random fields such as Ising models, Potts models, p-spins models, Gaussian random fields and restricted Boltzmann machine.

Hijazi, Hassan↗
DOE OSTI2017

Disordered Route to the Coulomb Quantum Spin Liquid: Random Transverse Fields on Spin Ice in Pr 2 Zr 2 O 7

Inelastic neutron scattering reveals a broad continuum of excitations in Pr 2 Zr 2 O 7 , the temperature and magnetic field dependence of which indicate a continuous distribution of quenched transverse fields ( Δ ) acting on the non-Kramers Pr 3 + crystal field ground state doublets. Spin-ice correlations are apparent within 0.2 meV of the Zeeman energy. In a random phase approximation an excellent account of the data is provided and contains a transverse field distribution ρ ( Δ ) ∝ ( Δ 2 + Γ 2 ) - 1 , where Γ = 0.27 ( 1 ) meV . Established during high temperature synthesis due to an underlying structural instability, it appears disorder in Pr 2 Zr 2 O 7 actually induces a quantum spin liquid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

End-to-end learning of multiple sequence alignments with differentiable Smith–Waterman

Abstract Motivation Multiple sequence alignments (MSAs) of homologous sequences contain information on structural and functional constraints and their evolutionary histories. Despite their importance for many downstream tasks, such as structure prediction, MSA generation is often treated as a separate pre-processing step, without any guidance from the application it will be used for. Results Here, we implement a smooth and differentiable version of the Smith–Waterman pairwise alignment algorithm that enables jointly learning an MSA and a downstream machine learning system in an end-to-end fashion. To demonstrate its utility, we introduce SMURF (Smooth Markov Unaligned Random Field), a new method that jointly learns an alignment and the parameters of a Markov Random Field for unsupervised contact prediction. We find that SMURF learns MSAs that mildly improve contact prediction on a diverse set of protein and RNA families. As a proof of concept, we demonstrate that by connecting our differentiable alignment module to AlphaFold2 and maximizing predicted confidence, we can learn MSAs that improve structure predictions over the initial MSAs. Interestingly, the alignments that improve AlphaFold predictions are self-inconsistent and can be viewed as adversarial. This work highlights the potential of differentiable dynamic programming to improve neural network pipelines that rely on an alignment and the potential dangers of optimizing predictions of protein sequences with methods that are not fully understood. Availability and implementation Our code and examples are available at: https://github.com/spetti/SMURF. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Summability of Connected Correlation Functions of Coupled Lattice Fields

We consider two nonindependent random fields ψ and ϕ defined on a countable set Z. For instance, Z=ℤ{sup d} or Z=ℤ{sup d}×I, where I denotes a finite set of possible “internal degrees of freedom” such as spin. We prove that, if the cumulants of ψ and ϕ enjoy a certain decay property, then all joint cumulants between ψ and ϕ are ℓ{sub 2}-summable in the precise sense described in the text. The decay assumption for the cumulants of ψ and ϕ is a restricted ℓ{sub 1} summability condition called ℓ{sub 1}-clustering property. One immediate application of the results is given by a stochastic process ψ{sub t}(x) whose state is ℓ{sub 1}-clustering at any time t: then the above estimates can be applied with ψ=ψ{sub t} and ϕ=ψ{sub 0} and we obtain uniform in t estimates for the summability of time-correlations of the field. The above clustering assumption is obviously satisfied by any ℓ{sub 1}-clustering stationary state of the process, and our original motivation for the control of the summability of time-correlations comes from a quest for a rigorous control of the Green–Kubo correlation function in such a system. A key role in the proof is played by the properties of non-Gaussian Wick polynomials and their connection to cumulants.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A scalable variational method for estimating the latent infection-rate field of an outbreak

In this paper, we explore whether the infection-rate of a disease can serve as a robust monitoring variable in epidemiological surveillance algorithms. The infection-rate is dependent on population mixing patterns that do not vary erratically day-to-day; in contrast, daily case-counts used in contemporary surveillance algorithms are corrupted by reporting errors. The technical challenge lies in estimating the latent infection-rate from case-counts. Here we devise a Bayesian method to estimate the infection-rate across multiple adjoining areal units, and then use it, via an anomaly detector, to discern a change in epidemiological dynamics. We extend an existing model for estimating the infection-rate in an areal unit by incorporating a Markov random field model, so that we may estimate infection-rates across multiple areal units, while preserving spatial correlations observed in the epidemiological dynamics. To carry out the high-dimensional Bayesian inverse problem, we develop an implementation of mean-field variational inference specific to the infection model and integrate it with the random field model to incorporate correlations across counties. The method is tested on estimating the COVID-19 infection-rates across all 33 counties in New Mexico using data from the summer of 2020, and then employing them to detect the arrival of the Fall 2020 COVID-19 wave. We perform the detection using a temporal algorithm that is applied county-by-county. We also show how the infection-rate field can be used to cluster counties with similar epidemiological dynamics.

60 APPLIED LIFE SCIENCES↗