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 1,099 records · Page 61

Less-Altered CI Lithology in the Kaidun Meteorite Breccia

One of the most interesting initial results of the characterization of the asteroid Ryugu regolith material returned by the Hayabusa2 spacecraft is the presence of a significant quantity of less-aqueously altered CI material, tentatively called “CI2”. This lithology is characterized by an approximately ten-fold increase in the amount of an hydrous ferromagnesian silicates over what is found n CI1 chondrites, with these being principally olivine, low- and high-calcium pyroxene. Remnant CAI phases are found as well as probable chondrule fragments. Magnetite, pyrrhotite and hydroxyapatite are abundant. The main carbonate is calcite. Matrix is mainly a mixture of saponite and serpentine phases. Assuming that this lithology should be a better representative of the more “primordial” mineralogy of early nebular solids than more processed CI1 material, we have been re-examining CI lithologies found as separate meteorites as well as clasts in other meteorites, seeking similar CI2 material. The Kaidun breccia is notorious for including almost every type of meteorite class except CR2 (which is ironic considering that Kaidun is still officially classified as CR2). We report here the discovery of a probable CI2 lithology in one carbonaceous chondrite-rich Kaidun sample.

Michael Zolensky↗

Precision Oxygen Isotope Measurements of Two C-Rich Hydrated Interplanetary Dust Particles

Introduction: Chondritic-smooth IDPs (Interplanetary Dust Particles) are low porosity objects whose mineralogy is dominated by aqueous alteration products such as Mg-rich phyllosilicates (smectite and serpentine group) and Mg-Fe carbonate minerals. Their hydrated mineralogy combined with low atmospheric entry velocities have been used to infer an origin largely from asteroidal sources. Spectroscopic studies show that the types and abundance of organic matter in CS IDPs is similar to that in CP IDPs. Although CS IDPs show broad similarities to primitive carbonaceous chondrites, only a few particles have been directly linked to specific meteorite groups such as CM and CI chondrites based on the presence of diagnostic minerals. Many CS IDPs however, have carbon contents that greatly exceed that of known meteorite groups suggesting that they either may derive from comets or represent samples of more primitive parent bodies than do meteorites. It is now recognized that many large, dark primitive asteroids in the outer main belt, as well as some trans-Neptunian objects, show spectroscopic evidence for aqueous alteration products on their surfaces. Some CS IDPs exhibit large bulk D enrichments similar to those observed in the cometary CP IDPs. While hydrated minerals in comets have not been unambiguously identified to date, the presence of the smectite group mineral nontronite has been inferred from infrared spectra obtained from the ejecta from comet 9P/Tempel 1 during the Deep Impact mission. Recent observations of low temperature sulfide minerals in Stardust mission samples suggest that limited aqueous activity occurred on comet Wild-2. All of these observations, taken together, suggest that the high-carbon hydrated IDPs are abundant and important samples of primitive solar system objects not represented in meteorite collections. Oxygen isotopic compositions of chondrites reflect mixing between a 16O-rich reservoir and a 17O,18O-rich reservoir produced via mass-independent fractionation. The composition of the 16O-rich reservoir is well constrained but material representing the 17O,18O-rich end-member is rare. Self-shielding models predict that cometary water, presumed to represent this reservoir, should be enriched in 17O and 18O by greater than 200 per mille. The high-carbon hydrated IDPs may be among the best materials available to search for preserved "cometary" H2O signatures. In order to better understand the origin and evolution of these particles, we have obtained 10 hydrated interplanetary dust particles for coordinated mineralogical, isotopic and organic analyses. We have previously reported the results of mineralogical and O isotopic measurements of two hydrated IDPs; here we present results of O isotopic measurements of three additional IDPs. Samples and Methods: Three interplanetary dust particles (L2079C35, L2083D46 and L2083E46) were embedded in S and partially ultramicrotomed into approximately 70 nanometer sections for analysis via transmission electron microscopy (TEM). The remainders of the unsliced particles were removed from S and pressed into high purity Au foil that was cleaned with HF acid and annealed at 800 degrees Centigrade. The pressed IDPs were analyzed via electron microprobe analysis (EPMA) for quantitative bulk chemical analysis. After EPMA analysis, the IDPs were subjected to precision O isotope analysis with the UCLA Cameca IMS-1270 ion probe. A 20 kiloelectronvolt, 0.5 nanoangstrom Cs+ primary beam of approximately 15 micrometers diameter was used for each measurement. Small particles of San Carlos olivine and Burma spinel were pressed into the Au foil for use as standards to correct for instrumental mass fractionation. The detection system was configured for multicollection, with 16O measured on a Faraday cup, and 17O and 18O measured on electron multipliers (EMs). Individual analyses consisted of 15 cycles of 10 seconds per cycle. Additionally, two microtome thin sections were measured for H isotopic compositions with the JSC NanoSIMS 50L ion probe. An 8 picoangstrom, 16 kiloelectronvolt Cs plus primary beam was used. Measurements consisted of H, C, 12C, 16O, and 18O collected with EMs in multicollection. Terrestrial biotite and kerogen were used for isotopic standards. A significant challenge in O isotope measurement of hydrated minerals is the interference from 16OH at mass 17O. We ensured that the 17O and 16OH peaks were fully resolved by using a mass resolution of greater than 7000 and by careful analyses of San Carlos olivine, Burma spinel and chlorite hydrated mineral standards. The hydride was further suppressed with a cold finger attached to an LN2 dewar to trap volatiles in the sample chamber. All sputtered ions were counted (i.e. presputtering was not used); after applying background, yield and deadtime corrections, we performed a change-point analysis on our data via R in order to determine when the sample reached sputtering equilibrium; data points collected prior to the change point were excluded. Change-point analysis was also used to determine whether the IDP had completely sputtered. Results: Mineralogy. IDPs C35 and E46 exhibited hydrated mineralogies, Fe-Ni sulfide grains, nanoglobules and occasional enstatite grains distributed throughout a fine-grained Mg-Fe saponite matrix. C35 also contained breunnerite (Mg,Fe)CO3; solar flare tracks were observed in enstatite, indicating minimal atmospheric entry heating. The mineralogy of D46 is dominated by a large FeS grain with a minor component of adhering silicate material. D46 was strongly heated during atmorpheric entry as evidenced by a well-developed magnetite rim. EPMA analyses show that both C35 and E46 have high carbon contents of 20 weight percentage (approximately 6X CI). D46 contains approximately 6 wt.weight percentage C. A significant amount of the carbon is present as carbon nanoglobules. Results: Hydrogen isotopes: Although the bulk delta D values of both sections of L2079C35 were within error of SMOW (minus 33 plus or minus 19 per mille, 1 plus or minus 14 per mille 1 sigma), several delta D-rich hotspots were also identified, reaching 2000 per mille. As shown in Fig. 1, these hotspots are clearly associated with discrete carbonaceous inclusions that are akin to nanoglobules found in many meteorites and other IDPs. Oxygen Isotopes. Results of the oxygen isotope measurements are shown in Figure 1. The oxygen isotope composition for L2079C35 was delta 18O equals plus11.6 plus or minus 1.9per mille, delta 17O equals plus 7.9 plus or minus 1.9per mille (2 standard errors). The oxygen isotope composition for L2083D46 was delta 18O equals minus 8.1 plus or minus 1.9 per mille, delta 17O equals minus 6.4 plus or minus 3.1 per mille (2 standard errors). The oxygen isotope composition for L2083E46 was delta 18O equals plus 12.0 plus or minus 1.9 per mille, delta 17O equals plus 9.2 plus or minus 2.0 per mille (2 standard errors). Discussion: Despite mineralogical similarities to highly aqueously altered carbonaceous chondrites, the hydrated IDPs we analyzed have oxygen isotopic compositions that are distinct from matrix materials in the CI, CM, and CR chondrites. The IDPs plot along the Young-Russell line, with delta 17O values for C35 and E46 suggestive of interaction with a 16O-poor reservoir. However, we have thus far not observed evidence of extreme 16O-poor reservoirs expected from self-shielding models and observed in Acfer 094 simplectite. The high carbon contents of the IDPs also set them apart from known meteoritic samples. The lack of atmospheric entry heating effects are consistent with low encounter velocities and suggest either an asteroidal source, or a low inclination, low eccentricity cometary origin. Conclusions: The unusual oxygen isotopic compositions, high carbon contents, and the abundance of Drich nanoglobules, together, suggest that the high-carbon, hydrated IDPs are derived from a primitive source that is not yet represented in meteorite collections.

Snead, C. J.↗

Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

Chemistry↗

Continuous Counter‐Current Microfluidic Liquid–Liquid Extraction Achieved Using a Pair of Wettable Screen Meshes

Continuous counter‐current microfluidic liquid–liquid extraction performs separations by flowing immiscible liquids in opposing directions within a single flow channel. In principle, this flow arrangement enables a large number of theoretical separation units in a small footprint, without using interstage valving, pumping, and phase separation. Despite its potential for excellent separation performance, this microfluidic scheme rarely appears in literature due to the requirement for capillary forces to be greater than hydrodynamic forces for stable flow. We present a novel microfluidic device and flow approaches that overcome this force‐balance challenge, enabling stable, long‐duration continuous counter‐current flow. Additionally, we cover a suite of methodologies for quantifying the performance of the microfluidic device, revealing the number of theoretical equilibrium stages achieved. The enabling technologies include a woven mesh screen‐based microfluidic device architecture that is easily fabricated outside of a clean room, surface functionalization strategies to promote conjugate (organic/aqueous) wettability, flow approaches to eliminate bubbles and carryover, and computer‐aided flow automation with optical measurement of extraction performance. The reported experiments lasted for over 36 h, terminated only at experiment conclusion, where the device still exhibited good performance. Automated Raman spectroscopy was used for solute quantitation of the ternary system tert‐butanol in a toluene/water matrix, a ternary system that was specifically chosen to analyze the device's performance with a small solute partition ratio and to enable in‐line Raman measurements of solute concentrations in both phases. The microfluidic device possessed a 55 mm contact length and a 38.5 µL internal volume. During counter‐current flow, we observed approximately 37 equilibrium stages (37 ± 13) based on a best‐fit of the solute fraction remaining in the aqueous phase using a Kremser Group Method analysis.

36 MATERIALS SCIENCE↗

Using computational singular perturbation as a diagnostic tool in ODE and DAE systems: a case study in heterogeneous catalysis

We have extended the computational singular perturbation (CSP) method to differential algebraic equation (DAE) systems and demonstrated its application in a heterogeneous-catalysis problem. The extended method obtains the CSP basis vectors for DAEs from a reduced Jacobian matrix that takes the algebraic constraints into account. Here we use a canonical problem in heterogeneous catalysis, the transient continuous stirred tank reactor (T-CSTR), for illustration. The T-CSTR problem is modelled fundamentally as an ordinary differential equation (ODE) system, but it can be transformed to a DAE system if one approximates typically fast surface processes using algebraic constraints for the surface species. We demonstrate the application of CSP analysis for both ODE and DAE constructions of a T-CSTR problem, illustrating the dynamical response of the system in each case. We also highlight the utility of the analysis in commenting on the quality of any particular DAE approximation built using the quasi-steady state approximation (QSSA), relative to the ODE reference case.

97 MATHEMATICS AND COMPUTING↗

Quasiparticle conductance in spin valve Josephson structures

We study the quasiparticle current in clean ferromagnetic Josephson structures of the form $S_1/F_1/N/F_2/S_2$, where $\textit{S, F,}$ and $\textit{N}$ denote superconducting, ferromagnetic, or normal layers, respectively. Our focus is on the structure of the conductance $\textit{G}$ as a function of bias $\textit{V$, emphasizing the subgap region. We use a fully self-consistent numerical method, coupled to a transfer matrix procedure to extract $\textit{G(V)}$. We choose material parameters appropriate to experimentally realized Co Cu Nb structures. In this work, we find a resonance peak structure as a function of the intermediate layer thickness and of the misalignment angle $\phi$ between $\textit{F}$ layers. To understand this resonance structure, we develop an approximate analytic method. For experimentally relevant thicknesses, the conductance has multiple subgap peaks, which oscillate in position between low and critical bias positions. These oscillations occur in both $\phi$ and the layer thicknesses. We compare our results with those obtained for the spin valve structures ($F_1/N/F_2/S_2$) and discuss the implications of our results for the fabrication of spin Josephson devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

powersqueeze

powersqueeze (psqz) is a truncated power iteration library intended for high-performance computing platforms. psqz efficiently produces low-dimensional, linear measurements of graph matrix spectra by combining classical power iteration with sparse Johnson-Lindenstrauss transforms. psqz is intended to produce high-quality, fast, data-oblivious low-dimensional representations of high-dimensional sparse data such as graphs and term-document matrices. psqz is intended to replace similar workflows that depend on directly approximating a truncated eigendecomposition (e.g., the first step of spectral clustering), which is a much more expensive operation.

Priest, BenjaminW [Lawrence Livermore National Lab↗

A study of strength evaluations of EDNi/EDCu/NARloy-Z bonded joints

Dissimilar material interfaces can be found in many materials and structural bonds such as composite materials, welded parts, inclusion in matrix, bond between metallic and ceramic materials, and so forth. One of such structural bonds can be seen in the main combustion chamber (MCC) of the Space Shuttle Main Engine (SSME). In this study, from a practical sense, the primary concern is to understand the systems response of EDNi/EDCu/NARloy-Z bonded joints using stress values approximated by the finite element method to determine an influence of the variation of structural bond parameters on the bonded joints, and consequently main combustion chamber of the Space Shuttle Main Engine. The results presented in this study could be an appropriate indicator for a good bond of EDNi/DECu/NARloy-Z layers with the desired thickness of copper deposition in the SSME MCC manufacturing process. Furthermore, the results from this study appear to be applicable to any bonded joints that can be characterized by the parameters and assumptions used in this analysis.

Min, J. B.↗

Coagulation of grains in static and collapsing protostellar clouds

The wavelength dependence of extinction in the diffuse interstellar medium implies that it is produced by particles of dominant size of approximately 10(exp -5) cm. There is some indication that in the cores of dense molecular clouds, sub-micron grains can coagulate to form larger particles; this process is probably driven by turbulence. The most primitive meteorites (carbonaceous chondrites) are composed of particles with a bimodal size distribution with peaks near 1 micron (matrix) and 1 mm (chondrules). Models for chondrule formation that involve processing of presolar material by chemical reactions or through an accretion shock during infall assume that aggregates of the requisite mass could form before or during collapse. The effectiveness of coagulation during collapse has been disputed; it appears to depend on specific assumptions. The first results of detailed numerical modeling of spatial and temporal variations of particle sizes in presolar clouds, both static and collapsing, is reported in this article.

Weidenschilling, S. J.↗

The Foreign Clast Populations of Anomalous Polymict Urelite Almahata Sitta (Asteroid 2008 TC(sub3) and Typical Polymict Ureilites: Implications for Asteroid-Meteorite Connections

Almahata Sitta (AhS) is the first meteorite to originate from an asteroid (2008 TC3) that had been studied in space before it hit Earth [1,2]. It is also unique because the fallen fragments comprise a variety of types: approximately 69% ureilites (achondrites) and 31% chondrites [3]. Two models have been proposed for the origin 2008 TC3: 1) an accretionary model [3,4]; or 2) a regolith model [5,6]. Typical polymict ureilites are interpreted to represent regolith, and contain a few % foreign clasts [7,8]. The most common are dark (CC matrix-like) clasts similar to those in many meteoritic breccias [9]. A variety of other chondrites, as well as achondrites (angrites), have also been reported [7,9,10]. We have been working to determine the full diversity of these clasts [10-13] for comparison with AhS. We discuss implications for mixing of materials in the early solar system and the origin of 2008 TC3.

Goodrich, C. A.↗

Vapor Grown Carbon Fiber/Phenolic Matrix Composites for Rocket Nozzles and Heat Shields

The ablation and mechanical and thermal properties of vapor grown carbon fiber (VGCF)/phenolic resin composites were evaluated to determine the potential of using this material in solid rocket motor nozzles. Composite specimens with varying VGCF loading (30%-50% wt) including one sample with ex-rayon carbon fiber plies were prepared and exposed to a plasma torch for 20 s with a heat flux of 16.5 MW/sq m at approximately 1650 C. Low erosion rates and little char formation were observed, confirming that these materials were promising for rocket motor nozzle materials. When fiber loadings increased, mechanical properties and ablative properties improved. The VGCF composites had low thermal conductivities (approximately 0.56 W/m-C) indicating they were good insulating materials. If a 65% fiber loading in VGCF composite can be achieved, then ablative properties are projected to be comparable to or better than the composite material currently used on the Space Shuttle Reusable Solid Rocket Motor (RSRM).

Patton, R. D.↗

Theoretical Investigations of Anharmonic Effects and Phonon Transport in the Cubic Phase of Crystalline Perovskite CsPbCl 3

The role of anharmonic effects on lattice dynamics and thermal transport has been investigated in the cubic phase of the CsPbCl 3 perovskite. Limitations of the harmonic approximation, which lead to phonon instabilities in the Brillouin zone, were addressed based on self-consistent phonon theory in combination with first-principles calculations and by incorporation of frequency renormalization effects from bubble and loop diagrams. This theoretical approach demonstrates significant improvements compared to results obtained using a self-consistent phonon dynamic matrix. Both the cubic-to-tetragonal phase transition temperature as well as the predicted lattice thermal conductivity values show good agreement with selected sets of available experimental data. In conclusion, the accuracy of the predicted thermal conductivity data is also tested against systems significantly larger than those allowed by quantum calculations, by using molecular dynamic simulations with machine-learning force fields.

Lattices↗

Accurate localization of Kosterlitz-Thouless-type quantum phase transitions for one-dimensional spinless fermions

We investigate the charge-density wave (CDW) transition for one-dimensional spinless fermions at half band filling with nearest-neighbor electron transfer amplitude t and interaction V. The model is equivalent to the anisotropic XXZ Heisenberg model for which the Bethe Ansatz provides an exact solution. For V>V c =2t, the CDW order parameter and the single-particle gap are finite but exponentially small, as is characteristic for a Kosterlitz-Thouless transition. It is notoriously difficult to locate such infinite-order phase transitions in the phase diagram using approximate analytical and numerical approaches. Second-order Hartree-Fock theory is qualitatively applicable for all interaction strengths, and predicts the CDW transition to occur at V$^{(2)}_{c,2}$≈1.5t. Second-order Hartree Fock theory is almost variational because the density of quasiparticle excitations is small. We apply the density-matrix renormalization group (DMRG) for periodic boundary conditions for system sizes up to 514 sites, which permits a reliable extrapolation of all physical quantities to the thermodynamic limit, apart from the critical region. We investigate the ground-state energy, the gap, the order parameter, the momentum distribution, the quasiparticle density, and the density-density correlation function to locate V c from the DMRG data. In conclusion, tracing the breakdown of the Luttinger liquid and the peak in the quasiparticle density at the band edge permits us to reproduce V c with an accuracy of one percent.

1-dimensional spin chains↗

Solution of free-boundary problems using finite-element/Newton methods and locally refined grids - Application to analysis of solidification microstructure

A new method is presented for the solution of free-boundary problems using Lagrangian finite element approximations defined on locally refined grids. The formulation allows for direct transition from coarse to fine grids without introducing non-conforming basis functions. The calculation of elemental stiffness matrices and residual vectors are unaffected by changes in the refinement level, which are accounted for in the loading of elemental data to the global stiffness matrix and residual vector. This technique for local mesh refinement is combined with recently developed mapping methods and Newton's method to form an efficient algorithm for the solution of free-boundary problems, as demonstrated here by sample calculations of cellular interfacial microstructure during directional solidification of a binary alloy.

Tsiveriotis, K.↗

On the Use of the Log-Normal Particle Size Distribution to Characterize Global Rain

Although most parameterizations of the drop size distributions (DSD) use the gamma function, there are several advantages to the log-normal form, particularly if we want to characterize the large scale space-time variability of the DSD and rain rate. The advantages of the distribution are twofold: the logarithm of any moment can be expressed as a linear combination of the individual parameters of the distribution; the parameters of the distribution are approximately normally distributed. Since all radar and rainfall-related parameters can be written approximately as a moment of the DSD, the first property allows us to express the logarithm of any radar/rainfall variable as a linear combination of the individual DSD parameters. Another consequence is that any power law relationship between rain rate, reflectivity factor, specific attenuation or water content can be expressed in terms of the covariance matrix of the DSD parameters. The joint-normal property of the DSD parameters has applications to the description of the space-time variation of rainfall in the sense that any radar-rainfall quantity can be specified by the covariance matrix associated with the DSD parameters at two arbitrary space-time points. As such, the parameterization provides a means by which we can use the spaceborne radar-derived DSD parameters to specify in part the covariance matrices globally. However, since satellite observations have coarse temporal sampling, the specification of the temporal covariance must be derived from ancillary measurements and models. Work is presently underway to determine whether the use of instantaneous rain rate data from the TRMM Precipitation Radar can provide good estimates of the spatial correlation in rain rate from data collected in 5(sup 0)x 5(sup 0) x 1 month space-time boxes. To characterize the temporal characteristics of the DSD parameters, disdrometer data are being used from the Wallops Flight Facility site where as many as 4 disdrometers have been used to acquire data over a 2 km path. These data should help quantify the temporal form of the covariance matrix at this site.

Meneghini, Robert↗

Diversity in C-Xanes Spectra Obtained from Carbonaceous Solid Inclusions from Monahans Halite

Monahans meteorite (H5) contains fluid inclusion- bearing halite (NaCl) crystals [1]. Microthermometry and Raman spectroscopy showed that the fluid in the inclusions is an aqueous brine and they were trapped near 25degC [1]. Their continued presence in the halite grains requires that their incorporation into the H chondrite asteroid was post metamorphism [2]. Abundant solid inclusions are also present in the halites. The solid inclusions include abundant and widely variable organics [2]. Analyses by Raman microprobe, SEM/EDX, synchrotron X-ray diffraction and TEM reveal that these grains include macromolecular carbon similar in structure to CV3 chondrite matrix carbon, aliphatic carbon compounds, olivine (Fo99−59), high- and low-Ca pyroxene, feldspars, magnetite, sulfides, lepidocrocite, carbonates, diamond, apatite and possibly the zeolite phillipsite [3]. Here we report organic analyses of these carbonaceous residues in Monahans halite using C-, N-, and O- X-ray absorption near edge structure (XANES). Samples and Methods: Approximately 100 nm-thick sections were extracted with a focused ion beam (FIB) at JSC from solid inclusions from Monahans halite. The sections were analyzed using the scanning transmission X-ray microscope (STXM) on beamline 5.3.2.2 at the Advanced Light Source, Lawrence Berkeley National Laboratory for XANES spectroscopy. Results and Discussion: C-XANES spectra of the solid inclusions show micrometer-scale heterogeneity, indicating that the macromolecular carbon in the inclusions have complex chemical variations. C-XANES features include 284.7 eV assigned to aromatic C=C, 288.4-288.8 eV assigned to carboxyl, and 290.6 eV assigned to carbonate. The carbonyl features obtained by CXANES might have been caused by the FIB used in sample preparation. No specific N-XANES features are observed. The CXANES spectra obtained from several areas in the FIB sections include type 1&2 chondritic IOM like, type 3 chondritic IOM like, and none of the above. The natures of the macromolecular carbon in the solid inclusions observed by C-XANES are consistent with the previous studies showing that the carbonaceous solid inclusions have not originated from Monahans parent body [1-3], and have various origins, including various chondritic meteorite parent bodies as well as other unknown source(s).

Kebukawa, Y.↗

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗

S-matrix method for the numerical determination of bound states.

A rapid numerical technique for the determination of bound states of a partial-wave-projected Schroedinger equation is presented. First, one needs to integrate the equation only outwards as in the scattering case, and second, the number of trials necessary to determine the eigenenergy and the corresponding eigenfunction is considerably less than in the usual method. As a nontrivial example of the technique, bound states are calculated in the exchange approximation for the e-/He+ system and l equals 1 partial wave.

Bhatia, A. K.↗