Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “test matrix”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The drivers and predictability of wildfire re-burns in the western United States (US)

Evidence is mounting that the effectiveness of using prescribed burns as a management tactic may be diminishing due to the higher incidence of wildfire re-burns. The development of predictive models of re-burns is thus essential to better understand their primary drivers so that forest management practices can be updated to account for these events. First, we assess the potential for human activity as a driver of re-burns by evaluating re-burn trends both within and outside of the wildland–urban interface (WUI) of the western US. Next, we investigate the predictability of re-burns through the application of both random forest and the explanatory machine learning non-negative matrix factorization using k-means clustering (NMFk) algorithms to predict re-burn occurrence over California based on a number of climate factors. Our findings indicate that while most states showed increasing trends within the WUI when trends were conducted over longer moving windows (e.g. 20 years), California was the only state where the rate of increase was consistently higher in the WUI, indicating a stronger potential for human activity as a driver in that location. Furthermore, we find model performance was found to be robust over most of California (Testing F1 scores = 0.688), although results were highly variable based on EPA level III Ecoregion (F1 scores = 0.0–0.778). Insights provided from this study will lead to a better understanding of climate and human activity drivers of re-burns and how these vary at broad spatial scales so that improvements in forest management practices can be tuned according to the level of change that is expected for a given region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

An ICME Modeling Framework for Titanium/Tungsten-Carbide Metal Matrix Composites

This report describes a collaborative project to develop a validated, predictive model for the high temperature mechanical properties of a titanium-matrix, tungsten-carbide/cobalt-reinforced metal matrix composite. The modeling approach was to first develop a detailed, microstructural model linking the material structure and the interfacial debonding properties to the effective properties of the material. The project then completed a throughput simulation campaign to generate a large number of simulations for discrete microstructures and different debonding parameters. Finally, the project trained a fast, Gaussian process surrogate model against this simulation database to provide a quick model linking the material compositions, structure, and processing parameters to the resulting material properties. This model was validated against high temperature tensile test data on a few particular composite compositions. The tests validate the model predictions for ultimate tensile strength and uniform elongation/ductility, meaning the final surrogate model can now be used to tune the material composition and processing parameters to identify optimal composite compositions for particular applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Microstructural and Micro-Chemical Evolutions in the Irradiated UCO Fuel Kernels of AGR-1 and AGR-2 TRISO Fuel Particles

AGR-1 and AGR-2 TRISO fuel particles were fabricated with slightly different fuel kernel chemical compositions, modified fabrication processes, different fuel kernel diameters, and changed 235U enrichments. To correlate those differences with the fuel kernel responses to neutron irradiations in terms of irradiated fuel microstructure, fission products chemical and physical states, and fission gas bubble evolutions, extensive microstructural and analytical characterizations were conducted. The studies used a state of art transmission electron microscopy (TEM) equipped with Energy-dispersive X-ray spectroscopy (EDS) of four silicon solid-state detectors which have super sensitivity and fast speed. The TEM specimens were prepared from selected AGR-1 and AGR-2 irradiated fuel kernels exposed to safety testing after irradiation. The particles were chosen to represent a representative irradiation conditions with a fuel burnup within the range from 10.8 to 18.6% FIMA, and the time-average volume-average temperatures vary from 1070 to 1287°C. The 235U enrichment was 19.74 wt.% for the AGR1 fuel kernels and 14.03 wt.% for the AGR-2 fuel kernels. The TEM results show that there were significant microstructural reconstructions in the irradiated fuel kernels for both the AGR-1 and AGR-2 fuels. There are four major phases including fuel matrix of UO2 and UC, U2RuC2, and UMoC2 in the irradiated AGR2 fuel kernel. Zr and Nb form solid solution in the UC phase. UMoC2 phase often shows a detectable concentration of Tc. Pd was found to mainly locate in the buffer layer or to be associated with fission gas bubble within the UMoC2 phase. The EDS maps qualitatively show that the rare-earth fission products (Nb, et al.) preferentially reside in the UO2 phase. In contrast, in the irradiated AGR1 fuel kernel, no U2RuC2 or UMoC2 precipitates were positively identified. Instead, there is a high number of rod-shape precipitates enriched with Ru, Tc, Rh, and Pd observed in the fuel kernel center and edge zone. The difference of microstructural and micro-chemical evolutions in irradiated fuel kernels between the AGR-1 and AGR-2 TRISO fuel particle may result from a combined factor of irradiation temperature, fuel geometry and chemical composition. However, the irradiation temperature probably play a more deterministic role. Limited electron energy loss spectroscopy (EELS) characterizations on the AGR2 fuel kernel show that there is nearly no carbon in the UO2 phase while a small fraction of oxygen was detected in the UC/UMoC2 phase.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Exploring viral particle, soil, and extraction buffer physicochemical characteristics and their impacts on extractable viral communities

Viruses are expected to be pivotal members of soil ecosystems, and recent advances in viral size fraction metagenomic (viromic) approaches have improved our ability to interrogate soil viral ecology. However, viromics relies on extraction buffers to effectively desorb viral particles from the soil matrix for downstream analysis, and viral extraction efficiency could be affected by the interplay between viral particles, soils, and extraction buffer chemistry. Here, we investigated whether extraction buffer chemistry affected extractable viral community composition measured by viromics from different soil types, for both biological (samples collected 1 m apart) and technical (subsamples from the same soil homogenate) replicates. We investigated pH manipulation of protein-supplemented phosphate-buffered saline buffer (PPBS, pHs 4.5, 5.5, 6.5, and 7.5) on forest, grassland and wetland soils that exhibited different soil edaphic properties, and we tested different buffer chemistries (PPBS, Carbonated Buffer, Glycine, and Saline Magnesium) on just the wetland soil. Spatial distance (i.e., biological replicate) was the primary driver of extractable viral community composition across all buffers and soils tested. Differences in viral community composition according to extraction buffer properties were only observed in the grassland technical replicates at PPBS buffer pH 4.5, and in both the wetland technical and biological replicates treated with different buffer chemistries, but the effects of buffer chemistry were secondary to spatial distance in the biological replicates. The lack of buffering capacity in the grassland soil technical replicates likely increased sorption of viral particles at pH 4.5, but neither protein composition nor isoelectric point explained this phenomenon. In conclusion, given that most soil viral ecological studies to date include sample collection over distances much farther apart than the 1-m distances considered here, results suggest that extraction buffer chemistry is likely of much lower importance than ecological considerations, such as spatial distance, in the design of future soil viral ecological investigations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Performance Restoration of a Tritium-Aged LaNi 4.25 Al 0.75 Sample

Hydride beds containing LaNi 4.25 Al 0.75 (LANA.75) are used to store significant quantities of tritium. These hydride beds have a limited service life due to radiolytic decay of tritium to 3He within the metal matrix. The crystal structure of the hydride is altered by trapped 3 He, which has a very low solubility in the metal. The altered structure induces the formation of a heel of trapped hydrogen isotopes and diminishes the reversible capacity of the hydride. With sufficient tritium exposure, the bed loses the ability to deliver 3He-free tritium, and replacement is needed. Demonstration of a means to regenerate tritium-aged LANA.75 in situ would delay or even eliminate the need to replace lanthanum nickel aluminum (LANA) hydride beds. This paper presents test results obtained during regeneration testing. The efficacy of regeneration testing was evaluated by comparing tritium desorption isotherms collected on the hydride before and after exposure to regeneration conditions. Testing was performed on a bench-scale tritium-aged LANA.75 sample that was previously isotopically exchanged (from tritium to deuterium), passivated, and recovered. Once transferred to a high-temperature test cell, the deuterium heel of the sample was isotopically exchanged with tritium, and a baseline desorption isotherm was collected for comparison purposes. The sample was then heated under vacuum, and comparative isotherms were gathered between regeneration evolutions. Shifts in isotherms show progressive improvements with higher-temperature exposure over the tritium-aged baseline. The heel was significantly reduced, and the reversible capacity of the hydride was essentially restored to near virgin values. For all tested conditions, the plateau pressure remained higher than virgin LANA.75.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Tuning the temperature range of superelastic Ni-Ti alloys for elastocaloric cooling via thermal processing

Abstract Caloric cooling enlisting solid-state refrigerants is potentially a promising eco-friendly alternative to conventional cooling based on vapor compression. The most common refrigerant materials for elastocaloric cooling to date are Ni-Ti based superelastic shape memory alloys. Here, we have explored tuning the operation temperature range of Ni 50.8 Ti 49.2 for elastocaloric cooling. In particular, we have studied the effect of thermal treatments (a.k.a. aging) on the transformation temperature, superelasticity, and elastocaloric effects of Ni 50.8 Ti 49.2 shape memory alloy tubes. The isothermal compressive test revealed that the residual strain of thermally-treated Ni-Ti tubes at room temperature approaches zero as aging time is increased. Short-time aging treatment at 400 °C resulted in good superelasticity and elastocaloric cooling performance with a large tunable austenite finish ( A f ) temperature range of 24.7 °C, as determined from the A f temperature of the samples that were aged 5–120 min. The main reason of the property change is the formation of a different amount of Ni 4 Ti 3 precipitates in the NiTi matrix. Our findings show that it is possible to tailor the A f temperature range for development of cascade elastocaloric cooling systems by thermally treating a starting single composition Ni-Ti alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spin Alignment of Vector Mesons in Heavy-Ion Collisions

Polarized quarks and antiquarks in high-energy heavy-ion collisions can lead to the spin alignment of vector mesons formed by quark coalescence. Here, using the relativistic spin Boltzmann equation for vector mesons derived from Kadanoff-Baym equations with an effective quark-meson model for strong interaction and quark coalescence model for hadronizaton, we calculate the spin density matrix element $ρ_{00}$ for Φ mesons and show that anisotropies of local field correlations with respect to the spin quantization direction lead to Φ meson's spin alignment. We propose that the local correlation or fluctuation of Φ fields is the dominant mechanism for the observed Φ meson's spin alignment and its strength can be extracted from experimental data as functions of collision energies. The calculated transverse momentum dependence of $ρ_{00}$ agrees with STAR's data. We further predict the azimuthal angle dependence of $ρ_{00}$ which can be tested in future experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Aging of Limestones and Silane–Siloxane-Based Protective Hydrophobics: The Impact of Heating–Cooling and Freeze–Thaw Cycles

Stones are traditionally used in construction and architectural applications as building elements due to their aesthetic and technical/structural performance. Like other environmental factors (rain, humidity, moisture, salt presence, biological activity, etc.), heating–cooling and freeze–thaw cycles significantly threaten the longevity of stone materials. Hence, considering the socio-economic and cultural value of stones, preventive actions such as hydrophobic coatings are applied to prevent or mitigate damage. The scope of this study is the performance assessment of limestones with different characteristics and the efficiency of various commercial silane/siloxane-based hydrophobic coatings when exposed to thermal variation and freeze–thaw. For that purpose, the standards EN 14066:2013 (determination of resistance to aging by thermal shock) and EN 12371:2010 (determination of frost resistance) were followed. Open porosity and static contact angles were estimated to assess the stone durability and water protection capabilities of the hydrophobics. Additionally, sound speed propagation velocity, quality of building material index, elastic modulus and flexural strength were measured to evaluate the variation of mechanical properties. Static contact angle revealed that the coatings maintained an efficient level of hydrophobicity even after thermal-shock and freeze–thaw weathering tests. The study also revealed a critical interaction between freeze–thaw cycles, hydrophobic coatings and structural integrity of the stones, mostly on more porous ones. When they are subjected to harsh environmental conditions, untreated porous limestones keep structural cohesion, allowing for the natural absorption and release of water during freezing and thawing. On the contrary, when limestones are treated, the hydrophobic coatings can moderately obstruct the water release due to the partial saturation of the porous framework by the products. It also probably resulted from the different mechanical behavior between the inner matrix and layer of stone coated, resulting in a premature breakout and mechanical damage of the stone.

Lisci, Carla (ORCID:000000032179130X)↗
DOE OSTI2022

Testing Lepton Flavor Universality with Pion, Kaon, Tau, and Beta Decays

We present an overview of searches for violation of lepton flavor universality with a focus on low energy precision probes using π, K, τ, and nuclear beta decays. We review the current experimental results, summarize the theoretical status within the context of the Standard Model, and discuss future prospects (both experimental and theoretical). We review the implications of these measurements for physics beyond the Standard Model by performing a global model-independent fit to modified W couplings to leptons and four-fermion operators. We also discuss new physics in the context of simplified models and review Standard Model extensions with a focus on those that can explain a possible deviation from unitarity of the Cabibbo–Kobayashi–Maskawa quark mixing matrix.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Wigner time delay in photoionization: a 1D model study

Abstract In scattering theory, the Wigner–Smith time delay, calculated through a phaseshift derivative or its multichannel generalization, has been demonstrated to measure the amount of delay or advance experienced by colliding particles during their interaction with the scattering potential. Fetic, Becker, and Milosevic argue that this concept cannot be extended to include photoionization, viewed as a half-scattering experiment. Their argument is based on the lack of information about scattering phaseshifts in the part of the wavefunction (satisfying the ingoing-wave boundary condition) going to the detector. This article aims to test this claim by examining a photodetachment process in a simple 1D model with a short-range symmetrical potential. Using time-dependent perturbation theory with a dipole interaction, the relevant wavepacket of the outgoing particle is analyzed and compared to the free wavepacket as a reference. Our findings confirm that, indeed, a time delay arises in the liberated fragmentation wavepacket, which is expressed as an energy derivative of the scattering phaseshift. Our study highlights that the source of the phaseshift content in the wavepacket arriving at the detector is the dipole matrix element, which is a direct consequence of imposing the ingoing-wave boundary condition. We illustrate our results through numerical simulations of both the non-free and free wavepackets. The amount of the observed time delay is found to be half of that appearing in a typical scattering experiment.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

MINE: maximally informative next experiment—toward a new GWAS experimental design and methodology

Abstract The computational methodology of Genome Wide Association Studies (GWAS) currently has several limitations: (i) the number of observations (rows) on a quantitative trait tends to be smaller than the number of single nucleotide polymorphisms (SNPs) (columns) in the design matrix; (ii) each SNP is usually modeled separately, failing to acknowledge interaction between each other (ie epistasis); (iii) there is implicit linkage disequilibrium (LD) between neighboring SNPs due to their linkage. To overcome these issues, we developed a tool that uses ensemble methods to fit mixed linear models to GWAS data, and these ensemble methods include the development of a new experimental design approach in GWAS, which uses the resultant models and data to select the next informative experiment over time. This new adaptive and staged approach for GWAS experimental design was developed and tested in a 3 yr adaptive model-guided discovery experiment against a fixed classical design. In Sorghum bicolor a total of 79, 86, and 78 accessions were tested in years 1, 2, and 3, respectively out of 343 accessions available in the Bioenergy Association Panel (BAP) each identified for 232,303 SNPs, 1 every 2–3 kb in the genomes. We demonstrated the feasibility of MINE enacted with 8 people in the field per year over 3 yr vs in 1 large classical design enacted with 20 people in 1 yr. The MINE results for chromosomal regions identified controlling dry weight were confirmed against results from previous sorghum GWAS experiments and 1 large classical design for the BAP panel.

Genetics & Heredity↗
DOE OSTI2024

Understanding EBC Lifetimes and Performance for Industrial Gas Turbines

Hydrogen or hydrogen blend fuels are expected to replace natural gas in land-based industrial gas turbines (IGTs) to support a greener power economy. Silicon carbide (SiC) base ceramic matrix composites (CMCs) are considered for replacement of Ni-based superalloys to facilitate future efficiency improvements. SiC CMCs require environmental barrier coatings (EBCs) to mitigate volatilization from high-temperature steam, thus making the EBC lifetime critical information for identifying CMC component lifetimes. The goal of this project is to determine the maximum bond coating temperature underneath the EBC for achieving an IGT component lifetime goal of 25,000 h, which is far greater than current CMC component lifetime requirements for aero-turbine applications. To provide data for the lifetime model, laboratory testing used plasma-sprayed rare-earth silicate EBCs on monolithic SiC substrates with an intermediate Si bond coating. Specimens exposed to 1-h thermal cycles in flowing air-steam environments and reaction kinetics were assessed from 700°-1350°C by measuring the thickness of the thermally grown silica scales. The silica growth and phase transformation appear critical in predicting EBC lifetime and several strategies have been explored to reduce the oxide growth rate and improve EBC durability at elevated temperatures. Advanced characterization using Raman spectroscopy has helped clarify this system.

Ridley, Mackenzie↗
DOE OSTI2020

A Rapid Sequential Separation Method for Determination of Actinides and TENORM in Fracking Wastes - 20223

With recent advances in unconventional drilling technology in the US and around the world, the risks of environmental contamination and exposure due to technologically enhanced naturally occurring radioactive materials (TENORM) from fracking wastes have also surfaced. Although the Permian Basin (southeastern New Mexico, western Texas) has long been a modest producer of oil, the advent of fracking technology a decade ago turned it into hot property for producers. A 2018 assessment of undiscovered, technically recoverable continuous oil and gas resources by the U.S. Geological Survey estimates an average 46.3 billion barrels of oil and 281 trillion cubic feet of gas. This is likely to increase the use of unconventional drilling and TENORM generation. TENORM concentrations in oil and gas exploration and production waste can be several hundred to several thousand times higher than background TENORM concentrations. Both the thorium and uranium decay chains contribute to airborne radionuclides arising from the radon gas escaping the ground and subsequently decaying as airborne particulates. Additionally, radon decay products {sup 210}Pb and {sup 210}Po can build up in scale on the internal surfaces of oil and gas handling pipes and in sludge in refineries, becoming potential inhalation and ingestion hazards for workers. Furthermore, southeastern NM is also home to world's only licensed and operating transuranic nuclear waste repository, the DoE's Waste Isolation Pilot Plant (WIPP). Plutonium isotopes ({sup 239+240}Pu) and {sup 241}Am, are expected to account for more than 99% of the total radioactivity scheduled for disposal in the WIPP repository. Thus, an improved understanding of the environmental fate and transport of TENORM, liberated by unconventional drilling, is essential to assess how best to protect individuals and the environment. In this context, accurate measurement of TENORM and actinides in environmental and biological samples is essential. In this presentation, a new sequential method for the separation and pre-concentration of actinides (Pu, Am, Np) and TENORM (Po, U and Th) derived from oil and gas exploration is proposed. The TEVA method involves a rapid co-precipitation step to remove matrix interferences from the samples, followed by plutonium oxidation state adjustment to Pu (IV) and an incubation period of ∼ 1 hour at 50-60 deg. C to allow the resultant Po (II) to oxidize into Po(IV). The polonium, neptunium, thorium and plutonium are then separated on a TEVA column, while americium and uranium are separated on a TRU column. After separation, the alpha counting source was prepared by micro-precipitation with copper sulfide (CuS) for polonium and neodymium fluoride (NdF{sub 3}) for actinides. The efficiency and reliability of the procedures were tested by analyzing filter, drinking water and frack sand samples. The method is simple, robust and can be performed quickly with excellent removal of interferences, high chemical recovery and very good alpha peak resolution. (authors)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Integrated TBC/EBC for SiC Fiber Reinforced SiC Matrix Composites for Next Generation Gas Turbines (Final Report)

In this ambitious project, a significant challenge preventing widespread use of ceramic matrix composites in hot section components of gas power generating turbines was addressed. Hot section components made from SiC fiber reinforced SiC matrix composites (SiC/SiC composites), have the potential to achieve a combined cycle energy efficiency of above 65%. To achieve this efficiency, the temperature of the combustion gases on the turbine component surface is expected to reach as high as ~1700ºC (3100ºF). However, SiC/SiC composites suffer from volatilization at high temperature under high velocity steam environment – typical in a turbine hot chamber. This integrated collaborative project, between two research groups at Clemson University and a team of engineers and scientists at GE Gas Power, focused on design, fabrication and evaluation of a new graded environmental barrier coating. Using the excellent compositional tailorability of polymer derived ceramics (PDCs), we developed an integrated EBC/BC that smoothly transitions from bottom SiC/SiC bond interface to hermetic sesquioxide (Y 2 O 3 ) top surface. This smooth compositional transition avoided the sharp CTE mismatch at interfaces to greatly reduce the thermal stress. The developed EBC is compatible with and bonds well with the bond coat that GE uses. The top hermetic oxide surface ensures low oxygen transportation and volatilization rates under the high velocity steam environment. The graded microstructure with increasing sesquioxide concentration also decreases the amount of Si and C elements near the top hermetic sesquioxide layer. Coatings made using a process suitable for manufacturing – atmospheric plasma spraying was demonstrated. The microstructure of coatings made using this process were fully characterized. Finally, the performance of these coatings was investigated under the most demanding conditions for this application. Specifically, high temperature, high velocity steam exposure; thermal shock test; and thermal cycling tests. The graded coatings performed exceptionally well in all these tests. Specific suggestions for future development in this area have also been identified. In addition to the excellent technical progress, this project helped initiate a robust collaboration between the research groups at Clemson and GE Gas Power. It also provided a high-quality learning and mentoring opportunity for the three post-doctoral research associates and two UG students at Clemson University who worked on this project. All proposed tasks, sub-tasks were completed, and all milestones achieved.

03 NATURAL GAS↗
DOE OSTI2024

Brittle failure analysis and modeling of high-burnup PWR fuel cladding alloys

The aim of this research is the development of methods for predicting mechanical behavior and identification of limiting conditions to prevent brittle failure of high-burnup (HBU) pressure water reactor (PWR) fuel cladding alloys. A finite element (FE) model of the ring compression test (RCT) was created to analyze the failure behavior of zirconium-based alloys with radial hydrides during the RCT. An elastic-plastic material model describes the zirconium alloy. The stress-strain curve needed for the elastic-plastic material model was derived by inverse finite element analyses. Cohesive zone modeling is used to reproduce sudden load drops during RCT loading. Based on the failure mechanism in non-irradiated ZIRLO (R) claddings, a micro-mechanical model was developed that distinguishes between brittle failure along hydrides and ductile failure of the zirconium matrix. Two different cohesive laws representing these types of failure are present in the same cohesive interface. The key differences between these constitutive laws are the cohesive strength, the stress at which damage initiates, and the cohesive energy, which is the damage energy dissipated by the cohesive zone. Statistically generated matrix-hydride distributions were mapped onto the cohesive elements and simulations with focus on the first load drop were performed. Computational results are in good agreement with the RCT results conducted on high-burnup M5 (R) samples. It could be shown that crack initiation and propagation strongly depend on the specific configuration of hydrides and matrix material in the fracture area.

Simbruner, Kai↗
DOE OSTI2025

Leveraging explainable AI to characterize floating-point exceptions in linear solvers

Linear solver packages are central to many scientific, engineering, and machine learning applications. When floating-point exceptions occur in these solvers, e.g., division by zero or overflow, numerical results are compromised and become unreliable. Existing static and dynamic analysis tools can detect such exceptions, but they do not explain why the exceptions occur in terms of the solver inputs. Here, we present a study to characterize the inputs that cause numerical exceptions in linear solver packages. Our approach uses explainable AI (XAI) to find the most relevant characteristics of input matrices that explain the occurrence of exceptions in the solvers. Since training data in this domain is scarce, we perform extensive data gathering and data augmentation to obtain exception-inducing inputs. Our approach uses a repair strategy on the features blamed by XAI to validate that such features indeed explain the exceptions. We compare the LIME and SHAP XAI techniques using a dozen matrix features with three classifiers. We evaluate the approach on three widely used linear solver packages and find that some input characteristics can explain the occurrence of exceptions 100% of the time, in specific solvers and preconditioners.

Explainable AI↗
DOE OSTI2023

Sparsity of the electron repulsion integral tensor using different localized virtual orbital representations in local second-order Møller–Plesset theory

Utilizing localized orbitals, local correlation theory can reduce the unphysically high system-size scaling of post-Hartree–Fock (post-HF) methods to linear scaling in insulating molecules. The sparsity of the four-index electron repulsion integral (ERI) tensor is central to achieving this reduction. For second-order Møller–Plesset theory (MP2), one of the simplest post-HF methods, only the (ia|jb) ERIs are needed, coupling occupied orbitals i, j and virtuals a, b. In this paper, we compare the numerical sparsity (called the “ragged list”) and two other approaches revealing the low-rank sparsity of the ERI. The ragged list requires only one set of (localized) virtual orbitals, and we find that the orthogonal valence virtual-hard virtual set of virtuals originally proposed by Subotnik et al. gives the sparsest ERI tensor. To further compress the ERI tensor, the pair natural orbital (PNO) type representation uses different sets of virtual orbitals for different occupied orbital pairs, while the occupied-specific virtual (OSV) approach uses different virtuals for each occupied orbital. Here, our results indicate that while the low-rank PNO representation achieves significant rank reduction, it also requires more memory than the ragged list. The OSV approach requires similar memory to that of the ragged list, but it involves greater algorithmic complexity. An approximation (called the “fixed sparsity pattern”) for solving the local MP2 equations using the numerically sparse ERI tensor is proposed and tested to be sufficiently accurate and to have highly controllable error. A low-scaling local MP2 algorithm based on the ragged list and the fixed sparsity pattern is therefore promising.

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

Coulomb excitation of the | T z | = 1 2 , A = 23 mirror pair

Background: Electric-quadrupole ( E 2 ) strengths relate to the underlying quadrupole deformation of a nucleus and present a challenge for many nuclear theories. Additionally, mirror nuclei in the vicinity of the line of N = Z represent a convenient laboratory for testing deficiencies in such models, making use of the isospin symmetry of the systems. Purpose: Uncertainties associated with literature E 2 strengths in Mg 23 are some of the largest in T z = | 1 2 | nuclei in the s d shell. The purpose of the present paper is to improve the precision with which these values are known, to enable better comparison with theoretical models. Methods: Coulomb-excitation measurements of Mg 23 and Na 23 were performed at the TRIUMF-ISAC facility using the TIGRESS spectrometer. They were used to determine the E 2 matrix elements of mixed E 2 / M 1 transitions. Results: Reduced E 2 transition strengths, B ( E 2 ) , were extracted for Mg 23 and Na 23 . Their precision was improved by factors of approximately 6 for both isotopes, while agreeing within uncertainties with previous measurements. Conclusions: A comparison was made with both shell-model and ab initio valence-space in-medium similarity renormalization group calculations. Valence-space in-medium similarity renormalization group calculations were found to underpredict the absolute E 2 strength, in agreement with previous studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗