SEARCH · Engineering Papers

Results for “MATRICES”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Simplified Models of the Bicarbonate Buffer for Scaled Simulations of CO 2 Electrolyzers

Bicarbonate electrolytes are used in a range of chemical processes; however, resolved simulation of these electrolytes is difficult, as disparate reaction time scales lead to numerical stiffness and the formation of fine boundary layers. Based on several physically motivated approximations, we reduce the full set of chemical reactions within a bicarbonate electrolyte to a simpler subset, eliminating the numerical stiffness. We supported this simplification via a two-variable singular perturbation expansion and demonstrated that under neutral conditions (6 < pH < 10) typically found within a fluid-fed CO 2 electrolyzer, the simplifications lead to negligible error. Here, we also discuss two alternative simplifications, one valid at high pH and another valid at arbitrary pH. These simplifications reduce the condition number of the matrices resulting from spatiotemporal discretization by up to 10 orders of magnitude and enable three-dimensional (3D) simulation of CO 2 electrolyzers containing carbonate solutions.

10 SYNTHETIC FUELS↗

DOE OSTI2024

A Rapid Microfluidic Neptunium Extraction Using a Supported Liquid Membrane Module

Extraction of neptunium from acidic matrices is important for its quantification, but its complex redox chemistry can cause variable yields. This study develops a microfluidic redox extraction for rapidly separating neptunium from submilliliter samples, achieving up to 90% process yield in less than 10 min for samples as small as 100 μL, with over 97% steady-state yield achieved after 20 min. It uses a supported liquid membrane module loaded with 30 vol % tributyl phosphate in n-dodecane, which performs forward- and back-extractions in a single, continuous step. Neptunium is first oxidized to +6 for extraction and then reduced during stripping. Bromate was selected as an oxidant over permanganate for its greater compatibility with the organic phase, achieving complete oxidation in under 30 s. Ascorbic acid and hydrogen peroxide were both effective reductants. Finally, the system’s high yield and rapid kinetics make it promising for future separations from complex mixtures.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2025

Atomic Orbital Energy Matching vs Overlap in Actinide-Ligand Dative Bonding

Prebonding conditions and dative bond formation in actinide(IV)hexachloride complexes with U, Np, and Pu are studied by density functional theory (DFT) calculations, focusing on the interplay of atomic orbital (AO) overlap and AO energy matching in the formation of molecular orbitals (MOs) and the subsequent identification of dative bonds. The extent of donation is tracked via population analysis, MO localization, and bond-order criteria. DFT Fock matrices are used to setup models in which varying numbers of valence AOs interact to form dative bonds. The results confirm, among other effects, a contribution to metal−ligand covalency from the actinide (An) 6p shell. Better matching of An(5f) and Cl(3p) energies occurs as the An(5f) level stabilizes with increasing An effective nuclear charge. A near-degeneracy occurs in the case of the α-spin orbitals of the Pu system, but it is inconsequential. Altogether, better AO energy matching for An(5f) and Cl(3p) along the series U, Np, Pu is counter-balanced by decreasing AO overlap and decreasing availability of 5f acceptor orbitals, leading to similar bond orders and extents of donation in the three systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Practical Guide to Chemometric Analysis of Optical Spectroscopic Data

The methodology and mathematical treatment of several classic multivariate methods for the analysis of spectroscopic data is demonstrated in a straightforward way that can be used as a basis for teaching an undergraduate introductory course on chemometric analysis. The multivariate techniques of classical least squares (CLS), principal component regression (PCR), and partial least squares (PLS), as well as the univariate Beer’s law method have been described and compared, building students’ understanding by starting with the univariate method and progressing step by step into the multivariate methods. Equations for the production of regression vectors from training set spectral data is described and their use demonstrated for the prediction of constituent concentrations on a separate validation set of spectra. Extreme care is taken to ensure consistency in variable formatting of data matrices. This provides a key foundation to understanding how spectral data are manipulated using these different mathematical approaches for building quantitative regression models. Each method is applied to a real-world data set, and the results are discussed to show students the types of information that can be gleaned from each method. A training set comprised of 20 infrared absorbance spectra containing 3 constituents (benzene, polystyrene, and gasoline) of known composition are used to demonstrate the matrix operations for each regression method. A separate set of 12 real-world napalm samples (containing benzene, polystyrene and gasoline) are used as a validation set to demonstrate the ability to utilize the regression models on an unknown dataset. A toolbox (PNNL Chemometric Toolbox) written in MATLAB language is supplied in the Supplemental Information file and can be used as a companion for understanding the development and deployment of the chemometric algorithms described in this paper. The datasets of the infrared spectra are also supplied, allowing users to build and inspect the chemometric models on their own. Finally, the Toolbox includes scripts to assist users in loading their own datasets into MATLAB and performing CLS, PCR, and PLS on their data.

Upper-Division Undergraduate, Analytical Chemistry↗

DOE OSTI2025

Nanopolysaccharide Builder: A User-Friendly Tool for Atomistic Models of Polysaccharide-Based Nanostructures

Here, we introduce Nanopolysaccharide Builder (NPB), a user-friendly software tool designed to construct polysaccharide nanostructures─mainly those based on cellulose, chitin, and chitosan─using experimental data or user-defined parameters. NPB enables the generation of cellulose and chitin allomorphs with customizable biochemical topologies and also facilitates the construction of large bundles that replicate nanostructures found in biological support systems, including plant cell walls and arthropod cuticles. The software outputs atomic Cartesian coordinates in Protein Data Bank (PDB) format and also provides atom connectivity files in PSF and PARM formats, ensuring seamless integration with major molecular dynamics (MD) engines such as NAMD, CHARMM, GROMACS, AMBER, OpenMM, and LAMMPS. Built on an interactive visualization framework, NPB features a graphical user interface (GUI) and supports both macOS and Linux operating systems. By enabling detailed atomic-scale studies of polysaccharide evolution in extracellular matrices and cell walls of algae, bacteria, fungi, and plants, NPB is poised to advance AI-guided research in sustainable chemical development and biomass utilization.

Wan, Zhangmin [Univ. of British Columbia, Vancouve↗

DOE OSTI2024

Robust Tensor Hypercontraction of the Particle–Particle Ladder Term in Equation-of-Motion Coupled Cluster Theory

One method of representing a high-rank tensor as a (hyper-)product of lower-rank tensors is the tensor hypercontraction (THC) method of Hohenstein et al. This strategy has been found to be useful for reducing the polynomial scaling of coupled-cluster methods by representation of a four-dimensional tensor of electron-repulsion integrals in terms of five two-dimensional matrices. Pierce et al. have already shown that the application of a robust form of THC to the particle–particle ladder (PPL) term reduces the cost of this term in couple-cluster singles and doubles (CCSD) from O(N 6 ) to O(N 5 ) with negligible errors in energy with respect to the density-fitted variant. In this work, we have implemented the least-squares variant of THC (LS-THC) which does not require a nonlinear tensor factorization, including the robust form (R-LS-THC), for the calculation of the excitation and electron attachment energies using equation-of-motion coupled cluster methods EOMEE-CCSD and EOMEA-CCSD, respectively. We have benchmarked the effect of the R-LS-THC-PPL approximation on excitation energies using the comprehensive QUEST database and the accuracy of electron attachment energies using the NAB22 database. Here, we find that errors on the order of 1 meV are achievable with a reduction in total calculation time of approximately 5x.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗

DOE OSTI2025

Linear-Scaling Local Natural Orbital-Based Full Triples Treatment in Coupled-Cluster Theory

We present an efficient, asymptotically linear-scaling implementation of the canonically O(N 8 ) coupled-cluster method with singles, doubles, and full triples excitations (CCSDT) method. We apply the domain-based local pair natural orbital (DLPNO) approach for computing CCSDT amplitudes. Our method, called DLPNO–CCSDT, uses the converged coupled-cluster amplitudes from a preceding DLPNO–CCSD(T) computation as a starting point for the solution of the CCSDT equations in the local natural orbital basis. To simplify the working equations, we t1-dress our two-electron integrals and Fock matrices, allowing our equations to take on the form of CCDT. With appropriate parameters, our method can recover more than 99.99% of the total canonical CCSDT correlation energy. In addition, we demonstrate that our method consistently yields sub-kJ mol –1 errors in relative energies when compared to canonical CCSDT, and, likewise, when computing the difference between CCSDT and CCSD(T). Finally, to highlight the low scaling of our algorithm, we present timings on linear alkanes (up to 30 carbons and 730 basis functions) and water clusters (up to 131 water molecules and 3144 basis functions).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗

DOE OSTI2025

Field-Driven Simulations to Probe the Impact of Ionic Correlations on Solution Transport Coefficients in Binary, Ternary, and Reciprocal Quaternary Aqueous Electrolytes

Ionic correlations play a critical role in governing transport properties of mixed-salt aqueous electrolyte solutions, yet their quantitative characterization remains challenging, particularly for multicomponent aqueous systems. Here, we develop a general nonequilibrium molecular dynamics framework to efficiently compute Onsager transport coefficients and ionic correlations in mixed-salt aqueous solutions. Using field-driven simulations, we obtain accurate Onsager matrices for LiCl/ KCl, KCl/KBr, and LiBr/KCl electrolyte solutions with significantly reduced computational cost relative to equilibrium Green−Kubo methods. The framework enables direct assessment of how attractive cation−anion and repulsive like-ion interactions contribute to conductivity and salt diffusivity across compositions. While static ion pairing exhibits strong composition dependence, dynamic ion correlations remain nearly invariant, leading to constant deviations from Nernst−Einstein predictions. These results highlight the disconnect between static ion association and dynamic transport correlations, and they establish a transferable approach for analyzing ion transport in complex electrolyte environments relevant to separation processes and electrochemical systems.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Speeding Up Hartree–Fock in JuliaChem with Density Fitting

In this work, the density fitting (DF) approximation is added to the restricted Hartree–Fock (RHF) implementation in the JuliaChem computational chemistry code. Utilizing a DF algorithm that uses symmetry and integral screening, a significant reduction in time to compute the Fock matrix is achieved. The symmetry and screening DF-RHF techniques were adapted to be performed on graphics processing units (GPUs), which are well suited to perform the matrix multiplications that comprise the bulk of the Fock build time in DF-RHF. The JuliaChem DF-RHF GPU algorithm employs a novel approach that automatically switches between two DF-RHF algorithms depending on the number of basis functions in the calculation. The JuliaChem GPU DF-RHF implementation demonstrates up to 2× speedup for Fock build times compared to the existing best-in-class GPU DF-RHF implementation by operating directly on screened intermediate matrices. Due to the high portability of the Julia language code, the JuliaChem CPU and GPU DF-RHF implementations could be benchmarked on a variety of CPU and GPU architectures from multiple hardware vendors.

Hayes, John J. [Ames Laboratory, and Iowa State Un↗

DOE OSTI2020

A Multireference Quantum Krylov Algorithm for Strongly Correlated Electrons

We introduce a multireference selected quantum Krylov (MRSQK) algorithm suitable for quantum simulation of many-body problems. MRSQK is a low-cost alternative to the quantum phase estimation algorithm that generates a target state as a linear combination of non-orthogonal Krylov basis states. This basis is constructed from a set of reference states via real-time evolution; thus, avoiding the numerical optimization of parameters. An efficient algorithm for the evaluation of the off-diagonal matrix elements of the overlap and Hamiltonian matrices is discussed and a selection procedure is introduced to identify a basis of orthogonal references that ameliorates the linear dependency problem. Preliminary benchmarks on linear H 6 , H 8 , and BeH 2 indicate that MRSQK can predict the energy of these systems accurately using very compact Krylov bases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Dynamics of Emim + in [Emim][TFSI]/LiTFSI Solutions as Bulk and under Confinement in a Quasi-liquid Solid Electrolyte

Quasi-liquid solid electrolytes are a promising alternative for next-generation Li batteries. Furthermore, these systems combine the safety of solid electrolytes with the desired properties of liquids and are typically formed by solutions of Li salts in ionic liquids incorporated into solid matrices. Here, we present a fundamental understanding of the transport properties in solutions of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([Emim][TFSI]), either in bulk form or incorporated in a boron nitride (BN) matrix. We performed a series of quasi-elastic neutron scattering experiments that, given the high incoherent neutron scattering cross section of hydrogen, allowed us to focus on the Emim + dynamics. First, [Emim][TFSI]/LiTFSI solutions (0.5 and 2.5 mol·kg –1 ) were investigated and we show how the increase in the concentration reduces the Emim+ mobility and increases the activation energy of their long-range motions. Then, the 0.5 mol·kg –1 solution was incorporated into the BN matrix and we report that the diffusivities of the Emim + cations that remain mobile under confinement are highly accelerated in comparison with the bulk sample and the activation energy of these motions is drastically reduced. We present the experimental evidence that this effect is related to the content of the Emim+ cations immobilized near the surfaces of the BN pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Numerical Simulation of Vibrational Sum Frequency Generation Intensity for Non-Centrosymmetric Domains Interspersed in an Amorphous Matrix: A Case Study for Cellulose in Plant Cell Wall

Vibrational sum frequency generation (SFG) spectroscopy can specifically probe molecular species non-centrosymmetrically arranged in a centrosymmetric or isotropic medium. This capability has been extensively utilized to detect and study molecular species present at the two-dimensional (2D) interface at which the centrosymmetry or isotropy of bulk phases is naturally broken. The same principle has been demonstrated to be very effective for the selective detection of non-centrosymmetric crystalline nanodomains interspersed in three-dimensional (3D) amorphous phases. However, the full spectral interpretation of SFG features has been difficult due to the complexity associated with the theoretical calculation of SFG responses of such 3D systems. This paper describes a numerical method to predict the relative SFG intensities of non-centrosymmetric nanodomains in 3D systems as functions of their size and concentration as well as their assembly patterns, i.e., the distributions of tilt, azimuth, and rotation angles with respect to the lab coordinate. We applied the developed method to predict changes in the CH and OH stretch modes characteristic to crystalline cellulose microfibrils distributed with various orders, which are relevant to plant cell wall structures. As a result, the same algorithm can also be applied to any SFG-active nanodomains interspersed in 3D amorphous matrices.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Adsorption and Catalytic Activity of Gold Nanoparticles in Mesoporous Silica: Effect of Pore Size and Dispersion Salinity

The assembled state of nanoparticles (NPs) within porous matrices plays a governing role in directing their biological, electronic, and catalytic properties. However, the effects of the spatial confinement and environmental factors, such as salinity, on the NP assemblies within the pores are poorly understood. In this study, we use adsorption isotherms, spectrophotometry, and small-angle neutron scattering to develop a better understanding of the effect of spatial confinement on the assembled state and catalytic performance of gold (Au) NPs in propylamine-functionalized SBA-15 and MCM-41 mesoporous silica materials (mSiO 2 ). We carry out a detailed investigation of the effect of pore diameter and ionic strength on the packing and spatial distribution of AuNPs within mSiO 2 to get a comprehensive insight into the structure, functioning, and activity of these NPs. We demonstrate the ability of the adsorbed AuNPs to withstand aggregation under high salinity conditions. We attribute the observed preservation of the adsorbed state of AuNPs to the strong electrostatic attraction between oppositely charged pore walls and AuNPs. The preservation of the structure allows the AuNPs to retain their catalytic activity for a model reaction in high salinity aqueous solution, here, the reduction of p-nitrophenol to p-aminophenol, which otherwise is significantly diminished due to bulk aggregation of the AuNPs. This fundamental study demonstrates the critical role of confinement and dispersion salinity on the adsorption and catalytic performance of NPs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Structural and Electronic Effect Driven Distortions in Visible Light Absorbing Polar Materials A Ta 2 V 2 O 11 ( A = Sr, Pb)

Complex vanadates of tantalum(V), such as ATa 2 V 2 O 11 (A = Sr, Pb), are rare and underrated materials, which have potential application domains that could be substantially expanded, mitigating the existing controversy on their atomic and electronic organization. Herein, we present a thorough structural examination combining synchrotron powder X-ray diffraction-aided distortion mode analysis with computational methods to study hettotypes of SrTa 2 V 2 O 11 (STVO) and PbTa 2 V 2 O 11 (PTVO). Being distinct from the perovskite family due to the presence of [VO 4 ] groups, both compounds are polar dielectric materials with certain similarities to SBT and PBT Aurivillius phases. Applying the model of anions of metallic matrices to the analysis of electron localization functions calculated on top of as-established equilibrium structures helps retrace the effects in the Sr and Pb surroundings on the respective crystal packings of STVO and PTVO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Facilitating Screening of MOFs for Mixed Matrix Membranes Using Machine Learning and the Maxwell Model

Metal organic framework (MOF)-based mixedmatrix membranes (MMMs), which embed MOF particles in polymer matrices, combine the advantages of polymeric and inorganic membranes. Multiple previous studies have used the Maxwell model together with molecular simulations and machine learning (ML) to predict the performance of MOF/polymer MMMs. However, the assumption of rigid MOF frameworks in molecular simulations limited the accuracy of the data used in the predictions, particularly in predicting molecular diffusivities. We developed a novel workflow integrating ML models with consideration of MOF flexibility to predict the permeability and selectivity of 131,722 MMMs for CO 2 /CH 4 , O 2 /N 2 and He/H 2 separations. The full range of achievable MMM performance within the Maxwell model was analyzed, and several promising MOFs were identified using this workflow. This approach offers an efficient tool for screening any polymer and MOF combination in gas separation applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Percolated Ionic Aggregate Morphologies and Decoupled Ion Transport in Precise Sulfonated Polymers Synthesized by Ring-Opening Metathesis Polymerization

Here, we describe a set of precise single-ion conducting polymers that form self-assembled percolated ionic aggregates in glassy polymer matrices and have decoupled transport of metal cations. These precise single-ion conductors (SICs), synthesized by a scalable ring-opening metathesis polymerization, consist of a polyethylene backbone with a sulfonated phenyl group pendant on every fifth carbon and are fully neutralized by a counterion X + (Li + , Na + , or Cs + ). Experimental X-ray scattering measurements and fully atomistic molecular dynamics (MD) simulations are in good agreement. The MD simulations show that the ionic groups nanophase separate from the polymer backbone to form percolating ionic aggregates. Using graph theory, we find that within the Li + - and Na + -neutralized polymers the percolated aggregates exhibit planar and ribbon-like configurations at intermediate length scales, while the percolated aggregates within the Cs + -neutralized polymers are more isotropic. Electrical impedance spectroscopy measurements show that the ionic conductivities exhibit Arrhenius behavior, with conductivities of 10 –7 to 10 –6 S/cm at 180 °C. In the MD simulations, the cations move between sulfonate groups in the percolated aggregates, larger ions travel further, and overall cations travel further than the polymer backbones, indicating a decoupled ion-transport mechanism. Thus, the percolated ionic aggregates in these polymers can serve as pathways to facilitate decoupled ion motion through a glassy polymer matrix.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗