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 73 records · Page 4

Direct ab initio calculation of the 4 He nuclear electric dipole polarizability

The calculation of nuclear electromagnetic sum rules by directly diagonalizing the nuclear Hamiltonian in a large basis is numerically challenging and has not been performed for A>2 nuclei. With the significant progress of high performance computing, we show that calculating sum rules using numerous discretized continuum states obtained by directly diagonalizing the ab initio no-core shell model Hamiltonian is achievable numerically. Specifically, we calculate the 4 He electric dipole (E1) polarizability, that is an inverse energy weighted sum rule, employing the Daejeon16 NN interaction. We demonstrate that the calculations are numerically tractable as the dimension of the basis increases and are convergent. Our results for the 4 He electric dipole polarizability are consistent with the most recent experimental data and are compared with those of other theoretical studies employing different techniques and various interactions.

Astronomy & Astrophysics↗

Intercalation leads to inverse layer dependence of friction on chemically doped MoS 2

Here, we present results of atomic-force-microscopy-based friction measurements on Re-doped molybdenum disulfide (MoS 2 ). In stark contrast to the widespread observation of decreasing friction with increasing number of layers on two-dimensional (2D) materials, friction on Re-doped MoS 2 exhibits an anomalous, i.e. inverse, dependence on the number of layers. Raman spectroscopy measurements combined with ab initio calculations reveal signatures of Re intercalation. Calculations suggest an increase in out-of-plane stiffness that inversely correlates with the number of layers as the physical mechanism behind this remarkable observation, revealing a distinctive regime of puckering for 2D materials.

42 ENGINEERING↗

Assessing mass balance-based inverse modeling methods via a pseudo-observation test to constrain NO x emissions over South Korea

This study constrains NO x emissions by utilizing mass balance–based inverse modeling methods and examines the optimization of inverse modeling conditions through a pseudo-observation test for South Korea, which has complex topography. Here, we applied mass balance–based inverse modeling methods such as Basic Mass Balance (BMB), Finite Difference Mass Balance (FDMB), and Iterative Finite Difference Mass Balance (IFDMB). We performed a numerical simulation of air quality using the Community Multi-scale Air Quality (CMAQ) model to calculate the NO 2 column density required for the inverse modeling. The pseudo-observation test was performed according to season, modeling resolution, and regridding methodology of satellite observation data to identify various conditions while applying the inverse modeling for South Korea. Comparing the inverse modeling results from the BMB, FDMB, and IFDMB methods, IFDMB was the most effective method in constraining NO x emissions in the South Korean region since it minimized smearing effects (i.e., transport-induced errors) through iterative calculations. The accuracy of the constrained NO x emissions using mass balance–based inversions in South Korea was the highest in the summer due to the minimized smearing effects, and the 9 km resolution modeling was the most efficient for inverse modeling of the region. In addition, the results of inverse modeling varied depending on the regridding methods, implying the importance of using a suitable regridding method and modeling resolution. This study used pseudo-observations, but the inversions are expected to be applied based on actual satellite data in the future.

54 ENVIRONMENTAL SCIENCES↗

CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.

Murdoch, Larry↗

Direct CP violation and the ΔI=1/2 rule in K→ππ decay from the standard model

We present a lattice QCD calculation of the Δ I = 1 / 2 , K → π π decay amplitude A 0 and ϵ ' , the measure of direct C P violation in K → π π decay, improving our 2015 calculation [1] of these quantities. Both calculations were performed with physical kinematics on a 32 3 × 64 lattice with an inverse lattice spacing of a - 1 = 1.3784 ( 68 ) GeV . However, the current calculation includes nearly 4 times the statistics and numerous technical improvements allowing us to more reliably isolate the π π ground state and more accurately relate the lattice operators to those defined in the standard model. We find Re ( A 0 ) = 2.99 ( 0.32 ) ( 0.59 ) × 10 - 7 GeV and Im ( A 0 ) = - 6.98 ( 0.62 ) ( 1.44 ) × 10 - 11 GeV , where the errors are statistical and systematic, respectively. The former agrees well with the experimental result Re ( A 0 ) = 3.3201 ( 18 ) × 10 - 7 GeV . These results for A 0 can be combined with our earlier lattice calculation of A 2 [2] to obtain Re ( ϵ ' / ϵ ) = 21.7 ( 2.6 ) ( 6.2 ) ( 5.0 ) × 10 - 4 , where the third error represents omitted isospin breaking effects, and Re ( A 0 ) / Re ( A 2 ) = 19.9 ( 2.3 ) ( 4.4 ) . The first agrees well with the experimental result of Re ( ϵ ' / ϵ ) = 16.6 ( 2.3 ) × 10 - 4 . A comparison of the second with the observed ratio Re ( A 0 ) / Re ( A 2 ) = 22.45 ( 6 ) , demonstrates the standard model origin of this “ Δ I = 1 / 2 rule” enhancement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Production of neutron-rich isotopes for 𝑍 ≥ 98 in the 238 U + 248 Cm reaction

Multinucleon transfer (MNT) reactions in actinide systems are a promising method to synthesize transuranium neutron-rich elements. Appropriate theoretical approaches are needed to understand the mechanisms behind MNT. We employ a microscopic approach to calculate neutron-rich isotope production in the reaction 238 U + 248 Cm system. Here, the stochastic mean-field (SMF) approach is used to calculate the primary cross sections in MNT reactions based on the quasifission and inverse quasifission processes, and a statistical de-excitation model with GEMINI ++ code to calculate the secondary fragment cross sections. The calculated cross sections using SMF and GEMINI ++ explain available experimental results for the 238 U + 248 Cm system at 𝐸 c.m. = 898.7 MeV energy. This shows the effectiveness and applicability of the quantal diffusion approach, based on the SMF theory, in heavy-ion collisions. Production of transuranium neutron-rich elements with a proton number up to 𝑍=101 is obtained with sizable cross sections. Theoretical results calculated for the 𝑍=102–105 region, for which there are no experimental data, show that the cross-section values would be lower than the microbarn level. SMF theory does not contain any adjustable parameters other than the standard parameters of the energy density functional used in the TDHF theory and is an important approach for the microscopic understanding of reaction mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Discrepancies between Theory and Experiment in Determining the Ionization Energy of NF 3

High-accuracy ab initio thermochemical predictions for the ionization energy of NF 3 , the barrier height (to inversion) of NF 3 + , and the dissociative ionization threshold of NF 3 to NF 2 + + F are presented and incorporated into Active Thermochemical Tables. The adiabatic ionization energy of the first ionization band of NF 3 , calculated at 12.647 ± 0.010 eV, is at odds with previous experimental interpretations by nearly 0.36 eV due to unfavorable Franck-Condon factors associated with this transition. The barrier (to inversion) height is calculated to be about 0.6 eV lower in energy than the prior interpretation, which instigates a discussion of the supposed vibrational structure of the first ionization band of NF 3 . Updated assignments of the photoelectron spectrum are proposed, and the loss in vibrational spacing on the high-energy side of the experimental ionization band is discussed. Rudimentary anharmonic Franck-Condon simulations qualitatively reproduce the broad spectral features observed in experiment.

Active Thermochemical Tables↗

Direct Characterization of Type-I Band Alignment in 2D Ruddlesden–Popper Perovskites

2D Ruddlesden–Popper halide perovskites have attracted considerable attention due to their desirable optoelectronic properties, high chemical and structural tunability, and improved environmental stability. However, the understanding of their structure–properties relationships is still limited. In particular, the energy level positions and band alignments at interfaces involving these materials, which are important features to control in the context of any applications, are still under debate. Here, the electronic structure of high-purity films of BA 2 MA n-1 Pb n I 3n+1 for n = 1–5 (where BA stands for butylammonium and MA for methylammonium) is investigated, using optical absorption, ultraviolet, and inverse photoemission spectroscopies, and density functional theory calculations. Further, this study determines the ionization energy and electron affinity of each compound and demonstrates a type-I band alignment for the BA 2 MA n-1 Pb n I 3n+1 series. This study further describes the evolution of the exciton binding energy as a function of the thickness of the inorganic layers.

2D Ruddlesden–Popper perovskites↗

Precision bootstrap for the $\mathcal{N}$ = 1 super-Ising model

In this note we report an improved determination of the scaling dimensions and OPE coefficients of the minimal supersymmetric extension of the 3d Ising model using the conformal bootstrap. We also show how this data can be used as input to the Lorentzian inversion formula, finding good agreement between analytic calculations and numerical extremal spectra once mixing effects are resolved.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Synthesis of [60]- and [70]Fullerene-Fused Tetrahydroquinoxaline Derivatives by Oxidative [4 + 2] Cycloaddition with Unusual Reactivity and Regioselectivity

The oxidative [4 + 2] reaction of o-phenylenediamine-derived disulfonamides with fullerene C 60 and C 70 is reported, in which electron-deficient reactants showed high reactivity. The reaction of C 70 exhibited unusual regioselectivity, yielding a [5,6]-adduct as the major product, which was characterized by 1 H, 13 C NMR and single-crystal X-ray diffraction. Here, DFT calculations revealed the reaction is an inverse-electron-demand Diels–Alder (IEDDA) reaction, and the [5,6]-adduct of C 70 is a kinetic product.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Radius of Contact During Friction Surfacing of Stainless Steel 304L: Effect of Spindle Speed and Rod Diameter

Abstract The “radius of contact” or the “real-rotational contact plane,” has been increasingly mentioned terminology in friction surfacing. However, the fundamental understanding of the flow dynamics behind this phenomenon is still very limited. The goal of this study was to understand the influence of spindle speed and consumable rod diameter on the flow dynamics and the radius of contact during friction surfacing of 304L stainless steel over a substrate of the same material. Friction surfacing was performed using consumable rods with diameters of 4.76 mm, 9.52 mm, and 12.7 mm while using spindle speeds from 1500 rpm to 20,000 rpm. The impact of spindle speed on deposition morphology, including the radius of contact, was studied. The radius of contact was calculated empirically and was found to be inversely proportional to the tangential velocity of the rod. The coupling between flow stresses and localized forces is hypothesized to be the key factor behind the variation of the radius of contact with processing conditions.

Engineering↗

Cationic Ordering and Its Influence on the Magnetic Properties of Co-Rich Cobalt Ferrite Thin Films Prepared by Reactive Solid Phase Epitaxy on Nb-Doped SrTiO 3 (001)

Here, we present the (element-specific) magnetic properties and cation ordering for ultrathin Co-rich cobalt ferrite films. Two Co-rich Co x Fe 3-x O 4 films with different stoichiometry (x=1.1 and x=1.4) have been formed by reactive solid phase epitaxy due to post-deposition annealing from epitaxial CoO/Fe 3 O 4 bilayers deposited before on Nb-doped SrTiO 3 (001). The electronic structure, stoichiometry and homogeneity of the cation distribution of the resulting cobalt ferrite films were verified by angle-resolved hard X-ray photoelectron spectroscopy. From X-ray magnetic circular dichroism measurements, the occupancies of the different sublattices were determined using charge-transfer multiplet calculations. For both ferrite films, a partially inverse spinel structure is found with increased amount of Co 3+ cations in the low-spin state on octahedral sites for the Co 1.4 Fe 1.6 O 4 film. These findings concur with the results obtained by superconducting quantum interference device measurements. Further, the latter measurements revealed the presence of an additional soft magnetic phase probably due to cobalt ferrite islands emerging from the surface, as suggested by atomic force microscope measurements.

36 MATERIALS SCIENCE↗

Spin-orbit-proximitized ferromagnetic metal by monolayer transition metal dichalcogenide: Atlas of spectral functions, spin textures, and spin-orbit torques in Co/MoSe 2 , Co/WSe 2 and Co/TaSe 2 heterostructures

In this work, the heterostructures composed of an ultrathin ferromagnetic metal (FM) and a material hosting strong spin-orbit (SO) coupling are the principal resource for SO torque and spin-to-charge conversion nonequilibrium effects in spintronics. The key quantity in theoretical description of these effects is nonequilibrium spin density, which can appear on any monolayer of the heterostructure through which the current is flowing as long as the monolayer bands are affected by the native or proximity induced SO coupling. Here we demonstrate how hybridization of wave functions of Co layer and a monolayer of transition metal dichalcogenides (TMDs)—such as semiconducting MoSe 2 and WSe 2 or metallic TaSe 2 —can lead to dramatic transmutation of electronic and spin structure of Co within some distance away from its interface with TMD, when compared to the bulk of Co or its surface in contact with a vacuum. This is due to proximity induced SO splitting of Co bands encoded in the spectral functions and spin textures on its monolayers, which we obtain using noncollinear density functional theory (ncDFT) combined with equilibrium Green function (GF) calculations. In fact, SO splitting is present due to structural inversion asymmetry of the bilayer—i.e., just the presence of the interface—even if SO coupling within TMD monolayer is artificially switched off in ncDFT calculations. However, switching it on makes the effects associated with proximity SO coupling within Co layer about five times larger. Injecting spin-unpolarized charge current through SO-proximitized monolayers of Co generates nonequilibrium spin density over them, so that its cross product with the magnetization of Co determines SO torque. The SO torque computed via first-principles quantum transport methodology, which combines ncDFT with nonequilibrium GF calculations, can be used as the screening parameter to identify optimal combination of materials and their interfaces for applications in spintronics. In particular, we identify heterostructure two-monolayer-Co/monolayer-WSe 2 as the most optimal one, at least in the clean limit.

2-dimensional systems↗

Influence of Polymorphs and Local Defect Structures on NMR Parameters of Graphite Fluorides

In this study, the role of local molecular structure on calculated 13 C and 19 F NMR chemical shifts for graphite fluoride materials was explored by using gauge-including projector augmented wave (GIPAW) computational methods for different periodic crystal polymorphs and density functional theory (DFT) gauge-including atomic orbital (GIAO) computational methods for individual graphite fluoride platelets, i.e., fluorinated graphene (FG). The impact of stacking sequences, d -spacing, and ring conformations on fully fluorinated graphite fluoride structures was investigated. A range of different defects including Stone–Wales, F and C vacancies, void formation, and F inversion were also evaluated using FG structures. These calculations show that distinct chemical shift signatures exist for many of these polymorphs and defects, therefore providing a basis for spectral assignment and development of models describing the mean local CF structure in disordered graphite fluoride materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Large-momentum effective theory’s asymptotic extrapolation vs the inverse problem

Large-momentum effective theory is a physics-guided systematic expansion to calculate light-cone parton distributions, including collinear (PDFs) and transverse-momentum-dependent ones, at any fixed momentum fraction 𝑥 within a range of [𝑥 min , 𝑥 max ]. It theoretically solves the ill-posed inverse problem that afflicts other theoretical approaches to collinear PDFs, such as short-distance factorizations. Recently, Dutrieux et al. raised practical concerns about whether current or even future lattice data will have sufficient precision in the subasymptotic correlation region to support an error-controlled extrapolation—and if not, whether it becomes an inverse problem where the relevant uncertainties cannot be properly quantified. While we agree that not all current lattice data have the desired precision to qualify for an asymptotic extrapolation, some calculations do, and more are expected in the future. We comment on the analysis and results in Dutrieux et al. and argue that a physics-based systematic extrapolation still provides the most reliable error estimates, even when the data quality is not ideal. In contrast, reframing the long-distance asymptotic extrapolation as a data-driven-only inverse problem with ad hoc mathematical conditioning could lead to unnecessarily conservative errors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scaling Ensembles of Data-Intensive Quantum Chemical Calculations for Millions of Molecules

Deep learning models are efficient computational tools that can accelerate the inverse design of molecules with desired functional properties by generating predictions at a fraction of the time required by traditional quantum chemical approaches. To ensure that a model maintains accuracy and transferability across broad regions of the chemical space explored during the inverse design, it must be trained on massively large volumes of simulation data. This requires running large-scale ensemble quantum chemical calculations on high-performance computing (HPC) systems for data collection. However, the efficient execution of such large ensemble calculations and the management of large volumes of output data require tools that can judiciously utilize computational resources and manage metadata overhead on the file system. Therefore, we present a high-performance, scalable, ensemble management framework for performing data-intensive quantum chemical electronic structure calculations for organic molecules. This framework provides abstractions to plug different ab initio, first principles, and first principles-based semi-empirical methods and executes them efficiently at large scale on HPC systems. It dynamically distributes tasks to resources and uses tiered storage for managing large collections of files. We employed this framework to process over ten million organic molecules and generate open-source datasets that provide UV-vis absorption spectra by running time-dependent density-functional tight-binding calculations. It is the largest database containing molecular optical spectra that were simulated with quantum chemical methods in a consistent manner.

Mehta, Kshitij↗

Tuning valley degeneracy with band inversion

Valley degeneracy is a key feature of the electronic structure that benefits the thermoelectric performance of a material. Despite recent studies which claim that high valley degeneracy can be achieved with inverted bands, our analysis of rock-salt IV–VI compounds using first-principles calculations and k•p perturbation theory demonstrates that mere band inversion is an insufficient condition for high valley degeneracy; rather, there is a critical degree to which the bands must be inverted to induce multiple carrier pockets. The so-called “band inversion parameter” is formalized as a chemically-tunable property, offering a design route to achieving high valley degeneracy in compounds with inverted bands. In conclusion, we predict that the valley degeneracy of rock-salt IV–VI compounds can be increased from N V = 4 to N V = 24, which could result in a corresponding increase in the thermoelectric figure of merit zT .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Divertor heat load estimates on NSTX and DIII-D using new and open-source 2D inversion analysis code

A thermography inversion algorithm has been developed in the open-source Python-based computer code, HYPERION, to calculate the heat flux incident on plasma-facing components (PFCs) in axisymmetric tokamaks. The chosen mesh size at the surface significantly affects the calculated transient heat flux results. The calculated transient heat flux will exceed the real value when the mesh size tends to zero but will underestimate the real value when the mesh size is large. A criterion for determining the appropriate mesh size for the transient heat flux calculation will be discussed. The numerical scheme for HYPERION uses a 2D fully implicit finite-difference approach, allowing temperature-dependent thermal properties of PFC materials. The inversion algorithm is benchmarked against established heat flux calculation codes, TACO and THEODOR, based on thermography data from NSTX and DIII-D respectively. The primary benefits of HYPERION compared to TACO and THEODOR are that it is open-source and it allows for the optimization of mesh thickness along the substrate. The algorithm also accounts for the thermal properties of thin surface layers that characteristically form on PFCs due to plasma-material interactions. The agreement between HYPERION and THEODOR is excellent, as the percent difference between the codes is ~5% on average in the case of the DIII-D data for moderate to high heat flux. Verification tests with TACO show slightly higher average percent differences of 8% and 12%. In using HYPERION to study filaments in heat flux, the initial results indicate that small ELMs filaments significantly broaden the divertor heat flux, and decrease divertor peak flux. Compared to the inter-ELM, the small ELM filaments decrease the divertor peak surface temperature. With intermittent divertor filaments, the divertor heat flux width is comparable with that found in L-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗