Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distance Based Approximation”

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 55 records · Page 3

TU-F-CAMPUS-T-03: Commissioning of a Custom Fetal Lead Shield for Radiotherapy

Purpose: To evaluate the peripheral dose (PD) to a fetus during radiation therapy of pregnant patients when using a newly designed fetal lead shield (FLS). Methods: A custom FLS has been designed and fabricated for our department. The FLS (1.1 TVLs for 6 MV) is mounted on a mobile frame and can be adjusted vertically with a motor actuator. PD measurements were acquired for multiple simple square fields and for a variety of potential treatment sites a pregnant patient may be treated for including brain, head and neck (H&N) and thorax. For measurements of the brain, H&N, and thorax, an ionization chamber and OSLDs were positioned on average at a distance of 48, 29 and 26 cm, respectively, from the edge of treatment fields to mimic the approximate position of the fundus. Results: Based on our measurements, applying a 90° collimator rotation and using tertiary MLCs to define the field aperture in combination with jaws resulted in an average dose reduction of 60%. When using these planning strategies in combination with the FLS, on average, the PD was reduced by additional 25% for simple square fields and 20% for clinical plans. Conclusion: The custom FLS is a safe, effective, and relatively easy system to position. Commissioning measurements have demonstrated that the PD to the fetus can be significantly reduced when using the FLS. The comprehensive dataset obviates the need for individual patient pre-treatment dose measurements as long as the geometry falls within the commissioning limits.

60 APPLIED LIFE SCIENCES↗

Gleipnir: toward practical error analysis for Quantum programs

Practical error analysis is essential for the design, optimization, and evaluation of Noisy Intermediate-Scale Quantum(NISQ) computing. However, bounding errors in quantum programs is a grand challenge, because the effects of quantum errors depend on exponentially large quantum states. In this work, we present Gleipnir, a novel methodology toward practically computing verified error bounds in quantum programs. Gleipnir introduces the (ρ,δ)-diamond norm, an error metric constrained by a quantum predicate consisting of the approximate state ρ and its distance δ to the ideal state ρ. This predicate (ρ,δ) can be computed adaptively using tensor networks based on the Matrix Product States. Gleipnir features a lightweight logic for reasoning about error bounds in noisy quantum programs, based on the (ρ,δ)-diamond norm metric. Furthermore, our experimental results show that Gleipnir is able to efficiently generate tight error bounds for real-world quantum programs with 10 to 100 qubits, and can be used to evaluate the error mitigation performance of quantum compiler transformations.

Tao, Runzhou↗

Dark Energy Survey Year 1 results: measurement of the baryon acoustic oscillation scale in the distribution of galaxies to redshift 1

ABSTRACT We present angular diameter distance measurements obtained by locating the baryon acoustic oscillations (BAO) scale in the distribution of galaxies selected from the first year of Dark Energy Survey data. We consider a sample of over 1.3 million galaxies distributed over a footprint of 1336 deg2 with 0.6 < $z$photo < 1 and a typical redshift uncertainty of 0.03(1 + $z$). This sample was selected, as fully described in a companion paper, using a colour/magnitude selection that optimizes trade-offs between number density and redshift uncertainty. We investigate the BAO signal in the projected clustering using three conventions, the angular separation, the comoving transverse separation, and spherical harmonics. Further, we compare results obtained from template-based and machine-learning photometric redshift determinations. We use 1800 simulations that approximate our sample in order to produce covariance matrices and allow us to validate our distance scale measurement methodology. We measure the angular diameter distance, DA, at the effective redshift of our sample divided by the true physical scale of the BAO feature, rd. We obtain close to a 4 per cent distance measurement of DA($z$eff = 0.81)/rd = 10.75 ± 0.43. These results are consistent with the flat Λ cold dark matter concordance cosmological model supported by numerous other recent experimental results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Information-theoretic discrepancy based iterative reconstructions (IDIR) for polychromatic x-ray tomography

Purpose: X-ray photons generated from a typical x-ray source for clinical applications exhibit a broad range of wavelengths, and the interactions between individual particles and biological substances depend on particles' energy levels. Most existing reconstruction methods for transmission tomography, however, neglect this polychromatic nature of measurements and rely on the monochromatic approximation. In this study, we developed a new family of iterative methods that incorporates the exact polychromatic model into tomographic image recovery, which improves the accuracy and quality of reconstruction.Methods: The generalized information-theoretic discrepancy (GID) was employed as a new metric for quantifying the distance between the measured and synthetic data. By using special features of the GID, the objective function for polychromatic reconstruction which contains a double integral over the wavelength and the trajectory of incident x-rays was simplified to a paraboloidal form without using the monochromatic approximation. More specifically, the original GID was replaced with a surrogate function with two auxiliary, energy-dependent variables. Subsequently, the alternating minimization technique was applied to solve the double minimization problem. Based on the optimization transfer principle, the objective function was further simplified to the paraboloidal equation, which leads to a closed-form update formula. Numerical experiments on the beam-hardening correction and material-selective reconstruction were conducted to compare and assess the performance of conventional methods and the proposed algorithms.Results: The authors found that the GID determines the distance between its two arguments in a flexible manner. In this study, three groups of GIDs with distinct data representations were considered. The authors demonstrated that one type of GIDs that comprises “raw” data can be viewed as an extension of existing statistical reconstructions; under a particular condition, the GID is equivalent to the Poisson log-likelihood function. The newly proposed GIDs of the other two categories consist of log-transformed measurements, which have the advantage of imposing linearized penalties over multiple discrepancies. For all proposed variants of the GID, the aforementioned strategy was used to obtain a closed-form update equation. Even though it is based on the exact polychromatic model, the derived algorithm bears a structural resemblance to conventional methods based on the monochromatic approximation. The authors named the proposed approach as information-theoretic discrepancy based iterative reconstructions (IDIR). In numerical experiments, IDIR with raw data converged faster than previously known statistical reconstruction methods. IDIR with log-transformed data exhibited superior reconstruction quality and faster convergence speed compared with conventional methods and their variants.Conclusions: The authors' new framework for tomographic reconstruction allows iterative inversion of the polychromatic data model. The primary departure from the traditional iterative reconstruction was the employment of the GID as a new metric for quantifying the inconsistency between the measured and synthetic data. The proposed methods outperformed not only conventional methods based on the monochromatic approximation but also those based on the polychromatic model. The authors have observed that the GID is a very flexible means to design an objective function for iterative reconstructions. Hence, the authors expect that the proposed IDIR framework will also be applicable to other challenging tasks.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Adaptive Hybridizable Discontinuous Galerkin discretization of the Grad–Shafranov equation by extension from polygonal subdomains

In this work, we propose a high-order adaptive numerical solver for the semilinear elliptic boundary value problem modeling magnetic plasma equilibrium in axisymmetric confinement devices. In the fixed boundary case, the equation is posed on curved domains with piecewise smooth curved boundaries that may present corners. The solution method we present is based on the hybridizable discontinuous Galerkin method and sidesteps the need for geometry-conforming triangulations thanks to a transfer technique that allows to approximate the solution using only a polygonal subset as computational domain. Moreover, the solver features automatic mesh refinement driven by a residual-based a posteriori error estimator. As the mesh is locally refined, the computational domain is automatically updated in order to always maintain the distance between the actual boundary and the computational boundary of the order of the local mesh diameter. Numerical evidence is presented of the suitability of the estimator as an approximate error measure for physically relevant equilibria with pressure pedestals, internal transport barriers, and current holes on realistic geometries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measurement of θ 13 in Double Chooz using neutron captures on hydrogen with novel background rejection techniques

We observed a measurement of the Double Chooz collaboration and the neutrino mixing angle θ 13 using reactor $\bar{v}$ e via the inverse beta decay reaction in which the neutron is captured on hydrogen. Our measurement is based on 462.72 live days data, approximately twice as much data as in the previous such analysis, collected with a detector positioned at an average distance of 1050 m from two reactor cores. Several novel techniques have been developed to achieve significant reductions of the backgrounds and systematic uncertainties. Accidental coincidences, the dominant background in this analysis, are suppressed by more than an order of magnitude with respect to our previous publication by a multi-variate analysis. Furthermore, these improvements demonstrate the capability of precise measurement of reactor $\bar{v}$ e without gadolinium loading. Spectral distortions from the $\bar{v}$ e reactor flux predictions previously reported with the neutron capture on gadolinium events are confirmed in the independent data sample presented here. A value of sin 2 2θ 13 = 0.095 0.039 +0.038 (stat+syst) is obtained from a fit to the observed event rate as a function of the reactor power, a method insensitive to the energy spectrum shape. A simultaneous fit of the hydrogen capture events and of the gadolinium capture events yields a measurement of sin 2 2θ 13 = 0.088 ± 0.033(stat+syst).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development of a machine-learning-based ionic-force correction model for quantum molecular dynamic simulations of warm dense matter

In this work Δ learning is used to map orbital-free density functional theory (OF-DFT) ionic forces to the corresponding Kohn-Sham (KS) DFT ionic forces. The development of the approximate force difference in terms of the ion positions is constructed and serves as a stand in for the ground truth force difference. Descriptor vectors for ion configurations are constructed using all distance between ions in conjunction with an indexing based on a nearest neighbor ranking. It is demonstrated that such a scheme of descriptors can uniquely describe an ionic configuration up to a rotation and reflection when no ambiguity in the nearest neighbor ranking exists. How to handle the case when an ambiguity exists in the nearest neighbor ranking is discussed. As a proof of principle, the model is trained and tested on warm dense hydrogen at temperatures between 1 and 15 eV. Once tested, the model was used to perform molecular dynamic simulations of warm dense hydrogen. Furthermore, the resulting energies and pressures are within 1% and 2% of their respective target KS values.

36 MATERIALS SCIENCE↗

Approximate expression for the electric potential around an absorbing particle in isotropic collisionless plasma

A new approximate expression for the potential distribution around an absorbing particle in isotropic collisionless plasma is proposed. The approximate expression is given by the sum of the Debye-Hückel potential with an effective screening length and the far-field asymptote obtained from the solution of the linearized Poisson equation. In contrast to analogous models, the effective screening length is not fixed but depends on the distance from the particle. This allows us to obtain a more accurate approximation for the potential distribution in the entire range of distances. The dependence of the screening length on the distance is predicted from the analysis of the charge density distribution function. This dependence contains two adjustable parameters, which are calculated by applying the procedure based on charge balance considerations. Using the obtained results, simple expressions for the parameters of the model are proposed. In addition, a simple expression for the characteristic screening length, which can be used to approximate the potential distribution near the particle, is obtained. The developed model potential is shown to be in excellent agreement with the solution of the nonlinear Poisson equation for typical conditions used in experiments with complex plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Conforming window functions for meshfree methods

Window functions provide a base for the construction of approximation functions in many meshfree methods. They control the smoothness and extent of the approximation functions and are commonly defined using Euclidean distances which helps eliminate the need for a meshed discretization, simplifying model development for some classes of problems. However, for problems with complicated geometries such as nonconvex or multi-body domains, poor solution accuracy and convergence can occur unless the extents of the window functions, and thus approximation functions, are carefully controlled, often a time consuming or intractable task. Here in this paper, we present a method to provide more control in window function design, allowing efficient and systematic handling of complex geometries. “Conforming” window functions are constructed using Bernstein–Bézier splines defined on local triangulations with constraints imposed to control smoothness. Graph distances are used in conjunction with Euclidean metrics to provide adequate information for shaping the window functions. The conforming window functions are demonstrated using the Reproducing Kernel Particle Method showing improved accuracy and convergence rates for problems with challenging geometries. Conforming window functions are also demonstrated as a means to simplify the imposition of essential boundary conditions.

, rapid design-to-analysis↗

Determination of lung segments in computed tomography images using the Euclidean distance to the pulmonary artery

Purpose: Computed tomography (CT) imaging is the modality of choice for lung cancer diagnostics. With the increasing number of lung interventions on sublobar level in recent years, determining and visualizing pulmonary segments in CT images and, in oncological cases, reliable segment-related information about the location of tumors has become increasingly desirable. Computer-assisted identification of lung segments in CT images is subject of this work.Methods: The authors present a new interactive approach for the segmentation of lung segments that uses the Euclidean distance of each point in the lung to the segmental branches of the pulmonary artery. The aim is to analyze the potential of the method. Detailed manual pulmonary artery segmentations are used to achieve the best possible segment approximation results. A detailed description of the method and its evaluation on 11 CT scans from clinical routine are given.Results: An accuracy of 2–3 mm is measured for the segment boundaries computed by the pulmonary artery-based method. On average, maximum deviations of 8 mm are observed. 135 intersegmental pulmonary veins detected in the 11 test CT scans serve as reference data. Furthermore, a comparison of the presented pulmonary artery-based approach to a similar approach that uses the Euclidean distance to the segmental branches of the bronchial tree is presented. It shows a significantly higher accuracy for the pulmonary artery-based approach in lung regions at least 30 mm distal to the lung hilum.Conclusions: A pulmonary artery-based determination of lung segments in CT images is promising. In the tests, the pulmonary artery-based determination has been shown to be superior to the bronchial tree-based determination. The suitability of the segment approximation method for application in the planning of segment resections in clinical practice has already been verified in experimental cases. However, automation of the method accompanied by an evaluation on a larger number of test cases is required before application in the daily clinical routine.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

An improved hyperbolic embedding algorithm

Because hyperbolic space has properties that make it amenable to graph representations, there is significant interest in scalable hyperbolic-space embedding methods. These embeddings enable constant-time approximation of shortest-path distances, and so are significantly more efficient than full shortest-path computations. In this article, we improve on existing landmark-based hyperbolic embedding algorithms for large-scale graphs. Whereas previous methods compute the embedding by using the derivative-free Nelder–Mead simplex optimization method, our approach uses the limited-memory BFGS (LBFGS) method, which is quasi-Newton optimization, with analytic gradients. Our method is not only significantly faster but also produces higher-quality embeddings. Moreover, we are able to include the hyperbolic curvature as a variable in the optimization. We compare our hyperbolic embedding method implementation in Python (called Hypy) against the best publicly available software, Rigel. Our method is an order of magnitude faster and shows significant improvements in the accuracy of the shortest-path distance calculations. Furthermore, tests are performed on a variety of real-world networks, and we show the scalability of our method by embedding a graph with 1.8 billion edges and 65 million nodes.

97 MATHEMATICS AND COMPUTING↗

Investigation on collisions of filament pairs in dielectric barrier discharge

Collisions of filament pairs in a hexagonal superlattice pattern in dielectric barrier discharge are investigated on different timescales. In the evolution of the pattern, the space scale of each hexagon cell decreases with the increasing voltage. The duration of one collision is seven half voltage cycles at least. Two stable orientations of a pair are approximately perpendicular to each other and the orientational changes occurring during the entire colliding process should be a multiple of 30°. The time interval between two consecutive collisions decreases with the increasing voltage. The distance between the paired spots decreases nonmonotonically. Based on the discharge order of the pattern, it is inferred that the collision should be the interaction between a discharging filament and the surface charges deposited by another discharged filament, and the nonmonotonic decrease of distance D is explained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding the many-body expansion for large systems. III. Critical role of four-body terms, counterpoise corrections, and cutoffs

Papers I and II in this series [R. M. Richard et al., J. Chem. Phys. 141, 014108 (2014); K. U. Lao et al., ibid. 144, 164105 (2016)] have attempted to shed light on precision and accuracy issues affecting the many-body expansion (MBE), which only manifest in larger systems and thus have received scant attention in the literature. Many-body counterpoise (CP) corrections are shown to accelerate convergence of the MBE, which otherwise suffers from a mismatch between how basis-set superposition error affects subsystem versus supersystem calculations. In water clusters ranging in size up to (H2O)37, four-body terms prove necessary to achieve accurate results for both total interaction energies and relative isomer energies, but the sheer number of tetramers makes the use of cutoff schemes essential. To predict relative energies of (H2O)20 isomers, two approximations based on a lower level of theory are introduced and an ONIOM-type procedure is found to be very well converged with respect to the appropriate MBE benchmark, namely, a CP-corrected supersystem calculation at the same level of theory. Results using an energy-based cutoff scheme suggest that if reasonable approximations to the subsystem energies are available (based on classical multipoles, say), then the number of requisite subsystem calculations can be reduced even more dramatically than when distance-based thresholds are employed. The end result is several accurate four-body methods that do not require charge embedding, and which are stable in large basis sets such as aug-cc-pVTZ that have sometimes proven problematic for fragment-based quantum chemistry methods. Even with aggressive thresholding, however, the four-body approach at the self-consistent field level still requires roughly ten times more processors to outmatch the performance of the corresponding supersystem calculation, in test cases involving 1500–1800 basis functions.

74 ATOMIC AND MOLECULAR PHYSICS↗

Nesting habitat use by river chubs in a hydrologically variable Appalachian tailwater

Hydrologic alteration continues to affect aquatic biodiversity asknowledge of the spawning requirements of fishes, especially keystone or foundation species, becomes more critical for conservation and management. Our objectives are to quantify the spawning micro- and mesohabitat use of river chub Nocomis micropogon, a gravel mound nesting minnow, in a hydrologically regulated river in North Carolina, USA. At the microhabitat scale, substrate sizes on nests were compared with pebble counts in 1-m2 adjacent quadrats. Average depths and current velocities at nests were compared with measurements from paired transects. At the mesohabitat scale, generalised linear mixed models (GLMMs) were used to identify the importance of average bed slope, average depth and percentages of rock outcrops (a measure of flow heterogeneity and velocity shelters) for predicting nest presence and abundance. To relate nesting activities to hydrologic alteration from dam operation, nest dimensions were measured before and after a scheduled discharge event approximately six times that of base flow. In addition, linear regression was used to predict changes in the use of flow refugia and overhead cover with increased fluvial distance from the dam. Microhabitats in which nests were placed had, on average, slower current velocities and shallower depths. Gravel diameters of nests were significantly smaller than substrate particles adjacent to nests. GLMMs revealed that mesohabitats with nests were shallower, had more moderate slopes and greater proportions of rock outcrops than mesohabitats without nests. Finally, the scheduled discharge event significantly flattened nests. Near the dam, nests were built in close proximity ( 2 m) to velocity shelters; this relationship diminished with distance from the dam. River chubs are spawning habitat specialists. Because multiple species rely on river chub nests for reproduction and food, the needs of this species should be considered when managing instream flows.

58 GEOSCIENCES↗

A Case Study on Pathogen Transport, Deposition, Evaporation and Transmission: Linking High-Fidelity Computational Fluid Dynamics Simulations to Probability of Infection

A high-fidelity, low-Mach computational fluid dynamics simulation tool that includes evaporating droplets and variable-density turbulent flow coupling is well-suited to ascertain transmission probability and supports risk mitigation methods development for airborne infectious diseases such as COVID-19. A multi-physics large-eddy simulation-based paradigm is used to explore droplet and aerosol pathogen transport from a synthetic cough emanating from a kneeling humanoid. For an outdoor configuration that mimics the recent open-space social distance strategy of San Francisco, maximum primary droplet deposition distances are shown to approach 8.1 m in a moderate wind configuration with the aerosol plume transported in excess of 15 m. In quiescent conditions, the aerosol plume extends to approximately 4 m before the emanating pulsed jet becomes neutrally buoyant. A dose–response model, which is based on previous SARS coronavirus (SARS-CoV) data, is exercised on the high-fidelity aerosol transport database to establish relative risk at eighteen virtual receptor probe locations.

59 BASIC BIOLOGICAL SCIENCES↗

Effect of Fe Doping and O Vacancies on the Magnetic Properties of Rutile TiO{sub 2}

Magnetism sources of Fe doping and O vacancies, which coexist in the presence of rutile TiO{sub 2}, are controversial and require a resolution. The effects of Fe doping and point defect on the magnetism of rutile TiO{sub 2} were studied using geometry optimization and energy calculation based on the first principle generalized gradient approximation + U method (GGA+U) of density functional theory. Fe/Vo ratio is 1:1 in rutile TiO{sub 2}. The next closest distance between Fe and Vo is 0.04082 nm. The formation energy of the doping system is the smallest. This doping system has the highest stability. The coexistence of Fe doping and O vacancy achieved ferromagnetic long-range order. The different ratios of Fe/Vo are 1:1, 1:2, 2:1, and 2:2. The magnetic properties of the doping system are significantly different. At a ratio of 2:1, the magnetic moment is the largest among all doping systems. Moreover, the Curie temperature of the doping system can be higher than room temperature. These findings can obtain dilute magnetic semiconductors with potential application value. At a Fe/V{sub O} ratio of 2:2, the Ti{sub 14}Fe{sub 2}O{sub 30} system produces a half-metallic property of 100% electron spin polarizability. This is the most valuable application for designing and preparing novel spintronic injection source dilute magnetic semiconductors. The magnetic source is mainly mediated by the holes generated by Fe doping and O vacancy complexes, thereby causing the spin polarization double exchange effect on the electron of Ti–3d orbital atoms near the O vacancy, O–2p electron orbital, and Fe–3d electron orbital. This feature is consistent with the following theories: the mean field approximation and the double exchange mechanism theories.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Locality and Digital Quantum Simulation of Power-Law Interactions

The propagation of information in nonrelativistic quantum systems obeys a speed limit known as a Lieb-Robinson bound. We derive a new Lieb-Robinson bound for systems with interactions that decay with distance r as a power law, 1 / r α . The bound implies an effective light cone tighter than all previous bounds. Our approach is based on a technique for approximating the time evolution of a system, which was first introduced as part of a quantum simulation algorithm by Haah et al ., FOCS’18. To bound the error of the approximation, we use a known Lieb-Robinson bound that is weaker than the bound we establish. This result brings the analysis full circle, suggesting a deep connection between Lieb-Robinson bounds and digital quantum simulation. In addition to the new Lieb-Robinson bound, our analysis also gives an error bound for the Haah et al. quantum simulation algorithm when used to simulate power-law decaying interactions. In particular, we show that the gate count of the algorithm scales with the system size better than existing algorithms when α > 3 D (where D is the number of dimensions).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interaction of macroparticles localized in Wigner–Seitz cells of various types of cubic lattices in an equilibrium plasma

The interaction of two charged point macroparticles located in Wigner–Seitz cells of simple cubic (SC), body-centered cubic (BCC), or face-centered cubic (FCC) lattices in an equilibrium plasma has been studied within the Debye approximation or, more specifically, based on the linearized Poisson–Boltzmann model. The shape of the outer boundary is shown to exert a strong influence on the pattern of electrostatic interaction between the two macroparticles, which transforms from repulsion at small interparticle distances to attraction as the interparticle distance approaches half the length of the computational cell. The macroparticle pair interaction potential in an equilibrium plasma is shown to be nevertheless the Debye one and purely repulsive for likely charged macroparticles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗