Results for “Core processes” — Page 87

DOE OSTI2020

Synthesis and Characterization of Pd-based Nanomaterials

Bimetallic nanoparticles (BNPs) consist of two different types of metals or alloys that are bonded together. Unique properties such as optical, electronic, thermal, and catalytic effects differ for each type of BNP. Important BNPs range from Au-Pd, Ag- Pt, Au-Pt, and Ag-Ni. Pd bimetallic nanoparticles are of interest due to their many applications such as catalysis and sensing. Bimetallic catalysts have increase reaction rates and have improved catalyst stability through the geometry and ligand distribution. Pd nanoparticles are considered to be a strong catalyst due to their high activity at low temperatures and high tolerance to moisture. The catalytic properties of bimetallic nanoparticles depend on the structural properties such as size and shape. Core-shell, hollow structure, and multi-shell alloy are three possible structures nanoparticles can form as bimetallic catalysts. BNPs can be synthesized through different methods to control the size, shape, and structure. To obtain different morphologies, a variety of methods can be performed. Different methods can range from the usage of the glancing angle deposition (GLAD) to the galvanic replacement reaction, but the methods all depend on the properties of the metals. The galvanic displacement reaction was the method used to obtain Pd-based nanoparticles. This reaction is best know for obtaining hollow shaped NPs. To determine what redox process was preformed, the activity series of metals was used. From the activity series of metals, silver (Ag) was selected to preform Pd-based nanoparticles. Objectives: Synthesize Ag nanoparticles and Ag-Pd nanoparticles to understand the morphology. Characterize the synthesized nanoparticles using scanning electron microscopy (SEM), phase analysis light scattering (PALS), dynamic light scattering (DLS), energy dispersive X-ray spectroscopy (EDS), and UV-Vis spectroscopy. Results: In the UV-Vis spectrum, the Ag-Pd bimetallic NP's plasmon band decreased as the volume of palladium increased. The surface charge increases as the concentration of palladium increases. The Pd{sup 2+} ions interact with the sodium citrate surface, and decrease the negative charge. Conclusion: Ag-Pd nanoparticles were successfully created and stabilized with sodium citrate. The addition of Pd decreased the prominent plasmon band of the Ag nanoparticles. The SEM analysis showed that Ag nanoparticles had a well-defined structure, while the Ag-Pd nanoparticles showed hollow and rough structure. The EDX analysis confirmed the presence of silver and palladium. This material can be used in many industrial and research fields such as organic synthesis, fuel cells, and environmental sensing and remediation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

SPC-70808 Rev 0 AL-30 Secondary Support Standoff Blocks Build to Print Fabrication for Use in MARVEL

A. Idaho National Laboratory (INL) is developing a microreactor to produce electrical power utilizing a small nuclear core and Stirling engines under the Microreactor Application Research Validation and Evaluation (MARVEL) program. This procurement specification defines the requirements for fabrication of the MARVEL AL-30 Secondary Support Standoff Blocks manufactured by Zircar. The design requires twelve blocks; however, the requested quantity shall be 16 to account for required testing after receipt of materials. B. Any conflict between this specification and referenced Codes and Standards, or any supplementary specifications in the procurement documents requires written clarification from the Contractor prior to proceeding with any work. Any deviation from the procurement documents requires approval by the Contractor with the change request process.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

DOE OSTI2025

SPC-70804 Rev 1 Zircar ZAL-45 Alumina Standoff Plates for Use in MARVEL

Idaho National Laboratory (INL) is developing a microreactor to produce electrical power utilizing a small nuclear core and Stirling engines under the Microreactor Application Research Validation and Evaluation (MARVEL) program. This procurement specification defines the requirements for fabrication of the MARVEL ZAL-45 standoff plates manufactured by Zircar. The design requires eight (8) ceramic standoff plates; however, the requested quantity shall be twelve to account for required testing after receipt of materials. Any conflict between this specification and referenced Codes and Standards, or any supplementary specifications in the procurement documents requires written clarification from the Contractor prior to proceeding with any work. Any deviation from the procurement documents requires approval by the Contractor with the change request process.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

DOE OSTI2022

Multiple ferroic orders and toroidal magnetoelectricity in the chiral magnet BaCoSiO 4

Discovering ferroic phase transitions and their consequential physical properties is at the core of condensed matter science due to rich physics and tremendous technological promises. BaCoSiO 4 , a chiral antiferromagnet, belongs to the tetrahedron-based chiral system, and exhibits diverse ferroic orders with coexisting chirality, polarity, trimerization, ferrorotational distortions, and magnetism. However, their mutual couplings remain to be explored. Here in this work, we used a comprehensive combination of several experimental tools—in situ x-ray, transmission electron microscopy, magnetization, and magnetoelectric measurements of single-crystalline BaCoSiO 4 —to investigate hierarchical phase transitions, their microscopic domain structures, and the resulting magnetoelectricity. We found that two different structural chiralities develop through distinct processes: global homochirality and local heterochirality induced by the ferrorotational distortions on top of existing polarization. In addition, magnetic chirality, with the simultaneous presence of net magnetic moment and magnetic toroidal moment, develops below 3.2 K due to the global chirality, which leads to magnetic field tunable toroidal magnetoelectricity. Thus, BaCoSiO 4 exhibits uniquely all four types of ferroic orders and provides an avenue to explore, for example, tunable or dynamic coupling of multiple ferroic degrees of freedom.

36 MATERIALS SCIENCE↗

DOE OSTI2026

PRISMA: PARALLEL REFINEMENT AND INTEGRATION SYSTEM FOR MULTI-AZIMUTHAL ANALYSIS

The Parallel Refinement and Integration System for Multi-azimuthal Analysis (PRISMA, version 1.1.0) is a Python application for processing X-ray diffraction (XRD) image data. PRISMA wraps GSAS-II to perform azimuthally-binned peak refinement, computes per-frame strain and d-spacing from those fits, and provides three PyQt5 graphical interfaces: (1) a Recipe Builder for selecting GSAS-II control (.imctrl) files, optional mask (.immask) files or threshold-ased masking, reference and experiment image sets, peaks, zimuthal range and bin size, and an optional ceria-based auto-calibration; (2) a Batch Processor that uses Dask on local workstations and pure MPI (mpi4py.futures.MPICommExecutor) on HPC to distribute GSAS-II refinement across cores or compute nodes and write results to a 4-dimensional (peaks x frames x azimuths x measurements) Zarr dataset; and (3) a Data Analyzer that renders heatmaps of fit parameters, strain, frame-to-frame deltas, and percent-change-vs-reference, and exports user-defined subsections to CSV or Excel. The peak-refinement algorithm is deterministic. Benchmark on ALCF Crux: a 20,000-image set, single-peak fit in frame mode with 44 azimuthal bins on 128 nodes x 128 workers, 48 seconds total wall time.

Lorenzo Martin, Maria De La Cinta [Argonne Nationa↗

DOE OSTI2019

Enabling 10mol/kg Swing Capacity via Heat Integrated Sub-ambient Pressure Swing Adsorption

This project focused on the development of a post-combustion CO 2 capture technology for coal-fired power stations. The core concept of the technology was the combination of an adsorption-based separation with combustion flue gas cooling, compression followed by heat and power recovery. The key features of the adsorption unit (which captures and purifies the CO 2 ) are the use of high flow fiber-based adsorption systems that incorporate materials that manage the significant amounts of heat released during the adsorption separation. The adsorbent material utilized was a hybrid organic-inorganic material (a metal-organic framework) that had the proper combination of separation selectivity and CO 2 capture capacity. In this case, record-setting, process-relevant CO 2 capture capacities (approximately 8-10 mol of CO 2 per kg of adsorbent “swing capacities”) were observed in experimental measurements of the adsorptive separation. The experimental work was accompanied by detailed multi-level simulation works, which investigated the molecular mechanism for the high CO 2 capture efficiencies, the heat and mass transport dynamics within the separation system, and the process-level integration of the various systems associated with the technology. Overall, the analysis conducted in this work indicates that this technology can potentially reach CO 2 capture costs as low as 40-50 dollars/tonne with capital costs of 250M-400M dollars for the 550 MW e basecase established by the National Energy Technology Laboratory.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2020

Verification and Validation of the STAT7 Code

This document presents the verification and validation of the STAT7 Version 1.0 code which has been developed to perform the steady-state thermal hydraulic safety calculation of the MIT Research Reactor required upon conversion of the reactor from high enrichment uranium fuel to low enrichment uranium fuel. The key capabilites of the code including the statistical processing, water property generation, and thermal hydraulic solution are verified independently. While the first two capabilites are checked using statistical tests and NIST data tables, each component of the thermal hydraulic solution is verified using simple geometry test problems. The comparison between results obtained by the code and hand calculations on these problems shows that each of these components work correctly. Additionally, a comprehensive comparison between STAT7 Version 1.0 and PLTEMP/ANL Version 4.2 calculations for a five-channel test problem is performed for the validation. The good agreement between the obtained results for all thermal hydraulic parameters belonging to the fuel core region is achieved. An issue in bypass flowrate calculation was detected when both fin-friction and bypass flow claculations are involved. This problem can be corrected in a later version of the code, however, does not affect the application of STAT7 code as it has been used to date where the flow distribution in the core region and bypass channel is specified as an input parameter.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Verification and Validation of the STAT7 Code

This document presents the verification and validation of the STAT7 Version 1.0 code which has been developed to perform the steady-state thermal hydraulic safety calculation of the MIT Research Reactor required upon conversion of the reactor from highly enriched uranium fuel to low- enriched uranium fuel. The key capabilities of the code including the statistical processing, water property generation, and thermal hydraulic solution are verified independently. While the first two capabilities are checked using statistical tests and NIST data tables, each component of the thermal hydraulic solution is verified using simple geometry test problems. The comparison between results obtained by the code and hand calculations on these problems shows that each of these components work correctly. Additionally, a comprehensive comparison between STAT7 Version 1.0 and PLTEMP/ANL Version 4.2 calculations for a five-channel test problem is performed for the validation. The good agreement between the obtained results for all thermal hydraulic parameters belonging to the fuel core region is achieved. An issue in bypass flowrate calculation was detected when both fin-friction and bypass flow calculations are involved. This problem can be corrected in a later version of the code. However, it does not affect the application of STAT7 as it has been used to date where the flow distribution in the core region and bypass channel is specified as an input parameter.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Dimensionality effects on multicomponent ionic transport and surface complexation in porous media

Coulombic interactions between charged species in pore water and at surface/solution interfaces are of pivotal importance for multicomponent reactive transport in porous media. In this study, we investigate the impact of domain dimensionality on electrostatically coupled dispersion and surface-solution reactions during transport of acidic plumes and major ions in porous media. Column and quasi two-dimensional flow-through experiments were performed, with identical silica porous media and under the same advection-dominated conditions. Equal mass fluxes of different electrolyte solutions (i.e., HCl - pH ~ 2.8, NaBr - 100 mM, HCl - pH ~ 2.8 plus NaBr - 100 mM) were continuously injected in the 1-D and 2-D experiments and breakthrough curves of pH and major ions were measured at the outlet of the domains. The presence of pronounced ionic strength gradients in the transverse direction in the 2-D setup caused distinct retardation and transport behaviors of protons and major ions which were not observed in the one-dimensional column experiments. Furthermore, in the cases of salt electrolytes injection, considerably enhanced release of H + (>61%) from the quartz surface was observed in the multidimensional system compared to the one-dimensional setup. Reactive transport modeling was performed to reproduce the experimental outcomes and to analyse the coupling between transport processes, based on the Nernst-Planck formulation of diffusive/dispersive fluxes and on surface complexation reactions at the solid-solution interface. Electrostatic interactions between Na + , B – , and H + , and deprotonation of the quartz surface upon the formation of sodium outer-sphere complexes, are the primary controllers of the spatial and temporal features displayed by the pH and major ions measurements. The reactive transport simulations allowed us to interpret the experimental observations, to visualize the distribution and spatio-temporal evolution of dissolved and solid species, to identify a spatially heterogeneous zonation of Coulombic interactions with distinct behavior at the fringe and core of the injected plumes in the multidimensional setup, and to quantify the different components of the Nerst-Planck fluxes of the charged solutes. This study demonstrates that the domain dimensionality directly affects electrostatic interactions between charged aqueous species in the pore water and surface complexation reactions at the solid-solution interface. The non-trivial effects of dimensionality on multicomponent ionic transport result in a significantly different behavior in 1-D and 2-D systems.

58 GEOSCIENCES↗

DOE OSTI2023

Structure of the antenna complex expressed during far-red light photoacclimation in Synechococcus sp. PCC 7335

Far-red light photoacclimation, or FaRLiP, is a facultative response exhibited by some cyanobacteria that allows them to absorb and utilize lower energy light (700–800 nm) than the wavelengths typically used for oxygenic photosynthesis (400–700 nm). During this process, three essential components of the photosynthetic apparatus are altered: photosystem I, photosystem II, and the phycobilisome. In all three cases, at least some of the chromophores found in these pigment–protein complexes are replaced by chromophores that have red-shifted absorbance relative to the analogous complexes produced in visible light. Recent structural and spectroscopic studies have elucidated important features of the two photosystems when altered to absorb and utilize far-red light, but much less is understood about the modified phycobiliproteins made during FaRLiP. We used single-particle, cryo-EM to determine the molecular structure of a phycobiliprotein core complex comprising allophycocyanin variants that absorb far-red light during FaRLiP in the marine cyanobacterium Synechococcus sp. PCC 7335. The structure reveals the arrangement of the numerous red-shifted allophycocyanin variants and the probable locations of the chromophores that serve as the terminal emitters in this complex. It also suggests how energy is transferred to the photosystem II complexes produced during FaRLiP. The structure additionally allows comparisons with other previously studied allophycocyanins to gain insights into how phycocyanobilin chromophores can be tuned to absorb far-red light. These studies provide new insights into how far-red light is harvested and utilized during FaRLiP, a widespread cyanobacterial photoacclimation mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment

This study presents a framework to quantify changes in soil saturation in response to flooding caused by extreme hydrologic perturbation on coastal ecosystems at the interfaces and transition between terrestrial and aquatic systems. Subsurface heterogeneity limits the use of in situ measurements to quantify subsurface flow during flooding due to the spatial discontinuity in the measured data. While geophysical methods, including time-lapse electrical resistivity imaging (ERI), are increasingly used to monitor soil hydrological processes, their abilities to parameterize flow models have been underutilized. This study combines background ERI, ground penetrating radar (GPR), time-lapse ERI, soil characterization, and a numerical flow model developed using an Advanced Terrestrial Simulator (ATS) code to quantify the infiltration pathway and describe the hydrological dynamics during a simulated flooding experiment. We assessed the use of two conceptual models developed using [1] ERI and GPR data that described the stratigraphic distribution, and time-lapse ERI that mapped permeability contrast, and [2] information from a national soil database for capturing changes in saturation. Combining the ERI and GPR results with soil core data revealed the stratigraphic heterogeneity at the site with a silty clay layer from 1 to 2 m between an overlying loamy topsoil and an underlying saturated silty sand. This silty clay layer could restrict deep infiltration. The time-lapse ERI showed up to a 35% decrease in resistivity, which correlated with soil moisture data (R 2 value > 0.53) and revealed preferential infiltration zones used to inform the flow model. Numerical simulation results from both the geophysics- and soil database-informed models quantified changes in soil saturation with calculated soil moistures that agreed with field data. The geophysics-informed model captured more of the system’s variability, reflective of shallow subsurface heterogeneities. The framework presented will serve as a precursor for a robust ecohydrological model that can describe the impacts of extreme events induced by climate change on coastal ecosystems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

A formal Fe III/V redox couple in an intercalation electrode

Iron redox cycling between low-valent oxidation states of Fe II and Fe III drives crucial processes in nature. The Fe II/III redox couple charge compensates the cycling of lithium iron phosphate, a positive electrode (cathode) for lithium-ion batteries. High-valent iron redox couples, involving formal oxidation higher than Fe III , could deliver higher electrochemical potentials and energy densities. However, because of the instability of high-valent Fe electrodes, they have proven difficult to probe and exploit in intercalation systems. Here, in this study, we report and characterize a formal Fe III/V redox couple by revisiting the charge compensation mechanism of (de)lithiation in Li 4 FeSbO 6 . Valence-sensitive experimental and computational core-level spectroscopy reveal a direct transition from Fe III (3d 5 ) to a negative-charge-transfer FeV (3d 5 L 2 ) ground state on delithiation, without forming Fe IV , or oxygen dimers. We identify that the cation ordering in Li 4 FeSbO 6 drives a templated phase transition to stabilize the unique Fe V species and demonstrate that disrupting cation ordering suppresses the Fe III/V redox couple. Exhibiting resistance to calendar aging, high operating potential and low voltage hysteresis, the Fe III/V redox couple in Li 4 FeSbO 6 provides a framework for developing sustainable, Fe-based intercalation cathodes for high-voltage applications.

Ramachandran, Hari [Stanford Univ., CA (United Sta↗

DOE OSTI2022

High Flux Isotope Reactor Low Enriched Uranium U-10Mo Fuel Design Parameters

Activities to convert the HFIR from HEU to LEU are ongoing as part of the US Department of Energy (DOE) National Nuclear Security Administration (NNSA) nuclear nonproliferation mission. Design activities to study the conversion of HFIR from HEU to LEU fuel explored different fuel design features and shapes with a uranium-molybdenum (U-10Mo) monolithic alloy fuel. This high-density alloy contains 90 wt % uranium and 10 wt % molybdenum and has a uranium density of 15.318gU/cm 3 . The goal of these studies is to generate several candidate HFIR LEU fuel designs of varying fuel fabrication complexity that meet the current HEU performance metrics and safety requirements. Recent advancements in modeling and simulation tools and design methods enabled a thorough analysis of the available design space with U-10Mo fuel. A surrogate model used this analysis as training data to quickly determine the performance of a design given specific design parameters. An optimization module used this surrogate model to quickly search this multidimensional search space given specific desired performance characteristics. This approach was made possible by the large available design space with U-10Mo fuel. Shift, a Monte Carlo tool optimized for high-performance computing (HPC) architectures, was used for faster calculation and better data management for reactor physics simulations. Once most of these design studies were complete, a new suite called the Python HFIR Analysis and Measurement Engine (PHAME) was developed to connect all fuel design analysis steps, making design studies more efficient and reproducible. The post-processing capabilities of these new tools are leveraged for the information provided herein. Leveraging these tools, several candidate fuel designs were selected with varying levels of feature complexity and reactor performance. This report provides design feature details for four selected HFIR LEU U-10Mo fuel designs and their corresponding performance and safety metrics. Nominal best-estimate design parameters and irradiation conditions, including fission rate densities, power densities, heat fluxes, and cumulative fission densities, are provided. Simulations show that the high uranium density of U-10Mo fuel provides a large potential design space that enables various LEU designs to meet HEU core performance metrics and safety requirements with a power increase from 85 MW (HEU) to 95 MW or 100 MW (LEU).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

UPC++ v1.0 Programmer’s Guide (Rev. 2023.9.0)

UPC++ is a C++ library that supports Partitioned Global Address Space (PGAS) programming. It is designed for writing efficient, scalable parallel programs on distributed-memory parallel computers. The key communication facilities in UPC++ are one-sided Remote Memory Access (RMA) and Remote Procedure Call (RPC). The UPC++ control model is single program, multiple-data (SPMD), with each separate constituent process having access to local memory as it would in C++. The PGAS memory model additionally provides one-sided RMA communication to a global address space, which is allocated in shared segments that are distributed over the processes. UPC++ also features Remote Procedure Call (RPC) communication, making it easy to move computation to operate on data that resides on remote processes. UPC++ was designed to support exascale high-performance computing, and the library interfaces and implementation are focused on maximizing scalability. In UPC++, all communication operations are syntactically explicit, which encourages programmers to consider the costs associated with communication and data movement. Moreover, all communication operations are asynchronous by default, encouraging programmers to seek opportunities for overlapping communication latencies with other useful work. UPC++ provides expressive and composable abstractions designed for efficiently managing aggressive use of asynchrony in programs. Together, these design principles are intended to enable programmers to write applications using UPC++ that perform well even on hundreds of thousands of cores.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

UPC++ v1.0 Programmer’s Guide (Revision 2022.3.0)

UPC++ is a C++ library that supports Partitioned Global Address Space (PGAS) programming. It is designed for writing efficient, scalable parallel programs on distributed-memory parallel computers. The key communication facilities in UPC++ are one-sided Remote Memory Access (RMA) and Remote Procedure Call (RPC). The UPC++ control model is single program, multiple-data (SPMD), with each separate constituent process having access to local memory as it would in C++. The PGAS memory model additionally provides one-sided RMA communication to a global address space, which is allocated in shared segments that are distributed over the processes. UPC++ also features Remote Procedure Call (RPC) communication, making it easy to move computation to operate on data that resides on remote processes. UPC++ was designed to support exascale high-performance computing, and the library interfaces and implementation are focused on maximizing scalability. In UPC++, all communication operations are syntactically explicit, which encourages programmers to consider the costs associated with communication and data movement. Moreover, all communication operations are asynchronous by default, encouraging programmers to seek opportunities for overlapping communication latencies with other useful work. UPC++ provides expressive and composable abstractions designed for efficiently managing aggressive use of asynchrony in programs. Together, these design principles are intended to enable programmers to write applications using UPC++ that perform well even on hundreds of thousands of cores.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

UPC++ v1.0 Programmer’s Guide, Revision 2023.3.0

UPC++ is a C++ library that supports Partitioned Global Address Space (PGAS) programming. It is designed for writing efficient, scalable parallel programs on distributed-memory parallel computers. The key communication facilities in UPC++ are one-sided Remote Memory Access (RMA) and Remote Procedure Call (RPC). The UPC++ control model is single program, multiple-data (SPMD), with each separate constituent process having access to local memory as it would in C++. The PGAS memory model additionally provides one-sided RMA communication to a global address space, which is allocated in shared segments that are distributed over the processes. UPC++ also features Remote Procedure Call (RPC) communication, making it easy to move computation to operate on data that resides on remote processes. UPC++ was designed to support exascale high-performance computing, and the library interfaces and implementation are focused on maximizing scalability. In UPC++, all communication operations are syntactically explicit, which encourages programmers to consider the costs associated with communication and data movement. Moreover, all communication operations are asynchronous by default, encouraging programmers to seek opportunities for overlapping communication latencies with other useful work. UPC++ provides expressive and composable abstractions designed for efficiently managing aggressive use of asynchrony in programs. Together, these design principles are intended to enable programmers to write applications using UPC++ that perform well even on hundreds of thousands of cores.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

UPC++ v1.0 Programmer’s Guide, Revision 2022.9.0

UPC++ is a C++ library that supports Partitioned Global Address Space (PGAS) programming. It is designed for writing efficient, scalable parallel programs on distributed-memory parallel computers. The key communication facilities in UPC++ are one-sided Remote Memory Access (RMA) and Remote Procedure Call (RPC). The UPC++ control model is single program, multiple-data (SPMD), with each separate constituent process having access to local memory as it would in C++. The PGAS memory model additionally provides one-sided RMA communication to a global address space, which is allocated in shared segments that are distributed over the processes. UPC++ also features Remote Procedure Call (RPC) communication, making it easy to move computation to operate on data that resides on remote processes. UPC++ was designed to support exascale high-performance computing, and the library interfaces and implementation are focused on maximizing scalability. In UPC++, all communication operations are syntactically explicit, which encourages programmers to consider the costs associated with communication and data movement. Moreover, all communication operations are asynchronous by default, encouraging programmers to seek opportunities for overlapping communication latencies with other useful work. UPC++ provides expressive and composable abstractions designed for efficiently managing aggressive use of asynchrony in programs. Together, these design principles are intended to enable programmers to write applications using UPC++ that perform well even on hundreds of thousands of cores.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Dynamical ejecta from binary neutron star mergers: Impact of a small residual eccentricity and of the equation of state implementation

Predicting the properties of the matter ejected during and after a neutron star merger is crucial to our ability to use electromagnetic observations of these mergers to constrain the masses of the neutron stars, the equation of state of dense matter, and the role of neutron star mergers in the enrichment of the Universe in heavy elements. Furthermore, our ability to reliably provide such predictions is however limited by a broad range of factors, including the finite resolution of numerical simulations, their treatment of magnetic fields, neutrinos, and neutrino-matter interactions, and the approximate modeling of the equation of state of dense matter. In this manuscript, we study specifically the role that a small residual eccentricity and different implementations of the same equation of state have on the matter ejected during the merger of a 1.3M ⊙ –1.4M ⊙ binary neutron star system. We find that a residual eccentricity e ~ 0.01, as measured ~ 4–6 orbits before merger, causes O(25%–30%) changes in the amount of ejected mass, mainly due to changes in the amount of matter ejected as a result of core bounces during merger. We note that O(1%) residual eccentricities have regularly been used in binary neutron star merger simulations as proxy for circular binaries, potentially creating an additional source of error in predictions for the mass of the dynamical ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗