Search NASA⌕ Search

SEARCH · Search NASA

Results for “Randomized methods”

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 433 records · Page 24

3D exemplar-based random walks for tooth segmentation from cone-beam computed tomography images

Purpose: Tooth segmentation is an essential step in acquiring patient-specific dental geometries from cone-beam computed tomography (CBCT) images. Tooth segmentation from CBCT images is still a challenging task considering the comparatively low image quality caused by the limited radiation dose, as well as structural ambiguities from intercuspation and nearby alveolar bones. The goal of this paper is to present and discuss the latest accomplishments in semisupervised tooth segmentation with adaptive 3D shape constraints. Methods: The authors propose a 3D exemplar-based random walk method of tooth segmentation from CBCT images. The proposed method integrates semisupervised label propagation and regularization by 3D exemplar registration. To begin with, the pure random walk method is to get an initial segmentation of the teeth, which tends to be erroneous because of the structural ambiguity of CBCT images. And then, as an iterative refinement, the authors conduct a regularization by using 3D exemplar registration, as well as label propagation by random walks with soft constraints, to improve the tooth segmentation. In the first stage of the iteration, 3D exemplars with well-defined topologies are adapted to fit the tooth contours, which are obtained from the random walks based segmentation. The soft constraints on voxel labeling are defined by shape-based foreground dentine probability acquired by the exemplar registration, as well as the appearance-based probability from a support vector machine (SVM) classifier. In the second stage, the labels of the volume-of-interest (VOI) are updated by the random walks with soft constraints. The two stages are optimized iteratively. Instead of the one-shot label propagation in the VOI, an iterative refinement process can achieve a reliable tooth segmentation by virtue of exemplar-based random walks with adaptive soft constraints. Results: The proposed method was applied for tooth segmentation of twenty clinically captured CBCT images. Three metrics, including the Dice similarity coefficient (DSC), the Jaccard similarity coefficient (JSC), and the mean surface deviation (MSD), were used to quantitatively analyze the segmentation of anterior teeth including incisors and canines, premolars, and molars. The segmentation of the anterior teeth achieved a DSC up to 98%, a JSC of 97%, and an MSD of 0.11 mm compared with manual segmentation. For the premolars, the average values of DSC, JSC, and MSD were 98%, 96%, and 0.12 mm, respectively. The proposed method yielded a DSC of 95%, a JSC of 89%, and an MSD of 0.26 mm for molars. Aside from the interactive definition of label priors by the user, automatic tooth segmentation can be achieved in an average of 1.18 min. Conclusions: The proposed technique enables an efficient and reliable tooth segmentation from CBCT images. This study makes it clinically practical to segment teeth from CBCT images, thus facilitating pre- and interoperative uses of dental morphologies in maxillofacial and orthodontic treatments.

60 APPLIED LIFE SCIENCES↗

Design Space Exploration of Emerging Memory Technologies for Machine Learning Applications

Memory design space exploration methods study memory systems’ performances and limitations before implementation. The computer memory design space has grown exponentially because of the enormous growth of memory types, memory controllers, and application software. Computer simulators are commonly used for memory design space exploration. However, complex memory simulations take an enormous amount of time. Hence, in this paper, we proposed a machine learning-based design space exploration method for dynamic random-access memory and non-volatile memory systems. We applied our method to the CosmoGAN and LeNet applications to predict the following six memory response parameters: (i) bandwidth, (ii) power, (iii) average latency, (iv) average total latency, (v) memory reads, and (vi) memory writes. Our experimental results show that machine learning models can predict memory response parameter values faster than simulations. We used support vector machine, random forest, and gradient boosting machine learning models. We observed that the support vector machine provides better performance for bandwidth, average latency, and average total latency. The random forest model works better for memory reads and writes. The gradient boosting model provides superior prediction performance for power. We provide a detailed discussion on learning curve characteristics, error analysis, and memory type recommendation.

Hasan, S M Shamimul↗

Shape evolution of neutron-rich 106,108,110 Mo isotopes in the triaxial degree of freedom

Background: Neutron-rich nuclei with mass number between 100 and 110 attract much attention, since several kinds of shapes, such as spherical, prolate, oblate, and triaxial shapes, are predicted. In particular, for neutron-rich Mo isotopes, different models predict different magnitudes and rigidity of triaxial deformation. Previous interpretations of experimental results based solely on low-lying $2^+_2$ states are insufficient to distinguish between the rigid triaxial shape, $\gamma$ vibration, or $\gamma$-soft rotor. Purpose: The shape evolution of 106 Mo, 108 Mo, and 110 Mo is investigated through their $2^+_1$-state lifetimes, decay-branching ratios of the $2^+_2$ state, and energies of the low-lying collective excited states with $K^π = 0^+, 2^+$, and $4^+$. Method: $\beta$-delayed $\gamma$-ray spectroscopy was employed for neutron-rich Nb and Zr isotopes produced at the RIKEN RI Beam Factory to populate excited states in 106 Mo, 108 Mo, and 110 Mo. The EUroball-RIKEN Cluster Array was used for high-resolution $\gamma$-ray detection and lifetimes of the $2^+_1$ states were determined using the UK fast-timing array of LaBr 3 (Ce) detectors. Results: New $\gamma$-ray transitions and levels are reported, including newly assigned $0^+_2$ states in 108,110 Mo. Quadrupole deformations were obtained for 106,108,110 Mo from their $2^+_1$ energies and lifetimes. The $\beta$-delayed neutron-emission probabilities of 108 Nb and 110 Nb were determined by examining the $\gamma$ rays of their respective daughter decays. Conclusions: In this work, the even-odd energy staggering in the $2^+_2$ band was compared with typical patterns of the $\gamma$-vibrational band, rigid triaxial rotor, and $\gamma$-soft rotor. The very small even-odd staggering of 106 Mo, 108 Mo, and 110 Mo favors a $\gamma$-vibrational band assignment. The kinematic moment of inertia for the $2^+_2$ band showed a trend similar to the ground-state band, which is as expected for the $\gamma$-vibrational band. Beyond-mean-field calculations employing the constrained Hartree-Fock-Bogoliubov and local quasiparticle-random-phase approximation method using the SLy5 + T interaction reproduced the ground and $2^+_2$ bands in 106 Mo and 108 Mo. The collective wave functions are consistent with the interpretation of the $2^+_2$ band as the $\gamma$-vibrational band of the prolate shape. However, the staggering pattern observed in 110 Mo differs from the one suggested in the calculations which predict a $\gamma$-soft rotor. There was no experimental indication of the oblate shape or the $\gamma$-soft rotor predicted in these Mo isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Use of adjoint methods in the probabilistic finite element approach to fracture mechanics

The adjoint method approach to probabilistic finite element methods (PFEM) is presented. When the number of objective functions is small compared to the number of random variables, the adjoint method is far superior to the direct method in evaluating the objective function derivatives with respect to the random variables. The PFEM is extended to probabilistic fracture mechanics (PFM) using an element which has the near crack-tip singular strain field embedded. Since only two objective functions (i.e., mode I and II stress intensity factors) are needed for PFM, the adjoint method is well suited.

Liu, Wing Kam↗

Stimulated resonance Raman scattering from organic dyes in a multiple-scattering medium as a potential method for determining their vibrational spectra

A method is described for deriving Raman spectra of organic dyes from their random lasing spectra. The method was tested using Rhodamine 6G. The Raman spectrum obtained for this dye agrees well with the spectra measured by standard techniques but is more structured, which allows unresolved features to be detected. The spectrum provides more detailed information owing to the interference between the Raman scattered light and amplified spontaneous emission of the dye molecules within a photon mean free path. One advantage of the method is that the luminescence of the dye helps to observe Raman lines, which allows one to work in the Stokes region and facilitates the measurement procedure. (nonlinear optical phenomena)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spin-crossover complexes: Self-interaction correction vs density correction

Complexes containing a transition metal atom with a 3d 4 –3d 7 electron configuration typically have two low-lying, high-spin (HS) and low-spin (LS) states. The adiabatic energy difference between these states, known as the spin-crossover energy, is small enough to pose a challenge even for electronic structure methods that are well known for their accuracy and reliability. In this work, we analyze the quality of electronic structure approximations for spin-crossover energies of iron complexes with four different ligands by comparing energies from self-consistent and post-self-consistent calculations for methods based on the random phase approximation and the Fermi–Löwdin self-interaction correction. Considering that Hartree–Fock densities were found by Song et al., J. Chem. Theory Comput. 14, 2304 (2018), to eliminate the density error to a large extent, and that the Hartree–Fock method and the Perdew–Zunger-type self-interaction correction share some physics, we compare the densities obtained with these methods to learn their resemblance. Here, we find that evaluating non-empirical exchange-correlation energy functionals on the corresponding self-interaction-corrected densities can mitigate the strong density errors and improves the accuracy of the adiabatic energy differences between HS and LS states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

A new insight into neutrino energy loss by electron capture of iron group nuclei in magnetar surfaces

Based on the relativistic mean-field effective interactions theory, and the Lai dong model, we discuss the influences of superstrong magnetic fields (SMFs) on electron Fermi energy, nuclear blinding energy, and single-particle level structure in magnetar surfaces. Using the Shell-Model Monte Carlo method and the Random Phase Approximation theory, we analyze the neutrino energy loss rates (NELRs) by electron capture for iron group nuclei in SMFs. First, when B {sub 12} < 100, we find that the SMFs have a slight influence on the NELRs for most nuclides at relativistic low temperatures (e.g., T {sub 9} = 0.233); nevertheless, the NELRs increase by more than four orders of magnitude at relativistic high temperatures (e.g., T {sub 9} = 15.53). When B {sub 12} > 100, the NELRs decrease by more than three orders of magnitude (e.g., at T {sub 9} = 15.53 for {sup 52–61}Fe, {sup 55–60}Co, and {sup 56–63}Ni). Second, for a certain value of magnetic field and temperature, the NELRs increase by more than four orders of magnitude when ρ{sub 7}⩽10{sup 3}, but as the density increases (i.e., when ρ{sub 7}>10{sup 3}), there is almost no influence on the density of NELRs. For the density around ρ{sub 7}=10{sup 2}, there is an abrupt increase in NELRs when B {sub 12} ≥ 10{sup 3.5}. Such jumps are an indication that the underlying shell structure has changed due to single-particle behavior by SMFs. Finally, we compare our NELRs with those of Fuller et al. (FFN) and Nabi and Klapdor-Kleingrothaus (NKK). For the case without SMFs, one finds that our rates for certain nuclei are close to about five orders of magnitude lower than FFN and NKK at relativistic low temperatures (e.g., T {sub 9} = 1). However, at a relativistic high temperature (e.g., T {sub 9} = 3), our results are in good agreement with NKK, but about one order of magnitude lower than FFN. For the case with SMFs, our NELRs for some iron group nuclei can be about five orders of magnitude higher than those of FFN and NKK. (Note that B {sub 12}, T {sub 9}, and ρ {sub 7} are in units of 10{sup 12} G, 10{sup 9} K, and 10{sup 7} g cm{sup −3}, respectively.)

79 ASTRONOMY AND ASTROPHYSICS↗

A partially reflecting random walk on spheres algorithm for electrical impedance tomography

In this work, we develop a probabilistic estimator for the voltage-to-current map arising in electrical impedance tomography. This novel so-called partially reflecting random walk on spheres estimator enables Monte Carlo methods to compute the voltage-to-current map in an embarrassingly parallel manner, which is an important issue with regard to the corresponding inverse problem. Our method uses the well-known random walk on spheres algorithm inside subdomains where the diffusion coefficient is constant and employs replacement techniques motivated by finite difference discretization to deal with both mixed boundary conditions and interface transmission conditions. We analyze the global bias and the variance of the new estimator both theoretically and experimentally. Subsequently, the variance of the new estimator is considerably reduced via a novel control variate conditional sampling technique which yields a highly efficient hybrid forward solver coupling probabilistic and deterministic algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Image processing of multiphase images obtained via X-ray microtomography: A review

Easier access to X-ray microtomography (μCT) facilities has provided much new insight from high-resolution imaging for various problems in porous media research. Pore space analysis with respect to functional properties usually requires segmentation of the intensity data into different classes. Image segmentation is a nontrivial problem that may have a profound impact on all subsequent image analyses. This review deals with two issues that are neglected in most of the recent studies on image segmentation: (i) focus on multiclass segmentation and (ii) detailed descriptions as to why a specific method may fail together with strategies for preventing the failure by applying suitable image enhancement prior to segmentation. In this way, the presented algorithms become very robust and are less prone to operator bias. Three different test images are examined: a synthetic image with ground-truth information, a synchrotron image of precision beads with three different fluids residing in the pore space, and a μCT image of a soil sample containing macropores, rocks, organic matter, and the soil matrix. Image blur is identified as the major cause for poor segmentation results. Other impairments of the raw data like noise, ring artifacts, and intensity variation can be removed with current image enhancement methods. Bayesian Markov random field segmentation, watershed segmentation, and converging active contours are well suited for multiclass segmentation, yet with different success to correct for partial volume effects and conserve small image features simultaneously.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Protein engineering for metabolic engineering: Current and next-generation tools

Protein engineering in the context of metabolic engineering is increasingly important to the field of industrial biotechnology. As the demand for biologically produced food, fuels, chemicals, food additives, and pharmaceuticals continues to grow, the ability to design and modify proteins to accomplish new functions will be required to meet the high productivity demands for the metabolism of engineered organisms. We review advances in selecting, modeling, and engineering proteins to improve or alter their activity. Some of the methods have only recently been developed for general use and are just beginning to find greater application in the metabolic engineering community. We also discuss methods of generating random and targeted diversity in proteins to generate mutant libraries for analysis. Recent uses of these techniques to alter cofactor use; produce non-natural amino acids, alcohols, and carboxylic acids; and alter organism phenotypes are presented and discussed as examples of the successful engineering of proteins for metabolic engineering purposes.

Marcheschi, RJ↗

Fullrmc, a rigid body reverse monte carlo modeling package enabled with machine learning and artificial intelligence

Here, a new Reverse Monte Carlo (RMC) package ‘fullrmc’ for atomic or rigid body and molecular, amorphous or crystalline materials is presented. fullrmc main purpose is to provide a fully modular, fast and flexible software, thoroughly documented, complex molecules enabled, written in a modern programming language (python, cython ,C and C++ when performance is needed) and complying to modern programming practices. fullrmc approach in solving an atomic or molecular structure is different from existing RMC algorithms and software. In a nutshell, traditional RMC methods and software randomly adjust atom positions until the whole system has the greatest consistency with a set of experimental data. In contrast, fullrmc applies smart moves endorsed with reinforcement machine learning to groups of atoms. While fullrmc allows running traditional RMC modelling, the uniqueness of this approach resides in its ability to customize grouping atoms in any convenient way with no additional programming efforts and to apply smart and more physically meaningful moves to the defined groups of atoms. Also fullrmc provides a unique way with almost no additional computational cost to recur a group’s selection, allowing the system to go out of local minimas by refining a group’s position or exploring through and beyond not allowed positions and energy barriers the unrestricted three dimensional space around a group.

97 MATHEMATICS AND COMPUTING↗

Comparison of Sedoanalgesia Versus Ultrasound-Guided Supraclavicular Brachial Plexus Block for the Prevention of the Pain During Endovascular Treatment of Dysfunctional Hemodialysis Fistulas

Purpose: Although intravenous sedation and analgesia have been widely used as a first choice to relieve pain during treatment of dysfunctional hemodialysis fistulas by interventional radiology, the sedoanalgesic drugs have a considerable risk of respiratory depression, especially in hemodialysis patients. In this study, we compared the utility and efficiency of ultrasound-guided supraclavicular brachial plexus block versus sedoanalgesia for the prevention of pain during endovascular treatment of dysfunctional hemodialysis fistulas Materials and Methods: Patients were randomized into two groups: ultrasound-guided supraclavicular brachial plexus block (n = 34) or sedoanalgesia group (n = 34). A visual analogue scale from no pain (= 0) to worst pain possible (= 10) was used to assess the pain intensity. Patient and operator satisfaction were graded from 0 to 2: 0, not satisfied at all; 1, partially satisfied; 2, satisfied (very well or complete satisfaction). Both groups were compared in terms of pain scores, patient and operator satisfaction as well as complications. Results: The median pain score was significantly lower in the block group compared to the sedoanalgesia group, 0 (0–4) versus 6 (2–10), p = 0.0001. Patient satisfaction and operator satisfaction were significantly higher in the block group than in the sedoanalgesia group (p = 0.0001). Severe oxygen desaturation occurred in five (14.7%) patients following the administration of sedoanalgesia. No side effects or complications related to block procedure occurred in any patient. Conclusion: Ultrasound-guided supraclavicular brachial plexus block has advantages over the sedoanalgesia during endovascular treatment of dysfunctional hemodialysis fistulas. It can provide safe and efficient analgesia with excellent procedural satisfaction in adult hemodialysis patients. Level of Evidence: Level 1 (randomized controlled trial).

62 RADIOLOGY AND NUCLEAR MEDICINE↗

UQpy: A general purpose Python package and development environment for uncertainty quantification

In this paper, we present the UQpy software toolbox, an open-source Python package for general uncertainty quantification (UQ) in mathematical and physical systems. The software serves as both a user-ready toolbox that includes many of the latest methods for UQ in computational modeling and a convenient development environment for Python programmers advancing the field of UQ. The paper presents an introduction to the software's architecture and existing capabilities, divided in the code in a set of modules centered around different UQ tasks such as sampling methods, generation of random processes and random fields, probabilistic inverse modeling, reliability analysis, surrogate modeling, and active learning. The paper also highlights the importance of the RunModel module, which is used to drive simulations in the uncertainty analyses performed in UQpy. This module conveniently allows the user to define computational models directly in Python, or to run simulations from a third-party software in serial or in parallel. To illustrate the various capabilities, two examples are tracked throughout the paper and analyzed repeatedly for various UQ tasks. The first is a Python model solving a nonlinear structural dynamics problem, used to illustrate UQpy's capabilities in sampling and forward propagation of high dimensional random vectors (stochastic processes), and probabilistic inference. The second model is a third-party Abaqus finite element model solving the thermomechanical response of a beam structure. This example is used to illustrate UQpy's capabilities in variance reduction sampling techniques, reliability analysis, surrogate modeling and active learning techniques.

97 MATHEMATICS AND COMPUTING↗

Analyzing the scale of the bainitic ferrite plates by XRD, SEM and TEM

Since the major strengthening mechanisms in nanocrystalline bainitic steels arise from the exceptionally small size of the bainitc ferrite plate, accurate determination of this parameter is fundamental for quantitative relating the microstructure to the mechanical properties. In this work, the thickness of the bainitic ferrite subunits obtained by different bainitic heat treatments was determined in two steels, with carbon contents of 0.3 and 0.7 wt.%, from SEM and TEM micrographs. As these measurements were made on 2D images taken from random sections, the method includes some stereological correction factors to obtain accurate information. Finally, the determined thicknesses of bainitic ferrite plates were compared with the crystallite size calculated from the analysis of X-ray diffraction peak broadening. Although in some case the values obtained for crystallite size and plate thickness can be similar, this study confirms that indeed they are two different parameters. - Highlights: •Bainitic microstructure in a nanostructured and sub-micron steel •Bainitic ferrite plate thickness measured by SEM and TEM •Crystallite size determined by X-ray analysis.

36 MATERIALS SCIENCE↗

Energy and mechanical properties predictions in Fe-Ni binary system by ab initio calculations

Here, the fcc random solid solutions in Fe-Ni binary system were investigated using ab initio calculations. Density functional theory (DFT) was applied along with supercells generated by the Special Quasi-Random Structure (SQS) method to predict the enthalpy of mixing and magnetic and elastic properties. The results were compared with relevant experimental and computational data produced by other groups. The effect of SQS supercell size and relaxation procedure on the predicted properties of the alloys were also studied. Our calculations show that 16-atoms SQS model provides adequate approach to the mixing energy of the alloy. However, calculation of mechanical properties requires larger supercell. Excellent agreement between experimental and calculated elastic properties was found in fcc region of Fe-Ni system.

36 MATERIALS SCIENCE↗

Probabilistic evolution of stochastic dynamical systems: A meso-scale perspective

Stochastic dynamical systems arise naturally across nearly all areas of science and engineering. Typically, a dynamical system model is based on some prior knowledge about the underlying dynamics of interest in which probabilistic features are used to quantify and propagate uncertainties associated with the initial conditions, external excitations, etc. From a probabilistic modeling standing point, two broad classes of methods exist, i.e. macro-scale methods and micro-scale methods. Classically, macro-scale methods such as statistical moments-based strategies are usually too coarse to capture the multi-mode shape or tails of a non-Gaussian distribution. Micro-scale methods such as random samples-based approaches, on the other hand, become computationally very challenging in dealing with high-dimensional stochastic systems. In view of these potential limitations, a meso-scale scheme is proposed here that utilizes a meso-scale statistical structure to describe the dynamical evolution from a probabilistic perspective. The significance of this statistical structure is twofold. First, it can be tailored to any arbitrary random space. Second, it not only maintains the probability evolution around sample trajectories but also requires fewer meso-scale components than the micro-scale samples. To demonstrate the efficacy of the proposed meso-scale scheme, a set of examples of increasing complexity are provided. Connections to the benchmark stochastic models as conservative and Markov models along with practical implementation guidelines are presented.

97 MATHEMATICS AND COMPUTING↗

Optimization of Random Phase Approximation Calculations for Improved Energies of Molecules, Solids, and Surfaces

We present an optimized random phase approximation method (optRPA26) that significantly improves upon conventional RPA through an optimized choice of reference orbitals and energy components, rather than a modification of the RPA correlation functional itself. The method employs an empirically constructed hybrid functional to generate DFT orbitals to evaluate the RPA correlation energy, which is then scaled by a constant. Comprehensive benchmarks across molecules, bulk solids, and surface systems demonstrate that optRPA26 consistently achieves high accuracy, with mean absolute errors of 0.05 eV for W4-11-RE reaction energies, 0.07 eV for cohesive energies, 0.09 eV for metal oxide formation energies, 0.11–0.12 eV for adsorption of small molecules on metals, and 0.06 eV for adsorption on oxides. In addition, optRPA26 correctly captures phase stability in metal oxides and magnetic metals. The optRPA26 approach can be run using standard RPA implementations, highlighting its potential as a general-purpose reference method that can accurately capture covalent, ionic, metallic, and van der Waals bonding in molecules, solids, and interfaces.

Adsorption↗