Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “random walks”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Distinct universality classes of diffusive transport from full counting statistics

The hydrodynamic transport of local conserved densities furnishes an effective coarse-grained description of the dynamics of a many-body quantum system. However, the full quantum dynamics contains much more structure beyond the simplified hydrodynamic description. Here we show that systems with the same hydrodynamics can nevertheless belong to distinct dynamical universality classes, as revealed by new classes of experimental observables accessible in synthetic quantum systems, which can, for instance, measure simultaneous site-resolved snapshots of all of the particles in a system. Specifically, we study the full counting statistics of spin transport, whose first moment is related to linear-response transport, but the higher moments go beyond. We present an analytic theory of the full counting statistics of spin transport in various integrable and nonintegrable anisotropic one-dimensional spin models, including the XXZ spin chain. We find that spin transport, while diffusive on average, is governed by a distinct non-Gaussian dynamical universality class in the models considered. We consider a setup in which the left and right half of the chain are initially created at different magnetization densities, and consider the probability distribution of the magnetization transferred between the two half-chains. We derive a closed-form expression for the probability distribution of the magnetization transfer, in terms of random walks on the half-line. We show that this distribution strongly violates the large-deviation form expected for diffusive chaotic systems, and explain the physical origin of this violation. Here, we discuss the crossovers that occur as the initial state is brought closer to global equilibrium. Our predictions can directly be tested in experiments using quantum gas microscopes or superconducting qubit arrays.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Primary fission fragment mass yields across the chart of nuclides

In this work, we have calculated a complete set of primary fission fragment mass yields, Y(A), for heavy nuclei across the chart of nuclides, including those of particular relevance to the rapid neutron capture process (r process) of nucleosynthesis. We assume that the nuclear shape dynamics are strongly damped, which allows for a description of the fission process via Brownian shape motion across nuclear potential-energy surfaces. The macroscopic energy of the potential was obtained with the Finite-Range Liquid-Drop Model (FRLDM), while the microscopic terms were extracted from the single-particle level spectra in the fissioning system by the Strutinsky procedure for the shell energies and the BCS treatment for the pairing energies. For each nucleus considered, the fission fragment mass yield, Y(A), is obtained from 50 000 to 500 000 random walks on the appropriate potential-energy surface. The full mass and charge yield, Y(Z,A), is then calculated by invoking the Wahl systematics. With this method, we have calculated a comprehensive set of fission-fragment yields from over 3800 nuclides bounded by 80 ≤ Z ≤ 130 and A ≤ 330; these yields are provided as an ASCII formatted database in the Supplemental Material. We compare our yields to known data and discuss general trends that emerge in low-energy fission yields across the chart of nuclides.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Formation of transfermium elements in reactions with Pb 208

Within the Langevin framework, we investigate the dynamics of the fusion process for production of transfermium elements in reactions of Ca 48 , Ti 50 , Cr 54 , and Fe 58 with Pb 208 . After the reacting nuclei have made contact, the early dynamical stage is dominated by the dissipation of the initial radial kinetic energy, while the subsequent shape evolution is diffusive. The probability for surmounting the inner barrier and forming a compound system is obtained by simulating the evolution as a Metropolis random walk in a five-dimensional potential-energy landscape. Good agreement with the available data is obtained, especially for the maximal formation probability. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Optimizing temperature distributions for training neural quantum states using parallel tempering

Parametrized artificial neural networks (ANNs) can be very expressive ansatzes for variational algorithms, reaching state-of-the-art energies on many quantum many-body Hamiltonians. Nevertheless, the training of the ANN can be slow and stymied by the presence of local minima in the parameter landscape. One approach to mitigate this issue is to use parallel tempering methods, and in this work, we focus on the role played by the temperature distribution of the parallel tempering replicas. Using an adaptive method that adjusts the temperatures in order to equate the exchange probability between neighboring replicas, we show that this temperature optimization can significantly increase the success rate of the variational algorithm with negligible computational cost by eliminating bottlenecks in the replicas' random walk. Furthermore, we demonstrate this using two different neural networks, a restricted Boltzmann machine and a feedforward network, which we use to study a toy problem based on a permutation invariant Hamiltonian with a pernicious local minimum and the 𝐽 1 −𝐽 2 model on a rectangular lattice.

Neural network simulations↗
DOE OSTI2020

Particlelike Phonon Propagation Dominates Ultralow Lattice Thermal Conductivity in Crystalline Tl 3 VSe 4

We explore the microscopic mechanisms of ultralow lattice thermal conductivity ($\kappa_{l}$) in Tl$_{3}$VSe$_{4}$~by combining a first-principles density-functional theory (DFT) based framework of anharmonic lattice dynamics with the Peierls-Boltzmann transport equation (PBTE) for phonons. We include contributions of the three- and four-phonon scattering processes to the phonon lifetimes as well as the temperature-dependent anharmonic renormalization of phonon energies arising from an unusually strong quartic anharmonicity in Tl$_{3}$VSe$_{4}$. In contrast to a recent report by Mukhopadhyay~\etal [\textcolor{blue}{Science 360, 1455 (2018)}] which suggested that a significant contribution to $\kappa_{l}$ arises from random walks among uncorrelated oscillators, we show that particle-like propagation of phonon excitations can successfully explain the experimentally observed ultralow $\kappa_{l}$. Our results are also supported by explicit calculations of the off-diagonal terms of the heat-current operator, which are found to be small and indicate that wave-like tunneling of heat-carrying vibrations is of minor importance. Our results (i) resolve the discrepancy between the theoretical and experimental $\kappa_{l}$, (ii) offer new insights into the minimum $\kappa_{l}$ achievable in \TlVSe, and (iii) highlight the importance of high-order anharmonicity in low-$\kappa_{l}$ systems. The methodology demonstrated here may be used to resolve the discrepancies between the experimentally measured and the theoretically calculated $\kappa_{l}$ in skutterides and perovskites, as well as to understand the glasslike $\kappa_{l}$ in complex crystals with strong anharmonicity, leading towards the goal of rational design of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Operator Lévy Flight: Light Cones in Chaotic Long-Range Interacting Systems

We argue that chaotic power-law interacting systems have emergent limits on information propagation, analogous to relativistic light cones, which depend on the spatial dimension d and the exponent α governing the decay of interactions. Using the dephasing nature of quantum chaos, we map the problem to a stochastic model with a known phase diagram. A linear light cone results for α ≥ $\textit{d}$ + 1/2. We also provide a Lévy flight (long-range random walk) interpretation of the results and show consistent numerical data for 1D long-range spin models with 200 sites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Shear strain alters the structure and migration mechanism of self-interstitial atoms in copper

Here, we use atomistic modeling to show that externally applied shear strain causes the lowest energy self-interstitial atom (SIA) structure in copper (Cu) to change from a <100>-type dumbbell to a <110>-type dumbbell. Concurrently, SIA migration switches from the 3-D random walk characteristic of <100>-type dumbbells to a 1-D mechanism analogous to that of crowdion SIAs. Furthermore, the relative energies of these two dumbbell structures as a function of strain are well predicted using elastic dipole tensors computed at zero strain, indicating that examination of these tensors may be used to assess the likelihood of strain-induced SIA structure transitions in other materials. Changes in lowest energy SIA structures and associated migration mechanisms stand to impact predictions of SIA behavior in irradiated solids.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microstructural and rheological training and memory of nanocolloidal soft glasses under cyclic shear

An intrinsic feature of disordered and out-of-equilibrium materials, such as glasses, is the dependence of their properties on their history. An important example is rheological memory, in which disordered solids obtain properties based on their deformation history. Here, in this study, we employ x-ray photon correlation spectroscopy with in situ rheometry to characterize memory formation in a nanocolloidal soft glass due to cyclic shear. During a cycle, particles undergo irreversible displacements composed of a combination of shear-induced diffusion and heterogeneous, residual strain fields. At lower shear amplitudes, the displacements resemble a random walk in which the directions in each cycle are independent of those in preceding cycles, while at high amplitude, the irreversible displacements in consecutive cycles become correlated. The magnitudes of the displacements decrease with each cycle before reaching a steady state where the microstructure has been trained to achieve enhanced reversibility even at shear amplitudes well above yielding and despite the presence of thermal fluctuations. At amplitudes below and near yielding, these decreases are monotonic, while well above yielding, they are nonmonotonic, suggesting evidence of shear banding. Accompanying this microstructural training are corresponding decreases in the dissipation during each cycle and the magnitude of the residual stress toward steady-state values. Memory of the training is revealed by measurements in which the amplitude of the shear is changed after steady state is reached. The magnitude of the particle displacements, as well as the dissipation and the change in residual stress, vary nonmonotonically with the new shear amplitude, having minima near the training amplitude, thereby revealing correlated microscopic and macroscopic signatures of memory.

Chen, Yihao [Johns Hopkins Univ., Baltimore, MD (U↗
DOE OSTI2026

Janus skyrmion: Interfacial quasiparticle with two-faced helicity

Janus particles are functional particles with at least two surfaces showing asymmetric properties. Here, we show at the interface between two magnetic regions with different antisymmetric exchange interactions, an alternative species of two-dimensional topological quasiparticles can emerge, in which different helicity structures can coexist. We name such an interfacial quasiparticle a “Janus skyrmion,” in analogy to the Janus particle. As the Janus skyrmion shows helicity asymmetry, its size could vary with both the in-plane and out-of-plane magnetic fields. A vertical spin current could drive the Janus skyrmion into one-dimensional motion along the interface without showing the skyrmion Hall effect, at a speed which depends on both the in-plane spin-polarization direction and current density. Thermal fluctuations could also lead to one-dimensional random walk of a Brownian Janus skyrmion. This work uncovers unique dynamics intrinsic to interfacial quasiparticles with exotic helicity, which may be realized in interface-engineered magnetic layers.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Fokker-Planck Equation Governing the Distribution of Walkers in Auxiliary-Field Quantum Monte Carlo

Auxiliary-field quantum Monte Carlo (AFQMC) is typically formulated as an open-ended random walk in an overcomplete space of Slater determinants, implemented through a Langevin equation. However, the explicit form of the underlying Fokker-Planck equation governing the walker population distribution has remained unknown. Here, in this Letter, we derive the Fokker-Planck equation for AFQMC and propose a novel numerical scheme to solve it. The solution of the Fokker-Planck equation reveals the wave function actually sampled by the AFQMC algorithm. Interestingly, we find that even when the exact ground state is used as a guiding wave function in constrained path AFQMC, contrary to the common assumption, the wave function sampled by AFQMC is not exact. Beyond clarifying several fundamental aspects of AFQMC, the availability of a Fokker-Planck equation formulation opens new avenues for systematically improving its accuracy, which we outline in this Letter.

Monte Carlo methods↗
DOE OSTI2020

Sub-Degree-Per-Hour MEMS Gyroscope for Measurement While Drilling at 300°C

The orientation module of MWD tool offers the critical drill bit orientation information, including azimuth, inclination and toolface in order to control the path of wellbore. Although high end navigation grade gyroscope can meet the accuracy requirement for azimuth finding, the requirements of cost, size and reliability under harsh environment have largely limited the deployment of gyroscopes in MWD tools. To overcome these constraints, GE Research has developed a low cost, MEMS based Multi-Ring Gyroscope (MRG) capable of azimuth-seeking in MWD applications at 300°C. The MRG prototype has been demonstrated to achieve angular random walk (ARW) better than 0.003 deg/rt(hr) and bias instability of 0.01 deg/hr, capable of meeting azimuth finding accuracy better than 0.25 deg. It has further been tested to remain fully operational at 300 °C.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Risk-Aware Measurement Synchronization and Recovery for DSSE With Heterogeneous Data Sources

Power distribution systems are increasingly integrating heterogeneous sensors with varying data reporting rates and types, which pose challenges to achieving observability at the desired temporal resolution of distribution system state estimation (DSSE). Multisensor failures caused by extreme events exacerbate these issues, introducing substantial uncertainties into DSSE. This article proposes a novel solution to these challenges by ensuring high-resolution system observability despite heterogeneous data sources and multisensor failures. First, a deep learning architecture combining long short-term memory (LSTM) and graph convolutional network (GCN) is employed to synchronize meters with different reporting rates, aiming to achieve system observability. A random-walk-model-based approach is introduced to generate pseudo-measurements while properly characterizing their uncertainties under multisensor failures. Finally, a disaster-risk-informed observability metric (RiOM) is defined to quantify the uncertainty associated with state estimation results. The proposed framework offers deeper insights into the system observability on the fly compared with conventional analysis. The effectiveness of the framework is demonstrated on an IEEE standard test case and a large-scale real-world distribution feeder in mid-Minnesota in the U.S.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Wave Dark Matter

Here we review the physics and phenomenology of wave dark matter: a bosonic dark matter candidate lighter than about 30 eV. Such particles have a de Broglie wavelength exceeding the average interparticle separation in a galaxy like the Milky Way and are, thus, well described as a set of classical waves. We outline the particle physics motivations for such particles, including the quantum chromodynamics axion as well as ultralight axion-like particles such as fuzzy dark matter. The wave nature of the dark matter implies a rich phenomenology: Wave interference gives rise to order unity density fluctuations on de Broglie scale in halos. One manifestation is vortices where the density vanishes and around which the velocity circulates. There is one vortex ring per de Broglie volume on average. For sufficiently low masses, soliton condensation occurs at centers of halos. The soliton oscillates and undergoes random walks, which is another manifestation of wave interference. The halo and subhalo abundance is expected to be suppressed at small masses, but the precise prediction from numerical wave simulations remains to be determined. For ultralight ∼10 −22 eV dark matter, the wave interference substructures can be probed by tidal streams or gravitational lensing. The signal can be distinguished from that due to subhalos by the dependence on stream orbital radius or image separation. Axion detection experiments are sensitive to interference substructures for wave dark matter that is moderately light. The stochastic nature of the waves affects the interpretation of experimental constraints and motivates the measurement of correlation functions. Current constraints and open questions, covering detection experiments and cosmological, galactic, and black hole observations, are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

dkainer/RWRtoolkit

a set of command-line and R tools for performing Random-walk with Restart analyses on multiplex networks in any species

Kainer, David↗
DOE OSTI2020

Heat Equation Neural Simulation

This code solves a steady state heat equation problem on a wire with random walks implemented using a spiking neural algorithm. SAND2020-12194 M

Reeder, Leah↗
DOE OSTI2024

Cylinders’ percolation: Decoupling and applications

In this paper we establish a strong decoupling inequality for the cylinder’s percolation process introduced by Tykesson and Windisch (Probab. Theory Related Fields 154 (2012) 165–191). This model features a very strong dependency structure, making it difficult to study, and this is why such decoupling inequalities are desirable. It is important to notice that the type of dependencies featured by cylinder’s percolation is particularly intricate, given that the cylinders have infinite range (unlike some models like Boolean percolation) while at the same time being rigid bodies (unlike processes such as random interlacements). Here our work introduces a new notion of fast decoupling, proves that it holds for the model in question and finishes with an application. More precisely, we prove that for a small enough density of cylinders, a random walk on a connected component of the vacant set is transient for all dimensions d≥3.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The Future of Computing: Integrating Scientific Computation on Neuromorphic Systems

Neuromorphic computing is known for its integration of algorithms and hardware elements that are inspired by the brain. Conventionally, this nontraditional method of computing is used for many neural or learning inspired applications. Unfortunately, this has resulted in the field of neuromorphic computing being relatively narrow in scope. In this paper we discuss two research areas actively trying to widen the impact of neuromorphic systems. The first is Fugu, a high-level programming interface designed to bridge the gap between general computer scientists and those who specialize in neuromorphic areas. The second aims to map classical scientific computing problems onto these frameworks through the example of random walks. This elucidates a class of scientific applications that are conducive to neuromorphic algorithms.

97 MATHEMATICS AND COMPUTING↗