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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Closo ‐Borate Gel Polymer Electrolyte with Remarkable Electrochemical Stability and a Wide Operating Temperature Window

Abstract A major challenge in the pursuit of higher‐energy‐density lithium batteries for carbon‐neutral‐mobility is electrolyte compatibility with a lithium metal electrode. This study demonstrates the robust and stable nature of a closo ‐borate based gel polymer electrolyte (GPE), which enables outstanding electrochemical stability and capacity retention upon extensive cycling. The GPE developed herein has an ionic conductivity of 7.3 × 10 −4 S cm −2 at room temperature and stability over a wide temperature range from −35 to 80 °C with a high lithium transference number ( = 0.51). Multinuclear nuclear magnetic resonance and Fourier transform infrared are used to understand the solvation environment and interaction between the GPE components. Density functional theory calculations are leveraged to gain additional insight into the coordination environment and support spectroscopic interpretations. The GPE is also established to be a suitable electrolyte for extended cycling with four different active electrode materials when paired with a lithium metal electrode. The GPE can also be incorporated into a flexible battery that is capable of being cut and still functional. The incorporation of a closo ‐borate into a gel polymer matrix represents a new direction for enhancing the electrochemical and physical properties of this class of materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Measurement of hadronic event shapes in high-p T multijet final states at $ \sqrt{s} $ = 13 TeV with the ATLAS detector

A measurement of event-shape variables in proton-proton collisions at large momentum transfer is presented using data collected at $ \sqrt{s} $ = 13 TeV with the ATLAS detector at the Large Hadron Collider. Six event-shape variables calculated using hadronic jets are studied in inclusive multijet events using data corresponding to an integrated luminosity of 139 fb -1 . Measurements are performed in bins of jet multiplicity and in different ranges of the scalar sum of the transverse momenta of the two leading jets, reaching scales beyond 2 TeV. These measurements are compared with predictions from Monte Carlo event generators containing leading-order or next-to-leading order matrix elements matched to parton showers simulated to leading-logarithm accuracy. At low jet multiplicities, shape discrepancies between the measurements and the Monte Carlo predictions are observed. At high jet multiplicities, the shapes are better described but discrepancies in the normalisation are observed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
DOE OSTI2023

Rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions for cometary and planetary applications, as predicted by mixed quantum-classical theory

We present new calculations of collision cross sections for state-to-state transitions between the rotational states in an H 2 O + H 2 O system, which are used to generate a new database of collisional rate coefficients for cometary and planetary applications. Calculations were carried out using a mixed quantum-classical theory approach that is implemented in the code MQCT. The large basis set of rotational states used in these calculations permits us to predict thermally averaged cross sections for 441 transitions in para- and ortho-H 2 O in a broad range of temperatures. It is found that all state-to-state transitions in the H 2 O + H 2 O system split into two well-defined groups, one with higher cross-section values and lower energy transfer, which corresponds to the dipole-dipole driven processes. The other group has smaller cross sections and higher energy transfer, driven by higher-order interaction terms. We present a detailed analysis of the theoretical error bars, and we symmetrized the state-to-state transition matrixes to ensure that excitation and quenching processes for each transition satisfy the principle of microscopic reversibility. We also compare our results with other data available from the literature for H 2 O + H 2 O collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗
DOE OSTI2020

Excited-state electronic structure of molecules using many-body Green’s functions: Quasiparticles and electron–hole excitations with VOTCA-XTP

Here, we present the open-source VOTCA-XTP software for the calculation of the excited-state electronic structure of molecules using many-body Green’s function theory in the GW approximation with the Bethe–Salpeter equation (BSE). This work provides a summary of the underlying theory and discusses the details of its implementation based on Gaussian orbitals, including resolution-of-identity techniques and different approaches to the frequency integration of the self-energy or acceleration by offloading compute-intensive matrix operations using graphics processing units in a hybrid OpenMP/Cuda scheme. A distinctive feature of VOTCA–XTP is the capability to couple the calculation of electronic excitations to a classical polarizable environment on an atomistic level in a coupled quantum- and molecular-mechanics (QM/MM) scheme, where a complex morphology can be imported from Molecular Dynamics simulations. The capabilities and limitations of the GW–BSE implementation are illustrated with two examples. First, we study the dependence of optically active electron–hole excitations in a series of diketopyrrolopyrrole-based oligomers on molecular-architecture modifications and the number of repeat units. Second, we use the GW–BSE/MM setup to investigate the effect of polarization on localized and intermolecular charge-transfer excited states in morphologies of low-donor content rubrene–fullerene mixtures. These showcases demonstrate that our implementation currently allows us to treat systems with up to 2500 basis functions on regular shared-memory workstations, providing accurate descriptions of quasiparticle and coupled electron–hole excited states of various characters on an equal footing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

𝐵 → 𝜌⁢ℓ⁢$\bar{v}$ Resonance Form Factors from 𝐵→ 𝜋⁢𝜋⁢ℓ⁢$\bar{v}$ in Lattice QCD

The decay 𝐵 → 𝜌⁢ℓ⁢$\bar{v}$ is an attractive process for determining the magnitude of the smallest Cabibbo-Kobayashi-Maskawa matrix element, |𝑉 𝑢⁢𝑏 |, and can provide new insights into the origin of the long-standing exclusive-inclusive discrepancy in determinations of this standard-model parameter. This requires a nonperturbative QCD calculation of the 𝐵 → 𝜌 form factors 𝑉, 𝐴 0 , 𝐴 1 , and 𝐴 12 . The unstable nature of the 𝜌 resonance has prevented precise lattice QCD calculations of these form factors to date. Here, we present the first lattice QCD calculation of the 𝐵 → 𝜌 form factors in which the 𝜌 is treated properly as a resonance in 𝑃-wave 𝜋⁢𝜋 scattering. To this end, we use the Lellouch-Lüscher finite-volume formalism to compute the 𝐵 → 𝜋⁢𝜋 form factors as a function of both momentum transfer and 𝜋⁢𝜋 invariant mass, and then analytically continue to the 𝜌 resonance pole. This calculation is performed with 2 + 1 dynamical quark flavors at a pion mass of approximately 320 MeV, and demonstrates a clear path toward results at the physical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

One-Body Properties and Their Perturbative Accuracy with Aufbau Suppressed Coupled Cluster Theory

In this work we derived and implemented the calculation of the one-body reduced density matrix for Aufbau suppressed coupled cluster theory, from which excited state natural orbitals and one-body properties, like atomic populations and dipole moments, are obtained. We utilized the natural orbitals to refine the ASCC solution for simple valence and Rydberg systems, exploring the process of repeatedly solving the ASCC equations in successive natural orbital bases to achieve independence from the starting molecular orbitals. For dipole moments in small molecules where high-level comparison data is available, we find that the accuracy of ASCC essentially matches that of linear response and equation-of-motion coupled cluster as long as care is taken to preserve the response's perturbative completeness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

First Simultaneous Determination of Inclusive and Exclusive | V u b |

The first simultaneous determination of the absolute value of the Cabibbo-Kobayashi-Maskawa matrix element V ub using inclusive and exclusive decays is performed with the full Belle data set at the $Υ$(4S) resonance, corresponding to an integrated luminosity of 711 fb -1 . We analyze collision events in which one B meson is fully reconstructed in hadronic modes. This allows for the reconstruction of the hadronic X u system of the semileptonic b → $πℓe_ℓ$ decay. We separate exclusive B → $πℓ\overline{ν}_ℓ$ decays from other inclusive B → $X_uℓ\overline{ν}_ℓ$ and backgrounds with a two-dimensional fit that utilizes the number of charged pions in the X u system and the four-momentum transfer q 2 between the B and X u systems. Combining our measurement with information from lattice QCD and QCD calculations of the inclusive partial rate as well as external experimental information on the shape of the B → $πℓe_ℓ$ form factor, we determine |$V^{excl}_{ub}$| = (3.78 ± 0.23 ± 0.16 ± 0.14) × 10 -3 and |$V^{incl}_{ub}$| = (3.88 ± 0.20 ± 0.31 ± 0.09) x 10 -3 , respectively, with the uncertainties being the statistical error, systematic errors, and theory errors. The ratio of |$V^{excl}_{ub}$|/|$V^{incl}_{ub}$| 0.97 ± 0.12 is compatible with unity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

First-Principles Quantum Transport Modeling of Spin-Transfer and Spin-Orbit Torques in Magnetic Multilayers

A unified approach for computing (i) spin-transfer torque in magnetic trilayers like spin valve and magnetic tunnel junction, where injected charge current flows perpendicularly to interfaces, and (ii) spin-orbit torque in magnetic bilayers of the type ferromagnet/spin-orbit-coupled material, where injected charge current flows parallel to the interface, was reviewed. The experimentally explored and technologically relevant spin-orbit-coupled materials include 5d heavy metals, topological insulators, Weyl semimetals, and transition metal dichalcogenides. This approach requires to construct the torque operator for a given Hamiltonian of the device and the steady-state nonequilibrium density matrix, where the latter is expressed in terms of the nonequilibrium Green’s functions and split into three contributions. Tracing these contributions with the torque operator automatically yields field-like and damping-like components of spin-transfer torque or spin-orbit torque vector, which is particularly advantageous for spin-orbit torque where the direction of these components depends on the unknown-in-advance orientation of the current-driven nonequilibrium spin density in the presence of spin-orbit coupling. Illustrative examples are provided by computing spin-transfer torque in a one-dimensional toy model of a magnetic tunnel junction and realistic Co/Cu/Co spin valve, both of which are described by first-principles Hamiltonians obtained from noncollinear density functional theory calculations, as well as by computing spin-orbit torque in a ferromagnetic layer described by a tight-binding Hamiltonian which includes spin-orbit proximity effect within ferromagnetic monolayers assumed to be generated by the adjacent monolayer transition metal dichalcogenide. In addition, it is shown here how spin-orbit proximity effect, quantified by computing (via first-principles retarded Green’s function) spectral functions and spin textures on monolayers of realistic ferromagnetic material like Co in contact with heavy metal or monolayer transition metal dichalcogenide, can be tailored to enhance the magnitude of spin-orbit torque. Furthermore, errors made in the calculation of spin-transfer torque are quantified when using Hamiltonian from collinear density functional theory, with rigidly rotated magnetic moments to create noncollinear magnetization configurations, instead of proper (but computationally more expensive) self-consistent Hamiltonian obtained from noncollinear density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Improved Evaluation of Large Network Matrices for Linear Power Flow Within Optimization Problems: Preprint

This work discusses methods for evaluating the Power Transfer Distribution Factor (PTDF) and Line Outage Distribution Factor (LODF) matrices by employing sparse linear algebra for large-scale computing applications. These matrices are critical in many power systems applications, such as the Unit Commitment Problem (UC), pre- and post-contingency power flow analysis, and transmission expansion. These matrices are typically dense, which means they require a significant amount of time and memory to be computed for large networks. However, by analyzing the structure of the matrices and their computation method, it is possible to use reduced memory methods based on sparse matrix operations. This paper shows that sparse linear algebra algorithms are faster and require less memory and time than traditional dense approaches. Additionally, we explore the effect of matrix sparsification by eliminating trailing digits on power flow calculations.

ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Observation of D s + → η ′ μ + ν μ , Precision Test of Lepton Flavor Universality with D s + → η ( ′ ) l + ν l , and First Measurements of D s + → η ( ′ ) μ + ν μ Decay Dynamics

By analyzing 7.33 fb − 1 of e + e − annihilation data collected at center-of-mass energies between 4.128 and 4.226 GeV with the BESIII detector, we report the observation of the semileptonic decay D s + → η ′ μ + ν μ , with a statistical significance larger than 10 σ , and the measurements of the D s + → η μ + ν μ and D s + → η ′ μ + ν μ decay dynamics for the first time. The branching fractions of D s + → η μ + ν μ and D s + → η ′ μ + ν μ are determined to be ( 2.235 ± 0.05 1 stat ± 0.05 2 syst ) % and ( 0.801 ± 0.05 5 stat ± 0.02 8 syst ) % , respectively, with precision improved by factors of 6.0 and 6.6 compared to the previous best measurements. Combined with the results for the decays D s + → η e + ν e and D s + → η ′ e + ν e , the ratios of the decay widths are examined both inclusively and in several ℓ + ν ℓ four-momentum transfer ranges. No evidence for lepton flavor universality violation is found within the current statistics. The products of the hadronic form factors f + , 0 η ( ′ ) ( 0 ) and the c → s Cabibbo-Kobayashi-Maskawa matrix element | V c s | are determined. The results based on the two-parameter series expansion are f + , 0 η ( 0 ) | V c s | = 0.452 ± 0.01 0 stat ± 0.00 7 syst and f + , 0 η ′ ( 0 ) | V c s | = 0.504 ± 0.03 7 stat ± 0.01 2 syst , which help to constrain present models on f + , 0 η ( ′ ) ( 0 ) . The forward-backward asymmetries are determined to be ⟨ A FB η ⟩ = − 0.059 ± 0.03 1 stat ± 0.00 5 syst and ⟨ A FB η ′ ⟩ = − 0.064 ± 0.07 9 stat ± 0.00 6 syst for the first time, which are consistent with the theoretical calculation. Published by the American Physical Society 2024

Physics↗
DOE OSTI2024

Improved Evaluation of Large Network Matrices for Linear Power Flow Within Optimization Problems

This work presents methods for evaluating the Power Transfer Distribution Factor (PTDF) and Line Outage Distribution Factor (LODF) matrices by employing sparse linear algebra for large-scale computing applications. These matrices play a critical role in many power system applications, such as the Unit Commitment Problem (UC), pre- and post-contingency power flow analysis, and transmission expansion. These matrices are typically dense, which means they require a significant amount of time and memory to be computed for large networks. However, by analyzing the structure of the matrices and their computation method, it is possible to use reduced memory methods based on sparse matrix operations. This paper shows that sparse linear algebra algorithms are faster and require less memory and time than traditional dense approaches. Additionally, we explore the effect of matrix sparsification by eliminating trailing digits on power flow calculations.

large scale↗
DOE OSTI2021

Shedding Light on Primary Donors in Photosynthetic Reaction Centers

Chlorophylls (Chl)s exist in a variety of flavors and are ubiquitous in both the energy and electron transfer processes of photosynthesis. The functions they perform often occur on the ultrafast (fs–ns) time scale and until recently, these have been difficult to measure in real time. Further, the complexity of the binding pockets and the resulting protein-matrix effects that alter the respective electronic properties have rendered theoretical modeling of these states difficult. Recent advances in experimental methodology, computational modeling, and emergence of new reaction center (RC) structures have renewed interest in these processes and allowed researchers to elucidate previously ambiguous functions of Chls and related pheophytins. This is complemented by a wealth of experimental data obtained from decades of prior research. Studying the electronic properties of Chl molecules has advanced our understanding of both the nature of the primary charge separation and subsequent electron transfer processes of RCs. In this review, we examine the structures of primary electron donors in Type I and Type II RCs in relation to the vast body of spectroscopic research that has been performed on them to date. Further, we present density functional theory calculations on each oxidized primary donor to study both their electronic properties and our ability to model experimental spectroscopic data. This allows us to directly compare the electronic properties of hetero- and homodimeric RCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Systematic improvement of x -dependent unpolarized nucleon generalized parton distributions in lattice-QCD calculation

We present a first study of the effects of renormalization-group resummation (RGR) and leading-renormalon resummation (LRR) on the systematic errors of the unpolarized isovector nucleon generalized parton distribution in the framework of large-momentum effective theory. This work is done using lattice gauge ensembles generated by the MILC Collaboration, consisting of 2 + 1 + 1 flavors of highly improved staggered quarks with a physical pion mass at lattice spacing a ≈ 0.09 fm and a box width L ≈ 5.76 fm . We present results for the nucleon H and E generalized parton distributions (GPDs) with average boost momentum P z ≈ 2 GeV at momentum transfers Q 2 = [ 0 , 0.97 ] GeV 2 at skewness ξ = 0 as well as Q 2 ∈ 0.23 GeV 2 at ξ = 0.1 , renormalized in the modified minimal subtraction ( MS ¯ ) scheme at scale μ = 2.0 GeV , with two- and one-loop matching, respectively. We demonstrate that the simultaneous application of RGR and LRR significantly reduces the systematic errors in renormalized matrix elements and distributions for both the zero and nonzero skewness GPDs, and that it is necessary to include both RGR and LRR at higher orders in the matching and renormalization processes. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2020

Quantum electron transport in degenerate donor–acceptor systems

In this work, we develop a mathematically rigorous theory for the quantum transfer processes in degenerate donor–acceptor dimers in contact with a thermal environment. We explicitly calculate the transfer rates and the acceptor population efficiency. The latter depends critically on the initial donor state. Additionally, we show that quantum coherence in the initial state enhances the transfer process. If the electron is initially shared coherently by the donor levels, then the efficiency can reach values close to 100%, while an incoherent initial donor state will significantly suppress the efficiency. The results are useful for a better understanding of the quantum electron transport in many chemical, solid state, and biological systems with complex degenerate and quasi-degenerate energy landscapes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Efficient Parallelization of Irregular Applications on GPU Architectures

With the enlarging computation capacity of general Graphics Processing Units (GPUs), leveraging GPUs to accelerate parallel applications has become a critical topic in academia and industry. However, a wide range of irregular applications with the computation-/memory-intensive nature cannot easily achieve high GPU utilization. The challenges mainly involve the following aspects: first, data dependence leads to coarse-grained kernel and inefficient parallelism; second, heavy GPU memory usage may cause frequent memory evictions and extra overhead of I/O; third, specific computation patterns produce memory redundancies; last, workload balance and data reusability conjunctly benefit the overall performance, but there may exist a dynamic trade-off between them. Targeting these challenges, this dissertation proposes multiple optimizations to accelerate two real-world applications: many-body correlation functions to simulate nuclear physics in a large-scale scientific system; the other is the eALS-based matrix factorization recommendation system. To accelerate the calculations of many-body correlation functions, this dissertation presents three frameworks in GPU memory management and multi-GPU scheduling. Firstly, an optimized systematic GPU memory management framework, MemHC, utilizes a series of new memory reduction designs in GPU memory allocation, CPU/GPU communications, and GPU memory oversubscription. Secondly, an enhanced multi-GPU scheduling framework, MICCO, particularly by taking both data dimension (e.g., data reuse and data eviction) and computation dimension into account. MICCO designs a heuristic scheduling algorithm and a machine learning-based regression model to generate the optimal settings of a proposed new concept to manage the trade-off. Thirdly, a locality-aware multi-GPU scheduling framework. This scheduler leverages pipeline batch generation with a looking-ahead strategy by building local dependency graphs for memory transfer reduction and better data reuse, achieving up to 79.92% memory cost reduction and 1.67x speedup. To parallelize the eALS-based recommendation system, this dissertation proposes an efficient CPU/GPU heterogeneous recommendation system, HEALS. HEALS employs newly designed architecture-adaptive data formats to achieve load balance and good data locality on CPU and GPU. To mitigate the data dependence, HEALS presents a CPU/GPU collaboration model for both task parallelism and data parallelism with multiple kernel computation optimizations. In summary, this dissertation efficiently accelerates two typical irregular applications on GPUs by building four frameworks, including CPU/GPU collaboration, GPU memory management, and multi-GPU scheduling.

Wang, Qihan↗