Search NASA⌕ Search

SEARCH · Search NASA

Results for “Generalized standard materials”

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 343 records · Page 19

LA17-V-InverseRadTransportProblem-PD2Kb. Final Project Report

In the inverse radiation transport problem, a radiation signal emanating from an unknown object is used to reconstruct geometric and material parameters of the object , e.g., the mass of special nuclear material (SNM), the isotopic composition of the SNM (e.g., the enrichment of 235U), the thickness of shielding, etc. The inverse radiation transport problem is one of the most fundamental and important in the emergency response mission, and has been addressed in previous NA-22 and NA-42/84 investments by using derivative-free optimization methods. Our project was aimed at (i) consolidating previous work by applying polynomial chaos expansion (PCE) for reducing significantly the number of forward model simulations required by derivative-free optimization methods; (ii) developing a general-purpose framework and stand-alone software for solving the inverse radiation transport problem by applying a new “predictive modeling” methodology (predictive modeling of coupled multi-physics systems, PM_CMPS) which enables the fusing of experimental and computational information; and (iii) developing several benchmarks (for unscattered gamma rays) for demonstrating the quantification of non-Gaussian features (e.g., skewness and kurtosis, which characterize asymmetries and long tails) of response distributions by deterministically and efficiently propagating parameter uncertainties (covariances and higher-order moments, where available) using high-order sensitivities computed by the recently developed comprehensive adjoint sensitivity analysis methodology (C-ASAM). With regard to (i), application of PCE, the most significant accomplishment of this project was an implementation of PCE into INVERSE that significantly reduced computation time for the DiffeRential Evolution Adaptive Metropolis (DREAM) and Differential Evolution (DE) methods. For inverse transport problems with measurements of neutron leakage, neutron multiplication, and neutron-induced gamma-rays, the solution times for the DREAM and DE methods were reduced by factors of 30-300, reducing run times from hours to minutes while maintaining the same accuracy as the standard (non-PCE) DREAM and DE. Speedups were further improved by exploring sparse grid quadrature sets using the Toolkit for Adaptive Stochastic Modeling and Non-Intrusive ApproximatioN (Tasmanian) toolset. An initial exploration of sparse grids indicated that further speedups factors of 5-10 are obtainable. The implementation of PCE in INVERSE is not in a form that users can easily access. With regard to (ii), solving the inverse problem using a predictive modeling paradigm, the most significant accomplishment of this project were the derivation and implementation of adjoint-based second derivatives in neutron and photon transport theory. There are two codes that compute these derivatives. SENSPG computes them for uncollided gamma rays and SENSMG computes them for neutron leakage. SENSMG is open-source code available on github. Other significant accomplishments included a stand-alone code, MULTI-PRED, that implements the predictive modeling paradigm, along with a detailed manual. The code was implemented in the existing INVERSE code for inverse problems, but a head-to-head comparison of inverse methods was not achieved. We also derived and implemented the capability to compute sensitivities for multiplicity measurements into the SENSMG code, and this capability was used in a data assimilation application, but it was not used for solving inverse problems. With regard to (iii), developing benchmarks to demonstrate the importance of high-order sensitivity and uncertainty analysis, the most significant accomplishment of this project were two test problems that showed how important second-order sensitivities are.

07 ISOTOPE AND RADIATION SOURCES↗

Continuum Damage Mechanics Modeling of High-Temperature Flaw Propagation: Application to Creep Crack Growth In 316H Standardized Specimens and Nuclear Reactor Components

Predicting creep crack growth (CCG) of flaws found during operation in high-temperature alloy components is essential for assessing the remaining lifetime of those components. While defect assessment procedures are available for this purpose in design codes, these are limited in their range of applicability. This study assesses the application of a local damage-based finite-element methodology as a more general technique for the prediction of CCG at high temperatures on a variety of structural configurations. Numerical results for stainless steel 316H, which are validated against experimental data, show the promise of this approach. This integration of continuum damage mechanics (CDM) based methodologies, together with adequate inelastic models, into assessment procedures can therefore inform the characterization of CCG under complex operating conditions, while avoiding excessive conservatism. This article shows that such modeling frameworks can be calibrated to experimental data and used to demonstrate that the degree of tri-axiality ahead of a growing creep crack affects its rate of growth. The framework is also successfully employed in characterizing CCG in a realistic reactor pressure vessel geometry under an arbitrary loading condition. These results are particularly relevant to the nuclear power industry for defect assessment and inspections as part of codified practices of structural components with flaws in high-temperature reactors.

42 ENGINEERING↗

Moiré metrology of energy landscapes in van der Waals heterostructures

The emerging field of twistronics, which harnesses the twist angle between two-dimensional materials, represents a promising route for the design of quantum materials, as the twist-angle-induced superlattices offer means to control topology and strong correlations. At the small twist limit, and particularly under strain, as atomic relaxation prevails, the emergent moiré superlattice encodes elusive insights into the local interlayer interaction. Here we introduce moiré metrology as a combined experiment-theory framework to probe the stacking energy landscape of bilayer structures at the 0.1 meV/atom scale, outperforming the gold-standard of quantum chemistry. Through studying the shapes of moiré domains with numerous nano-imaging techniques, and correlating with multi-scale modelling, we assess and refine first-principle models for the interlayer interaction. We document the prowess of moiré metrology for three representative twisted systems: bilayer graphene, double bilayer graphene and H-stacked MoSe 2 /WSe 2 . Moiré metrology establishes sought after experimental benchmarks for interlayer interaction, thus enabling accurate modelling of twisted multilayers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Machine-learning accelerated geometry optimization in molecular simulation

Geometry optimization is an important part of both computational materials and surface science because it is the path to finding ground state atomic structures and reaction pathways. These properties are used in the estimation of thermodynamic and kinetic properties of molecular and crystal structures. This process is slow at the quantum level of theory because it involves an iterative calculation of forces using quantum chemical codes such as density functional theory (DFT), which are computationally expensive and which limit the speed of the optimization algorithms. It would be highly advantageous to accelerate this process because then one could do either the same amount of work in less time or more work in the same time. Here, we provide a neural network (NN) ensemble based active learning method to accelerate the local geometry optimization for multiple configurations simultaneously. We illustrate the acceleration on several case studies including bare metal surfaces, surfaces with adsorbates, and nudged elastic band for two reactions. In all cases, the accelerated method requires fewer DFT calculations than the standard method. In addition, we provide an Atomic Simulation Environment (ASE)-optimizer Python package to make the usage of the NN ensemble active learning for geometry optimization easier.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Progress Report on Graphite-Salt Interaction Studies (FY21)

This report summarizes the activities performed in FY21 to investigate the interactions between nuclear graphite and molten salts, such as FLiNaK. First, building on last year’s achievements, the team improved the procedure for measurements on pressurized salt intrusion in nuclear graphite and performed new experiments at variable pressures with specimens of different sizes and shapes. Additionally, the team actively provided suggestions and comments for ASTM D8091, Standard Guide for Impregnation of Graphite with Molten Salts , which came up for periodic revision in 2021. Second, the team expanded its capabilities to further characterize salt-impregnated graphite by procuring and installing a laser-induced breakdown spectroscopy (LIBS) instrument capable of producing a 3D chemical composition of a surface and subsurface region of specimens. A methodology for accurately analyzing information was developed, and a manuscript was submitted for publication. The LIBS instrument sample cell is sealed in controlled inert atmosphere, which is an advantage for using humidity-sensitive samples and Be-contaminated samples. This aligns with the team’s effort to expand the capabilities of handling FLiBe and Becontaining materials in two newly installed glovebox units. Third, the team performed preliminary tribology tests for pebble graphite wear in contact with stainless steel in dry state and in molten salt at high temperature. These tests are needed to establish the baseline for the wear of pebble graphite by friction against metallic walls. The friction and wear of pebbles in the fluoride salt-cooled hightemperature reactor are expected to generate dust, which is a concern for other components’ properties and safe reactor operation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Mechanistic Elucidation of Electronically Conductive PEDOT:PSS Tailored Binder for a Potassium‐Ion Battery Graphite Anode: Electrochemical, Mechanical, and Thermal Safety Aspects

Potassium-ion batteries (KIBs) are considered more appropriate for grid-scale storage than lithium-ion batteries (LIBs) due to similar operating chemistry, abundant precursors, and compatibility with low-cost graphite anodes. However, a larger ion reduces rate capabilities and exacerbates capacity fading from volumetric expansion. In this report, conductive polymer, poly(3,4-ethylenedioxythiophene) polystyrene sulfonate (PEDOT:PSS), is substituted for standard insulating polyvinylidene fluoride (PVDF). Half-cells using carbon black (CB) in continuously conductive PEDOT:PSS/CB binder outperforms PVDF/CB by mitigating electrically isolated “dead” graphite, improving 100 cycle capacity retention at C/10 from 63 to 80%. Enhanced electrical contact with PEDOT:PSS/CB also reduces ion impedance and improves rate capabilities. Without CB however, PEDOT:PSS binder performs poorly in electrochemical studies despite promising ex situ electronic conductivity. This discrepancy is mechanistically elucidated through identification of redox activity between PEDOT:PSS and K + which results in high impedances in the anode operating voltage window. Additionally, the impact of conducting binder on mechanical properties and thermal safety of the anode is investigated. Brittleness and poor wettability of PEDOT:PSS are identified as issues, but greater stability against reactive KC 8 reduces overall heat generation. Binder substitution offers a promising means of mitigating issues with current KIB anodes regardless of active material, and the work herein addresses issues towards further improvement.

electrochemical impedance spectroscopy↗

Techno-economic analysis of advanced small modular nuclear reactors

Here, small modular nuclear reactors (SMRs) represent a robust opportunity to develop low-carbon and reliable power with the potential to meet cost parity with conventional power systems. This study presents a detailed, bottom-up economic evaluation of a 12 × 77 MW e (924 MW e total) light-water SMR (LW-SMR) plant, a 4 × 262 MW e (1,048 MW e ) gas-cooled SMR (GC-SMR) plant, and a 5 × 200 MW e (1,000 MW e total) molten salt SMR (MS-SMR) plant. Cost estimates are derived from equipment costs, labor hours, material inputs, and process-engineering models. The advanced SMRs are compared to natural gas combined cycle plants with and without post-combustion carbon capture and a conventional large nuclear reactor. Overnight capital cost (OCC) and levelized cost of energy (LCOE) estimates are developed. The OCC of the LW-SMR, GC-SMR, and MS-SMR are found to be $\$4,844$/kW, $\$4,355$/kW, and $\$3,985$/kW respectively. The LCOE of the LW-SMR, GC-SMR, and MS-SMR are found to be $\$89.6$/MWh, $\$81.5$/MWh, and $\$80.6$/MWh respectively. A Monte Carlo analysis is performed, for which the OCC and construction time of the LW-SMR is found to have a lower mean and standard deviation than a conventional large reactor. The LW-SMR OCC is found to have a mean of $\$5,233$/kW with a standard deviation of $\$658$/kW and a 90 % probability of remaining between $\$4,254$/kW and $\$6,399$/kW, while the construction duration is found to have a mean of 4.5 years with a standard deviation of 0.8 years and a 90 % probability of remaining between 3.4 and 6.0 years. The economic impact of economies of scale, simplification, modularization, and construction time for SMRs are discussed. Additionally, policy implications for direct SMR capital subsidies and the impact of a carbon tax on natural gas emissions are explored.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Recipe for coating ceramic blades for ion trapping

The first batches of ion traps patterned and coated were processed per the standard 3-step clean, air fire, and metallization processes. The third or fourth lot using this process resulted in poorly adhering metallization. Up until this point, the standard process was used to metallize and pattern ceramic ion traps without fail. At about the 4th batch of parts something changed. After the 5th batch, the ceramic ion traps received generally came with some unknown contamination that does not come off in a standard 3-step clean (Lenium Vapor Degreaser, Acetone, IPA) and air fire (860C for 1 hour) for which this process removes the vast majority of all contamination for most ceramic metallization. This is highly unusual. Using HF + Boiling H 2 O 2 is extreme for cleaning the ceramic ion traps. The contamination was never identified and is stubborn to effectively clean. Standard as-fired ceramic should be very easy to clean as if s fired at temperatures greater than 1400°C and not much in terms of contamination should exist at these temperatures, so there must be an intermediate step/process which is imparting this contamination. It is likely a polishing compound or previous polishing contaminant, but also not easily visually distinguishable until after metallization. The halo marks observed on parts might be fingerprints (less likely) or potential polishing marks (more likely) as metallization typically doesn't cover/hide any damage or contamination, but rather quite clearly the opposite, it accentuates it. Blotchy appearances in the metallization usually indicated an adhesion issue. As a result of the fragility of the parts (yield loss due to handling) and difficulty in identifying the contamination during cleaning, we have taken a conservative approach of HF + H 2 O 2 cleaning for all batches after the contamination and adhesion issues were identified.

36 MATERIALS SCIENCE↗

Safety-Related Instrumentation and Control Pilot Upgrade (Initial Scoping Phase Implementation Report)

This research report (1) describes the process followed and products developed during the SR I&C Pilot Project Initial Scoping Phase, and (2) captures lessons learned. Exelon Generation and LWRS collaborated to develop a Digital Transformation Strategy as part of a larger Advanced Concept of Operations. The proposed SR I&C Pilot Upgrade provides a foundation stone for this Digital Transformation that will improve plant safety, reliability, and operational performance while lowering plant Total Cost of Ownership (TCO). Initial Scoping Phase activities for this Pilot Project have been performed in accordance with industry processes that have been adapted to better support digital upgrades. These processes include IP-ENG-001, Standard Design Process (SDP), NISP-EN-04, Standard Digital Engineering Process (SDEP), and Electric Power Research Institute (EPRI) Report 3002011816, Digital Engineering Guide (DEG). Completing Initial Scoping Phase Engineering and Operations, Licensing, and Project Management Activities was necessary to sufficiently bound the scope, schedule, and estimated cost of the Project to enable utility management to authorize moving into the Conceptual Design Phase. A significant finding of the Business Case Analysis (BCA) methodology developed and applied as part of this effort was that the growth rate of material costs for sustaining the operation of obsolete SR I&C equipment is accelerating. This directly contributed to the Project Economic Analysis created to justify continuing the Project. Project Initial Scoping Phase lessons learned have also been captured to assist the larger industry in understanding the Digital Transformation Strategy and SR I&C Pilot Project Initial Scoping Phase efforts. This is in keeping with the public/private partnership that has been established between the Department of Energy (DOE) and Exelon for this effort with engagement from the NRC. By addressing first-of-a-kind (FOAK) risks and capturing lessons learned, the SR I&C Pilot Upgrade Project addresses technical, regulatory, and business risks to enable subsequent implementers of similar upgrade.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A high strength Al-2Ni-0.5Zr conductor alloy fabricated via laser powder bed fusion

There is a current need for new aluminum alloy design strategies to target applications requiring high strength and conductivity with reductions in mass. A new lightweight Al-2Ni-0.5Zr (wt. %) conductor alloy was fabricated using laser powder bed fusion. A design of experiments probed the alloy's solidification cracking susceptibility. It was observed that solidification cracking was generally reduced with fast scan speeds, above 1500 mm/s, and smaller hatch spacings. The different cooling rates throughout the melt pool produced a heterogeneous distribution of cellular and equiaxed Al 3 Ni precipitates in the as-printed alloy. Additionally, the rapid solidification characteristic of laser powder bed fusion created a super-saturated Zr solid solution. An aging heat treatment at 375 °C for 24 h imparted strengthening through the precipitation of L1 2 -Al 3 Zr nanoprecipitates, which counteracted the softening caused by the fragmentation and coarsening of Al 3 Ni precipitates. The yield strength increased from 138 MPa in the as-printed condition to 168 MPa after aging, while the ductility remained constant at ∼21%. The aging treatment simultaneously increased the electrical conductivity from 40.8% IACS (International Annealed Copper Standard) to 53.5% IACS. Modeling of the strengthening mechanisms and electrical conductivity contributions rationalized the simultaneous increase in strength and conductivity upon aging. Furthermore, the strengthening efficacy of the Al 3 Ni and L1 2 -Al 3 Zr precipitates, combined with the low Ni and Zr solubility in the FCC Al matrix, facilitated both high strength and electrical conductivity. Overall, the combination of strength and electrical conductivity positions this alloy as a suitable choice for additively manufactured lightweight conductors.

Additive manufacturing↗

A simplified approach for creep damage and remaining life calculation for materials surveillance in advanced reactors

High temperature reactor materials will experience a combination of mechanical degradation caused by creep and fatigue and environmental degradation caused by neutron irradiation and, in some cases, exposure to corrosive coolants. A materials surveillance program is one option for ensuring the safe, reliable operation of key reactor components under these conditions. Such a program would monitor the degradation in key material properties over time, using this data to predict how changes in the material properties affect components performance. This could then be used to make plant operational decisions and support eventual plant license extensions. Past work de-scribes such a material surveillance program based on passively actuated mechanical test articles that impose creep-fatigue type loading on test materials driven only by changes in temperature, for example those experienced by components under standard operating cycles. This report focuses on two aspects of the proposed material surveillance program: how to design the test articles to mimic the mechanical response of the corresponding component and how to use ex-situ test data to deter-mine the amount of damaged experienced by the test article (and hence the corresponding compo-nent) in service. Specifically, this report derives and validated simplified methods for both tasks, replacing earlier, more complicated approaches based on simulating the test article response and finding the best specimen design or current damage via complex numerical optimization. The sim-plified methods for each task developed here can be implemented in spreadsheet software and are simple enough for practical use in future operating plants. This report derives the methods (one for sizing and two options for damage inference), verifies the simplified approaches versus more so-phisticated methods, and compares the results of applying the simplified methods to previous re-sults using the complex numerical optimization approach. The report also provides a worked sam-pled problem applying the simplified techniques to a realistic high temperature reactor component.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Stable Higgs mode in anisotropic quantum magnets

Low-energy excitations associated with the amplitude fluctuation of an order parameter in condensed-matter systems can mimic the Higgs boson, an elementary particle in the standard model, and are dubbed as Higgs modes. Identifying the condensed-matter Higgs mode is challenging because it is known in many cases to decay rapidly into other low-energy bosonic modes, which renders the Higgs mode invisible. Therefore, it is desirable to find a way to stabilize the Higgs mode, which can offer an insight into the stabilization mechanism of the Higgs mode in condensed-matter physics. In quantum magnets, magnetic order caused by spontaneous symmetry breaking supports transverse (magnons) and longitudinal (Higgs modes) fluctuations. Additionally, when a continuous symmetry is broken, the Goldstone magnon mode generally has an excitation energy lower than that of the Higgs mode, causing a rapid decay of the latter. In this work, we show that a stable Higgs mode exists in anisotropic quantum magnets near the quantum critical point between the dimerized and magnetically ordered phases. We find that an easy axis anisotropy increases the magnon gap such that the magnon mode is above the Higgs mode near the quantum critical point, and the decay of the Higgs mode into the magnon mode is forbidden kinematically. Our results suggest that the anisotropic quantum magnets provide ideal platforms to explore the Higgs physics in condensed-matter systems.

36 MATERIALS SCIENCE↗

Approach to Quality Assurance for Complex Environmental Modeling - 20407

A strong Quality Assurance (QA) program for complex environmental modeling is essential for regulatory and public acceptance and trust, but it does not have to be onerous. Neptune and Company, Inc. (Neptune) has developed a strong QA program that improves transparency, traceability, reproducibility, and therefore, defensibility and trust. Neptune's QA program has evolved over the past 27 years, transitioning from an ad hoc QA program, to a program that is currently Nuclear Quality Assurance-1 (NQA-1) compliant, and will soon be NQA-1/DOE approved for DOE EM modeling work. In addition, Neptune is a qualified laboratory assessor and qualified auditor for reference materials, a proficiency testing provider, and fully complies with the Environmental Protection Agency (EPA) QA program. Neptune's President and CEO, Kelly Black, was recently named the Chairperson of the International Organization for Standardization (ISO) Technical Committee 69, Application of Statistical Methods. A strong QA program has been developed for Neptune's radiological performance assessment (PA) program. Although pieces of their QA program are currently in development, their current program includes document control using Subversion, issue tracking and work flow tracking using JIRA, and transparency and traceability via a system of 'calc sheets' for documentation of all data analysis and modeling combined with 'check print' documentation of QA checking, and rigorous model testing and configuration control. All work is reviewed by an independent subject matter expert who is not associated with the collection and assembly of information, for an internal peer review. The program is enforced with a handful of Standard Operation Procedures (SOPs), Work Instructions, QA Project Plans (QAPPs), and Quality Management Plans (QMPs) that are updated frequently, with required annual training and acknowledgment. In addition, effective communication of modeling approaches and results to clients and stakeholders are integral to their QA program. Many of their models are built using the GoldSim modeling platform, for which their models are well-known for their level of transparency, documentation, and QA traceability. This level of QA is also applied to their process-level models. Neptune has taken some lessons learned from the extremely rigorous QA program required for the Yucca Mountain Project (YMP), and imposed the YMP requirements of complete traceability, transparency, and reproducibility, but avoided the inflexibility of a QA program that likely contributed to the suspect e-mails in 2005 that resulted in loss of trust in the YMP, and in nuclear waste disposal (and nuclear power) in general. There are far too many examples of loss of public trust due to poor QA that could have been easily avoided with a simple and straight-forward QA program combined with stakeholder engagement. The purpose of this paper is to share features of Neptune's Radiological Performance Assessment program QA program, and some of their QA related lessons learned over the past 27 years of QA for complex environmental modeling. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Hybrid Tandem Photovoltaics

Tandem solar cell structures are the only strategy demonstrated to surpass the detailed balance efficiency limit of high-quality single-junction solar cells. To continue to improve the efficiencies of cost-effective terrestrial solar power, hybrid tandems of dissimilar subcells are being considered by many around the world, especially designs that incorporate silicon solar cells as a bottom subcell. In this project, we studied a wide variety of tandem design possibilities including those with three-terminal (3T) and four-terminal (4T) configurations. The use of 3T and 4T designs could be useful for efficient and economical hybrid tandem designs that utilize the best available subcell materials such as emerging perovskite materials. Three-terminal configurations, in particular, have not been sufficiently studied previously. We have laid the foundational groundwork in this project for understanding the operation of 3T tandems: developing a taxonomy for naming, a methodology for measuring and interconnecting, and models for simply characterizing 3T tandems. Electrical and optical subcell coupling between the subcells was also measured and modeled. An important part of this work was the fabrication of novel example tandem structures, including 4T GaAs/Si, 3T GaInP/Si, 3T GaAs/Si, and 3T GaInP/GaAs devices. Using these high-quality tandem cells, we have been able to clearly demonstrate the achievability of high-efficiencies, and subtle physical effects such as photon recycling and luminescent coupling. We have developed and demonstrated essential building-block tools such as transparent conductive adhesives (TCA) and 3T silicon bottom cells with interdigitated back contacts (IBC) that can also be used in many other tandem designs. We have tested the reliability of these tools and devices under standardized testing and outdoor measurements. We have found 4T GaAs/Si tandems to be relatively straightforward to fabricate and robust in real-world outdoor conditions. While we have demonstrated working hybrid 3T III-V/TCA/Si IBC tandems, we experienced low yields even with our best process flows yet. Further work is still needed to improve the processing yield of these devices. We therefore also created tandem cells using an all-III-V 3T tandem process which was very robust with high yields, allowing for the creation of voltage-matched strings in many different configurations using 8 nearly identical 3T tandems. Using these robust 3T tandem examples, we were able measure and precisely characterize 3T tandem behaviors to predict their operation under changing spectrum and temperature. The optoelectronic equivalent-circuit model was shown to be very general and applicable to hybrid tandems, and encompassed the operation 3T Si IBC cells. This general model has been distributed to the public in as open-source Python-based software called PVcircuit. We have calculated the implications of these new tandem device designs on the real-world energy production and shown how the relative performance of different tandem configurations is situational and can be engineered using the tools developed here.

14 SOLAR ENERGY↗

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

BSE versus StarTrack: Implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

As a result of their formation via massive single and binary stellar evolution, the masses of stellar-remnant black holes (BH) are subjects of great interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star merger events. In this work, we present new developments in the stellar-remnant formation and related schemes of the current N-body evolution program NBODY7. We demonstrate that the newly implemented stellar-wind and remnant-formation schemes in the stellar-evolutionary sector or BSE of the NBODY7 code, such as the “rapid” and the “delayed” supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produce neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main-sequence (ZAMS) mass and metallicity, with those from the widely recognised StarTrack population-synthesis program. We also demonstrate the new, recipe-based implementations of various widely debated mechanisms of natal kicks on NSs and BHs, such as “convection-asymmetry-driven”, “collapse-asymmetry-driven”, and “neutrino-emission-driven” kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. All the above newly implemented schemes are also shared with the standalone versions of SSE and BSE. Furthermore, all these demonstrations are performed with both the updated standalone BSE and the updated NBODY7/BSE. When convolved with stellar and primordial-binary populations as observed in young massive clusters, such remnant-formation and natal-kick mechanisms crucially determine the accumulated number, mass, and mass distribution of the BHs retained in young massive, open, and globular clusters (GCs); these BHs would eventually become available for long-term dynamical processing. Among other conclusions, we find that although the newer, delayed SN remnant formation model gives birth to the largest number (mass) of BHs, the older remnant-formation schemes cause the largest number (mass) of BHs to survive in clusters, when incorporating SN material fallback onto the BHs. The SN material fallback also causes the convection-asymmetry-driven SN kick to effectively retain similar numbers and masses of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry-driven SN kick would cause nearly all BHs to be retained in clusters irrespective of their mass, remnant-formation model, and metallicity, whereas the inference of a large population of BHs in GCs would potentially rule out the neutrino-driven SN kick mechanism. Pre-SN mergers of massive primordial binaries would potentially cause BH masses to deviate from the theoretical, single-star ZAMS to mass-remnant mass relation unless a substantial of the total merging stellar mass of up to ≈40% is lost during a merger process. In particular, such mergers, at low metallicities, have the potential to produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster and be available for subsequent dynamical interactions. As recent studies indicate, the new remnant-formation modelling reassures us that young massive and open clusters would potentially contribute to the dynamical BBH merger detection rate to a similar extent as their more massive GC counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗

Raman spectroscopy of thermally perturbed carbon fibers: Discriminating spectral responses of modulus classes and defect types

The mechanical strength properties of carbon fibers are generally thought to be correlated with the presence of underlying atomic-level defects. These defects serve to break the underlying translational symmetry of the graphitic or graphene subunits of the fiber construction, resulting in the emergence of new Raman-active spectral features. However, historical attempts to classify the precise origin of defect contributions to the Raman spectra have been challenging because of indistinct and overlapping features in the carbon fiber Raman spectra. Further, while substantial research exists on high-temperature exposure in inert atmospheres for carbon fiber composites, comparatively less addresses microscale behaviors and Raman spectral alterations to monofilament carbon fibers exposed to high-temperature atmospheric environments. Here, we report Raman spectral responses of nine commercially available high-performance, polyacrylonitrile-based carbon fibers exposed to various atmospheric heat treatments. In this work, we introduce a model-independent characterization method, the integrated absolute difference, to quantify the spectral responses to heat treatment of different fiber modulus classes. We combine this new method with a newly reported strategy to standardize spectral fitting for carbon fibers. With this approach, we show that atomic-scale defects in the underlying fiber microstructure manifest in measurably distinct manners and have distinct responses to thermal perturbation. These combined approaches may lay the foundation for disentangling contributions from specific defects in the Raman spectra of carbon fibers.

36 MATERIALS SCIENCE↗