Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Markovian transport processes in a heterogeneous, variably saturated watershed: A multi-domain spatial Markov model

Transport processes in watersheds remain challenging to describe accurately with upscaled models because of their intricate complexities at multiple scales that can lead to multi-modal, non-Fickian breakthrough curves. One of the recent advances in solute transport modeling has been the Lagrangian spatial Markov model (SMM), which describes transport using velocity increments and models the correlated transitions between them. The method has been applied previously to saturated porous media and fractured media, but these successful applications suggest it may be useful for describing watershed scale processes as well. An existing, 3-D, heterogeneous variably saturated model of a headwaters catchment and classical random walk particle tracking were used to generate a forward simulation and the resulting transport was analyzed with a spatial Markov framework. Existing SMM methods were used to analyze the transport paths through the complex, unconfined domain and we found that the vadose zone plays an important role in the correlation structure. As such, a new kind of SMM is proposed, termed a multi-domain SMM (MD-SMM), that uses separate correlation models for auto-transitions in the vadose zone and the saturated zone, and for the cross-transitions between them. The MD-SMM improves representation of transport relative to the single domain model and can naturally delineate times when a particle is in the vadose zone or the saturated zone, adding an additional degree of realism to the upscaled model. Lastly, we show that the transition probability matrix for transport in the saturated portion of this unconfined system is approximately Markovian, which is an important validation of the proposed framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

R –8 Dispersion Interaction: Derivation and Application to the Effective Fragment Potential Method

The anisotropic and isotropic R –8 dispersion contributions (disp8) are derived and implemented within the framework of the effective fragment potential (EFP) method formulated with imaginary frequency-dependent Cartesian polarizability tensors distributed at the centroids of the localized molecular orbitals (LMOs). Two forms of damping functions, intermolecular overlap-based and Tang–Toennies, are extended for disp8. To obtain LMO polarizability tensors centered at LMO centroids, an origin-shifting transformation is derived and implemented for the dipole–octopole polarizability tensor and the quadrupole–quadrupole polarizability tensor. The analytic gradient is derived and implemented for the isotropic disp8 contribution. Relative to the previously implemented empirical EFP disp8 energy, the isotropic disp8 component of the interaction energy improves the overall agreement of the EFP dispersion energies with the symmetry-adapted perturbation theory (SAPT) benchmarks, reducing the mean absolute errors (MAEs) and mean absolute percentage errors for most of the databases examined in this work. Further, while the anisotropic disp8 can further enhance the accuracy of the EFP dispersion energy and yield smaller MAEs, significantly overbound dispersion energies are predicted by the anisotropic disp8 when the maximum element in the intermolecular overlap matrix is greater than 0.1, possibly due to the breakdown of the approximations made in the EFP dispersion derivation at a short range. For potential energy scan databases, the newly developed EFP dispersion model with isotropic disp8 yields the overall correct curvature and good agreement with SAPT benchmarks around equilibrium and longer but overestimates the dispersion interactions at a short range. While the overlap-based dispersion-damping functions produce better MAEs than Tang–Toennies damping functions, further improvement is needed to better screen the large attractive dispersion energies at a short range (overlap >0.1).

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

Classical Benchmarks for Variational Quantum Eigensolver Simulations of the Hubbard Model

Simulating the Hubbard model is of great interest to a wide range of applications within condensed matter physics, however its solution on classical computers remains challenging in dimensions larger than one. The relative simplicity of this model, embodied by the sparseness of the Hamiltonian matrix, allows for its efficient implementation on quantum computers, and for its approximate solution using variational algorithms such as the variational quantum eigensolver. While these algorithms have been shown to reproduce the qualitative features of the Hubbard model, their quantitative accuracy in terms of producing true ground state energies and other properties, and the dependence of this accuracy on the system size and interaction strength, the choice of variational ansatz, and the degree of spatial inhomogeneity in the model, remains unknown. Here we present a rigorous classical benchmarking study, demonstrating the potential impact of these factors on the accuracy of the variational solution of the Hubbard model on quantum hardware, for systems with up to 32 qubits. We find that even when using the most accurate wavefunction ansätze for the Hubbard model, the error in its ground state energy and wavefunction plateaus for larger lattices, while stronger electronic correlations magnify this issue. Concurrently, spatially inhomogeneous parameters and the presence of off-site Coulomb interactions only have a small effect on the accuracy of the computed ground state energies. Our study highlights the capabilities and limitations of current approaches for solving the Hubbard model on quantum hardware, and we discuss potential future avenues of research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Processing Compressed Expanded Natural Graphite for Phase Change Material Composites

Phase change materials (PCMs) are used in various thermal energy storage applications but are limited by their low thermal conductivity. One method to increase conductivity involves impregnating organic PCMs into highly porous conductive matrix materials. Of these materials, compressed expanded natural graphite (CENG) matrices have received the most attention. Despite this attention, the effect that CENG processing has on PCM saturation and overall matrix thermal conductivity has not been fully investigated. Therefore, the effect of the heat treatment process used to expand intercalated graphite flakes is evaluated here. Higher heat treatment temperatures yielded higher saturation rates and overall saturation at similar matrix porosities. For example, increasing temperature from 300 °C to 700 °C resulted in approximately 60%-70% increase in pore saturation after 100 minutes of soaking. The exposure time to heat treatment had less of an effect on PCM saturation. The exposure time had negligible effect above 30 min and above 500 °C heating temperatures. However, because the expanded graphite was found to oxidize around 700 °C, the use of longer exposure time in manufacturing applications can be beneficial if a shortened impregnation time is needed. Heat treatment conditions did not impact thermal conductivity. The composite latent heat of fusion was also reduced approximately proportionally to the PCM mass fraction. A local maximum in axial thermal conductivity was observed at around 83% porosity, which is similar to previous studies. The observed conductivity at this maximum was a factor of 81 times greater than the conductivity of the PCM.

36 MATERIALS SCIENCE↗
DOE OSTI2020

BCS-BEC crossover in a (t 2g ) 4 excitonic magnet

The condensation of spin-orbit-induced excitons in t 2 g 4 electronic systems is attracting considerable attention. At large Hubbard U , antiferromagnetism was proposed to emerge from the Bose-Einstein Condensation (BEC) of triplons ( J eff = 1 ). Furthermore, we show that even at intermediate U regimes, the spin-orbit exciton condensation is possible leading also to staggered magnetic order. The canonical electron-hole excitations (excitons) transform into local triplon excitations at large U , and this BEC strong coupling regime is smoothly connected to the intermediate U excitonic insulator region. We solved the degenerate three-orbital Hubbard model with spin-orbit coupling ( λ ) in one dimension using the density matrix renormalization group, while in two dimensions we use the Hartree-Fock approximation (HFA). Employing these techniques, we provide the full λ versus U phase diagrams for both one- and two-dimensional lattices. Our main result is that at intermediate Hubbard U , increasing λ at fixed U the system transitions from an incommensurate spin-density-wave metal to a Bardeen-Cooper-Schrieffer (BCS) excitonic insulator, with coherence length r coh of O ( a ) and O ( 10 a ) in 1 d and 2 d , respectively, with a being the lattice spacing. Further increasing λ , the system eventually crosses over to the BEC limit (with r coh ≪ a ).

1-dimensional systems↗
DOE OSTI2022

A linear surrogate for optimising functions of an orthogonal matrix with applications in wave function theory

The technique of surrogate optimisation is to use a simpler function to approximate a complex function that is time-consuming to evaluate. Here we show that the maximum of a special type of surrogate function f(U)=Tr(AU),UϵO(n) is at A T (AA T ) 1/2 , and that there is one and only one local maximum both in SO(n) and O(n)–SO(n). This function f(U) has been found to be useful in various aspects of electronic structure theory, including proving the Carlson-Keller theorem, and localising orbitals. As one other example, we apply it here to optimise the ground state of molecules using the Generalised Valence Bond wavefunction.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Statistics of Green's functions on a disordered Cayley tree and the validity of forward scattering approximation

The accuracy of the forward scattering approximation for two-point Green's functions of the Anderson localization model on the Cayley tree is studied. A relationship between the moments of the Green's function and the largest eigenvalue of the linearized transfer-matrix equation is proved in the framework of the supersymmetric functional-integral method. The new large-disorder approximation for this eigenvalue is derived and its accuracy is established. Using this approximation the probability distribution of the two-point Green's function is found and compared with that in the forward scattering approximation (FSA). It is shown that FSA overestimates the role of resonances and thus the probability for the Green's function to be significantly larger than its typical value. The error of FSA increases with increasing the distance between points in a two-point Green's function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Uncertainty quantification of an empirical shell-model interaction using principal component analysis

Recent investigations have emphasized the importance of uncertainty quantification (UQ) in nuclear theory. Here, we carry out UQ for configuration-interaction shell-model calculations in the 1$\textit{s}$–0$\textit{d}$ valence space, investigating the sensitivity of observables to perturbations in the 66 parameters (matrix elements) of a high-quality empirical interaction. The large parameter space makes computing the corresponding Hessian numerically costly, so we compare a cost-effective approximation, using the Feynman-Hellmann theorem, to the full Hessian and find it works well. Diagonalizing the Hessian yields the principal components of the interaction: linear combinations of parameters ordered by sensitivity. This approximately decoupled distribution of parameters facilitates theoretical uncertainty propagation onto structure observables: electromagnetic transitions, Gamow-Teller decays, and dark matter-nucleus scattering matrix elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Low and moderate x gluon contribution to exclusive Compton scattering processes

We revisit the high energy semi-classical description of the exclusive processes DVCS, TCS, and Double DVCS by explicitly keeping track of the Feynman x dependence in both the hard and the hadronic matrix elements. This is achieved by a modification of the standard shock wave approximation to derive the effective Feynman rules, which leads to a generic expression on which we then perform a partial twist expansion to get rid of quantities suppressed by the proper physical scales. We obtain a compact factorized master formula that can be used to investigate the Bjorken limit at leading twist. In particular, we recover the full one-loop result in the collinear limit for pure gluon exchange with the target. Finally, we discuss the subtleties in taking the simultaneous collinear and small x limit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Uncovering the role of nanoscale precipitates on martensitic transformation and superelasticity

Here we characterize the role of coherent nanoscale B2 Ni 50 Al 50 precipitates on the temperature- and stress-induced martensitic phase transformation in nanocrystalline Ni 63 Al 37 shape memory alloys using multi-million-atoms molecular dynamics (MD) simulations. We studied two types of precipitates: one with single crystal precipitates (SXP) and a second where grain boundaries cut through precipitates (PXP). Simulations reveal that the presence of B2 precipitates stabilizes the cyclic flag-shaped stress-strain response, characteristic of superelasticity, and reduces remnant strain. In contrast, single-phase nanocrystalline Ni 63 Al 37 exhibits degradation of the reverse transformation during cyclic loading and, eventually, incomplete reversible transformation within a few cycles. This is consistent with previous experimental findings of ultra-low fatigue in Ni-Ti-Cu alloys with Ti 2 Cu precipitates. The simulations reveal that the presence of precipitates significantly improves the reversibility of the transformation by acting as elastic zones that partially shield the martensitic transformation and drive the reverse transformation. A detailed analysis of the MD trajectories reveals that the martensitic transformation of the matrix induces ultra-large elastic deformation in some of the B2 precipitates (approximately 12%) to the point of resulting in a martensite-like atomic structure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Direct ink writing of shear exfoliated two-dimensional nanomaterial- elastomeric multifunctional nanocomposite

Direct ink writing (DIW) of polymer nanocomposites with high loadings of two-dimensional (2D) nanofillers (graphene and hexagonal boron nitride (hBN)) is challenging because of potential clogging, use of hazardous solvents, and agglomeration. Here, in this work, a shear exfoliation and sieving method to prepare DIW ink with high loading of nanofillers produced from low-cost bulk layered materials such as graphite and bulk hBN powder for successful DIW printing without the use of any solvents, binders, or plasticizers. The single-step exfoliation technique resulted in a composite with substantial layer reduction along the c-axis, as confirmed by SEM, TEM, XRD, and Raman analysis. Incorporating exfoliated graphene (40 wt%) increased viscosity by ∼6 orders of magnitude due to enhanced particle–matrix interactions, leading to pronounced yield stress behavior and a yield stress of approximately 1598 Pa, which enabled excellent shape retention during extrusion. Using the DIW technique, porous structures such as desalination membranes, self-sensing bone scaffolds, thermal management coating, and serpentine strain sensors were fabricated. When tested in a direct contact membrane distillation setup, the fabricated membrane demonstrated a promising permeate flux of 21.85 Lm −2 h −1 and a salt rejection of 74.3 %. The fabricated serpentine sensor exhibited stable signal variations under cyclic tensile loading, with a working range of 0–200 % strain and a maximum gauge factor of 43,735. A cell culture test using the printed bone scaffold demonstrated promising cell attachment and proliferation. The DIW printed hBN nanocomposite exhibited reversible shape change under heat, demonstrating potential 4D printing capability and efficient thermal management when exposed to high heat or flame.

Desalination↗
DOE OSTI2022

CRISPR-based engineering of phages for in situ bacterial base editing

Investigation of microbial gene function is essential to the elucidation of ecological roles and complex genetic interactions that take place in microbial communities. While microbiome studies have increased in prevalence, the lack of viable in situ editing strategies impedes experimental progress, rendering genetic knowledge and manipulation of microbial communities largely inaccessible. Here, we demonstrate the utility of phage-delivered CRISPR-Cas payloads to perform targeted genetic manipulation within a community context, deploying a fabricated ecosystem (EcoFAB) as an analog for the soil microbiome. First, we detail the engineering of two classical phages for community editing using recombination to replace nonessential genes through Cas9-based selection. We show efficient engineering of T7, then demonstrate the expression of antibiotic resistance and fluorescent genes from an engineered λ prophage within an Escherichia coli host. Next, we modify λ to express an APOBEC-1-based cytosine base editor (CBE), which we leverage to perform C-to-T point mutations guided by a modified Cas9 containing only a single active nucleolytic domain (nCas9). We strategically introduce these base substitutions to create premature stop codons in-frame, inactivating both chromosomal ( lacZ ) and plasmid-encoded genes (mCherry and ampicillin resistance) without perturbation of the surrounding genomic regions. Furthermore, using a multigenera synthetic soil community, we employ phage-assisted base editing to induce host-specific phenotypic alterations in a community context both in vitro and within the EcoFAB, observing editing efficiencies from 10 to 28% across the bacterial population. The concurrent use of a synthetic microbial community, soil matrix, and EcoFAB device provides a controlled and reproducible model to more closely approximate in situ editing of the soil microbiome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Photoinduced anomalous Hall effect in two-dimensional transition metal dichalcogenides

A circularly polarized ac pump field illuminated near resonance on two-dimensional transition metal dichalcogenides (TMDs) produces an anomalous Hall effect in response to a dc bias field. In this work, we develop a theory for this photoinduced anomalous Hall effect in undoped TMDs irradiated by a strong coherent laser field. The strong field renormalizes the equilibrium bands and opens up a dynamical energy gap where single-photon resonance occurs. The resulting photon dressed states, or Floquet states, are treated within the rotating-wave approximation. A quantum kinetic equation approach is developed to study the nonequilibrium density matrix and time-averaged transport currents under the simultaneous influence of the strong ac pump field and the weak dc probe field. Dissipative effects are taken into account in the kinetic equation that captures relaxation and dephasing. The photoinduced longitudinal and Hall conductivities display notable resonant signatures when the pump field frequency reaches the spin-split interband transition energies. Rather than valley polarization, we find that the anomalous Hall current is mainly driven by the intraband response of photon-dressed electron populations near the dynamical gap at both valleys, accompanied by a smaller contribution due to the interband response. In conclusion, these findings highlight the importance of photon-dressed bands and nonequilibrium distribution functions in achieving a proper understanding of the photoinduced anomalous Hall effect in a strong pump field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Nanoscale elemental mapping reveals the effect of light and temperature on Gephyrocapsa huxleyi coccoliths

Coccolithophores produce vast quantities of stable calcium carbonate coccoliths, occupying the distinctive position as one of the most productive calcifying organisms on Earth. Gephyrocapsa huxleyi, the most abundant species of coccolithophores across modern oceans, plays an important role in the biological carbon pump and may serve as inspiration for nature-based carbon capture and sequestration technologies. Despite intensive study into G. huxleyi, the biological mechanisms of coccolith synthesis remain unclear. Here, we employ synchrotron X-ray fluorescence to quantitatively map elemental incorporation (calcium, Ca; strontium, Sr; sulfur, S; iron, Fe; and silicon, Si) into G. huxleyi coccoliths at 40-nanometer resolution. Trace elements play a role in determining the thermodynamic stability of coccolith carbonate sequestration and may provide insight into unknown mechanisms of coccolithogenesis. Although the coccolith Sr/Ca ratio has been widely applied as a paleoceanographic proxy, understanding of the Sr/Ca phenomena is limited. Our results confirm that Sr/Ca is temperature dependent but challenge the hypothesis that growth rate is a primary determinant. We find that increased light irradiance led to an ~30% decrease in calcium incorporated into the coccoliths. Surprisingly, X-ray fluorescence spectra and mapping revealed evidence of sulfur, at approximately three times the concentration of Sr, embedded in and around the calcite matrix, suggesting that G. huxleyi coccoliths may play an expanded role in the global sulfur cycle. Here, based on our elemental maps, we postulate that sulfur is active in the intracellular coccolith vesicle and plays a direct role in formation of G. huxleyi coccoliths.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Advancing the Theory of Nuclear Data Evaluations [Abstract]

We present recent advances in the R-matrix formalism as well as the Bayesian evaluation framework for improved nuclear data evaluations. The advances in the R matrix formalism include: 1) direct processes, 2) doorway, as well as multistep, processes, and 3) various forms of the Reich-Moore approximation for eliminated capture channels. Furthermore, to address unreasonably small posterior uncertainties often encountered in nuclear data evaluations of large data sets using the conventional form of the Bayes’ theorem, we introduce imperfections (of the data or the model) as a formal evaluation tool for taming the evaluated uncertainties in harmony with Bayes’ theorem. These theoretical advances were motivated by the nuclear data evaluations of differential resolved resonance cross section data using the code SAMMY, as well as the integral benchmark experiments using the SCALE code system, being performed at Oak Ridge National Laboratory for the Nuclear Criticality Safety Program. Some pedagogical applications of the new formalism, as well as a snapshot of the SAMMY modernization efforts, will be presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

R-matrix calculations for opacities: I. Methodology and computations

Abstract An extended version of the R -matrix methodology is presented for calculation of radiative parameters for improved plasma opacities. Contrast and comparisons with existing methods primarily relying on the distorted wave approximation are discussed to verify accuracy and resolve outstanding issues, particularly with reference to the opacity project (OP). Among the improvements incorporated are: (i) large-scale Breit–Pauli R -matrix calculations for complex atomic systems including fine structure, (ii) convergent close coupling wave function expansions for the ( e + ion) system to compute oscillator strengths and photoionization cross sections, (iii) open and closed shell iron ions of interest in astrophysics and experiments, (iv) a treatment for plasma broadening of autoionizing resonances as function of energy-temperature-density dependent cross sections, (v) a ‘top-up’ procedure to compare convergence with R -matrix calculations for highly excited levels, and (vi) spectroscopic identification of resonances and bound ( e + ion) levels. The present R -matrix monochromatic opacity spectra are fundamentally different from OP and lead to enhanced Rosseland and Planck mean opacities. An outline of the work reported in other papers in this series and those in progress is presented. Based on the present re-examination of the OP work, opacities of heavy elements might require revisions in high temperature-density plasma sources.

Pradhan, A. K. (ORCID:0000000187753643)↗
DOE OSTI2020

Spin–Orbit Matrix Elements for a Combined Spin-Flip and IP/EA approach

We introduce a practical approach for computing the Breit–Pauli spin–orbit matrix elements of multiconfigurational systems with both spin and spatial degeneracies based on our recently developed RAS-nSF-IP/EA method (Houck, S. E.; et al. J. Chem. Theory Comput. 2019, 15, 2278). The spin–orbit matrix elements over all the multiplet components are computed using a single one-particle reduced density matrix as a result of the Wigner–Eckart theorem. A mean field spin–orbit approximation was used to account for the two-electron contributions. Basis set dependence as well as the effect of including additional excitations is presented. The effect of correlating the core and semicore orbitals is also examined. Surprisingly accurate results are obtained for spin–orbit coupling constants, despite the fact that the efficient wave function approximations we research neglect the bulk of dynamical correlation.

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

Toward more accurate adiabatic connection approach for multireference wavefunctions

A multiconfigurational adiabatic connection (AC) formalism is an attractive approach to compute the dynamic correlation within the complete active space self-consistent field and density matrix renormalization group (DMRG) models. Practical realizations of AC have been based on two approximations: (i) fixing one- and two-electron reduced density matrices (1- and 2-RDMs) at the zero-coupling constant limit and (ii) extended random phase approximation (ERPA). This work investigates the effect of removing the “fixed-RDM” approximation in AC. The analysis is carried out for two electronic Hamiltonian partitionings: the group product function- and the Dyall Hamiltonians. Exact reference AC integrands are generated from the DMRG full configuration interaction solver. Two AC models are investigated, employing either exact 1- and 2-RDMs or their second-order expansions in the coupling constant in the ERPA equations. Calculations for model molecules indicate that lifting the fixed-RDM approximation is a viable way toward improving the accuracy of existing AC approximations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗