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 217 records · Page 12

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↗

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↗

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↗

Diffusive Escape of a Nanoparticle from a Porous Cavity

Narrow escape from confinement through a nano-channel is the critical step of complex transport processes including size exclusion-based separations, oil and gas extraction from the microporous subsurface environment, and RNA translocation through nuclear pore complex channels. While narrow escape has been studied using theoretical and computational methods, experimental quantification is rare because of the difficulty in confining a particle into a microscopic space through a nanoscale hole. In this paper, we studied narrow escape in the context of continuous nanoparticle diffusion within the liquid-filled void space of an ordered porous material. Specifically, we quantified the spatial dependence of nanoparticle motion and the sojourn times of individual particles in the interconnected confined cavities of a liquid-filled inverse opal film. We found that nanoparticle motion was inhibited near cavity walls and cavity escape was slower than predicted by existing theories and random-walk simulations. Here, a combined computational-experimental analysis indicated that translocation through a nano-channel is barrier-controlled rather than diffusion-controlled.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Accumulation-Driven Surfactant-Free Synthesis of Architectured Immiscible Metallic Nanoalloys with Enhanced Catalysis

Accumulation-mediated chemical reactions are a ubiquitous phenomenon in nature. In this work, we explore microbubble-induced accumulation of precursor ions to achieve surfactant-free synthesis of immiscible metallic nanoalloys and simultaneously pattern the nanoalloys into targeted architectures for their enhanced catalytic applications. Additionally, our unified spatiotemporal synthesis and structuring (US3) strategy, whereby millisecond-scale accumulation of the ions takes place in a highly confined laser-mediated microbubble trap (MBT), drives ultrafast alloy synthesis in sync with the structuring process. As a case in point, we employ the US3 strategy for the in situ surfactant-free synthesis and patterning of traditionally immiscible rhodium-gold (RhAu) nanoalloys. Stochastic random walk simulations justify the millisecond-scale accumulation process, leading to a 3-order reduction in synthesis time. The catalytic activity and structure-property relationship were evaluated using the reduction of p-nitrophenol with NaBH4. Our in situ synthesis and structuring strategy can be translated for high-throughput production and screening of multimetallic systems with tailored catalytic, optoelectronic, and magnetic functions.

36 MATERIALS SCIENCE↗

Rigorous Data Fusion for Computationally Expensive Simulations

This manuscript comprises the final report for the 1-year, FY19 LDRD project "Rigorous Data Fusion for Computationally Expensive Simulations," wherein an alternative approach to Bayesian calibration was developed based a new sampling technique called VoroSpokes. Vorospokes is a novel quadrature and sampling framework defined with respect to Voronoi tessellations of bounded domains in $R^d$ developed within this project. In this work, we first establish local quadrature and sampling results on convex polytopes using randomly directed rays, or spokes, to approximate the quantities of interest for a specified target function. A theoretical justification for both procedures is provided along with empirical results demonstrating the unbiased convergence in the resulting estimates/samples. The local quadrature and sampling procedures are then extended to global procedures defined on more general domains by applying the local results to the cells of a Voronoi tessellation covering the domain in consideration. We then demonstrate how the proposed global sampling procedure can be used to define a natural framework for adaptively constructing Voronoi Piecewise Surrogate (VPS) approximations based on local error estimates. Finally, we show that the adaptive VPS procedure can be used to form a surrogate model approximation to a specified, potentially unnormalized, density function, and that the global sampling procedure can be used to efficiently draw independent samples from the surrogate density in parallel. The performance of the resulting VoroSpokes sampling framework is assessed on a collection of Bayesian inference problems and is shown to provide highly accurate posterior predictions which align with the results obtained using traditional methods such as Gibbs sampling and random-walk Markov Chain Monte Carlo (MCMC). Importantly, the proposed framework provides a foundation for performing Bayesian inference tasks which is entirely independent from the theory of Markov chains.

97 MATHEMATICS AND COMPUTING↗