Search NASA⌕ Search

SEARCH · Search NASA

Results for “inverse calculation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Improved Earthquake Source Parameters with 3D Wavespeed Models in California and Nevada

Seismic tomography harnesses earthquake data to explore the inaccessible structure of the Earth. Adjoint waveform tomography (AWT), a method of seismic tomography, updates the tomographic model by optimizing the fit between observed earthquake data and synthetic waveforms. The synthetic data are calculated by solving the wave equation through a given 3D model. An important requirement to calculating synthetics is the source information (location, centroid time, depth, and moment tensor). Errors in source information affect the quality of the synthetics produced, which in turn can limit how structure can be inferred in the AWT workflow. Here, to test the effect of updating source information, we used MTTime (Chiang, 2020), a time-domain full-waveform moment tensor inversion code, to calculate the moment tensors and depths of 118 earthquakes that occurred in California and Nevada over a 20-yr period. We calculated 3D Green’s functions using a 3D seismic wavespeed model of California and Nevada (Doody et al., 2023b). We show that the inverted solutions provide better waveform fits than the Global Centroid Moment Tensor catalog and increase usable, well-correlated data by up to 7%. Therefore, we argue that recalculating source parameters should be considered in AWT workflows, particularly for smaller magnitude events (⁠M w > 5.0).

58 GEOSCIENCES↗
DOE OSTI2021

Efficient Precoding for Single Carrier Modulation in Multi-User Massive MIMO Networks

Abstract—Frequency domain (FD) multi-user detection (MUD) and precoding has been shown to be an effective means of approaching the theoretical per-user capacity for single carrier modulation (SCM) schemes in massive MIMO scenarios with highly dispersive channels. When a cyclic prefix is added to the SCM waveform, the circulant structures of the resulting convolutional channel matrix allows for relatively simple expressions for the FD detection and precoding. In this paper, we further refine the computationally efficient minimum mean squared error (MMSE) FD-MUD technique in a time-division duplexing (TDD) massive MIMO setup by pre-computing the scale factor needed to produce an unbiased estimate. The resulting scale factor and the matrix inverse that was calculated for the uplink is then reused for the downlink to form a computationally efficient precoding scheme with superior performance compared to zero-forcing (ZF) precoding. Performance analysis is provided and confirmed through simulation results.

5G and Beyond Communications↗
DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Performance-Portable GPU Acceleration of the EFIT Tokamak Plasma Equilibrium Reconstruction Code

This paper presents the steps followed to GPU-offload parts of the core solver of EFIT-AI, an equilibrium reconstruction code suitable for tokamak experiments and burning plasmas. For this work, we will focus on the fitting procedure that consists of a Grad–Shafranov (GS) equation inverse solver that calculates equilibrium reconstructions on a grid. We will show profiling results of the original code (CPU-baseline), as well as the directives used to GPU-offload the most time-consuming function, initially to compare OpenACC and OpenMP on NVIDIA and AMD GPUs and later on to assess OpenMP performance portability on NVIDIA, AMD and Intel GPUs. We will make a performance comparison for different spatial grid sizes and show the speedup achieved on NVIDIA A100 (Perlmutter-NERSC), AMD MI250X (Frontier-OLCF) and Intel PVC GPUs (Sunspot-ALCF). Finally, we will draw some conclusions and recommendations to achieve high-performance portability for an equilibrium reconstruction code on the new HPC architectures

GPU↗
DOE OSTI2021

Removal of spurious data in Bragg coherent diffraction imaging: an algorithm for automated data preprocessing

Bragg coherent diffraction imaging (BCDI) provides a powerful tool for obtaining high-resolution structural information from nanocrystalline materials. Here a BCDI sample consisting of a large number of randomly oriented nanoscale crystals is considered. Ideally, only one crystal is oriented to produce a Bragg peak on the detector. However, diffraction from other crystals often produces additional signals on the detector. Before the measured diffraction patterns can be processed into structural images, scientists routinely need to manually identify and remove the `alien' intensities from sources other than the intended crystal. With the development of modern high-coherence storage rings, such as the upgraded Advanced Photon Source (APS), the already slow process of manual preprocessing will be untenable for the large volumes of data that will be produced. An automated method of identifying and deleting alien intensities is proposed. This method exploits the fact that BCDI of a perfect crystal produces diffraction data with inversion symmetry around the Bragg peak. This approach uses the machine learning clustering method DBSCAN to distinguish between diffraction from multiple sources, and then calculates cluster size and inversion symmetry to assess whether clusters of intensity belong to desired data or alien signals. This approach can dramatically reduce the amount of time spent manually processing data, allowing BCDI data processing capabilities to keep pace with the technological advances of fourth-generation synchrotron light sources.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Code for Calculating Regional Seismic Travel Time

The RSTT software computes predictions of the travel time of seismic energy traveling from a source to a receiver through 2.5D models of the seismic velocity distribution within the Earth. The two primary applications for the RSTT library are tomographic inversion studies and seismic event location calculations. In tomographic inversions studies, a seismologist begins with number of source-receiver travel time observations and an initial starting model of the velocity distribution within the Earth. A forward travel time calculator, such as the RSTT library, is used to compute predictions of each observed travel time and all of the residuals (observed minus predicted travel time) are calculated. The Earth model is then modified in some systematic way with the goal of minimizing the residuals. The Earth model obtained in this way is assumed to be a better model than the starting model if it has lower residuals. The other major application for the RSTT library is seismic event location. Given an Earth model, an initial estimate of the location of a seismic event, and some number of observations of seismic travel time thought to have originated from that event, location codes systematically modify the estimate of the location of the event with the goal of minimizing the difference between the observed and predicted travel times. The second application, seismic event location, is routinely implemented by the military as part of its effort to monitor the Earth for nuclear tests conducted by foreign countries.

BALLARD, SANFORD↗
DOE OSTI2025

Artificial intelligence driven laser parameter search: Inverse design of photonic surfaces using greedy surrogate-based optimization

Photonic surfaces designed with specific optical characteristics are becoming increasingly crucial for novel energy harvesting and storage systems. The design of these surfaces can be achieved by texturing materials using lasers. The optimal adjustment of laser fabrication parameters to achieve target surface optical properties is an open challenge. Thus, we develop a surrogate-based optimization approach. Our framework employs the Random Forest algorithm to model the forward relationship between the laser fabrication parameters and the resulting optical characteristics. During the optimization process, we use a greedy, prediction-based exploration strategy that iteratively selects batches of laser parameters to be used in experimentation by minimizing the predicted discrepancy between the surrogate model’s outputs and the user-defined target optical characteristics. This strategy allows for efficient identification of optimal fabrication parameters without the need to model the error landscape directly. We demonstrate the efficiency and effectiveness of our approach on two synthetic benchmarks and two specific experimental applications of photonic surface inverse design targets. By calculating the average performance of our algorithm compared to other state of the art optimization methods, we show that our algorithm performs, on average, twice as well across all benchmarks. Additionally, a warm starting inverse design technique for changed target optical characteristics enhances the performance of the introduced approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Size resolved particle volatility for 30, 60, and 90 nm particles collected at the EPCAPE Mount Soledad site from 04/23/2023 to 06/15/2023

VTDMA GF-PDF: Growth factor probability density functions from the Gysel inversion are located here, where probability is out of 200. The instrument measured three different diameters 30 nm, 60nm and 90nm each diameter was selected for an hour. During that hour three different temperatures where scanned (40°C, 80°C, 160°C) in between each temperature a bypass scan was taken at room temperature. Each scan took 10-minutes and scanned growth factors from 0.2-1.2. Calibrations of the instrument with NaCl were performed on 4/26/2023 from 11:20 to 14:20, on 5/16/2023 from 11:40 to 16:00, and on 5/30/2023 from 10:20 to 14:20. Calibrations with (NH4)2SO4 were performed on 5/8/2023 from 11:45 to 14:40, on 5/23/2023 from 11:20 to 3:00, and on 6/14/2023 from 13:00 to 16:00. Time is recorded in seconds since 1/1/1904. GF_Vol_avg: Growth factor volume average calculated by the Gysel inversion. The instrument measured three different diameters 30 nm, 60nm and 90nm each diameter was selected for an hour. During that hour three different temperatures where scanned (40°C, 80°C, 160°C) in between each temperature a bypass scan was taken at room temperature. Each scan took 10-minutes and scanned growth factors from 0.2-1.2. Calibrations of the instrument with NaCl were performed on 4/26/2023 from 11:20 to 14:20, on 5/16/2023 from 11:40 to 16:00, and on 5/30/2023 from 10:20 to 14:20. Calibrations with (NH4)2SO4 were performed on 5/8/2023 from 11:45 to 14:40, on 5/23/2023 from 11:20 to 3:00, and on 6/14/2023 from 13:00 to 16:00. Time is recorded in seconds since 1/1/1904.

30nm growth factor volume average at 160C (GF_vol_↗
DOE OSTI2025

Probing the Fine Symmetry Breaking in High-Temperature Superconductor Bi 2 Sr 2 CaCu 2 O 8+δ with Angle-Resolved Nonreciprocal Transport

After decades of research, symmetry breaking in high-temperature cuprate superconductors remains a key issue to resolve and is relevant to understanding their exotic quantum phases. In the prototypical cuprate superconductor, Bi 2 Sr 2 CaCu 2 O 8+δ (Bi2212), the possible symmetry breaking has been mostly examined microscopically with scanning tunneling microscopy and photoemission spectroscopy. However, macroscopic evidence and the direct implications for electronic transport have remained elusive. Using superconductivity-enhanced nonreciprocal transport, we report macroscopic evidence of inversion symmetry breaking in Bi2212. While the inversion symmetry breaking is subtle, its effect on nonreciprocal transport is significantly enhanced by the vortex motion during the superconducting transition, leading to a robust manifestation of inversion symmetry breaking in macroscopic transport. Combining angle-resolved nonreciprocal transport and 3D tight-binding model calculations, we derive that the inversion symmetry breaking is due to subtle crystal distortions that give rise to both in-plane and out-of-plane polar axes. Our work not only establishes nonreciprocal transport as a sensitive macroscopic probe to provide electrical transport-based evidence of fine symmetry breaking in Bi2212, but also paves the way for novel device applications such as high-temperature superconducting diodes and superconducting spintronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Peptide nucleic acid Hoogsteen strand linker design for major groove recognition of DNA thymine bases

Sequence-specific targeting of double-stranded DNA and non-coding RNA via triple-helix-forming peptide nucleic acids (PNAs) has attracted considerable attention in therapeutic, diagnostic and nanotechnological fields. An E-base (3-oxo-2,3-dihydropyridazine), attached to the polyamide backbone of a PNA Hoogsteen strand by a side-chain linker molecule, is typically used in the hydrogen bond recognition of the 4-oxo group of thymine and uracil nucleic acid bases in the major groove. We report on the application of quantum chemical computational methods, in conjunction with spatial constraints derived from the experimental structure of a homopyrimidine PNA·DNA-PNA hetero-triplex, to investigate the influence of linker flexibility on binding interactions of the E-base with thymine and uracil bases in geometry-optimised model systems. Hydrogen bond formation between the N2 E-base atom and target pyrimidine base 4-oxo groups in model systems containing a β-alanine linker (J Am Chem Soc 119:11116, 1997) was found to incur significant internal strain energy and the potential disruption of intra-stand aromatic base stacking interactions in an oligomeric context. Additionally, in geometry-optimised model systems containing a 3-trans olefin linker (Bioorg Med Chem Lett 14:1551, 2004) the E-base swung out away from the target pyrimidine bases into the solvent. These findings are in qualitative agreement with calorimetric measurements in hybridisation experiments at T–A and U–A inversion sites. In contrast, calculations on a novel 2-cis olefin linker design indicate that it could permit simultaneous E-base hydrogen bonding with the thymine 4-oxo group, circumvention and solvent screening of the thymine 5-methyl group, and maintenance of triplex intra-stand base stacking interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Residual Stress Measurements in Extreme Environments for Hazardous, Layered Specimens

In nuclear fuel plates of low-enriched U-10Mo (LEU) clad with aluminum by hot isostatic pressing (HIP), post-irradiation stresses arising during reactor shutdown are a major concern for safe reactor operations. Measurement of those residual stresses has not previously been possible because the high radioactivity of the plates requires handling only by remote manipulation in a hot cell. The incremental slitting method for measuring through-thickness stress profiles was modified, and a system for automated, remote operation was built and tested. Experimental modifications consisted of replacing electric-discharge machining (EDM) with a small end mill and strain-gauge measurements with cantilever displacement measurements. The inverse method used to calculate stresses was the pulse-regularization method modified to allow discontinuities across material interfaces. The new system was validated by comparing with conventional slitting on a depleted U-10Mo (DU) fuel plate. The new system was applied to two measurements each on six as-fabricated (pre-irradiation) LEU miniature fuel plates. Variations between the measurements at two locations in the same plate were strongly correlated with measured geometrical heterogeneity in the plate—a tilt in the fuel foil. Compressive stresses in the U-10Mo were shown to increase from 20 to 250 MPa as the ratio of aluminum thickness to U-10Mo thickness increased causing increased constraint during cooling. Faster cooling rates during processing also increased stress magnitudes. The measurements trends agreed with data in the literature from similar plates made with DU, which further validates the method. Because other methods are impractical in a hot cell, the modified slitting method is now poised for the first measurements of post-irradiation stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Charge Density Waves and Electronic Properties of Superconducting Kagome Metals

We report that kagome metals AV 3 Sb 5 (A = K, Rb, and Cs) exhibit intriguing superconductivity below 0.9 ~ 2.5 K, a charge density wave (CDW) transition around 80 ~ 100 K, and $\mathbb{Z}$ 2 topological surface states. The nature of the CDW phase and its relation to superconductivity remains elusive. In this work, we investigate the electronic and structural properties of CDW by first-principles calculations. We reveal an inverse Star of David deformation as the 2 x 2 x 2 CDW ground state of the kagome lattice. The kagome lattice shows softening breathing-phonon modes, indicating the structural instability. However, electrons play an essential role in the CDW transition via Fermi surface nesting and van Hove singularity. The inverse Star of David structure agrees with recent experiments by scanning tunneling microscopy (STM). The CDW phase inherits the nontrivial $\mathbb{Z}$ 2 -type topological band structure. Further, we find that the electron-phonon coupling is too weak to account for the superconductivity T c in all three materials. It implies the existence of unconventional pairing of these kagome metals. Our results provide essential knowledge toward understanding the superconductivity and topology in kagome metals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Size resolved particle hygroscopcity for 30, 60, and 90 nm particles collected at the EPCAPE Mount Soledad site from 04/20/2023 to 06/15/2023

Contains both the growth factor probability density function (probability is out of 200) and the average growth factors calculated by the Gysel inversion. The instrument measured three different diameters 30 nm, 60nm and 90nm each diameter was selected for an hour. During that hour, five 12-minute scans were taken scanning growth factors from 0.5-3.0. RH was consistently kept at 85%. Calibrations of the instrument with NaCl were performed on 4/26/2023 from 11:20 to 14:20, on 5/16/2023 from 11:40 to 16:00, and on 5/30/2023 from 10:20 to 14:20. Calibrations with (NH4)2SO4 were performed on 5/8/2023 from 11:45 to 14:40, on 5/23/2023 from 11:20 to 3:00, and on 6/14/2023 from 13:00 to 16:00. Time is recorded in seconds since 1/1/1904.

30 nm particle hygroscopicity↗
DOE OSTI2023

Current rectification based on noncentrosymmetric quantum materials

Rectification is a process that converts electromagnetic fields into direct current (DC). Such a process underlies a wide range of technologies, including wireless communication, wireless charging, energy harvesting, and infrared detection. Existing rectifiers are mostly based on semiconductor diodes, with limited applicability to small voltages or high frequency inputs. Here, we present an alternative approach to current rectification that uses the electronic properties of quantum crystals without semiconductor junctions. We identify a new mechanism for rectification from skew scattering due to the chirality of itinerant electrons in time-reversal-invariant but inversion-breaking materials. Our calculations reveal large, tunable rectification effects in graphene multilayers and transition metal dichalcogenides. These effects can be used in high-frequency rectifiers by rational material design and quantum wavefunction engineering.

Isobe, Hiroki↗
DOE OSTI2021

High-throughput search for magnetic topological materials using spin-orbit spillage, machine learning, and experiments

Magnetic topological insulators and semi-metals have a variety of properties that make them attractive for applications including spintronics and quantum computation. Here, we use systematic high-throughput density functional theory calculations to identify magnetic topological materials from the ≈ 40000 three-dimensional materials in the JARVIS-DFT database. First, we screen materials with net magnetic moment > 0.5 μB and spin-orbit spillage > 0.25, resulting in 25 insulating and 564 metallic candidates. The spillage acts as a signature of spin-orbit induced band-inversion. Then, we carry out calculations of Wannier charge centers, Chern numbers, anomalous Hall conductivities, surface bandstructures, and Fermi-surfaces to determine interesting topological characteristics of the screened compounds. We also train machine learning models for predicting the spillage, bandgaps, and magnetic moments of new compounds, to further accelerate the screening process. We experimentally synthesize and characterize a few candidate materials to support our theoretical predictions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Topological Circular Dichroism in Chiral Multifold Semimetals

Uncovering the physical contents of the nontrivial topology of quantum states is a critical problem in condensed matter physics. Here, we study the topological circular dichroism in chiral semimetals using linear response theory and first-principles calculations. We show that, when the low-energy spectrum respects emergent SO(3) rotational symmetry, topological circular dichroism is forbidden for Weyl fermions, and thus is unique to chiral multifold fermions. This is a result of the selection rule that is imposed by the emergent symmetry under the combination of particle-hole conjugation and spatial inversion. Using first-principles calculations, we predict that topological circular dichroism occurs in CoSi for photon energy below about 0.2 eV. In conclusion, our Letter demonstrates the existence of a response property of unconventional fermions that is fundamentally different from the response of Dirac and Weyl fermions, motivating further study to uncover other unique responses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗