Search NASA⌕ Search

SEARCH · Search NASA

Results for “Soft computing”

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 199 records · Page 11

Fast equilibration protocol for million atom systems of highly entangled linear polyethylene chains

Equilibrated systems of entangled polymer melts cannot be produced using direct brute force equilibration due to the slow reptation dynamics exhibited by high molecular weight chains. Instead, these dense systems are produced using computational techniques such as Monte Carlo-Molecular Dynamics hybrid algorithms, though the use of soft potentials has also shown promise mainly for coarse-grained polymeric systems. Through the use of soft-potentials, the melt can be equilibrated via molecular dynamics at intermediate and long length scales prior to switching to a Lennard-Jones potential. We will outline two different equilibration protocols, which use various degrees of information to produce the starting configurations. In one protocol, we use only the equilibrium bond angle, bond length, and target density during the construction of the simulation cell, where the information is obtained from available experimental data and extracted from the force field without performing any prior simulation. In the second protocol, we moreover utilize the equilibrium radial distribution function and dihedral angle distribution. This information can be obtained from experimental data or from a simulation of short unentangled chains. Both methods can be used to prepare equilibrated and highly entangled systems, but the second protocol is much more computationally efficient. These systems can be strictly monodisperse or optionally polydisperse depending on the starting chain distribution. Our protocols, which utilize a soft-core harmonic potential, will be applied for the first time to equilibrate a million particle system of polyethylene chains consisting of 1000 united atoms at various temperatures. Calculations of structural and entanglement properties demonstrate that this method can be used as an alternative towards the generation of entangled equilibrium structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New-Sum: A Novel Online ABFT Scheme For General Iterative Methods

Emerging high-performance computing platforms, with large component counts and lower power margins, are anticipated to be more susceptible to soft errors in both logic circuits and memory subsystems. We present an online algorithm-based fault tolerance (ABFT) approach to efficiently detect and recover soft errors for general iterative methods. We design a novel checksum-based encoding scheme for matrix-vector multiplication that is resilient to both arithmetic and memory errors. Our design decouples the checksum updating process from the actual computation, and allows adaptive checksum overhead control. Building on this new encoding mechanism, we propose two online ABFT designs that can effectively recover from errors when combined with a checkpoint/rollback scheme.

Instantaneous Footprint↗

Probing many-body localization on a noisy quantum computer

A disordered quantum system of interacting particles exhibits localized behavior when the disorder is large compared to the interaction strength. Studying this phenomenon on a quantum computer with no, or limited, error correction is challenging because even weak coupling to a thermal environment destroys most signatures of localization. Fortunately, spectral functions of local operators are known to contain features that can survive the presence of noise. In these spectra, discrete peaks and a soft gap at low frequencies compared to the thermal phase indicate localization. Here, we present the computation of spectral functions on a trapped-ion quantum computer for a one-dimensional Heisenberg model with disorder. Further, we design an error-mitigation technique which is effective at removing the noise from the measurement allowing clear signatures of localization to emerge as the disorder increases. Furthermore, we show that spectral functions can serve as a robust and scalable diagnostic of many-body localization on current and future generations of quantum computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Color-octet scalars in Dirac gaugino models with broken R symmetry

In this work we study the collider phenomenology of color-octet scalars ( sgluons ) in supersymmetric models with Dirac gaugino masses that feature an explicitly broken R symmetry ( R - broken models ). We construct such models by augmenting minimal R -symmetric models with a fairly general set of supersymmetric and softly supersymmetry-breaking operators that explicitly break R symmetry. We then compute the rates of all significant two-body decays and highlight new features that appear as a result of R symmetry breaking, including enhancements to extant decay rates, novel tree- and loop-level decays, and improved cross sections of single sgluon production. We demonstrate in some detail how the familiar results from minimal R -symmetric models can be obtained by restoring R symmetry. In parallel to this discussion, we explore constraints on these models from the Large Hadron Collider. We find that, in general, R symmetry breaking quantitatively affects existing limits on color-octet scalars, perhaps closing loopholes for light CP-odd (pseudoscalar) sgluons while opening one for a light CP-even (scalar) particle. Qualitatively, however, we find that — much as for minimal R -symmetric models, despite stark differences in phenomenology — scenarios with broken R symmetry and two sgluons below the TeV scale can be accommodated by existing searches.

79 ASTRONOMY AND ASTROPHYSICS↗

Organic Matrix Derived from Host–Microbe Interplay Contributes to Pathological Renal Biomineralization

Matrix stones are a rare form of kidney stones. They feature a high percentage of hydrogel-like organic matter, and their formation is closely associated with urinary tract infections. Herein, comprehensive materials and biochemical approaches were taken to map the organic-inorganic interface and gather insights into the host-microbe interplay in pathological renal biomineralization. Surgically extracted soft and slimy matrix stones were examined using micro-X-ray computed tomography and various microspectroscopy techniques. Higher-mineral-density laminae were positive for calcium-bound Alizarin red. Lower-mineral-density laminae revealed periodic acid-Schiff-positive organic filamentous networks of varied thickness. These organic filamentous networks, which featured a high polysaccharide content, were enriched with zinc, carbon, and sulfur elements. Neutrophil extracellular traps (NETs) along with immune response-related proteins, including calprotectin, myeloperoxidase, CD63, and CD86, also were identified in the filamentous networks. Expressions of NETs and upregulation of polysaccharide-rich mucin secretion are proposed as a part of the host immune defense to "trap" pathogens. These host-microbe derived organic matrices can facilitate heterogeneous nucleation and precipitation of inorganic particulates, resulting in macroscale aggregates known as "matrix stones". These insights into the plausible aggregation of constituents through host-microbe interplay underscore the unique "double-edged sword" effect of the host immune response to pathogens and the resulting renal biominerals.

59 BASIC BIOLOGICAL SCIENCES↗

Recommended direct simulation Monte Carlo collision model parameters for modeling ionized air transport processes

A systematic approach for calibrating the direct simulation Monte Carlo (DSMC) collision model parameters to achieve consistency in the transport processes is presented. The DSMC collision cross section model parameters are calibrated for high temperature atmospheric conditions by matching the collision integrals from DSMC against ab initio based collision integrals that are currently employed in the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) and Data Parallel Line Relaxation (DPLR) high temperature computational fluid dynamics solvers. The DSMC parameter values are computed for the widely used Variable Hard Sphere (VHS) and the Variable Soft Sphere (VSS) models using the collision-specific pairing approach. The recommended best-fit VHS/VSS parameter values are provided over a temperature range of 1000-20 000 K for a thirteen-species ionized air mixture. Use of the VSS model is necessary to achieve consistency in transport processes of ionized gases. The agreement of the VSS model transport properties with the transport properties as determined by the ab initio collision integral fits was found to be within 6% in the entire temperature range, regardless of the composition of the mixture. The recommended model parameter values can be readily applied to any gas mixture involving binary collisional interactions between the chemical species presented for the specified temperature range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Code-to-code comparison and validation of the radiation-hydrodynamics capabilities of the FLASH code using a laboratory astrophysical jet

The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and behavior of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities in the publicly available FLASH radiation-hydrodynamics code, we present detailed code-to-code comparisons between FLASH and the HYDRA code developed at Lawrence Livermore National Laboratory using previously published HYDRA simulations from Grava et al. [Phys. Rev. E 78, 016403 (2008)]. That study describes a laser experiment that produced a jet-like feature that the authors compare to astrophysical jets. Importantly, the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment included detailed x-ray interferometric measurements of electron number densities and a time-integrated measurement of the soft x-ray spectrum. Despite markedly different methods for treating the computational mesh, and different equations of state and opacity models, the FLASH results resemble the results from HYDRA and, most importantly, the experimental measurements of electron density. Having validated the FLASH code in this way, we use the code to further investigate and understand the formation of the jet seen in the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment and discuss its relation to the Wan et al. [Phys. Rev. E 55, 6293 (1997)] experiment at the NOVA laser.

Orban, Chris (ORCID:0000000312004538)↗

Bosonic Superfluid on the Lowest Landau Level

We develop a low-energy effective field theory of a two-dimensional bosonic superfluid on the lowest Landau level at zero temperature and identify a Berry term that governs the dynamics of coarse-grained superfluid degrees of freedom. For an infinite vortex crystal we compute how the Berry term affects the low-energy spectrum of soft collective Tkachenko oscillations and nondissipative Hall responses of the particle number current and stress tensor. This term gives rise to a quadratic in momentum term in the Hall conductivity, but does not generate a nondissipative Hall viscosity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

BonVoision: Leveraging Spatial Data Smoothness For Recovery From Memory Soft Errors

The increasing soft error rates in memory systems raise an emerging concern for modern computing systems. As a result, detectable but uncorrectable errors (DUEs) become potentially more frequent and affect HPC applications. Today, upon encountering a DUE, HPC applications crash, incurring significant performance, storage, and energy overheads. In this paper, we propose a technique to continue application execution past a DUE through the repair of the corrupted memory data by leveraging spatial data smoothness. We present BonVoision, a run-time system that intercepts DUE events, analyzes the binary to identify data elements in the structural neighborhood of the event, and fixes the corrupted data elements by interpolating from the values in their neighborhood. Our evaluation demonstrates that BonVoision incurs negligible overhead and outperforms other recovery strategies by a factor of 2×, on average. We demonstrate that BonVoision also improves the efficiency of existing checkpointing/restart schemes by approximately increasing the optimal checkpoint interval by 23%.

Fang, Bo↗

Pileup and Infrared Radiation Annihilation (PIRANHA): a paradigm for continuous jet grooming

Jet grooming is an important strategy for analyzing relativistic particle collisions in the presence of contaminating radiation. Most jet grooming techniques introduce hard cutoffs to remove soft radiation, leading to discontinuous behavior and associated experimental and theoretical challenges. In this paper, we introduce Pileup and Infrared Radiation Annihilation (Piranha), a paradigm for continuous jet grooming that overcomes the discontinuity and infrared sensitivity of hard-cutoff grooming procedures. We motivate Piranha from the perspective of optimal transport and the Energy Mover’s Distance and review Apollonius Subtraction and Iterated Voronoi Subtraction as examples of Piranha-style grooming. We then introduce a new tree-based implementation of Piranha, Recursive Subtraction, with reduced computational costs. Finally, we demonstrate the performance of Recursive Subtraction in mitigating sensitivity to soft distortions from hadronization and detector effects, and additive contamination from pileup and the underlying event.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Physics-based shape matching for intraoperative image guidance

Purpose: Soft-tissue deformations can severely degrade the validity of preoperative planning data during computer assisted interventions. Intraoperative imaging such as stereo endoscopic, time-of-flight or, laser range scanner data can be used to compensate these movements. In this context, the intraoperative surface has to be matched to the preoperative model. The shape matching is especially challenging in the intraoperative setting due to noisy sensor data, only partially visible surfaces, ambiguous shape descriptors, and real-time requirements. Methods: A novel physics-based shape matching (PBSM) approach to register intraoperatively acquired surface meshes to preoperative planning data is proposed. The key idea of the method is to describe the nonrigid registration process as an electrostatic–elastic problem, where an elastic body (preoperative model) that is electrically charged slides into an oppositely charged rigid shape (intraoperative surface). It is shown that the corresponding energy functional can be efficiently solved using the finite element (FE) method. It is also demonstrated how PBSM can be combined with rigid registration schemes for robust nonrigid registration of arbitrarily aligned surfaces. Furthermore, it is shown how the approach can be combined with landmark based methods and outline its application to image guidance in laparoscopic interventions. Results: A profound analysis of the PBSM scheme based on in silico and phantom data is presented. Simulation studies on several liver models show that the approach is robust to the initial rigid registration and to parameter variations. The studies also reveal that the method achieves submillimeter registration accuracy (mean error between 0.32 and 0.46 mm). An unoptimized, single core implementation of the approach achieves near real-time performance (2 TPS, 7–19 s total registration time). It outperforms established methods in terms of speed and accuracy. Furthermore, it is shown that the method is able to accurately match partial surfaces. Finally, a phantom experiment demonstrates how the method can be combined with stereo endoscopic imaging to provide nonrigid registration during laparoscopic interventions. Conclusions: The PBSM approach for surface matching is fast, robust, and accurate. As the technique is based on a preoperative volumetric FE model, it naturally recovers the position of volumetric structures (e.g., tumors and vessels). It cannot only be used to recover soft-tissue deformations from intraoperative surface models but can also be combined with landmark data from volumetric imaging. In addition to applications in laparoscopic surgery, the method might prove useful in other areas that require soft-tissue registration from sparse intraoperative sensor data (e.g., radiation therapy)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Anticipating gelation and vitrification with medium amplitude parallel superposition (MAPS) rheology and artificial neural networks

Abstract Anticipating qualitative changes in the rheological response of complex fluids (e.g., a gelation or vitrification transition) is an important capability for processing operations that utilize such materials in real-world environments. One class of complex fluids that exhibits distinct rheological states are soft glassy materials such as colloidal gels and clay dispersions, which can be well characterized by the soft glassy rheology (SGR) model. We first solve the model equations for the time-dependent, weakly nonlinear response of the SGR model. With this analytical solution, we show that the weak nonlinearities measured via medium amplitude parallel superposition (MAPS) rheology can be used to anticipate the rheological aging transitions in the linear response of soft glassy materials. This is a rheological version of a technique called structural health monitoring used widely in civil and aerospace engineering. We design and train artificial neural networks (ANNs) that are capable of quickly inferring the parameters of the SGR model from the results of sequential MAPS experiments. The combination of these data-rich experiments and machine learning tools to provide a surrogate for computationally expensive viscoelastic constitutive equations allows for rapid experimental characterization of the rheological state of soft glassy materials. We apply this technique to an aging dispersion of Laponite ® clay particles approaching the gel point and demonstrate that a trained ANN can provide real-time detection of transitions in the nonlinear response well in advance of incipient changes in the linear viscoelastic response of the system.

Lennon, Kyle R. (ORCID:0000000212515461)↗

Bulk viscous corrections to screening and damping in QCD at high temperatures

Non-equilibrium corrections to the distribution functions of quarks and gluons in a hot and dense QCD medium modify the \hard thermal loops" (HTL). The HTLs determine the retarded, advanced, and symmetric (time-ordered) propagators for gluons with soft momenta as well as the Debye screening and Landau damping mass scales. Here, we compute such corrections to a thermal as well as to a non-thermal fixed point. The screening and damping mass scales are sensitive to the bulk pressure and hence to (pseudo-) critical dynamical scaling of the bulk viscosity in the vicinity of a second-order critical point. This could be reected in the properties of quarkonium bound states in the deconfined phase and in the dynamics of soft gluon fields.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Explicit soft supersymmetry breaking in the heterotic M-theory B - L MSSM

The strongly coupled heterotic M-theory vacuum for both the observable and hidden sectors of the B - L MSSM theory is reviewed, including a discussion of the “bundle” constraints that both the observable sector SU(4) vector bundle and the hidden sector bundle induced from a single line bundle must satisfy. Gaugino condensation is then introduced within this context, and the hidden sector bundles that exhibit gaugino condensation are presented. The condensation scale is computed, singling out one line bundle whose associated condensation scale is low enough to be compatible with the energy scales available at the LHC. The corresponding region of Kähler moduli space where all bundle constraints are satisfied is presented. The generic form of the moduli dependent F-terms due to a gaugino superpotential — which spontaneously break N = 1 supersymmetry in this sector — is presented and then given explicitly for the unique line bundle associated with the low condensation scale. The moduli-dependent coefficients for each of the gaugino and scalar field soft supersymmetry breaking terms are computed leading to a low-energy effective Lagrangian for the observable sector matter fields. We then show that at a large number of points in Kähler moduli space that satisfy all “bundle” constraints, these coefficients are initial conditions for the renormalization group equations which, at low energy, lead to completely realistic physics satisfying all phenomenological constraints. Finally, we show that a substantial number of these initial points also satisfy a final constraint arising from the quadratic Higgs-Higgs conjugate soft supersymmetry breaking term.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Algorithm-Based Fault Tolerance for Convolutional Neural Networks

Convolutional neural networks (CNNs) are becoming more and more important for solving challenging and critical problems in many fields. CNN inference applications have been deployed in safety-critical systems, which may suffer from soft errors caused by high-energy particles, high temperature, or abnormal voltage. Of critical importance is ensuring the stability of the CNN inference process against soft errors. Traditional fault tolerance methods are not suitable for CNN inference because error-correcting code is unable to protect computational components, instruction duplication techniques incur high overhead, and existing algorithm-based fault tolerance (ABFT) techniques cannot protect all convolution implementations. In this paper, we focus on how to protect the CNN inference process against soft errors as efficiently as possible, with the following three contributions. (1) We propose several systematic ABFT schemes based on checksum techniques and analyze their fault protection ability and runtime thoroughly. Unlike traditional ABFT based on matrix-matrix multiplication, our schemes support any convolution implementations. (2) We design a novel workflow integrating all the proposed schemes to obtain a high detection/correction ability with limited total runtime overhead. (3) We perform our evaluation using ImageNet with well-known CNN models including AlexNet, VGG-19, ResNet-18, and YOLOv2. Here, experimental results demonstrate that our implementation can handle soft errors with very limited runtime overhead (4%~8% in both error-free and error-injected situations).

97 MATHEMATICS AND COMPUTING↗

On proton CT reconstruction using MVCT-converted virtual proton projections

Purpose: To describe a novel methodology of converting megavoltage x-ray projections into virtual proton projections that are otherwise missing due to the proton range limit. These converted virtual proton projections can be used in the reconstruction of proton computed tomography (pCT). Methods: Relations exist between proton projections and multispectral megavoltage x-ray projections for human tissue. Based on these relations, these tissues can be categorized into: (a) adipose tissue; (b) nonadipose soft tissues; and (c) bone. These three tissue categories can be visibly identified on a regular megavoltage x-ray computed tomography (MVCT) image. With an MVCT image and its projection data available, the x-ray projections through heterogeneous anatomy can be converted to the corresponding proton projections using predetermined calibration curves for individual materials, aided by a coarse segmentation on the x-ray CT image. To show the feasibility of this approach, mathematical simulations were carried out. The converted proton projections, plotted on a proton sinogram, were compared to the simulated ground truth. Proton stopping power images were reconstructed using either the virtual proton projections only or a blend of physically available proton projections and virtual proton projections that make up for those missing due to the range limit. These images were compared to a reference image reconstructed from theoretically calculated proton projections. Results: The converted virtual projections had an uncertainty of {+-}0.8% compared to the calculated ground truth. Proton stopping power images reconstructed using a blend of converted virtual projections (48%) and physically available projections (52%) had an uncertainty of {+-}0.86% compared with that reconstructed from theoretically calculated projections. Reconstruction solely from converted virtual proton projections had an uncertainty of {+-}1.1% compared with that reconstructed from theoretical projections. If these images are used for treatment planning, the average proton range uncertainty is estimated to be less than 1.5% for an imaging dose in the milligray range. Conclusions: The proposed method can be used to convert x-ray projections into virtual proton projections. The converted proton projections can be blended with existing proton projections or can be used solely for pCT reconstruction, addressing the range limit problem of pCT using current therapeutic proton machines.

60 APPLIED LIFE SCIENCES↗

Brownian Dynamics Simulation of Protein Solutions: Structural and Dynamical Properties

The study of solutions of biomacromolecules provides an important basis for understanding the behavior of many fundamental cellular processes, such as protein folding, self-assembly, biochemical reactions, and signal transduction. Here, we describe a Brownian dynamics simulation procedure and its validation for the study of the dynamic and structural properties of protein solutions. In the model used, the proteins are treated as atomically detailed rigid bodies moving in a continuum solvent. The protein-protein interaction forces are described by the sum of electrostatic interaction, electrostatic desolvation, nonpolar desolvation, and soft-core repulsion terms. The linearized Poisson-Boltzmann equation is solved to compute electrostatic terms. Simulations of homogeneous solutions of three different proteins with varying concentrations, pH, and ionic strength were performed. The results were compared to experimental data and theoretical values in terms of long-time self-diffusion coefficients, second virial coefficients, and structure factors. The results agree with the experimental trends and, in many cases, experimental values are reproduced quantitatively. There are no parameters specific to certain protein types in the interaction model, and hence the model should be applicable to the simulation of the behavior of mixtures of macromolecules in cell-like crowded environments.

59 BASIC BIOLOGICAL SCIENCES↗

Quantitative Bioluminescence Tomography-Guided Conformal Irradiation for Preclinical Radiation Research

Widely used cone beam computed tomography (CBCT)-guided irradiators in preclinical radiation research are limited to localize soft tissue target because of low imaging contrast. Knowledge of target volume is a fundamental need for radiation therapy (RT). Without such information to guide radiation, normal tissue can be overirradiated, introducing experimental uncertainties. This led us to develop high-contrast quantitative bioluminescence tomography (QBLT) for guidance. The use of a 3-dimensional bioluminescence signal, related to cell viability, for preclinical radiation research is one step toward biology-guided RT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗