Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting methods”

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

Ab Initio Methods for L-edge X-ray Absorption Spectroscopy

The theoretical prediction of X-ray absorption spectra (XAS) has become common- place in electronic structure theory. The ability to better model and understand L-edge spectra is of great interest in the study of transition metal complexes and a wide variety of solid state materials. However, until recently few rst-principles works have mod- eled L-edge XAS due to the presence of strong spin-orbit coupling in the 2p orbitals which splits the observed peaks into multiple groups of features. Therefore, a proper description of spin-orbit coupling is vital for the successful prediction of L-edge spectra. A number of new approaches that incorporate spin-orbit coupling have recently made advances in the computation of L-edge spectra. In this review, we describe recent work in computational L-edge XAS and how these methods may continue to improve in the future. Comparison of the advantages and disadvantages of the various approaches are considered, with special attention to not only the computational cost of the level of theory, but also the various approaches that can be used to compute the absorption spectra with a large number of high energy excited states.

Kasper, Joseph M.↗
DOE OSTI2023

Direct Laser Writing of Multimetal Bifunctional Catalysts for Overall Water Splitting

Water electrolysis is of interest as a sustainable way to produce clean hydrogen and oxygen fuel and help mitigate the rising problems of climate change while meeting global energy demands. High-efficiency, stable, and earth-abundant bifunctional catalysts are needed to enable more effective electrochemical cells for the hydrogen evolution reaction (HER) and oxygen evolution reaction (OER). Here, we investigate the synthesis, composition, performance, and mechanism of multimetal catalysts serving dual functionality in both OER and HER of water electrolysis. Through a laser synthesis method, we synthesized heterogeneous catalysts of nanocrystalline multimetallic alloy pockets embedded within an amorphous oxide matrix. We evaluated the performance and composition of a range of mixed transition-metal oxide materials for both OER and HER, ultimately synthesizing a Cr 0.01 Fe 0.27 Co 0.34 Ni 0.38 O x /C y catalyst that has a stable, high-rate, and competitive overall water splitting performance of 1.76 V at 100 mA cm –2 in an alkaline medium. Using density functional theory to gain insight as the active site and mechanism, we propose that the inclusion of a minor amount of Cr increases the degeneracy of energetic states that lowers the cost of forming the O 2 p–d bond and H 1 s–d bond due to the hybridization of s, p, and d orbitals from Cr. Using a two-electrode water electrolysis cell with a constant potential of 1.636 V to mimic the setup for fuel production, we found the catalyst to be stable at 14–15 mA cm –2 for 40 h. Furthermore, this laser synthesis method allowing for facile and rapid synthesis of complex multimetal systems demonstrates how doping a Fe, Co, and Ni heterogeneous amorphous/nanocrystalline structure with small amounts of Cr is important for bifunctional catalytic behavior, particularly for increasing HER functionality in advancing our understanding for future electrocatalytic design.

08 HYDROGEN↗
DOE OSTI2023

Design Considerations and Analysis of Multi-Level Erasure Coding in Large-Scale Data Centers

Multi-level erasure coding (MLEC) has seen large deployments in the field, but there is no in-depth study of design considerations for MLEC at scale. In this paper, we provide comprehensive design considerations and analysis of MLEC at scale. We introduce the design space of MLEC in multiple dimensions, including various code parameter selections, chunk placement schemes, and various repair methods. We quantify their performance and durability, and show which MLEC schemes and repair methods can provide the best tolerance against independent/correlated failures and reduce repair network traffic by orders of magnitude. To achieve this, we use various evaluation strategies including simulation, splitting, dynamic programming, and mathematical modeling. We also compare the performance and durability of MLEC with other EC schemes such as SLEC and LRC and show that MLEC can provide high durability with higher encoding throughput and less repair network traffic over both SLEC and LRC.

Wang, Meng↗
DOE OSTI2022

Low-Cost Preform and Molding Processes

The entry cost for prototyping a composite component for manufacture using automated, high rate processes is prohibitively expensive in many cases, especially for small business, where tooling costs may be several $100k. Discussions with industry also indicate that many small companies, tier 1 and 2 suppliers, have an interest to mold composite parts but do not want to deal with the capital cost, material handling issues, and labor associated with dry fiber preforming operations. While the molders may locate near the end user for logistics reasons, it may be more cost effective for the performer to remain regional and invest in capital equipment to support preform automation, thus keeping costs to a minimum. This project was designed to explore and demonstrate several options to meet these industry needs. Dry fiber preforming approaches were evaluated which allow for low pressure resin infusion, single sided tooling options such a vacuum assisted resin transfer molding (VARTM) or low pressure resin transfer molding (RTM-light). Unlike sheet molding compound, SMC compression molding where typical molding pressures of 1000 psi are required to push material into the desired location; positioning of a dry fiber preform into the desired location on the tool allows for low molding pressures of 10-50 psi. Lower molding pressures allow for use of low cost, additive fabrication of polymeric tooling. Polymeric tooling is suitable for rapid part prototyping and limited production. Dry fiber preforming approaches evaluated included use of commercial chopped strand mat, robotic chopper gun deposition, and continuous fiber preform augmentation using tailored fiber placement (TFP). Use of chopped strand mat does not require a robotic deposition method, however a cutting table is generally required and there is typically 20-30% scrap generation. While various fiber areal weights are available, the preform is not readily optimized for minimal fiber use or weight savings. In contrast, a robotic chopper gun approach allows for localized deposition where fiber is required to meet structural requirements. The robotic method is highly automated and minimizes fiber scrap, however the capital cost of the equipment and engineering labor for programming can result in higher preform cost compared to chopped strand mat in certain cases depending on preform complexity. Dry fiber preforming using the robotic chopper gun method allows for creation of three dimensional forms. This approach may be ideal for molding in-house, or if the preforms stack together densely to allow for efficient shipping. Applications evaluated for this program considered trade-off between fabrication of a fully 3D preform versus production of a flat preform which is designed to readily drape into the final desired shape. Such a preform design greatly simplifies robotic programming and requires no specialized tooling. The flat preforms are easily stacked and shipped to the final molding location. Flat preforms are much easier to augment with TFP continuous fiber to provide local reinforcement. The demonstration and evaluation of these preforming and tooling methods were completed on three component applications. The first application was a battery box cover for an electric vehicle which was highly three dimensional. The second demonstrator article was comprised of complex contours and was used to demonstrate the use of TFP and RTM-light molding process. The third demonstration article was the roof of an operator’s cab for large construction equipment. The roof is relatively flat however it is comprised of complex changes in thickness which clearly demonstrate the advantage of robotic chopper gun approach as compared to using numerous preform layers of chopped strand mat. The cost trades for the various preforming methods are summarized to help guide the reader as to preforming method considerations. Finally, these demonstrations all used glass fiber roving. A fourth, exploratory task was added to evaluate the ability to make preforms using Zoltek’s carbon fiber split tow roving. We were able to adapt the chopper gun to make flat preforms for laminate testing, but further development effort would be required to make suitable preforms.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sustainable Oxide Electrocatalyst for Hydrogen- and Oxygen-Evolution Reactions

In this work, an electrocatalyst for water splitting based on earth-abundant metals is reported. This perovskite-oxide catalyst, CaSrFe 0.75 Co 0.75 Mn 0.5 O 6-δ (CSFCM), is examined using both experimental and computational methods. It demonstrates a combination of properties, which include (a) very high activity for the oxygen-evolution reaction with an overpotential of eta = 0.19 V at 10 mA/cm 2 , (b) high stability over 1000 cycles of catalysis, (c) the ability to catalyze the hydrogen-evolution reaction effectively in both acidic and basic conditions, and (d) catalytic activity as a single-phase bulk material without the need for any additional processing, multicomponent composite preparation, or nanofabrication. Therefore, the catalytic activity of CSFCM is intrinsic, making it a good benchmark compound for future studies of electrocatalytic parameters. This work also highlights the impact of systematic structural design on electrocatalytic activity. Results from density functional theory calculations indicate that in addition to an optimal e g occupancy of similar to 1, an additional descriptor, i.e., maximizing the number of free e g carriers, correlates with the electrocatalytic activity.

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

Band structures and Z 2 invariants of two-dimensional transition metal dichalcogenide monolayers from fully relativistic Dirac-Kohn-Sham theory using Gaussian-type orbitals

Two-dimensional (2D) materials exhibit a wide range of remarkable phenomena, many of which owe their existence to the relativistic spin-orbit coupling (SOC) effects. To understand and predict properties of materials containing heavy elements, such as the transition-metal dichalcogenides (TMDs), relativistic effects must be taken into account in first-principles calculations. We present an all-electron method based on the four-component Dirac Hamiltonian and Gaussian-type orbitals (GTOs) that overcomes complications associated with linear dependencies and ill-conditioned matrices that arise when diffuse functions are included in the basis. Until now, there has been no systematic study of the convergence of GTO basis sets for periodic solids either at the nonrelativistic or the relativistic level. Here, in this study, we provide such a study of relativistic band structures of the 2D TMDs in the hexagonal (2H), tetragonal (1T), and distorted tetragonal (1T') structures, along with a discussion of their SOC-driven properties (Rashba splitting and Z 2 topological invariants). We demonstrate the viability of our approach even when large basis sets with multiple basis functions involving various valence orbitals (denoted triple- and quadruple-ζ) are used in the relativistic regime. Our method does not require the use of pseudopotentials and provides access to all electronic states within the same framework. Our study paves the way for direct studies of material properties, such as the parameters in spin Hamiltonians, that depend heavily on the electron density near atomic nuclei where relativistic and SOC effects are the strongest.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A high-order finite difference method for moving immersed domain boundaries and material interfaces

Here, we present a high-order sharp treatment of immersed moving domain boundaries and material interfaces, and apply it to the advection-diffusion equation in two and three dimensions. The spatial discretization combines dimension-split finite difference schemes with an immersed boundary treatment based on a weighted least-squares reconstruction of the solution, providing stable discretizations with up to sixth order accuracy for diffusion terms and third order accuracy for advection terms. The temporal discretization relies on a novel strategy for maintaining high-order temporal accuracy in problems with moving boundaries that minimizes implementation complexity and allows arbitrary explicit or diagonally-implicit Runge-Kutta schemes. The approach is broadly compatible with popular PDE-specialized Runge-Kutta time integrators, including low-storage, strong stability preserving, and diagonally implicit schemes. Through numerical experiments we demonstrate that the full discretization maintains high-order spatial and temporal accuracy in the presence of complex 3D geometries and for a range of boundary conditions, including Dirichlet, Neumann, and flux conditions with large jumps in coefficients.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Investigations on the multiple-sector hard-copper X-band accelerating structures

The development of advanced, high gradient accelerating structures is one of the leading activity of the particle accelerator community. In the technological research of new construction methods for these devices, high-power testing is a critical step for the verification of their viability. Recent experiments showed that accelerating cavities made out of hard copper, fabricated without high-temperature processes, can achieve better performance as compared with soft copper ones. Recently, we have built cavities using Tungsten Inert Gas welding and the high-power experiments confirmed that this joining process is a robust and low-cost alternative to brazing or diffusion bonding. This is a good solution for high-gradient operation, with a gradient of about 150 MV/m in X-band, at a breakdown rate of 10 -3 /pulse/meter using a shaped RF pulse with a 150 ns flat part. We continue the design, construction and high power tests of three-cells standing-wave X-band accelerating cavities fabricated with a split-open geometry, made of hard copper and vacuum-sealed by welding. Our aim is the fabrication of accelerating structures made out of hard copper alloys by using innovative cost-effective technologies. Moreover, our method of multiple-sector structures opens the way to new technological approaches and design methods, for example providing the possibility of placing dampers of the parasitic higher-order modes. Such structures with dampers allow for acceleration of multi-bunch beams, both in standing and in travelling-wave accelerating modality. Here, this paper describes the design of two cavities made of multiple sectors that were fabricated at the Italian Institute for Nuclear Physics in Frascati, Italy (INFN-LNF), assembled by means of clamping and then joined with the Tungsten Inert Gas welding in order to preserve the hardness of the metal. These are three-cell standing-wave X-band accelerating cavities, to be operated in the π-mode. We report the low-power RF tests of the first one, made of two halves, and the RF design of the second one, made of four quadrants. Both structures were built for high-gradient tests at SLAC National Accelerator Laboratory.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Photon energy-resolved velocity map imaging from spectral domain ghost imaging

Abstract We present an approach that combines photon spectrum correlation analysis with the reconstruction of three-dimensional momentum distribution from velocity map images in an efficient, single-step procedure. We demonstrate its efficacy with the results from the photoionization of the 2 p -shell of argon using the Free-electron LASer in Hamburg free-electron laser (FEL). Distinct spectral features due to the spin-orbit splitting of Ar + ( 2 p − 1 ) are resolved, despite the large average bandwidth of the ionizing pulses from the FEL. This demonstrates a clear advantage over the conventional analysis method, and it will be broadly beneficial for velocity map imaging experiments with FEL sources. The retrieved linewidth of the binding energy spectrum approaches the resolution limitation prescribed by the spectrometers used to collect the data. Our approach presents a path to extend spectral-domain ghost imaging to the case where the photoproduct observable is high-dimensional.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Deformation and Collectivity in Doubly Magic 208 Pb

Lead-208 is the heaviest known doubly magic nucleus and its structure is therefore of special interest. Despite this magicity, which acts to provide a strong restorative force toward sphericity, it is known to exhibit both strong octupole correlations and some of the strongest quadrupole collectivity observed in doubly magic systems. In this Letter, we employ state-of-the-art experimental equipment to conclusively demonstrate, through four Coulomb-excitation measurements, the presence of a large, negative, spectroscopic quadrupole moment for both the vibrational octupole 3$^{−}_{1}$ and quadrupole 2$^{+}_{1}$ state, indicative of a preference for prolate deformation of the states. The observed quadrupole moment is discussed in the context of the expected splitting of the 3 − ⊗3 − two-phonon states, due to the coupling of the quadrupole and octupole motion. These results are compared with theoretical values from three different methods, which are unable to reproduce both the sign and magnitude of this deformation. Thus, in spite of its well-studied nature, 208 Pb remains a puzzle for our understanding of nuclear structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Sequential formulation of all‐way coupled finite strain thermoporomechanics for largely deformable gas hydrate deposits

We develop a numerically stable sequential formulation of thermoporomechanics for largely deformable gas hydrate deposits, extended from the fixed stress split of infinitesimal transformation. Constitutive equations are based on the total Lagrangian approach for both flow and geomechanics, including dynamic full tensor permeability and thermal conductivity updated from the deformation gradient. For space discretization, we take the cell-centered finite volume and node-based finite element method for flow and geomechanics, respectively. Then, we propose a sequential implicit method for all-way coupled thermoporomechanics, where the nonisothermal multiphase flow problem of gas hydrates is solved implicitly first and then the geomechanics problem is solved implicitly at the next step. During solution of the flow problem, we fix the rate of first Pioal total stress for numerical stability as well as apply porosity correction and entropy correction to account for geomechanical effects. We test numerical examples where flow and geomechanics parameters are based on deep oceanic gas hydrate deposits. When applying depressurization, even though the results between the infinitesimal transformation and finite strain geomechanics are similar in the early stages due to small deformation, we find differences between them in the late times as deformation becomes large. Accordingly, permeability and thermal conductivity tensors become nonisotropic full tensors although they are initially isotropic. Furthermore, we identify numerical stability of the developed sequential method from the test cases that exhibit the highly complex coupled gas hydrate systems with large deformation. Thus, the proposed sequential formulation can be applied in largely deformable gas hydrate systems.

42 ENGINEERING↗
DOE OSTI2026

Mixed-precision numerics in scientific applications: survey and perspectives

The explosive demand for artificial intelligence (AI) workloads has led to a significant increase in silicon area dedicated to lower-precision computations on recent high-performance computing hardware designs. However, mixed-precision capabilities, which can achieve performance improvements of up to 8x compared to double-precision in extreme compute-intensive workloads, remain largely untapped in most scientific applications. A growing number of efforts have shown that mixed-precision algorithmic innovations can deliver superior performance without sacrificing accuracy. These developments should prompt computational scientists to seriously consider whether their scientific modeling and simulation applications could benefit from the acceleration offered by new hardware and mixed-precision algorithms. In this survey, we (1) review progress across diverse scientific domains—fluid dynamics, weather and climate, quantum chemistry, and computational genomics—that have begun adopting mixed-precision strategies; (2) examine state-of-the-art algorithmic techniques such as iterative refinement, splitting and emulation schemes, and adaptive precision solvers; (3) assess their implications for accuracy, performance, and resource utilization; and (4) survey the emerging software ecosystem that enables mixed-precision methods at scale. We conclude with perspectives and recommendations on cross-cutting opportunities, domain-specific challenges, and the role of co-design between application scientists, numerical analysts, and computer scientists. Collectively, this survey underscores that mixed-precision numerics can reshape computational science by aligning algorithms with the evolving landscape of hardware capabilities.

Graphics processing units↗
DOE OSTI2022

Main and Satellite Features in the Ni 2p XPS of NiO

The origin and assignment of the complex main and satellite XPS features of the cations in ionic compounds has been the subject of extensive theoretical studies using different methods. There is agreement that within a molecular orbital model one needs to take into account different types of configurations. Specifically, those where a core electron is removed but no other configuration changes are made and those where in addition to ionization there are also shake or charge transfer changes to the ionic configuration. However, there are strong disagreements about the assignment of XPS features to these configurations. The present work is directed toward resolving the origin of main and satellite features for the Ni 2p XPS of NiO based on ab initio molecular orbital wavefunctions for a cluster model of NiO. A major problem in earlier ab initio XPS studies of ionic compounds has been the use of a common set of orbitals that was not able to properly describe all the ionic configurations that contribute to the full XPS spectra. This is resolved in the present work by using orbitals that are optimized for averages of the occupations of the different configurations that contribute to the XPS. The approach of using State Averaged orbitals is validated through comparisons between different averages and through use of higher order excitations in the wavefunctions for the ionic states. It represents a major extension of our earlier work on the main and satellite features of the Fe 2p XPS of Fe 2 O 3 and proves the reliability and the generality of the assignments of the character and origin of the different features of the XPS obtained with orbitals optimized for State Averages. These molecular orbital methods permit the characterization of the ionic states in terms of the importance of shake excitations and of the coupling of ionization of 2p 1/2 and 2p 3/2 spin-orbit split sub shells. Further, the work lays the foundation for definitive assignments of the character of main and satellite XPS features and points to their origin in the electronic structure of the material.

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

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗
DOE OSTI2022

Accelerated Discovery of Solar Thermochemical Hydrogen Production Materials via High-Throughput Computational and Experimental Methods

In this project, combinatorial synthesis and testing methods were combined with high-throughput materials theory calculations to greatly accelerate the discovery of thermodynamically suitable candidates for green hydrogen production via a two-stage solar thermochemical water splitting (STCH) process. Over the course of the project, more than 8000 quinary and higher oxide compositions were computationally screened for STCH viability, and detailed stability calculations were performed for more than 30 of the most promising identified compositional archetypes. As a result, three new STCH capable compositional families were discovered and experimentally verified. The first, Ce x Sr 2-x MnO 4 (CSM), represents the first known Ruddlesden-Popper compound to show STCH activity, and thus demonstrates that perovskite-related structures may hold promise for this application. The second family, Sr 1-x Ce x MnO 3 (SCM), is the simple perovskite sister-analog to CSM. Sr 0.7 Ce 0.3 MnO 3 (SCM30), a member of this compositional family, was found to produce the highest hydrogen yields of any compound tested in this project, exceeding the end of project milestone target of > 150 μmol H 2 /gram oxide at a reduction temperature of 1350 °C, although only at steam-to-hydrogen ratios greater than 1000:1. Finally, we proved that a third novel Sr-and Mn-containing family, Sr 1-x Ca x Ti 1-y Mn y O 3 (SCTM), which was identified by Materials Project tools, also splits water. The behavior of the SCTM system was found to be similar to the previously discovered Sr 1-x La x Al 1-y Mn y O 3 (SLMA) family, albeit with lower H 2 yields. Across the three thrusts of the project (computational, combinatorial, and bulk testing), five journal articles were published. As part of Program End Analysis and Data Dissemination, relevant data used for the publications was uploaded to the HydroGEN Data Hub for public access, and in certain cases, results were added to public materials databases.

08 HYDROGEN↗
DOE OSTI2025

HydroGEN STCH node support: FY25 Q1 quarterly report – LLNL's input

For reducing the cost of H 2 production by enhancing the efficiency and cyclability of water splitting reactions, we aim to understand the relationship between microstructure and surface reaction performance on perovskite oxides. In this quarter, mesoscale simulations were carried with various 3D porous microstructures to characterize the microstructure-dependent H 2 production. During the water splitting process, incoming gas including steam (H 2 O) passes through porous regions, which can induce local pressure variations according to the Bernoulli’s principle. Accordingly, we developed the mesoscale simulation method to characterize the reaction kinetics with inhomogeneous local pressures within porous microstructures.

08 HYDROGEN↗
DOE OSTI2020

Millimeter-wave Point Sources from the 2500 Square Degree SPT-SZ Survey: Catalog and Population Statistics

In this work, we present a catalog of emissive point sources detected in the SPT-SZ survey, a contiguous 2530- square-degree area surveyed with the South Pole Telescope (SPT) from 2008 - 2011 in three bands centered at 95, 150, and 220GHz. The catalog contains 4845 sources measured at a significance of 4.5σ or greater in at least one band, corresponding to detections above approximately 9.8, 5.8, and 20.4 mJy in 95, 150, and 220 GHz, respectively. Spectral behavior in the SPT bands is used for source classification into two populations based on the underlying physical mechanisms of compact, emissive sources that are bright at millimeter wavelengths: synchrotron radiation from active galactic nuclei and thermal emission from dust. The latter population includes a component of high-redshift sources often referred to as submillimeter galaxies (SMGs). In the relatively bright flux ranges probed by the survey, these sources are expected to be magnified by strong gravitational lensing. The survey also contains sources consistent with protoclusters, groups of dusty galaxies at high redshift undergoing collapse. We cross-match the SPT-SZ catalog with external catalogs at radio, infrared, and X-ray wavelengths and identify available redshift information. The catalog splits into 3980 synchrotron-dominated and 865 dust-dominated sources and we determine a list of 506 SMGs. Ten sources in the catalog are identified as stars. We calculate number counts for the full catalog, and synchrotron and dusty components, using a bootstrap method and compare our measured counts with models. This paper represents the third and final catalog of point sources in the SPT-SZ survey.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Split representation of adaptively compressed polarizability operator

The polarizability operator plays a central role in density functional perturbation theory and other perturbative treatment of first principle electronic structure theories. The cost of computing the polarizability operator generally scales as O(Ne4) where Ne is the number of electrons in the system. The recently developed adaptively compressed polarizability operator (ACP) formulation [L. Lin, Z. Xu and L. Ying, Multiscale Model. Simul. 2017] reduces such complexity to O(Ne 3 ) in the context of phonon calculations with a large basis set for the first time, and demonstrates its effectiveness for model problems. In this paper, we improve the performance of the ACP formulation by splitting the polarizability into a near singular component that is statically compressed, and a smooth component that is adaptively compressed. The new split representation maintains the O(Ne 3 ) complexity, and accelerates nearly all components of the ACP formulation, including Chebyshev interpolation of energy levels, iterative solution of Sternheimer equations, and convergence of the Dyson equations. For simulation of real materials, we discuss how to incorporate nonlocal pseudopotentials and finite temperature effects. In this work, we demonstrate the effectiveness of our method using one-dimensional model problem in insulating and metallic regimes, as well as its accuracy for real molecules and solids.

97 MATHEMATICS AND COMPUTING↗