Search NASA⌕ Search

SEARCH · Search NASA

Results for “transfer matrix calculations”

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

In-pile thermal conductivity of uranium dioxide at low burnup

In-pile thermal conductivity of uranium dioxide (UO 2 ) was investigated at low burnup levels (<0.2 kWd/kg-UO 2 ) to elucidate the concurrent effects of the fission damage and the thermal recovery. In-pile experiments were performed in the Halden reactor. The fuel centerline and the cladding temperatures were measured during the experiments, and the power level of specially designed test rodlets were also monitored. The uncertainty of the pellet-cladding gap's thermal resistivity was minimized with a liquid-metal bond, while comparable measurements were made in reference test rodlets without liquid metal bonding. The experimental data were analyzed using the inverse heat transfer approach. The heat conduction equation was solved by applying the measured temperature of the cladding as a boundary condition, and the calculated fuel centerline temperature was compared with the measured fuel centerline temperature to determine the relative matrix resistivity. Results showed that the in-pile matrix resistivity, the phonon-lattice interaction terms of the thermal conductivity, was ~1.5–2.5 times higher than its unirradiated value due to fission-induced damage. Finally, the in-pile annealing caused a significant recovery of the matrix resistivity, and the amount of recovery increased from ~20 to ~33%, while the annealing temperature increased from 700 to 1,000 °C.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Anisotropic quantum Hall states in the presence of interactions with fourfold rotational symmetry

We study the effects of anisotropic interactions in the quantum Hall effect in the presence of a fourfold discrete rotational (C 4 ) symmetry. Employing the density matrix renormalization group technique on an infinite cylinder geometry (iDMRG), we calculate the anisotropy response of the Laughlin state at ν = 1/3 and the composite-Fermi liquid (CFL) state at ν = 1/2. We find that the anisotropy transferred from the interaction potential to the ν = 1/3 state is stronger when compared with the complementary case of an anisotropic band. Further, the strength of anisotropy reduces as the interaction is made shorter ranged. Quite surprisingly, at ν = 1/2, the deformation in the CF Fermi surface changes sign as the interaction range is reduced. Our results imply that the short-distance and long-distance parts of the interaction potential have opposite effects on the CFL state in the presence of C 4 -symmetric anisotropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Photon-odderon interference in exclusive χ c charmonium production at the Electron-Ion Collider

Exclusive C = + 1 scalar, axial-vector, and tensor quarkonium production in high-energy electron-proton scattering requires a C -odd t -channel exchange of a photon or a three-gluon ladder. We derive the expressions for the corresponding amplitudes. The relative phase of the photon vs three-gluon exchange amplitudes is determined by the sign of the light-front matrix element of the eikonal color current operator d a b c J + a J + b J + c at moderate x , and is not affected by small- x QCD evolution. Model calculations predict constructive interference, which is particularly strong for momentum transfer | t | ∼ 1 GeV 2 where the cross section for χ c J production exceeds that for pure photon exchange by up to a factor of 4. Exclusive χ c J electroproduction at the high luminosity Electron-Ion Collider occurs with well measurable rates and measurements of these processes should find evidence for the perturbative odderon exchange. We also compute the total electroproduction cross section as a function of energy and provide first estimates of the number of χ c J events per month at the Electron-Ion Collider design luminosity. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Characterization of shale using Helium and Argon at high pressures

In order to estimate the shale gas in place and the eventual recovery during shale gas operations, one of the key requirements is to accurately characterize the shale’s petrophysical and transport properties such as porosity, permeability, diffusivity, and storage capacity. Despite the many efforts reported in the technical literature aiming to provide an improved understanding of the complex pore structures and the associated fluid flow in gas shales, complete characterization of organic-rich shale samples still poses a big challenge. Here, we have characterized mass transfer and sorption in shale at different length scales using Helium (He) and Argon (Ar) as probe gases. Thermogravimetric analysis (TGA) with a shale cube of ~1 cm 3 in size and gas expansion experiments with a full-diameter core (3.5” in diameter) were used to measure sorption kinetics/isotherms and mass transfer, respectively. Both samples are from the same depth/location in the Marcellus shale formation. The TGA steady-state technique was initially used to generate excess sorption isotherms for Ar, while dynamic TGA experiments were used to study its sorption kinetics. The TGA experiments demonstrate that Ar, which has a similar sorption potential as Methane, but is generally assumed to be inert, adsorbs onto the surfaces of the mesoporous and microporous regions of the shale samples according to a Langmuir-type behavior. Helium expansion experiments, on the full-diameter core, were used to measure the overall porosity, on the basis that He is a non-sorbing and inert gas as compared to Ar. The He expansion experiments, furthermore, allow us to delineate the mass transfer of gas across the inherent hierarchy of pore sizes, including macropores (macro- and microcracks), mesopores and micropores. Similar expansion experiments were also performed with Ar to study the combined impact of mass transfer and sorption. A triple-porosity model (TPM) was utilized to interpret the He expansion experiments with the shale core and to extract (estimate) relevant transport parameters. We report and compare here the diffusivities and permeabilities of the whole core for both He and Ar, as calculated from the modeling and fitting of the experimental data. On the premise that the shale cube is representative of the matrix region of the core, the Ar sorption kinetics from the cube experiments were subsequently combined with the transport parameters extracted from the He experiments to predict the behavior of the Ar expansion test with the full-diameter core. An excellent agreement is observed between the model predictions and the experimental data. The experimental observations and their interpretation indicate that one must be cautious when using Ar to estimate the true porosity and permeability of shales. In addition, we demonstrate that He and Ar probe gases, when used in tandem, can be employed effectively as a tool to characterize shales in terms of mass transfer and sorption dynamics across scales.

04 OIL SHALES AND TAR SANDS↗

Nature of polariton transport in a Fabry-Perot cavity

Fabry-Perot microcavities can strongly enhance interactions between light and molecules, leading to the formation of hybrid light-matter states known as polaritons. Polaritons possess much smaller effective masses and much larger group velocities when the molecules are resonant with cavity modes that have finite (nonzero) in-plane wave vectors, giving rise to the possibilities of long-range and ultrafast ballistic transport. In this paper, we present the results of numerical simulations of the ultrafast ballistic transport phenomenon in real space and time during and after initialization with a short, spatially localized pulse. Here, we address this problem with two approaches: the standard transfer-matrix method (TMM) for planar structures and a second simulation based on the numerical solution of the Maxwell-Bloch equations that can be used for general configurations (complex metasurfaces, for instance) and excitation modes. The agreement between the TMM and the full numerical calculation with the Maxwell-Bloch equations when applied to multilayer planar configurations provides proof of the validity of both approaches for the present analyses. Overall, we find that the transport of the molecular excitons inside the cavity synchronize with the evolution of the enhanced electromagnetic field inside the cavity. Moreover, the synchronized transport rate is in good agreement with the group velocities predicted from a calculated dispersion relation across a wide range of frequencies. Finally, we relate the group velocity to the Hopfield coefficient obtained from quantum modeling and suggest that the dependence of light-matter coupling on the in-plane wave vector can be an important but overlooked factor for understanding the transport behavior of polaritons. These simulations provide an intuitive tool for understanding the collective motion of light and excitons and helps us to better understand how experimental observations of polaritons should be interpreted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Radiative Width of K*(892) from Lattice Quantum Chromodynamics

In this dissertation, we use lattice quantum chromodynamics to explore the radiative transitions of ?K to K, to calculate the radiative-width of the resonant K?(892) which appears in the P-wave ?K ? ?K transition amplitude. The matrix elements are extracted from three-point functions calculated in a finite-volume discretized lattice with a pion mass of 284 MeV. The finite-volume amplitudes, which are constrained over a large number of ?K energy-points and four-momentum transfers, are mapped to the infinite volume transition amplitude by using the Lellouch-Lüscher formalism. The radiative width is determined to be ? = 35 ± 8 keV by analytically continuing the amplitude into the complex energy plane and calculating the residue at the K? pole.

Radhakrishnan, Archana↗

Full Counting Statistics of Charge in Chaotic Many-Body Quantum Systems

We investigate the full counting statistics of charge transport in U(1)-symmetric random unitary circuits. We consider an initial mixed state prepared with a chemical potential imbalance between the left and right halves of the system and study the fluctuations of the charge transferred across the central bond in typical circuits. Using an effective replica statistical mechanics model and a mapping onto an emergent classical stochastic process valid at large on-site Hilbert space dimension, we show that charge transfer fluctuations approach those of the symmetric exclusion process at long times, with subleading t –1/2 quantum corrections. Here, we discuss our results in the context of fluctuating hydrodynamics and macroscopic fluctuation theory of classical nonequilibrium systems and check our predictions against direct matrix-product state calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tracking of Marangoni driven motion during laser powder bed fusion

Laser Powder Bed Fusion (LPBF) can be used to create Functionally Graded Materials (FGMs) by selectively distributing property augmenting particles throughout the build matrix. With this method parts can be created with nonhomogeneous properties tuned to spatially varied, pre-determined design constraints. To unlock this capability, the fluid motion that occurs in the molten liquid phase must be thoroughly understood. More importantly, the actual tracking of motion during the LPBF process is needed. Further, a numerical model of heat transfer and fluid flow in LPBF process was created. Using results from the model the motion of particles starting at many different locations was tracked. The average change in position for these particles based on starting location was calculated. The effects of non-dimensional parameters on melt pool size were examined. These results were validated experimentally and compared to what is available in the literature.

42 ENGINEERING↗

On the novel 3-D neutron transport kinetic tRAPID algorithm and its validation

The Real-time Analysis for Particle-transport and In-situ Detection (RAPID) Code System, based on the Multi-stage Response-function Transport (MRT) methodology, allows for real-time simulation of nuclear systems based on 3-D continuous-energy particle transport. RAPID's steady-state (criticality) neutron transport algorithm is based on the Fission Matrix (FM) method, and has been extensively verified and validated against computational benchmarks and experiments. This paper introduces the novel 3-D time-dependent transport algorithm that has been implemented into the code, tRAPID, and its validation using the JSI TRIGA Mark-II reactor. tRAPID accurately and efficiently calculates neutron kinetics parameters (such as β{sub eff}, l{sub eff} , Λ, α{sub Rossi}) and 3-D time-dependent neutron fission source distribution and neutron importances for both prompt and delayed neutrons. tRAPID is used to simulate a rod insertion experiment performed at the JSI TRIGA Mark-II reactor, during which signals from four fission chambers at four different locations in the core were collected. The results demonstrate how tRAPID is capable of calculating detailed and accurate results with only a minimal use computational resources and time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An extended moments model of quantum efficiency for metals and semiconductors

The complexity of photocathode designs and detector materials, and the need to model their performance for short pulse durations, response to high frequency photons, presence of coatings and/or thinness of the absorptive layer, necessitates modifications to three-step and moments models of photoemission that are used in simulation codes. Here, methods to include input from computationally intensive approaches, such as density functional theory to model optical properties and transfer matrix approaches to treat emission from the surface or transport past coatings, by means of parametric models are demonstrated. First, a technique to accurately represent optical behavior so as to model reflectivity and penetration depth is given. Second, modifications to bulk models arising from the usage of thin film architectures, and a means to rapidly calculate them, is provided. Third, a parameterization to model the impact of wells associated with coatings and surface layers on the transmission probably is given. In all cases, the methods are computationally efficient and designed to allow for including input from numerically intensive approaches that would otherwise be unavailable for simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Subnanometer Scale Mapping of Hydrogen Doping in Vanadium Dioxide

Hydrogen donor doping of correlated electron systems such as vanadium dioxide (VO 2 ) profoundly modifies the ground state properties. The electrical behavior of H x VO 2 is strongly dependent on the hydrogen concentration; hence, atomic scale control of the doping process is necessary. It is however a nontrivial problem to quantitatively probe the hydrogen distribution in a solid matrix. As hydrogen transfers its sole electron to the material, the ionization mechanism is suppressed. Here in this study, a methodology mapping the doping distribution at subnanometer length scale is demonstrated across a H x VO 2 thin film focusing on the oxygen-hydrogen bonds using electron energy loss spectroscopy (EELS) coupled with first-principles EELS calculations. The hydrogen distribution was revealed to be nonuniform along the growth direction and between different VO 2 grains, calling for intricate hydrogenation mechanisms. Our study points to a powerful approach to quantitatively map dopant distribution in quantum materials relevant to energy and information science

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Porting fragmentation methods to GPUs using an OpenMP API: Offloading the resolution-of-the-identity second-order Møller–Plesset perturbation method

Here, using an OpenMP Application Programming Interface, the resolution-of-the-identity second-order Møller–Plesset perturbation (RI-MP2) method has been off-loaded onto graphical processing units (GPUs), both as a standalone method in the GAMESS electronic structure program and as an electron correlation energy component in the effective fragment molecular orbital (EFMO) framework. First, a new scheme has been proposed to maximize data digestion on GPUs that subsequently linearizes data transfer from central processing units (CPUs) to GPUs. Second, the GAMESS Fortran code has been interfaced with GPU numerical libraries (e.g., NVIDIA cuBLAS and cuSOLVER) for efficient matrix operations (e.g., matrix multiplication, matrix decomposition, and matrix inversion). The standalone GPU RI-MP2 code shows an increasing speedup of up to 7.5× using one NVIDIA V100 GPU with one IBM 42-core P9 CPU for calculations on fullerenes of increasing size from 40 to 260 carbon atoms using the 6-31G(d)/cc-pVDZ-RI basis sets. A single Summit node with six V100s can compute the RI-MP2 correlation energy of a cluster of 175 water molecules using the correlation consistent basis sets cc-pVDZ/cc-pVDZ-RI containing 4375 atomic orbitals and 14 700 auxiliary basis functions in ~0.85 h. In the EFMO framework, the GPU RI-MP2 component shows near linear scaling for a large number of V100s when computing the energy of an 1800-atom mesoporous silica nanoparticle in a bath of 4000 water molecules. The parallel efficiencies of the GPU RI-MP2 component with 2304 and 4608 V100s are 98.0% and 96.1%, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Skewness-dependent moments of the pion GPD from nonlocal quark-bilinear correlators

We present lattice QCD calculations of the odd Mellin moments of pion valence-quark generalized parton distribution up to fifth order ⟨𝑥 4 ⟩ and for the skewness range [−0.33, 0] using operator product expansion of bilocal quark-bilinear operators. The calculations are performed on an ensemble with lattice spacing 𝑎 = 0.04 fm and valence pion mass 300 MeV, employing boosted pion states with momenta up to 2.428 GeV and momentum transfers reaching 2.748 GeV 2 . We employ ratio-scheme renormalization and next-to-leading logarithmic resummed perturbative matching. At zero skewness, our results are consistent with previous lattice studies. By combining matrix elements at multiple values of skewness and momentum transfer, skewness-dependent moments are obtained through simultaneous polynomiality-constrained fits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evaluation of partial volume correction and analysis of longitudinal [ 18 F]GTP1 tau PET imaging in Alzheimer's disease using linear mixed-effects models

Purpose: We evaluated the impact of partial volume correction (PVC) methods on the quantification of longitudinal [ 18 F]GTP1 tau positron-emission tomography (PET) in Alzheimer's disease and the suitability of describing the tau pathology burden temporal trajectories using linear mixed-effects models (LMEM). Methods: We applied van Cittert iterative deconvolution (VC), 2-compartment, and 3-compartment, and the geometric transfer matrix plus region-based voxelwise methods to data acquired in an Alzheimer's disease natural history study over 18 months at a single imaging site. We determined the optimal PVC method by comparing the standardized uptake value ratio change (%ΔSUVR) between diagnostic and tau burden–level groups and the longitudinal repeatability derived from the LMEM. The performance of LMEM analysis for calculating %ΔSUVR was evaluated in a natural history study and in a multisite clinical trial of semorinemab in prodromal to mild Alzheimer's disease by comparing results to traditional per-visit estimates. Results: The VC, 2-compartment, and 3-compartment PVC methods had similar performance, whereas region-based voxelwise overcorrected regions with a higher tau burden. The lowest within-subject variability and acceptable group separation scores were observed without PVC. The LMEM-derived %ΔSUVR values were similar to the per-visit estimates with lower variability. Conclusion: The results indicate that the tested PVC methods do not offer a clear advantage or improvement over non-PVC images for the quantification of longitudinal [ 18 F]GTP1 PET data. LMEM offers a robust framework for the longitudinal tau PET quantification with low longitudinal test–retest variability. Clinical trial registration: NCT02640092 and NCT03289143.

60 APPLIED LIFE SCIENCES↗

A matrix completion algorithm for efficient calculation of quantum and variational effects in chemical reactions

This work examines the viability of matrix completion methods as cost-effective alternatives to full nuclear Hessians for calculating quantum and variational effects in chemical reactions. The harmonic variety-based matrix completion (HVMC) algorithm, developed in a previous study (https://doi.org/10.1063/5.0018326), exploits the low-rank character of the polynomial expansion of potential energy to recover, using a small sample, vibrational frequencies (square roots of nuclear Hessian eigenvalues) constituting the reaction path. Furthermore, these frequencies are essential for calculating rate coefficients using variational transition state theory with multidimensional tunneling (VTST-MT). HVMC performance is examined for four SN2 reactions and five hydrogen transfer reactions, with each H-transfer reaction consisting of at least one vibrational mode strongly coupled to the reaction coordinate. HVMC is robust and captures zero-point energies, vibrational free energies, zero-curvature tunneling, and adiabatic ground state and free energy barriers as well as their positions on the reaction coordinate. For medium to large reactions involving H-transfer, with the exception of the most complex Ir catalysis system, less than 35% of total eigenvalue information is necessary for accurate recovery of key VTST-MT observables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural insight into host plasma membrane association and assembly of HIV-1 matrix protein

Oligomerization of Pr55 Gag is a critical step of the late stage of the HIV life cycle. It has been known that the binding of IP6, an abundant endogenous cyclitol molecule at the MA domain, has been linked to the oligomerization of Pr55 Gag . However, the exact binding site of IP6 on MA remains unknown and the structural details of this interaction are missing. Here, we present three high-resolution crystal structures of the MA domain in complex with IP6 molecules to reveal its binding mode. Additionally, extensive Differential Scanning Fluorimetry analysis combined with cryo- and ambient-temperature X-ray crystallography and GNM-based transfer entropy calculations identify the key residues that participate in IP6 binding. Our data provide novel insights about the multilayered HIV-1 virion assembly process that involves the interplay of IP6 with PIP2, a phosphoinositide essential for the binding of Pr55 Gag to membrane. IP6 and PIP2 have neighboring alternate binding sites within the same highly basic region (residues 18–33). This indicates that IP6 and PIP2 bindings are not mutually exclusive and may play a key role in coordinating virion particles’ membrane localization. Based on our three different IP6-MA complex crystal structures, we propose a new model that involves IP6 coordination of the oligomerization of outer MA and inner CA domain’s 2D layers during assembly and budding.

60 APPLIED LIFE SCIENCES↗

Unifying Combinatorial and Graphical Methods in Artificial Intelligence

Recently, a new graph Laplacian, called the inner product Laplacian, was introduced which generalizes many existing Laplacians, including the normalized and combinatorial Laplacian and their weighted variants. The key observation behind the inner product Laplacian is that by defining appropriate inner product spaces on the vertices and edges, the standard Laplacians can be recovered as Hodge Laplacians over the simplicial complex formed by the edges and vertices. These inner product spaces form a natural way to incorporate non-combinatorial information into the definition of a domain-specific Laplacian. In particular, in contrast to current domain-specific weighting schemes which rely solely on edge weights, information regarding the similarity of non-adjacent vertices and arbitrary pairs of edges can be effectively incorporated into the Laplacian. In order to illustrate this approach we consider the problem of calculating the potential energy of an atomistic configuration using Graph Neural Networks. In comparison with start-of-the-art approaches, such as SchNet, our approach replaces a learned (via auto-encoder) representation of the atom types with an inner product space on atoms based on scientific knowledge (e.g., electronegativity). We will illustrate how this approach captures key chemical properties of the molecules and compare the energy calculations with state-of-the-art neural network approaches. However, to compute the resulting Laplacian involves a mixture of sparse and dense matrix computation and yields a dense matrix as the basis for the graph convolution. This dense convolutional kernel necessitates moving away from the standard message passing framework for graph neural networks and increases the computational cost of applying the kernel. In order to mitigate these costs we investigate means of leveraging the mixed sparse and dense computations to reduce the overall computational cost and how these approaches can be automatically transferred to energy efficient hardware (e.g., field programmable gate arrays (FPGAs)).

97 MATHEMATICS AND COMPUTING↗

Mechanistic Insights into Visible-Light-Induced ATRA Reactions Powered by the Symbiotic Relationship between Cu(II)/Cu(I) Phenanthroline Complexes

As photoredox catalysis continues to yield promising chemical transformations, there is an increased need to understand how specific photocatalysts function to improve reaction efficiencies while expanding their scope. Copper-phenanthroline-based photocatalysts such as Cu II (dap)Cl 2 (dap = 2,9-di(p-anisyl)-1,10-phenanthroline) and [Cu I (dap) 2 ]Cl were both found to be equally capable of olefin activation through electrophilic atom transfer radical addition (ATRA) reactions. Although these molecular catalysts have proven successful, many intermediates suggested in the proposed catalytic cycles have never been detected. One undetermined aspect in this chemistry is related to how one equivalent of Cu II (dap)Cl 2 generates half of an equivalent of [Cu I (dap) 2 ] + during the photocatalytic sequence. To this end, we initially used more synthetically accessible model systems, namely, [Cu I (dpp) 2 ]Cl and Cu II (dpp)Cl 2 (dpp = 2,9-diphenyl-1,10-phenanthroline), to glean detailed mechanistic insights into this unusual symbiotic relationship. We directly detected several intermediates involved in the ATRA photocatalytic cycle using these model chromophores in conjunction with electronic spectroscopy, infrared matrix-assisted laser desorption electrospray ionization (IR-MALDESI) mass spectrometry, electronic structure calculations, EPR spin-trap experiments, and 1 H NMR spectroscopy. We found that the unique ligand lability and coordinating properties of acetonitrile enable both the in situ oxidation of [Cu I (dpp) 2 ] + by tosyl chloride into Cu II (dpp)Cl 2 and the visible-light-induced homolysis of the Cu II –Cl bond, which initiates the conversion to the Cu I species [Cu I (dpp) 2 ][Cu I Cl 2 ]. Here, the combined findings from the present study of the catalytic cycle demonstrate that the symbiotic relationship between Cu II (dpp)Cl 2 and [Cu I (dpp) 2 ] + , as well as between Cu II (dap)Cl 2 and [Cu I (dap) 2 ] + , is the critical factor enabling the ATRA photoreaction by departing from either photocatalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗