Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Material Point Method”

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

Hand Calculation Methods for Nuclear Criticality Safety

This primer provides an overview of the most common hand calculation methods used for criticality safety calculations. The most widely used tools available to a nuclear criticality safety (NCS) practitioner are probably the common Monte Carlo or deterministic criticality safety codes, which can be used to model very complex systems. However, use of these codes can obscure the parameters to which a particular fissile system may be sensitive, whereas the hand calculation methods can be used to delve into the ways each parameter may affect the reactivity of a fissile material system. Furthermore, practitioners must avoid using computer codes as devices that take inputs and simply provide outputs (i.e., a “black box”). Many years ago, pioneers such as Joe Thomas, David Smith, and Hugh Paxton, among others in the field of nuclear criticality safety, took the time before the advent of high-speed desktop computers to create simple hand methods for criticality safety analyses. Some of the methods can be used for single fissile units; others are applicable to fissile units arranged into simple array configurations. This primer discusses the applicability of the various methods, illustrates how they are used, and provides an interpretation of the various results. The NCS practitioner will need to spend time to master the methods that could be most useful; however, they can provide the practitioner with fast and accurate answers to criticality safety problems if they are used correctly and if critical data exist for the problem at hand. Hand calculation methods can be used as a starting point for more advanced calculations, and in many circumstances, they can provide sensitivity and perturbation information more quickly than using a criticality code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Upcycling Polyethylene Waste to Advanced Carbon Materials for Energy Storage Applications

Polyethylene is notoriously difficult to upcycle because it decomposes into light gases at approximately 350-400 °C which prevents processing it at higher temperatures to convert it into high-value carbon materials, such as graphene or graphite. We address this challenge by using an air-based, thermo-oxidation process, which heats polyethylene (PE) just below the decomposition point to initiate oxidation and cross-linking of PE alkyl chains. These molecular transformations allow PE to be further graphenized or graphitized at higher temperatures without decomposing. The PE-derived graphene has a specific surface area up to 1,800 m 2/g and Raman ID/IG ratio of 0.85, which enables it to be used as an electrode material for a symmetric supercapacitor with the 1 M H2SO4 electrolyte. The PE-derived graphene material has a comparable electrochemical capacitive performance, such as power density, specific capacitance, and long cycle stability, to the commercial state of art porous carbon electrode. The stabilized PE could also be converted into highly crystalline flake graphite via low-temperature catalytic graphitization. The PE-derived graphite is used for the lithium-ion battery anode, which showed comparable electrochemical battery performance, such as reversible rate performance and long-term cyclic stability, to the current use battery-grade graphite, thus providing a scalable method to upcycle PE plastic waste into high-value carbon material.

Gao, Yuan↗
DOE OSTI2023

Solution and sensitivity analysis of nonlinear equations using a hypercomplex-variable Newton-Raphson method

Here, the classical Newton-Raphson (NR) method for solving nonlinear equations is enhanced in two ways through the use of hypercomplex variables and algebra. In particular, i) the Jacobian is computed in a highly accurate and automated way, and ii) the derivative of the solution to the nonlinear equations is computed with respect to any parameter contained within the system of equations. These advances provide two significant enhancements in that it is straightforward to provide an accurate Jacobian and to construct a reduced order model (ROM) of arbitrary order with respect to any parameter of the system. The ROM can then be used to approximate the solution for other parameter values without requiring additional solutions of the nonlinear equations. Several case studies are presented including 1D and 2D academic examples with fully functioning Python code provided. Additionally, a case of study of the catenary of an elastic cable subject to its own weight and a vertical point load. Derivatives up to 10th order were computed with respect to material, loading, and geometrical parameters. The derivatives were used to generate reduced order models of the cable deformation and reaction forces at its ends with respect to multiple input parameters. Results show that from a single hypercomplex evaluation of the cable under a single vertical point load, it is possible to generate an accurate reduced order model capable of predicting the cable deformation with 1.5 times the load in the opposite direction and with 3.5 times the load in the same direction without resolving the system of equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 2. Quantitative Moisture Dynamics Estimation Model

Abstract The long‐term containment of high‐level radioactive waste in geological disposal repositories relies on Engineered Barrier Systems (EBS), with bentonite clay emerging as a candidate material due to its unique properties. Understanding moisture dynamics within bentonite buffers is crucial for EBS performance, as it directly influences the material's swelling capacity, thermal and hydraulic conductivity, mechanical properties, and long‐term evolution under complex thermal‐hydrological‐mechanical (THM) processes. This study develops an advanced Electrical Resistivity Tomography (ERT)‐based framework to quantitatively monitor moisture dynamics under THM conditions. Our framework extends the Waxman‐Smits model to incorporate the coupled effects of temperature, water content, fluid chemistry, and mechanical changes on bentonite's electrical properties. Utilizing HotBENT‐Lab data from our companion paper, which includes electrical conductivity, CT density, and thermocouple measurements, this study offers a novel methodological framework bridging different scales of the model. Our results show that the extended model can estimate water content from ERT data, capturing spatial and temporal variations in moisture distribution within bentonite columns. However, the model tends to overestimate water content compared to CT density‐derived measurements. We address this discrepancy by incorporating a simplified swelling effect model, which improves agreement between ERT and CT density‐based water content estimates. We also discuss model limitations, including simplified treatment of swelling and micropore effects, and propose a conceptual framework for transitioning from laboratory to field applications, addressing challenges such as parameter scalability, field validation methods, and integration of diverse data sources. This ERT‐based framework can potentially advance real‐world moisture monitoring of bentonite‐based EBS in nuclear waste repositories. Plain Language Summary Safely containing high‐level radioactive waste depends on barriers made from materials like bentonite clay, which is effective because it swells and seals in the waste. To ensure these barriers work well over time, it's important to understand how moisture moves through the clay. Our study developed a new method using ERT to monitor moisture levels in bentonite under conditions that mimic those in actual storage sites, including changes in temperature, water content, and mechanical stress. This study improved an existing model to better account for how these factors affect the clay, allowing us to create more accurate moisture maps. Initially, the proposed model overestimated the amount of water in the clay, but its accuracy was improved by factoring in how the clay swells when wet. This study also identified some limitations of the model and suggested ways to adapt it for use in real‐world waste storage sites. This new approach could lead to better monitoring and safety checks for nuclear waste storage systems, helping to ensure long‐term containment. Key Points This work develops an ERT‐based framework extending the Waxman‐Smits model to monitor bentonite moisture dynamics during coupled THM processes The extended model accurately estimates water content from Electrical Resistivity Tomography data, incorporating swelling effects to improve precision This work proposes a conceptual framework for transitioning from laboratory to field applications, advancing EBS monitoring in nuclear waste repositories

Chen, Hang↗
DOE OSTI2021

Machine learning for neutron reflectometry data analysis of two-layer thin films

Neutron reflectometry (NR) is a powerful tool for probing thin films at length scales down to nanometers. We investigated the use of a neural network to predict a two-layer thin film structure to model a given measured reflectivity curve. Application of this neural network to predict a thin film structure revealed that it was accurate and could provide an excellent starting point for traditional fitting methods. Employing prediction-guided fitting has considerable potential for more rapidly producing a result compared to the labor-intensive but commonly-used approach of trial and error searches prior to refinement. A deeper look at the stability of the predictive power of the neural network against statistical fluctuations of measured reflectivity profiles showed that the predictions are stable. We conclude that the approach presented here can provide valuable assistance to users of NR and should be further extended for use in studies of more complex n-layer thin film systems. This result also opens up the possibility of developing adaptive measurement systems in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Impact of Thermal Scattering Law on Similarity Assessment in Light-Water or Polyethylene-Moderated Systems

A collaborative effort between Pacific Northwest National Laboratory (PNNL) and Oak Ridge National Laboratory (ORNL) is underway to provide a technical basis and methodology for the criticality safety community to use the sum of fractions (SoF) method for generating limits for mixtures of “selected actinide nuclides” included in the ANSI/ANS 8.15 standard. The PNNL scope in this project is to define a range of mixtures of 233 U, 235 U, and 239 Pu moderated with either light water or polyethylene and to examine the critical masses for these mixtures. The ORNL scope is primarily to provide validation support for these studies. More complete discussion of the project and its validation aspects will be presented at the upcoming International Conference on Nuclear Criticality Safety (ICNC) this Fall in Sendai, Japan. Clear differences in benchmark similarity to application systems as assessed by the integral parameter ck are noted in the validation studies performed as part of this project as a function of moderator. The c k value is a correlation coefficient that represents that amount of shared uncertainty in k eff due to cross sections between two systems. Individual nuclide-reaction contributions between the two systems can be simply summed to arrive at the total c k value. Specifically, the c k values for light-water–moderated solution experiments are higher for a water-moderated application than for a polyethylene-moderated application. This result is neither totally unexpected nor surprising, but the magnitude of the difference was difficult to anticipate. The TSUNAMI sequence, in the SCALE 6.2.4 code package developed by ORNL, was used to generate eigenvalues and reactivity effects with perturbation-theory based approach through sensitivity coefficients for all nuclides in the system with all reactions and energy groups. The TSUNAMI-Indices and Parameters (IP) sequence then uses the sensitivity data generated through TSUNAMI to generate relational parameters (i.e., c k ) to determine the degree of similarity between systems. One detail of the SCALE material and data implementation must be discussed at this point. Several thermal scattering laws (TSLs) are available for 1 H. SCALE uses a different nuclide ID number for each TSL; essentially, each version of 1 H is treated as a unique nuclide. For example, 1 H bound in water ( 1 H-H2O) is assigned the nuclide ID 1001, whereas 1 H bound in polyethylene (h-poly) is assigned the nuclide ID 9001001. The same cross section data are used for all reactions in 1 H, regardless of TSL, except for scattering below the TSL cutoff energy. TSUNAMI-IP treats different nuclide IDs as different nuclides; thus, no uncertainty is shared between 1 H-H2O and h-poly, despite much of the same data, including covariance data, being used for both nuclides. This presents a question: how much of the difference in assessed similarity between water- and polyethylene-moderated systems is due to the moderators, and how much is caused by the treatment of 1 H-H2O and h poly with cross section and covariance data. The extended edits generated by TSUANMI-IP allow for an investigation of this issue specifically, as well as a demonstration of the general techniques available within TSUNAMI to understand the results of the similarity assessment. This paper presents and analyzes the similarity assessment of both water- and polyethylene-moderated systems for a single benchmark: PST-002-001.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Bulk and Nanocrystalline Cesium Lead-Halide Perovskites as Seen by Halide Magnetic Resonance

Lead-halide perovskites increasingly mesmerize researchers because they exhibit a high degree of structural defects and dynamics yet nonetheless offer an outstanding (opto)electronic performance on par with the best examples of structurally stable and defect-free semiconductors. This highly unusual feature necessitates the adoption of an experimental and theoretical mindset and the reexamination of techniques that may be uniquely suited to understand these materials. Surprisingly, the suite of methods for the structural characterization of these materials does not commonly include nuclear magnetic resonance (NMR) spectroscopy. The present study showcases both the utility and versatility of halide NMR and NQR (nuclear quadrupole resonance) for probing the structure and structural dynamics of CsPbX 3 (X = Cl, Br, I), in both bulk and nanocrystalline forms. The strong quadrupole couplings, which originate from the interaction between the large quadrupole moments of, e.g., the 35 Cl, 79 Br, and 127 I nuclei, and the local electric-field gradients, are highly sensitive to subtle structural variations, both static and dynamic. The quadrupole interaction can resolve structural changes with accuracies commensurate with synchrotron X-ray diffraction and scattering. It is shown that space-averaged site-disorder is greatly enhanced in the nanocrystals compared to the bulk, while the dynamics of nuclear spin relaxation indicates enhanced structural dynamics in the nanocrystals. The findings from NMR and NQR were corroborated by ab initio molecular dynamics, which point to the role of the surface in causing the radial strain distribution and disorder. These findings showcase a great synergy between solid-state NMR or NQR and molecular dynamics simulations in shedding light on the structure of soft lead-halide semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Completing the Ba–As Compositional Space: Synthesis and Characterization of Three New Binary Zintl Arsenides, Ba3As4, Ba5As4, and Ba16As11

Three novel binary barium arsenides, Ba3As4, Ba5As4, and Ba16As11, were synthesized and their crystal and electronic structures were investigated. Structural data collected via the single-crystal X-ray diffraction method indicate that the anionic substructures of all three novel compounds are composed of structural motifs based on the homoatomic As–As contacts, with [As2]4− dimers found in Ba5As4 and Ba16As11, and an [As4]6− tetramer found in Ba3As4. Ba3As4 and Ba5As4 crystallize in the orthorhombic crystal system—with the non-centrosymmetric space group Fdd2 (a = 15.3680(20) Å, b = 18.7550(30) Å, c = 6.2816(10) Å) for the former, and the centrosymmetric space group Cmce (a = 16.8820(30) Å, b = 8.5391(16) Å, and c = 8.6127(16) Å) for the latter—adopting Eu3As4 and Eu5As4 structure types, respectively. The heavily disordered Ba16As11 structure was solved in the tetragonal crystal system with the space group P4¯21m (a = 12.8944(12) Å and c = 11.8141(17) Å). The Zintl concept can be applied to each of these materials as follows: Ba3As4 = (Ba2+)3[As4]6−, Ba5As4 = (Ba2+)5(As3−)2[As2]4−, and 2 × Ba16As11 = (Ba2+)32(As3−) ≈ 20[As2]4− ≈ 1, pointing to the charge-balanced nature of these compounds. Electronic structure calculations indicate narrow bandgap semiconducting behavior, with calculated bandgaps of 0.47 eV for Ba3As4, 0.34 eV for Ba5As4, and 0.33 eV for Ba16As11.

Crystallography↗
DOE OSTI2022

Sustainable Aviation Fuel from Hydrothermal Liquefaction of Wet Wastes

Hydrothermal liquefaction (HTL) uses heat and pressure to liquefy the organic matter in biomass/waste feedstocks to produce biocrude. When hydrotreated the biocrude is converted into transportation fuels including sustainable aviation fuel (SAF). Further, by liquifying the organic matter in wet wastes such as sewage sludge, manure, and food waste, HTL can prevent landfilling or other disposal methods such as anerobic digestion, or incineration. A significant roadblock to the development of a new route for SAF is the strict approval process, and the large volumes required (>400 L) for testing. Tier α and β testing can predict some of the properties required for ASTM testing with <400 mL samples. The current study is the first to investigate the potential for utilizing wet-waste HTL biocrude (WWHTLB) as an SAF feedstock. Herein, several WWHTLB samples were produced from food waste, sewage sludge, and fats, oils, and grease, and subsequently hydrotreated and distilled to produce SAF samples. The fuels (both undistilled and distilled samples) were analyzed via elemental and 2D-GC-MS. Herein, we report the Tier α and β analysis of an SAF sample derived originally from a WWHTLB. The results of this work indicate that the upgraded WWHTLB material exhibits key fuel properties, including carbon number distribution, distillation profile, surface tension, density, viscosity, heat of combustion, and flash point, which all fall within the required range for aviation fuel. WWHTLB has therefore been shown to be a promising candidate feedstock for the production of SAF.

09 BIOMASS FUELS↗
DOE OSTI2024

Development of a Colorimetric Loop-Mediated Isothermal Amplification Assay for the Detection of Trypanosoma cruzi in Low-Resource Settings

Chagas disease is an inflammatory parasitic infection caused by Trypanosoma cruzi (T. cruzi). Early diagnosis is crucial in guiding treatment and slowing disease progression; however, current diagnostic methods have insufficient detection limits and often require skilled technicians. Molecular tests, especially isothermal nucleic acid assays, are advantageous due to their excellent sensitivity, specificity, speed, and simplicity. Here, we optimized a colorimetric loop-mediated isothermal amplification (LAMP) assay for T. cruzi. We can detect as few as 2 genomic copies/reaction using three different T. cruzi strains. We examined selectivity using other parasitic protozoans and successfully detected T. cruzi DNA extracted from parasites in human whole blood down to 1.2 parasite equivalents/reaction. We also performed a blinded study using canine blood samples and established a 100% sensitivity, specificity, and accuracy for the colorimetric LAMP assay. Finally, we used a heated 3D printer bed and an insulated thermos cup to demonstrate that the LAMP incubation step could be performed with accessible, low-cost materials. Altogether, we have developed a high-performing assay for T. cruzi with a simple colorimetric output that would be ideal for rapid, low-cost screening at the point of use.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Application of the variational autoencoder to detect the critical points of the anisotropic Ising model

We generalize the previous study on the application of variational autoencoders to the two-dimensional Ising model to a system with anisotropy. Due to the self-duality property of the system, the critical points can be located exactly for the entire range of anisotropic coupling. This presents an excellent test bed for the validity of using a variational autoencoder to characterize an anisotropic classical model. Furthermore, we reproduce the phase diagram for a wide range of anisotropic couplings and temperatures via a variational autoencoder without the explicit construction of an order parameter. Considering that the partition function of ($d$ + 1)-dimensional anisotropic models can be mapped to that of the $d$-dimensional quantum spin models, the present study provides numerical evidence that a variational autoencoder can be applied to analyze quantum systems via the quantum Monte Carlo method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Azimuthal Anchoring Strength in Photopatterned Alignment of a Nematic

Spatially-varying director fields have become an important part of research and development in liquid crystals. Characterization of the anchoring strength associated with a spatially-varying director is difficult, since the methods developed for a uniform alignment are seldom applicable. Here we characterize the strength of azimuthal surface anchoring produced by the photoalignment technique based on plasmonic metamsaks. The measurements used photopatterned arrays of topological point defects of strength +1 and −1 in thin layers of a nematic liquid crystal. The integer-strength defects split into pairs of half-integer defects with lower elastic energy. The separation distance between the split pair is limited by the azimuthal surface anchoring, which allows one to determine the strength of the latter. The strength of the azimuthal anchoring is proportional to the UV exposure time during the photoalignment of the azobenzene layer.

36 MATERIALS SCIENCE↗
DOE OSTI2020

GIS-Based Modeling of Contaminated Soil Volumes at Multiple Sites in the Formerly Utilized Sites Remedial Action Program - 20149

The remediation of hazardous, toxic, and radioactive waste (HTRW) sites produces cost-related risks associated with the estimation of contaminated soil or debris volumes. Historical risk-management techniques include cost contingencies to cover volume uncertainties that affect project budgeting and decision-making. The Buffalo District teamed with project partners to lessen volume uncertainty and reduce project risks at multiple HTRW sites managed under the Formerly Utilized Sites Remedial Action Program (FUSRAP). Historical remedial investigations under FUSRAP commonly identified the presence of radiological material in site media, the associated human health risk, and then areas of remediation. To manage remedial execution and reduce risk, pre-design or remediation-phase sampling essentially 'chased' contamination, which was not conducive to efficient predictive budgeting derived from Feasibility Study (FS) cost analyses. The Buffalo District first optimized their approach to better understand volume uncertainty by utilizing the Argonne National Laboratory's Bayesian Approaches for Adaptive Spatial Sampling (BAASS) software [1]. BAASS processed soft data (e.g., gamma walk-over data) and spatial sampling data to estimate the lateral extent of contaminated soil irrespective of depth (i.e., gross contamination extent) and define areas of contaminant uncertainty. The software performed a binary transformation of contaminant concentrations at all sampling points based upon remedial action goals or a sum of ratios approach (i.e., clean, impacted, or range of impacts in soil). The model produced two-dimensional (horizontal) contaminant probability contours and statistical uncertainty in the sampling coverage and resulting contaminant extents. This method was translated vertically by partitioning the sampling data into depth brackets that produced a stacked representation of contaminant extents and uncertainty in the subsurface (i.e., similar to construction lifts). The results commonly led to a better understanding of project uncertainty and the need for sampling strategies that produce high-confidence soil volumes, which control costs. The BAASS-based delineations were eventually replaced by Empirical Bayesian Kriging (EBK) methods available in ArcGIS Spatial or 3D Analysts [2]. The EBK method calculates contaminant probability zones derived from user-controlled semivariograms of the spatial datasets. The resulting probability zones (e.g., 50% or 80% of contaminant probability) represent the two-dimensional surface delineation of the overall horizontal remedial area, similarly to BAASS. However, unlike BAASS, the vertical sampling data within these probability zones became vertical control points to contour a subterranean surface that connects subsurface points to the land-surface delineations of contamination. The resulting representation of horizontal and vertical impacts within an enclosed envelop (volume) of soil included uncertainty distributions that are used to plan uncertainty-reduction sampling. These data-driven and math-based models of three-dimensional sampling results produced well-bounded remedial volumes for project planning and better uncertainty predictions during project budgeting. The EBK method was applied to several FUSRAP sites managed by the Buffalo District and compared to less rigorously modeled sites previously remediated by the District. The comparison of modeled to actual remediated volumes provide a basis for validating the volume-estimation method. This comparison is important to ensure modeled volumes match physical boundaries of site remediation. FUSRAP sites with denser investigative sampling and lesser volume uncertainty proved useful in remedial planning and contracting. The Buffalo District noted that historical sites with sparser sampling arrays had greater disparity between estimated volumes and final remedial volumes. The benefit achieved over the cost of detailed soil sampling appears positive for FUSRAP projects, especially where impacts vary widely and appear unbounded by investigation-phase sampling. The subsequent Empirical Bayesian Kriging of contamination coupled with vertical contouring for soil estimations reduces uncertainty in soil volumes or indicates where sampling is required to reduce uncertainty, which together optimize remedial planning and budgeting. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Reproducibility in G 0 W 0 calculations for solids

Ab initio many-body perturbation theory within the GW approximation is a Green's function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn–Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange–correlation functional “starting point”, are used to construct G and W, and then perturbatively corrected by the resultant GW self-energy. In practice, there are several ways to construct the GW self-energy, and these can lead to variations in predicted quasiparticle energies. For example, for ZnO and TiO 2 , the GW fundamental gaps reported in the literature can vary by more than 1 eV depending on the GW code used. In this work, we calculate and analyze GW quasiparticle (QP) energies of these and other systems with three different GW codes: BERKELEYGW, ABINIT and YAMBO. Through a systematic analysis of the GW implementation of these three codes, we identify the primary origin of major discrepancies between codes reported in prior literature to be the different implementations the Coulomb divergence in the Fock exchange term and the frequency integration scheme of the GW self-energy. We then eliminate these discrepancies by using common numerical methods and algorithms, demonstrating that the same quasiparticle energies for a given material can be obtained with different codes, within numerical differences ascribable to the technical details of the underling implementations. This work will be important for users and developers in assessing the precision of future GW applications and methods.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Li interaction-induced phase transition from black to blue phosphorene

A comprehensive first-principle calculation has been carried out and revealed that at sufficiently high Li concentration and certain well-defined configurations a phase transition from black to blue phosphorene can take place. Blue phosphorene, a newly predicted allotrope of phosphorus, possesses unique crystalline and electronic structure and is a promising candidate, not only for fundamental research but also for electronic and optoelectronic applications. Methods to growth high quality blue phosphorene layers are highly desirable but challenging. Here, a novel kinetic pathway to grow blue phosphorene layers from black phosphorene layers via Li intercalation is proposed based on first principle study. This study pointed out that Li atoms intercalated in black phosphorene could act as ‘catalyst’ in the ‘reactive region’ of the lone pair of P atoms, leading to a P-P bond breaking and subsequently, a local structural transformation from orthorhombic lattice to an assembly of parallel narrow nanoribbons with rhombohedra-like symmetry. During Li deintercalation, these nanoribbons are self-mended and form blue phosphorene layers. The interlayer distance was found 4.60 Å for double layer with AA stacking and 4.13 Å for multilayer with ABC stacking, respectively, indicating a monolayer blue phosphorene can be mechanically exfoliated. Furthermore, this study also points out the possibility of new phases in other systems, where intercalation can lead to an unexpected structural phase transition and even a discovery of novel materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Steel Equivalency Workbook: An X-Ray Transmission Calculator

The Steel Equivalency Spreadsheet was created as a robust, user-friendly method for the following: 1. Determining if it is possible to image an object. 2. Reducing time needed when selecting equipment and preparing for field imaging activities. 3. Reducing the amount of equipment taken for field imaging activities. 4. Determining a starting point for exposure settings prior to imaging an object. 5. Reducing the amount of dose deposited to complete an imaging operation (in the spirit of As Low As Reasonably Achievable, ALARA). The spreadsheet was developed using various physical models to account for different phenomena. Future implementations aim to expand beyond the Digital Radiography and Computed Tomography Single Munition Scanner (DRCT SMS) in standard configuration to include high energy XRGs (Betatrons) and sub-MeV pulsed XRGs (XRS4). Alternative computational methods that compensate for incoherent scattering when calculating relative transmission are also being pursued. Ultimately the spreadsheet exceeded the developmental goal of having less than 10% average error when comparing calculations to real-world images.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Classifying photonic topology using the spectral localizer and numerical K -theory

Recently, the spectral localizer framework has emerged as an efficient approach for classifying topology in photonic systems featuring local nonlinearities and radiative environments. In nonlinear systems, this framework provides rigorous definitions for concepts such as topological solitons and topological dynamics, where a system’s occupation induces a local change in its topology due to nonlinearity. For systems embedded in radiative environments that do not possess a shared bulk spectral gap, this framework enables the identification of local topology and shows that local topological protection is preserved despite the lack of a common gap. However, as the spectral localizer framework is rooted in the mathematics of C*-algebras, and not vector bundles, understanding and using this framework requires developing intuition for a somewhat different set of underlying concepts than those that appear in traditional approaches for classifying material topology. In this tutorial, we introduce the spectral localizer framework from a ground-up perspective and provide physically motivated arguments for understanding its local topological markers and associated local measure of topological protection. In doing so, we provide numerous examples of the framework’s application to a variety of topological classes, including crystalline and higher-order topology. We then show how Maxwell’s equations can be reformulated to be compatible with the spectral localizer framework, including the possibility of radiative boundary conditions. To aid in this introduction, we also provide a physics-oriented introduction to multi-operator pseudospectral methods and numerical K-theory, two mathematical concepts that form the foundation for the spectral localizer framework. Finally, we provide some mathematically oriented comments on the C*-algebraic origins of this framework, including a discussion of real C*-algebras and graded C*-algebras that are necessary for incorporating physical symmetries. Looking forward, we hope that this tutorial will serve as an approachable starting point for learning the foundations of the spectral localizer framework.

97 MATHEMATICS AND COMPUTING↗