Search NASA⌕ Search

SEARCH · Search NASA

Results for “force 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Evaluation of Gamma Ray Production - Quarterly Report FY22Q2

The LLNL team has performed the first systematic study of the effect of two-body weak currents in β decay rates across the entire mass table. In light and medium-mass nuclei, two-body weak currents from chiral effective field theory account for a significant portion of the phenomenological quenching of Gamow-Teller transition matrix elements. We examined the systematic effects of two-body axial currents on Gamow-Teller strength and β-decay rates in heavy nuclei within energy-density functional theory. We extended the charge-changing finiteamplitude pnFAM code to incorporate the contributions of two-body currents to the usual onebody linear response in the Gamow-Teller channel. We found that the two-body currents, as expected, usually quench both summed Gamow-Teller strength and decay rates, but by an amount that decreases as the neutron excess grows; see Fig. 1. In addition, they can enhance individual low-lying transitions, leading to decay rates that are quite different from those that an energy-independent quenching would produce, particularly in neutron-rich nuclei. This new physics capability is unique and provides the highest-fidelity description of β decay on the market. The underlying computational implementation can be reused to calculate γ-strength functions for finite-range nuclear potentials such as the Gogny force, which will allow us to benchmark some of our future developments against published results. The paper is in print in Phys. Rev. C.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Utilizing Novel Field and Data Exploration Methods to Explore Hot Moments in High-Frequency Soil Nitrous Oxide Emissions Data: Opportunities and Challenges

Soil nitrous oxide (N 2 O) emissions are an important driver of climate change and are a major mechanism of labile nitrogen (N) loss from terrestrial ecosystems. Evidence increasingly suggests that locations on the landscape that experience biogeochemical fluxes disproportionate to the surrounding matrix (hot spots) and time periods that show disproportionately high fluxes relative to the background (hot moments) strongly influence landscape-scale soil N 2 O emissions. However, substantial uncertainties remain regarding how to measure and model where and when these extreme soil N 2 O fluxes occur. High-frequency datasets of soil N 2 O fluxes are newly possible due to advancements in field-ready instrumentation that uses cavity ring-down spectroscopy (CRDS). Here, we outline the opportunities and challenges that are provided by the deployment of this field-based instrumentation and the collection of high-frequency soil N 2 O flux datasets. While there are substantial challenges associated with automated CRDS systems, there are also opportunities to utilize these near-continuous data to constrain our understanding of dynamics of the terrestrial N cycle across space and time. Finally, we propose future research directions exploring the influence of hot moments of N 2 O emissions on the N cycle, particularly considering the gaps surrounding how global change forces are likely to alter N dynamics in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Lowest-energy broad 𝛼-cluster resonances in 19 F

There is strong astrophysical interest in the structure of 19 F states near the α-decay threshold, as their properties are important for understanding the development of α clustering in the 20 Ne region. The emergence of clustered states, and more generally of states that couple strongly to reaction channels near their decay thresholds, is also a topic of current interest in theoretical nuclear physics. In this work, we determine the parameters of broad low-spin states in 19 F near the α-decay threshold and present a theoretical study of these states. The analysis is restricted to ℓ = 0 and 1 resonances in the α + 15 N system close to the α threshold in 19 F. Excitation functions for 15 N(α,α) elastic scattering were measured using the Thick Target Inverse Kinematics (TTIK) method, and these new data, together with older high–energy-resolution measurements, were analyzed within the R-matrix framework. The nuclear structure of 19 F was calculated using configuration-interaction methods with a recently developed effective interaction Hamiltonian. As a result, we identify a sequence of α-cluster resonances in 19 F and map the distribution of clustering strength, which is relevant for astrophysical reaction modeling and for the theoretical understanding of many-body dynamics and the emergence of clustering in loosely bound or unstable nuclei. Furthermore, the work advances theoretical insight into the origins of clustering and highlights open questions for future theoretical and experimental studies.

19F↗
DOE OSTI2024

Effect of parallel flow on resonant layer responses in high beta plasmas

Abstract Resonant layers in a tokamak respond to non-axisymmetric magnetic perturbations by amplifying the mode amplitude and balancing the plasma rotation through magnetic reconnection and force balance, respectively. This resonant response can be characterized by local layer parameters and especially by a single quantity in the linear regime, the so-called inner-layer Δ. The computation of Δ under two-fluid drift-MHD formalism has been progressed by reducing the order of the system in the phase space, where the shielding current is approximated as being only carried by electrons, a posteriori . In this study, we relax the approximation and compute Δ accounted for by the parallel flow associated with the ion shielding current. The posteriori is numerically verified in great agreement with the original SLAYER developed in a previous paper (J.-K. Park 2022 Phys. Plasmas 29 072506). Extending the resonant layer response theory to high β plasmas, our research findings answer two important questions: how the parallel flow influences the resonant layer response and why the parallel flow effect appears in high β plasmas. The complicated plasma compression in high β regime allows the parallel flow response to give rise to the ion shielding current, which not only shifts the zero-crossing condition of the ExB flow but also enhances the field penetration threshold. Technically, the Riccati matrix transformation method is adapted to handle the numerical stiffness due to the increased order of the system. The high fidelity of this numerical method makes use of further extension of the model to higher-order systems to take other physical phenomena into account. This work is envisaged to predict the resonant layer response under high β fusion reactor conditions.

Lee, Yeongsun (ORCID:000000034474416X)↗
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Isospin mixing and the cubic isobaric multiplet mass equation in the lowest T=2, A=32 quintet

The isobaric multiplet mass equation (IMME) is known to break down in the first T=2,A=32 isospin quintet. In this work we combine high-resolution experimental data with state-of-the-art shell-model calculations to investigate isospin mixing as a possible cause for this violation. The experimental data are used to validate isospin-mixing matrix elements calculated with newly developed shell-model Hamiltonians. Our analysis shows that isospin mixing with nonanalog T=1 states contributes to the IMME breakdown, making the requirement of an anomalous cubic term inevitable for the multiplet.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Enhancing membrane performance for CO 2 capture from flue gas with ultrahigh MW polyvinylamine

Membranes for post-combustion CO 2 capture are required to have a high CO 2 permeance due to the limited driving force. For facilitated transport membranes synthesized with polyvinylamine (PVAm), a defect-free selective layer of <200 nm is usually required to render sufficient permeance with high CO 2 /N 2 selectivity. In order to meet such a demand through the knife-coating process, the coating solution needs to have a high viscosity at a relatively low concentration to minimize its penetration into the substrate. The demand was met by synthesizing PVAm with an ultrahigh molecular weight (MW) via inverse emulsion polymerization (IEP) in this study. Compared to solution polymerization, IEP isolates the reaction in inverse micelles suspended in a continuous organic phase, which allows excellent dissipation of the heat generated by the reaction and reduces gel formation drastically. The polymerization parameters, including monomer concentration, initiator concentration, and reaction temperature, were investigated to obtain the optimal MW of PVAm for membrane performance. As compared to solution polymerization, IEP enhanced the MW of PVAm from 1.2 to 12.7 MDa, which minimized the penetration of the coating solution into the substrate and hence the extra mass transfer resistance. The effect of MW and degree of hydrolysis of PVAm on the membrane transport properties were studied. Furthermore, by employing PVAm with a MW of 12.7 MDa to strengthen the polymer matrix, the loading of piperazine glycinate in the membrane was increased up to 85 wt.%. In conclusion, the resultant membrane achieved a CO 2 permeance of 839 GPU and a CO 2 /N 2 selectivity of 161 at the typical flue gas temperature of 57°C.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Quasi-deuteron model at low renormalization group resolution

The quasi-deuteron model introduced by Levinger is used to explain cross sections for knocking out high-momentum protons in photoabsorption on nuclei. This is within a framework we characterize as exhibiting high renormalization group (RG) resolution. Assuming a one-body reaction operator, the nuclear wave function must include two-body short-range correlations (SRCs) with deuteronlike quantum numbers. In Phys. Rev. C 104, 034311 (2021), we showed that SRC physics can be naturally accounted for at low RG resolution. We describe the quasi-deuteron model at low RG resolution and determine the Levinger constant, which is proportional to the ratio of nuclear photoabsorption to that for photodisintegration of a deuteron. We extract the Levinger constant based on the ratio of momentum distributions at high relative momentum. We compute momentum distributions evolved under similarity RG (SRG) transformations where the SRC physics is shifted into the operator as a universal two-body term. The short-range nature of this operator motivates using local-density approximations with uncorrelated wave functions in evaluating nuclear matrix elements, which greatly simplifies the analysis. The operator must be consistently matched to the RG scale and scheme of the interaction for a reliable extraction. We apply SRG transformations to different nucleon-nucleon (NN) interactions and use the deuteron wave functions and Weinberg eigenvalues to determine approximate matching scales. We predict the Levinger constant for several NN interactions and a wide range of nuclei comparing to experimental extractions. The predictions at low RG resolution are in good agreement with experiment when starting with a hard NN interaction and the initial operator. Similar agreement is found using soft NN interactions when the additional two-body operator induced by evolution from hard to soft is included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Machine learning Frenkel Hamiltonian parameters to accelerate simulations of exciton dynamics

In this manuscript, we develop multiple machine learning (ML) models to accelerate a scheme for parameterizing site-based models of exciton dynamics from all-atom configurations of condensed phase sexithiophene systems. This scheme encodes the details of a system’s specific molecular morphology in the correlated distributions of model parameters through the analysis of many single-molecule excited-state electronic-structure calculations. These calculations yield excitation energies for each molecule in the system and the network of pair-wise intermolecular electronic couplings. In this work, we demonstrate that the excitation energies can be accurately predicted using a kernel ridge regression (KRR) model with Coulomb matrix featurization. We present two ML models for predicting intermolecular couplings. The first one utilizes a deep neural network and bi-molecular featurization to predict the coupling directly, which we find to perform poorly. The second one utilizes a KRR model to predict unimolecular transition densities, which can subsequently be analyzed to compute the coupling. We find that the latter approach performs excellently, indicating that an effective, generalizable strategy for predicting simple bimolecular properties is through the indirect application of ML to predict higher-order unimolecular properties. Such an approach necessitates a much smaller feature space and can incorporate the insight of well-established molecular physics.

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

Machine Learning–Augmented Laser-Induced Breakdown Spectroscopy for Spectral Discrimination of Iron Oxalates

Enhanced characterization and phase identification of post-PUREX Pu Oxalates (PuOXA) are pivotal for nonproliferation and pre-detonation nuclear forensics. Despite significant advances in the characterization of PuO 2 samples, little is known about the impact of both the chemical structure and oxidation states of PuOXA (i.e., Pu(III) and Pu(IV)) have on optical emission signatures. Here, we demonstrate the analytical capabilities of laser-induced breakdown spectroscopy (LIBS) applied to Fe(II) and Fe(III) oxalate samples as surrogates for PuOXA, highlighting the discriminating features in the LIBS emission spectra arising from differences in the oxidation states within mixed FeOXA samples. We report the enhancement of spectral feature selection using Principal Component Analysis (PCA), which enables the analytical superiority of machine learning algorithms such as Linear Discriminant Analysis (LDA), Quadratic Discriminant Analysis (QDA), Partial Least Squares Regression (PLSR), Support Vector Regression (SVR), and Random Forest Regression (RFR) over conventional univariate techniques for phase discrimination and chemometric analysis. Cluster analysis revealed how both matrix effects and laser ablation influence cluster separability by introducing spectral artifacts that misdirect the maximization of variance. PCA-selected emission lines were used in the regression models, demonstrating that both univariate and multivariate linear regression models (i.e., PLSR and SVR) can achieve acceptable performance, with machine learning models outperforming conventional calibration regressions. Furthermore, the application of non-linearly activated PCA-selected emission lines illustrates how simplifying the data while retaining captured variance enables the use of less complex and more computationally efficient models. Furthermore, this is particularly evident in the underperformance of RFR, which suffers from increased computational costs and overfitting owing to its high complexity.

Oxalates↗
DOE OSTI2022

Ab initio studies of the double–Gamow-Teller transition and its correlation with neutrinoless double- β decay

Here, we use chiral interactions and several ab initio methods to compute the nuclear matrix elements (NMEs) for ground-state-to-ground-state double Gamow-Teller transitions in a range of isotopes and explore the correlation of these NMEs with those for neutrinoless double beta decay produced by the exchange of a light Majorana neutrino. When all the NMEs of both isospin-conserving and isospin-changing transitions from the ab initio calculations are considered, the correlation is strong. For the experimentally relevant isospin-changing transitions by themselves, however, the correlation is weaker and may not be helpful for reducing the uncertainty in the NMEs for neutrinoless double beta decay.

20 ≤ A ≤ 38↗
DOE OSTI2024

Bayesian mixture model approach to quantifying the empirical nuclear saturation point

The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the nuclear saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), providing tight constraints for microscopic interactions derived from chiral effective field theory (EFT). However, when considering several DFT predictions for $(n_0,E_0)$ from Skyrme and Relativistic Mean Field (RMF) models together, a discrepancy between these model classes emerges at high confidence levels that each model prediction's uncertainty cannot explain. How can we leverage these DFT constraints to rigorously benchmark nuclear saturation properties of chiral interactions? To address this question, we present a Bayesian mixture model that combines multiple DFT predictions for $(n_0,E_0)$ using an efficient conjugate prior approach. The inferred posterior distribution for the saturation point's mean and covariance matrix follows a Normal-inverse-Wishart class, resulting in posterior predictives in the form of correlated, bivariate $t$-distributions. The DFT uncertainty reports are then used to mix these posteriors using an ordinary Monte Carlo approach. At the 95\% credibility level, we estimate $n_0 \approx 0.157 \pm 0.010 \, \mathrm{fm}^{-3}$ and $E_0 \approx -15.97 \pm 0.40 \, \mathrm{MeV}$ for the marginal (univariate) $t$-distributions. Combined with chiral EFT calculations of the pure neutron matter EOS, we obtain bivariate normal distributions for the nuclear symmetry energy and its slope parameter evaluated at $n_0$: $S_v \approx 32.0 \pm 1.1 \, \mathrm{MeV}$ and $L\approx 52.6\pm 8.1 \, \mathrm{MeV}$ (95\%), respectively. Furthermore, our Bayesian framework is publicly available, so practitioners can readily use and extend our results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Source apportionment of soot particles and aqueous-phase processing of black carbon coatings in an urban environment

Abstract. The impacts of soot particles on climate and human health depend on the concentration of black carbon (BC) as well as the thickness and composition of the coating material, i.e., organic and inorganic compounds internally mixed with BC. In this study, the size-resolved chemical composition of BC-containing aerosol was measured using a high-resolution soot particle aerosol mass spectrometer (SP-AMS) during wintertime in Fresno, California, a location influenced by abundant combustion emissions and frequent fog events. Concurrently, particle optical properties were measured to investigate the BC absorption enhancement. Positive matrix factorization (PMF) analysis was performed on the SP-AMS mass spectral measurements to explore the sources of soot particles and the atmospheric processes affecting the properties of BC coatings. The analysis revealed that residential woodburning and traffic are the dominant sources of soot particles. Alongside primary soot particles originating from biomass burning (BBOABC) and vehicles (HOABC) two distinct types of processed BC-containing aerosol were identified: fog-related oxidized organic aerosol (FOOABC) and winter-background OOABC (WOOABC). Both types of OOABC showed evidence of having undergone aqueous processing, albeit with differences. The concentration of FOOABC was substantially elevated during fog events, indicating the formation of aqueous secondary organic aerosol (aqSOA) within fog droplets. On the other hand, WOOABC was present at a relatively consistent concentration throughout the winter and was likely related to the formation of secondary organic aerosol (SOA) in both the gas phase and aerosol liquid water. By comparing the chemical properties and temporal variations in FOOABC and WOOABC, we gain insights into the key aging processes of BC aerosol. It was found that aqueous-phase reactions facilitated by fog droplets had a significant impact on the thickness and chemical composition of BC coatings, thereby affecting the light absorption and hygroscopic properties of soot particles. These findings underscore the important role of chemical reactions occurring within clouds and fogs and influencing the climate forcing of BC aerosol in the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

SAS4A/SASSYS-1 Commercial Grade Dedication Example Report for a Generic Sodium Pool-Type Fast Reactor Application

In the U.S., a key component of the commercialization of advanced reactors is completion of a license application, which must ultimately be approved by the Nuclear Regulatory Commission (NRC). The approval of the license application by the NRC is contingent on satisfactory demonstration of the design basis and the response of the advanced reactor design to transient and accident scenarios using accepted codes and methods. This report describes the qualification and dedication requirements that the advanced reactor safety analysis system software SAS4A/SASSYS-1 are expected to need to fulfill to be used for sodium-cooled pool-type fast reactor licensing. The qualification and dedication requirements are identified through performance critical characteristics and evaluation model acceptance criteria representative of the advanced reactor design considered for licensing. This document captures, additionally, the verification process developed to demonstrate that the software fulfills the qualification and dedication requirements for a generic sodium-cooled pool-type fast reactor as part of the commercial grade dedication process. Like most software that has primarily existed in the research and development space, the most significant challenge facing SAS4A/SASSYS-1 for use in a licensing framework is the availability of a documentation basis describing the code pedigree. SAS4A/SASSYS-1 has been used for licensing of the fast flux test facility (FFTF) and the JOYO sodium-cooled fast reactor in Japan, as well as the design of the CRBR Plant. However, the historical verification and validation (V&V) activities supporting SAS4A/SASSYS-1 development do not align with modern software quality assurance (SQA) and V&V requirements. Two approaches to use of SAS4A/SASSYS-1 in a commercial licensing framework have been identified: commercial-grade dedication (CGD) and software qualification. The methods and requirements prescribed in the ASME NQA-1-2008/2009 Standard and Regulatory Guide 1.203 on the evaluation model development and assessment process (EMDAP) have been used as guidance to define the CGD and qualification processes, respectively. A qualification and dedication requirements matrix has been developed which utilizes fundamental software verification. In this process, software verification is defined as a software quality process aimed at defining software requirement specifications, developing software design documentation, and performing and documenting acceptance testing of the code against requirements. A key element of software qualification and dedication includes determination of software acceptance with respect to critical characteristics relevant to the functional requirements of the software. To assist with identification of cross-cutting transient phenomena and functional requirements, domestic SFR vendor designs have been reviewed to identify a reference SFR design. For this report, the reference design is defined as a pool-type reactor with metal alloy fuel, a liquid-metal intermediate heat transport system, and passive decay heat rejection systems. Given this reference, a series of high-level cross-cutting phenomena was identified for a general class of single-fault undercooling or reactivity insertion transients that scopes the design basis space, with the goal of assisting with prioritization of documentation development efforts for key transient models in SAS: 1) Reactivity feedback response prior to scram; 2) System-wide thermal inertia; 3) Transition in natural circulation flow regime in heat removal systems; 4) Decay heat generation; 5) Steady-state fuel characterization; 5) Clad/fuel behavior at elevated temperatures; 6) Point kinetics and decay heat; 7) Pump coastdown behavior; 8) Core flow redistribution in loss of forced convection; 9) Pool stratification. As a demonstration of CGD of SAS4A/SASSYS-1 for a sodium pool reactor, a software qualification and dedication gap analysis as it relates to code documentation has been performed. This effort leverages the framework established as part of the SAS4A/SASSYS-1 SQA Program. This CGD demonstration provides a framework that vendors can build upon to demonstrate the applicability of the SAS4A/SASSYS-1 software for licensing a sodium-cooled pool-type fast reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Mapping Thermal Conductivity at the Atomic Scale: A Step toward the Thermal Design of Materials

We describe a spatial decomposition of the thermal conductivity, termed site-projected thermal conductivity, which quantifies the thermal conduction activity at each atomic site—a critical parameter for the thermal design of materials. The method is based on the Green–Kubo formula and the harmonic approximation and requires the force-constant and dynamical matrices, as well as a relaxed structural model. Throughout the study, it uses high quality models previously tested and compared to many experiments. It discusses the method and underlying approximations for amorphous silicon, carries the detailed analysis for amorphous silicon, and then examines an amorphous-crystal silicon interface and representative carbon materials. This study identifies the sites and local structures that reduce heat transport, and quantify these (estimate the spatial range) over which these “thermal defects” are effective. It identifies filamentary structures in the amorphous silicon network which impact heat transport and electronic structure (the Urbach edge) and electronic transport.

36 MATERIALS SCIENCE↗
DOE OSTI2022

7.2 kV Three-Port SiC Single-Stage Current-Source Solid-State Transformer With 90 kV Lightning Protection

This article proposes a multiport modular single-stage current-source solid-state transformer (SST) for applications like photovoltaic, energy storage integration, electric vehicle fast charging, data center, etc. The 7.2 kV 50 kVA current-source SST consists of five input-series output-parallel modules, each based on 3.3 kV SiC reverse-blocking MOSFET-plus-diode modules. The proposed SST has some unique features. First, compared to the voltage-source or matrix converter-based SSTs, the current-source SST has a unique advantage of single-stage isolated AC/DC or AC/AC conversion with an inductive DC link, but no medium-voltage (MV) AC experiments have been reported. This article for the first time demonstrates MV AC current-source SST up to 7.5 kV peak. Second, the multiport SST has a buffer port for active power decoupling (APD) or energy storage integration. The double-line-frequency power ripple from single-phase AC grid normally results in a large capacitor size in MV SSTs. The APD scheme is proposed in MV applications for the first time to enable a reduced DC link and the electrolytic capacitor-less SST with high reliability. Third, as a direct grid-connected converter without line-frequency transformer, insulation and protection are critical. A medium-frequency transformer design passes 55 kV basic-insulation level (BIL) and 60 kV high potential dielectrics withstand test with only 0.09% leakage inductance. Importantly, a lightning protection scheme is presented to protect the SST itself from 90 kV BIL impulse. Fourth, the proposed current-source SST topology is a modular soft-switching solid-state transformer (M-S4T) with full-range zero-voltage switching and controlled dv/dt for low electromagnetic interference. Furthermore, these concepts are verified in a three-port M-S4T prototype with forced oil cooling under single-module, stacked-module, steady-state, and dynamic operations.

14 SOLAR ENERGY↗
DOE OSTI2021

Comparing different density-matrix expansions for long-range pion exchange

Empirical energy density functionals (EDFs) are generally successful in describing nuclear properties across the table of nuclides. But their limitations motivate using the density-matrix expansion (DME) to embed long-range pion interactions into a Skyrme functional. Recent results on the impact of the pion were both encouraging and puzzling, necessitating a careful re-examination of the DME implementation. In this study, we take the first steps, focusing on two-body scalar terms in the DME. Exchange energies with long-range one-pion contributions are well approximated by all DME implementations considered, with preference for variants that do not truncate at two derivatives in every EDF term. The use of the DME for chiral pion contributions is therefore supported by this investigation. For scalar-isovector energies it is important to treat neutrons and protons separately. The results are found to apply under broad conditions, although self-consistency is not yet tested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Present Status of Nuclear Shell-Model Calculations of 0 νββ Decay Matrix Elements

Neutrinoless double beta (0νββ) decay searches are currently among the major foci of experimental physics. The observation of such a decay will have important implications in our understanding of the intrinsic nature of neutrinos and shed light on the limitations of the Standard Model. The rate of this process depends on both the unknown neutrino effective mass and the nuclear matrix element (M 0ν ) associated with the given 0νββ transition. The latter can only be provided by theoretical calculations, hence the need of accurate theoretical predictions of M 0ν for the success of the experimental programs. This need drives the theoretical nuclear physics community to provide the most reliable calculations of M 0ν . Among the various computational models adopted to solve the many-body nuclear problem, the shell model is widely considered as the basic framework of the microscopic description of the nucleus. Here, we review the most recent and advanced shell-model calculations of M 0ν considering the light-neutrino-exchange channel for nuclei of experimental interest. We report the sensitivity of the theoretical calculations with respect to variations in the model spaces and the shell-model nuclear Hamiltonians.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗