Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic sampling”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Design and Optimization of a Gas-Cooled, Airfoil Fin Microchannel Heat Exchanger

High-performance microchannel heat exchangers are needed to supply heat for power conversion for nuclear microreactors. An airfoil fin microchannel design, constructed of Alloy 617 with helium as the working fluid, was analyzed and optimized using a design of experiments with artificial intelligence and machine learning techniques. The use of airfoil fins offers the potential to reduce pressure drop across the heat exchanger, as compared to other types of channel configurations. A framework for topology optimization of airfoil fin PCHEs has been developed that can be readily extended to different fin sizes and shapes, as well as different inlet and operating conditions, materials of construction, and working fluids. An optimization procedure was developed that employs computational fluid dynamics for a set of design points identified using Latin hypercube sampling. STAR-CCM+ was used to analyze a simplified two-channel configuration where five parameters were varied – inlet angle, fin scale, extent of staggering, transverse and longitudinal pitches. Two methods were compared for generating surrogate models – a 5D polynomial and a regression neural network. A response surface approximation was created from the surrogate models and input to a genetic algorithm. The genetic algorithm identified a set of optimal points on the Pareto front. The optimal geometry was found across six channel Reynolds numbers ranging from 1000 to 5000 to analyze how varying inlet conditions affects the optimal design. A set of optimal designs that maximizes heat transfer and minimizes pressure drop was identified, and a thermal stress analysis was performed on the optimal design. This work has developed a digital framework for the expedient topology design and evaluation of PCHE designs for gas-cooled microreactor applications. Correlations for the Nusselt number and Darcy friction factor were developed that can be useful for thermal hydraulic analyses using system codes. A thermal stress analysis was conducted and a brief discussion of the status of code cases of PCHEs for nuclear applications is given. Testing and thermomechanical modeling is needed to facilitate future code compliance of PCHEs for high pressure and high temperature applications.

42 ENGINEERING↗
DOE OSTI2022

Quiescent Galaxy Size, Velocity Dispersion, and Dynamical Mass Evolution

Abstract We use surveys covering the redshift range 0.05 < z < 3.8 to explore quiescent galaxy scaling relations and the redshift evolution of the velocity dispersion, size, and dynamical mass at fixed stellar mass. For redshift z < 0.6, we derive mass-limited samples and demonstrate that these large samples enhance constraints on the evolution of the quiescent population. The constraints include 2985 new velocity dispersions from the SHELS F2 survey. In contrast with the known substantial evolution of size with redshift, evolution in the velocity dispersion is negligible. The dynamical-to-stellar-mass ratio increases significantly as the universe ages, in agreement with recent results that combine high-redshift data with the Sloan Digital Sky Survey. Like other investigators, we interpret this result as an indication that the dark matter fraction within the effective radius increases as a result of the impact of the minor mergers that are responsible for size growth. We emphasize that dense redshift surveys covering the range 0.07 < z < 1 along with strong and weak lensing measurements could remove many ambiguities in evolutionary studies of the quiescent population.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Quantitative Prediction of Uranium Speciation and Amidoxime Binding in Seawater from Advanced Simulation Techniques

The goal of this project was to quantitatively predict stability of uranium species and thermochemistry of uranium binding with amidoxime ligands in realistic seawater conditions. The specific aims of the project were: (i) determine the stable species of uranium and its solvation structure in realistic seawater conditions; also speciation and solvation of uranium's major competing species in seawater, vanadium; (ii) accurately predict the free energies (and thereby stability constants) of binding between uranium and amidoxime; (iii) accurately predict the reaction enthalpies and entropies of binding between uranium and amidoxime as well as between vanadium and amidoxime. Main results and conclusions are as follows. First-principles molecular dynamic simulations found that the structure of the Ca 2 UO 2 (CO 3 ) 3 complex is very stable and that one Ca ion binds to the center UO 2 (CO 3 ) 3 4- anion stronger than the other Ca ion. This finding suggests that using time-resolved EXAFS spectra may confirm the asymmetry in binding of the two Ca ions in the aqueous Ca 2 UO 2 (CO 3 ) 3 complex. To consider the common ions, solvation of the Ca 2 UO 2 (CO 3 ) 3 complex in seawater was simulated by classical molecular dynamics simulation. It was found that the structure of the Ca 2 UO 2 (CO 3 ) 3 complex is very stable in the model seawater. A Na+ ion was found to be closely associated with the Ca 2 UO 2 (CO 3 ) 3 complex by indirectly interacting with one axial oxygen atom of the UO 2 group bridged by a water molecule. In addition, the Na + ion interacts closely with one Ca 2+ ion than the other. The present simulations revealed the key role of common ions such as Na + in impacting the solvation, structure, and apparent charge of the Ca 2 UO 2 (CO 3 ) 3 complex in seawater. The sequential processes of UO 2 2+ binding with the three CO 3 2- groups and the two Ca 2+ ions in pure water and in 0.1 M NaCl were simulated by classical molecular dynamics with both the non-constrained brute-force approach and umbrella sampling. The simulated free energies show excellent agreement with the experiment. The coordination structures of UO 2 (CO 3 ) 3 4- and CaUO 2 (CO 3 ) 3 2- were found to be significantly affected by the presence of Na + ions, leading to a monodentate binding of a carbonate group to U and a water molecule entering the first coordination shell of U. To shed light on the binding between uranyl and glutardiamidoxime (H 2 B), an important model and state-of-the-art ligand for seawater uranium extraction, a suite of computational methods was employed. From molecular dynamics (MD) simulations, it was found that the binding configurations of B 2- with uranyl favor the twofold distorted η 2 binding between the oximate ends (C=N-O-) and U, while HB - prefers the chelating mode for the oximate end with the neutral end being solvated by water. The free energies of sequential ligand binding to form UO 2 B, [UO 2 B 2 ] 2- , and [UO 2 (HB)B] - were simulated with umbrella sampling and very good agreement with the experimental values was achieved, which corroborates the structural insights into the binding mode. State-of-the-art polymeric sorbents employ both amidoximate and carboxylate groups on the side chains to achieve optimal U uptake and selectivity, so we simulated the binding of a model amidoximate–carboxylate bifunctional ligand with uranyl. Classical MD and free-energy simulations in 0.5 M NaCl showed that the carboxylate group binds first to uranyl, leading to a loose intermediate state, and then, the amidoximate group binds, resulting in a more stable and tight chelate state. Binding of the second bifunctional ligand follows a similar process, and the two ligands prefer a trans configuration around the uranyl group. The simulated free energies indicate that the two bifunctional ligands bind with uranyl 14 kcal/mol stronger than the two ligands with only amidoximate groups, confirming an important synergy between amidoximate and carboxylate groups in binding uranyl.

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

Linear Viscoelastic Properties of Putative Cyclic Polymers Synthesized by Reversible Radical Recombination Polymerization (R3P)

Linear viscoelastic properties in both melt and solution states are reported for a series of poly(3,6-dioxa-1,8-octanedithiol) (polyDODT) made by reversible radical recombination polymerization (R 3 P) under conditions designed to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT). PolyDODT is amorphous (T g < -50 °C) and highly flexible (entanglement molecular weight M e,lin ≈ 1850 g/mol for LDODT). PolyDODT’s low T g and low M e,lin enable characterization over a wide dynamic range and a wide range of dimensionless weight-average molecular weight Z w = M w /M e,lin . Measurements at temperatures from -57 °C to 100 °C provide up to 18 deades of reduced frequency, which is necessary to characterize RDODT melts with Z w from 23 to 300. The two highest molecular weight polymers in the present RDODT series have such high M w (406k and 556k g/mol) that mass spectrometry, NMR spectroscopy, and even chemical assays for chain ends are unable to rule out up to 2 mol% linear contaminant. By studying the samples in solution (using dilution to reduce Z w ) we could compare their dynamics with those of previously established high-purity polystyrene (PS) rings (limited to Z w ≤ 13.6). RDODT solutions with Z w < 15 (concentrations <5wt% for RDODT-406k and 556k) have dynamic moduli G^* that accord with LCCC-purified PS rings in terms of the frequency dependence (including the absence of a plateau), the progression of shapes of G* as a function of Z w , and the linear scaling of their zero-shear viscosity η0 with M w . The shape of G* as a function of Z w for solutions of RDODT-406k and -556k also accord with lower M w RDODT melts (which have ≤ 1.3mol% linear contaminant). Thus, measurement of the linear viscoelastic properties of appropriate concentrations of high M w (>200k g/mol) putative cyclic polymers, in which linear chains evade spectroscopic detection, may provide an alternative means (though not fully proven) of validation of sample purity. When Z w > 15 (including all seven RDODT melts and eight of their solutions), G* has a rubbery plateau. This suggests that the onset of entanglement-like behavior in rings requires 4-5-fold greater Z w than is required for linear chains. Further, the plateau moduli of RDODT samples are indistinguishable from G N o of the corresponding LDODT (melt or matched-concentration solutions). In entangled linear polymers, the observation that G N o is independent of Z w follows from limitations on lateral fluctuations due to neighboring chains becoming independent of position along a given chain. The present results for RDODT suggest that this holds for sufficiently long endless chains, too. While the RDODT have the same G N o as entangled LDODT, when Z w > 60 the terminal relaxation, if reached at all, of RDODT extends to orders of magnitude lower frequency than an entangled linear polymer of the same Z w . Consequently, the viscosity of RDODT with Z w > 60 increases with Z w much more strongly than the 3.4-power observed for entangled linear polymers. Lastly, these novel polymers, with disulfide-linked backbone and broad relaxation time distribution may prove important in relation to biodegradable elastomers and materials with exceptional low-frequency dissipation, extending at least 12 decades below the onset of the rubbery plateau.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Uncovering extreme nonlinear dynamics in solids through time-domain field analysis

Time-domain analysis of harmonic fields with sub-cycle resolution is now experimentally viable due to the emergence of sensitive, on-chip techniques for petahertz-scale optical-field sampling. We demonstrate how such a time-domain, field-resolved analysis uncovers the extreme nonlinear electron dynamics responsible for high-harmonic generation within solids. Time-dependent density functional theory was used to simulate harmonic generation from a solid-state band-gap system driven by near- to mid-infrared waveforms. Particular attention was paid to regimes where both intraband and interband emission mechanisms play a critical role in shaping the nonlinear response. We show that a time-domain analysis of the harmonic radiation fields identifies the interplay between intra- and interband dynamical processes underlying the nonlinear light generation. With further analysis, we show that changes to the dominant emission regime can occur after only slight changes to the peak driving intensity and central driving wavelength. Furthermore, time-domain analysis of harmonic fields also reveals, for the first time, the possibility of rapid changes in the dominant emission mechanism within the temporal window of the driving pulse envelope. Finally, we examine the experimental viability of performing time-domain analysis of harmonic fields with sub-cycle resolution using realistic parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Dynamic reflectance changes in high-purity zirconium when it is shocked from the α to ω phase

Flash lamp-illuminated integrating spheres have been used to measure the dynamic reflectance at the interface between a shocked metal target and a window maintaining the shock stress. In conjunction with radiance measurements, this method was recently used to study the pressure versus temperature phase boundary of the beta phase to body-centered-tetragonal transition of tin at around 13 GPa and, in separate experiments, to show a relationship between the reflectance and the phase fraction. We have used this technique to measure the ratio, R, of the dynamic reflectance relative to the ambient value for a high-purity zirconium sample attached to a sapphire window to study the alpha to omega phase transition. Nine experiments cover the range of stress, σ, from 4 to 16 GPa at six wavelengths from 500 to 1550 nm. The two phases in zirconium differ in volume by <1%, so a phase change between them can be difficult to identify from anomalies in the particle velocity measured on the shock Hugoniot. The phase boundary is believed to cross the Hugoniot around 3 GPa, but dynamically it is observed at much higher stresses >6 GPa. Our results show that the relative reflectance decreases upon initial shock in the alpha state but increases as the material transforms to its omega state. The slope dR/dσ changes the sign around the phase transition. For shocks below 6 GPa, dR/dσ < 0, and for σ above 10 GPa, dR/dσ > 0. In the stress region between 6 and 10 GPa, dR/dσ varies somewhat, perhaps due to the band structure or because the transition is incomplete there. I. INTRODUCTION

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Effect of Variations in the Gas-Dynamic Unloading on the Laser Initiation of the PETN–Aluminum Composite

Thresholds of explosive decomposition H{sub cr} under irradiation using the first harmonic of a pulsed neodymium laser (14 ns) are experimentally determined for pentaerythritol tetranitrate (PETN) with aluminum nanoparticles (100 nm) at concentrations ranging from 0.05 to 1 wt %. Amplitudes of optoacoustic signals are measured versus concentrations of impurities in the sample at a fixed laser fluence. The experiments are performed in two regimes. In the first one, the irradiated surface of the sample is covered with a finite-weight glass plate that impedes the gas-dynamic unloading. In the second regime, an external pressure of no less than 10{sup 7} Pa is applied to the glass plate to prevent the gas-dynamic unloading. The interpretation of the experimental effects is proven by the results of theoretical simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

Using Importance Sampling Monte Carlo to Analyze Aircraft Takeoff and Landing

Aircraft have achieved high levels of reliability, so the probability of undesirable dynamics is very low. The low probability of a bad takeoff or landing challenges nondeterministic methods to quantify the uncertainty of these improbable events. In the present work, a Monte Carlo method is defined to more effectively quantify these unlikely events that is suitable for uncertainty analysis of a mature design with many uncertainty parameters. The proposed Monte Carlo method has been applied to the takeoff and landing of the X-59 Quiet SuperSonic Technology (QueSST). Takeoff and landing simulation profiles were defined using a combination of industry standards and piloted simulation experience. The proposed method is shown to be able to quantify events with significantly fewer samples than would be required with conventional Monte Carlo. These results improve understanding of the sensitivity of the takeoff and landing dynamics to inform further studies and test planning.

Jeffrey Ouellette↗
NASA NTRSDate not supplied

Using Importance Sampling Monte Carlo to Analyze Aircraft Takeoff and Landing

Aircraft have achieved high levels of reliability, so the probability of undesirable dynamics is very low. The low probability of a bad takeoff or landing challenges nondeterministic methods to quantify the uncertainty of these improbable events. In the present work, a Monte Carlo method is defined to more effectively quantify these unlikely events that is suitable for uncertainty analysis of a mature design with many uncertainty parameters. The proposed Monte Carlo method has been applied to the takeoff and landing of the X-59 Quiet SuperSonic Technology (QueSST). Takeoff and landing simulation profiles were defined using a combination of industry standards and piloted simulation experience. The proposed method is shown to be able to quantify events with significantly fewer samples than would be required with conventional Monte Carlo. These results improve understanding of the sensitivity of the takeoff and landing dynamics to inform further studies and test planning.

Jeffrey Ouellette↗
DOE OSTI2019

Considerations for imaging thick, low contrast, and beam sensitive samples with liquid cell transmission electron microscopy

Transmission electron microscopy of whole cells is hindered by the inherently large thickness and low atomic contrast intrinsic of cellular material. Liquid cell transmission electron microscopy, allows samples to remain in their native hydrated state and may permit visualizing cellular dynamics in-situ. However, imaging biological cells with this approach remains challenging and identifying an optimal imaging regime using empirical data would help foster new advancements in the field. Recent questions about the role of the electron beam inducing morphological changes or damaging cellular structure and function necessitates further investigation of electron beam-cell interactions, but is complicated by variability in imaging techniques used across various studies currently present in literature. The necessity for using low electron fluxes for imaging biological samples requires finding an imaging strategy which produces the strongest contrast and signal to noise ratio for the electron flux used. Here, we experimentally measure and evaluate signal to noise ratios and damage mechanisms between liquid and cryogenic samples for cells using multiple electron imaging modalities all on the same instrument and with equivalent beam parameters to standardize the comparison. We also discuss considerations for optimal electron microscopy imaging conditions for future studies on whole cells within liquid environments.

Moser, Trevor H. (ORCID:0000000193990514)↗
DOE OSTI2022

Atomistic simulations of nanoindentation on nanoglasses: Effects of grain size and gradient microstructure on the mechanical properties

The use of nanoglass (NG) microstructures has been shown as an effective strategy to improve the ductility of rather brittle metallic glass alloys. To evaluate the effects of grain size and gradient design on the mechanical properties of NGs, we perform molecular dynamics simulations of nanoindentation on Cu 64 Zr 36 NGs. We consider samples with uniform 3 and 7 nm grain sizes as well as gradient microstructures with gradient grain sizes varying from 3 to 7 nm. Here, the results show that the deformation mechanism in NGs with small grain sizes is dominated by the activation and evolution of multiple shear transformation zones. Increasing grain sizes enhances the elastic modulus and hardness at the cost of reduced abrasion resistance. Though the average grain size at the indentation surface plays a crucial role in the deformation behavior of gradient NG models, their strain localization and plastic deformation states are affected by the grain sizes far from the indenter. Desired combinations of mechanical properties can be realized by different NG microstructure designs.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Neutron and x-ray scattering study of phonon dispersion and diffuse scattering in ( Na , Bi ) Ti O 3 - x BaTi O 3 single crystals near the morphotropic phase boundary

Neutron and x-ray scattering measurements were performed on (Na 1/2 Bi 1/2 )TiO 3-x at % BaTiO 3 (NBT-xBT) single crystals (x = 4, 5, 6.5, and 7.5) across the morphotropic phase boundary (MPB), as a function of both composition and temperature, and probing both structural and dynamical aspects. In addition to the known diffuse scattering pattern near the gamma points, our measurements revealed new, faint superlattice peaks, as well as an extensive diffuse scattering network, revealing a short-range ordering of polar nanoregions (PNR) with a static stacking morphology. Furthermore, in samples with compositions closest to the MPB, our inelastic neutron scattering investigations of the phonon dynamics showed two unusual features in the acoustic phonon branches, between the superlattice points, and between the superlattice points and gamma points, respectively. Finally, these critical elements are not present in the other compositions away from the MPB, which suggests that these features may be related to the tilt modes coupling behavior near the MPB.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

MemGaze: Rapid and Effective Load-Level Memory Trace Analysis

A major challenge of memory analysis tools is combining high-resolution analysis and low overhead measurement. Currently, hardware/software-based analysis of load-level sequences incurs time slowdowns of O(100×). We present MemGaze, a tool for low-overhead, high-resolution memory analysis. MemGaze uses Intel’s Processor Tracing (PT) instruction ptwrite to collect sampled and compressed memory address traces for load-level, sequence-aware analysis of data reuse. We describe multi-resolution analysis for locations vs. operations, accesses vs. spatio-temporal reuse, and reuse (distance, rate, volume) vs. access patterns. Both trace size and resolution are controllable. We use MemGaze to elucidate the memory effects of different data structures and algorithms. For sampled traces that are ˜1% of a full one, analysis metrics have 1-25% MAPE for histograms of varying dynamic sequence lengths. With current suboptimal kernel support (PT runs continuously), MemGaze’s time overhead is typically 10–95%; 7× at worst. However, when PT runs only during samples, overhead is 10–35% on memory intensive regions and correlates with executed ptwrites.

Kilic, Ozgur O.↗
NASA NTRS1987

Identification of structural interface characteristics using component mode synthesis

The inability to adequately model connections has limited the ability to predict overall system dynamic response. Connections between structural components are often mechanically complex and difficult to accurately model analytically. Improved analytical models for connections are needed to improve system dynamic predictions. This study explores combining Component Mode synthesis methods for coupling structural components with Parameter Identification procedures for improving the analytical modeling of the connections. Improvements in the connection properties are computed in terms of physical parameters so the physical characteristics of the connections can be better understood, in addition to providing improved input for the system model. Two sample problems, one utilizing simulated data, the other using experimental data from a rotor dynamic test rig are presented.

Huckelbridge, A. A.↗
NASA NTRS1987

Identification of structural interface characteristics using component mode synthesis

The inability to adequately model connections has limited the ability to predict overall system dynamic response. Connections between structural components are often mechanically complex and difficult to accurately model analytically. Improved analytical models for connections are needed to improve system dynamic predictions. This study explores combining Component Mode synthesis methods for coupling structural components with Parameter Identification procedures for improving the analytical modeling of the connections. Improvements in the connection properties are computed in terms of physical parameters so the physical characteristics of the connections can be better understood, in addition to providing improved input for the system model. Two sample problems, one utilizing simulated data, the other using experimental data from a rotor dynamic test rig, are presented.

Huckelbridge, A. A.↗
NASA NTRS1989

Identification of structural interface characteristics using component mode synthesis

The inability to adequately model connections has limited the ability to predict overall system dynamic response. Connections between structural components are often mechanically complex and difficult to accurately model analytically. Improved analytical models for connections are needed to improve system dynamic predictions. This study explores combining Component Mode synthesis methods for coupling structural components with Parameter Identification procedures for improving the analytical modeling of the connections. Improvements in the connection properties are computed in terms of physical parameters so the physical characteristics of the connections can be better understood, in addition to providing improved input for the system model. Two sample problems, one utilizing simulated data, the other using experimental data from a rotor dynamic test rig, are presented.

Huckelbridge, A. A.↗
DOE OSTI2022

Spatial turnover of soil viral populations and genotypes overlain by cohesive responses to moisture in grasslands

Viruses shape microbial communities, food web dynamics, and carbon and nutrient cycling in diverse ecosystems. However, little is known about the patterns and drivers of viral community composition, particularly in soil, precluding a predictive understanding of viral impacts on terrestrial habitats. To investigate soil viral community assembly processes, here we analyzed 43 soil viromes from a rainfall manipulation experiment in a Mediterranean grassland in California. We identified 5,315 viral populations (viral operational taxonomic units [vOTUs] with a representative sequence ≥10 kbp) and found that viral community composition exhibited a highly significant distance–decay relationship within the 200-m 2 field site. This pattern was recapitulated by the intrapopulation microheterogeneity trends of prevalent vOTUs (detected in ≥90% of the viromes), which tended to exhibit negative correlations between spatial distance and the genomic similarity of their predominant allelic variants. Although significant spatial structuring was also observed in the bacterial and archaeal communities, the signal was dampened relative to the viromes, suggesting differences in local assembly drivers for viruses and prokaryotes and/or differences in the temporal scales captured by viromes and total DNA. Despite the overwhelming spatial signal, evidence for environmental filtering was revealed in a protein-sharing network analysis, wherein a group of related vOTUs predicted to infect actinobacteria was shown to be significantly enriched in low-moisture samples distributed throughout the field. Overall, our results indicate a highly diverse, dynamic, active, and spatially structured soil virosphere capable of rapid responses to changing environmental conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Magnetic tunnel junction random number generators applied to dynamically tuned probability trees driven by spin orbit torque

Abstract Perpendicular magnetic tunnel junction (pMTJ)-based true-random number generators (RNGs) can consume orders of magnitude less energy per bit than CMOS pseudo-RNGs. Here, we numerically investigate with a macrospin Landau–Lifshitz-Gilbert equation solver the use of pMTJs driven by spin–orbit torque to directly sample numbers from arbitrary probability distributions with the help of a tunable probability tree. The tree operates by dynamically biasing sequences of pMTJ relaxation events, called ‘coinflips’, via an additional applied spin-transfer-torque current. Specifically, using a single, ideal pMTJ device we successfully draw integer samples on the interval [0, 255] from an exponential distribution based on p -value distribution analysis. In order to investigate device-to-device variations, the thermal stability of the pMTJs are varied based on manufactured device data. It is found that while repeatedly using a varied device inhibits ability to recover the probability distribution, the device variations average out when considering the entire set of devices as a ‘bucket’ to agnostically draw random numbers from. Further, it is noted that the device variations most significantly impact the highest level of the probability tree, with diminishing errors at lower levels. The devices are then used to draw both uniformly and exponentially distributed numbers for the Monte Carlo computation of a problem from particle transport, showing excellent data fit with the analytical solution. Finally, the devices are benchmarked against CMOS and memristor RNGs, showing faster bit generation and significantly lower energy use.

77 NANOSCIENCE AND NANOTECHNOLOGY↗