Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large data sets”

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 775 records · Page 43

Distribution Substation Planning Toolkit (dsp-toolkit) v1.0

The Distribution Substation Planning Toolkit (DSP Toolkit) is a software suite designed to streamline the planning and optimization of distribution substations. This toolkit offers a comprehensive set of tools and APIs for data curation, short-term electric load forecasting, and weather-sensitive load adjustment, making it an essential resource for utility companies, engineers, and researchers. Features • Data Preprocessing and Curation: Efficiently manage and preprocess large datasets to ensure high-quality input for analysis. • Short-Term Load Forecasting: Utilize data-driven models to predict short-term electric loads accurately. • Weather-Sensitive Modeling: Automatically adjust load forecasts based on weather data to predict future peak demands more precisely. Uses The DSP Toolkit is ideal for planning and optimizing distribution substations, providing a user-friendly interface and comprehensive documentation. It is suitable for both novice and experienced users, facilitating efficient and accurate planning processes. Advantages • Efficiency: Automates complex planning tasks, reducing manual effort and minimizing errors. • Scalability: Handles large datasets and complex models, making it suitable for large-scale projects. • Community and Support: Open-source with active community contributions, ensuring continuous improvement and support. • Extensibility: Easily extendable with custom modules and plugins, allowing users to tailor the toolkit to their specific needs. The DSP Toolkit stands out by offering a robust, flexible, and user-friendly solution for distribution substation planning. Public Abstract

Li, Han [Lawrence Berkeley National Laboratory (LB↗

Population Structure of White Sturgeon (Acipenser transmontanus) in the Columbia River Inferred from Single-Nucleotide Polymorphisms

White sturgeon (Acipenser transmontanus) are the largest freshwater fish in North America, with reproducing populations in the Sacramento-San Joaquin, Fraser, and Columbia River Basins. Of these, the Columbia River is the largest, but it is also highly fragmented by hydroelectric dams, and many segments are characterized by declining abundance and persistent recruitment failure. Efforts to conserve and supplement these fish requires an understanding of their spatial genetic structure. Here, we assembled a large set of samples from throughout the Columbia River Basin, along with representative collections from adjacent basins, and genotyped them using a panel of 325 single-nucleotide markers. Results from individual- and group-based analyses of these data indicate that white sturgeon in the uppermost Columbia River Basin, in the Kootenai and upper Snake Rivers, are the most distinct, while the remaining populations downstream in the basin can be described as a genetic gradient consistent with an isolation-by-distance effect. Notably, the population in the lowest reaches of the Columbia River is more distinct from the middle or upper reaches than from outside basins, and suggests historically a higher or more recent gene exchange through coastal routes than with populations in the interior Columbia Basin. Nonetheless, proximal reaches were generally only marginally or non-significantly divergent, suggesting that transplanting larvae or juveniles from nearby sources poses relatively little risk of outbreeding depression. Indeed, we inferred examples of dispersal between reaches via close-kin mark-recapture and genetic mark-recapture that indicate movement between nearby reaches is not unusual. Samples from the Kootenai and upper Snake Rivers exhibited notably lower genetic diversity than the remaining samples as a result of population bottlenecks, genetic drift, and/or historical divergence. Conservation actions, such as supplementation, are underway to maintain population viability and will require balanced efforts to increase demographic abundance while maintaining genetic diversity.

Willis, Stuart C.↗

DeepZipper: A Novel Deep-learning Architecture for Lensed Supernovae Identification

Large-scale astronomical surveys have the potential to capture data on large numbers of strongly gravitationally lensed supernovae (LSNe). To facilitate timely analysis and spectroscopic follow-up before the supernova fades, an LSN needs to be identified soon after it begins. To quickly identify LSNe in optical survey data sets, we designed ZipperNet, a multibranch deep neural network that combines convolutional layers (traditionally used for images) with long short-term memory layers (traditionally used for time series). We tested ZipperNet on the task of classifying objects from four categories—no lens, galaxy-galaxy lens, lensed Type-Ia supernova, lensed core-collapse supernova—within high-fidelity simulations of three cosmic survey data sets: the Dark Energy Survey, Rubin Observatory’s Legacy Survey of Space and Time (LSST), and a Dark Energy Spectroscopic Instrument (DESI) imaging survey. Among our results, we find that for the LSST-like data set, ZipperNet classifies LSNe with a receiver operating characteristic area under the curve of 0.97, predicts the spectroscopic type of the lensed supernovae with 79% accuracy, and demonstrates similarly high performance for LSNe 1–2 epochs after first detection. We anticipate that a model like ZipperNet, which simultaneously incorporates spatial and temporal information, can play a significant role in the rapid identification of lensed transient systems in cosmic survey experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A differentiable approach to the maximum independent set problem using dataless neural networks

The success of machine learning solutions for reasoning about discrete structures has brought attention to its adoption within combinatorial optimization algorithms. Such approaches generally rely on supervised learning by leveraging datasets of the combinatorial structures of interest drawn from some distribution of problem instances. Reinforcement learning has also been employed to find such structures. Here, in this paper, we propose a different approach in that no data is required for training the neural networks that produce the solution. In this sense, what we present is not a machine learning solution, but rather one that is dependent on neural networks and where backpropagation is applied to a loss function defined by the structure of the neural network architecture as opposed to a training dataset. In particular, we reduce the popular combinatorial optimization problem of finding a maximum independent set to a neural network and employ a dataless training scheme to refine the parameters of the network such that those parameters yield the structure of interest. Additionally, we propose a universal graph reduction procedure to handle large-scale graphs. The reduction exploits community detection for graph partitioning and is applicable to any graph type and/or density. Experimental results on both real and synthetic graphs demonstrate that our proposed method performs on par or outperforms state-of-the-art learning-based methods in terms of the size of the found set without requiring any training data.

97 MATHEMATICS AND COMPUTING↗

Pattern of global spin alignment of ϕ and K*0 mesons in heavy-ion collisions

The strong force, as one of the four fundamental forces at work in the universe, governs interactions of quarks and gluons, and binds together the atomic nucleus. Notwithstanding decades of progress since Yukawa first developed a description of the force between nucleons in terms of meson exchange 1 , a full understanding of the strong interaction remains a major challenge in modern science. One remaining difficulty arises from the non-perturbative nature of the strong force, which leads to the phenomenon of quark confinement at distance scales on the order of the size of the proton. Here we show that in relativistic heavy-ion collisions, where quarks and gluons are set free over an extended volume, two species of produced vector (spin-1) mesons, namely φ and K∗0 , emerge with a surprising pattern of global spin alignment. In particular, the global spin alignment for φ is unexpectedly large, while that for K∗0 is consistent with zero. The observed spin-alignment pattern and magnitude for the φ cannot be explained by conventional mechanisms, while a model with strong force fields2, 3 accommodates the current data. This is the first time that the strong force field is experimentally supported as a key mechanism that leads to global spin alignment. We extract a quantity proportional to the intensity of the field of the strong force. Within the framework of the Standard Model, where the strong force is typically described in the quark and gluon language of Quantum Chromodynamics, the field being considered here is an effective proxy description. This is a qualitatively new class of measurement, which opens a new avenue for studying the behaviour of strong force fields via their imprint on spin alignment.

Abdallah, M. S.↗

Structure-aware Initialization via Numerical Continuation and Informed Priors

Scientific machine learning (SciML) often operates in ill-conditioned, weakly identifiable regimes due to limited data or indirect observations. In such settings, optimization and inference are highly sensitive to the starting point, making initialization--often under-reported--a consequential degree of freedom. Random initialization is not a neutral default as it induces an implicit prior over candidate solutions and can systematically bias the result, producing large run-to-run variability. Here, we formalize this view by treating initialization as a hidden confounder in SciML and develop a unifying theory for structure-aware initialization via numerical continuation, constructing warm starts from related problem instances. Across representative tasks, including physics-informed neural networks, maximum likelihood estimation, and variational inference, warm starts have been shown to consistently reduce optimization effort and improve reliability.

Data integrity↗

SECARB Regional Project Assessment

The SECARB Cranfield Early Test Project was designed to safely demonstrate large-scale, long-term CO 2 injection and storage in a CO 2 -enhanced oil recovery (EOR) operation and associated saline reservoir. Storage during EOR operations has the potential to offset commercial costs, holding a significant promise for future large-scale deployment in the southeastern United States. Denbury Onshore LLC began CO 2 -EOR operations in July 2008 at the Cranfield oilfield located east of Natchez, Mississippi (Figure 1). The location offered an ideal setting to investigate storage monitoring at a commercial CO 2 -EOR operation. The field team, led by the Gulf Coast Carbon Center of the Bureau of Economic Geology at the University of Texas at Austin, deployed the largest suite of monitoring, verification, accounting (MVA) technologies to collect data to assess long-term CO 2 utilization and storage. The CO 2 injection interval is within the 3,000-meter-deep fluvial lower Tuscaloosa Formation, a large and regionally extensive saline reservoir, at a broad four-way structural closure that defines the Cranfield oilfield. Tests were conducted in the oil-producing area as well as in the non-oil-producing brine aquifer. The CO 2 used during EOR operation was transported via pipeline from the Jackson Dome, a natural source of CO 2 . The study focused on four integrated research areas: (1) the High-Volume Injection Test area (HiVIT); (2) the Detailed Area of Study (DAS); (3) the Geomechanical area; and (4) the near-surface observatory, also called the “P-site.” CO 2 injection occurred at the HiVIT and the DAS (Figure 2 and Figure 3). The DAS was located down-dip of the oil-water interface of the Cranfield oilfield. The well layout included one injection well and two down-dip observation wells. The team was the first of the regional partnerships to monitor more than 1 million metric tons of CO 2 . At the end of the project, the team monitored more than 11 million metric tons of CO 2 injected and 5 million metric tons of CO 2 stored. The Cranfield Project was the fifth project worldwide, and the first such project in the United States, to reach this CO 2 injection volume while being monitored to demonstrate storage effectiveness.

01 COAL, LIGNITE, AND PEAT↗

New physics versus quenching factors in Coherent Neutrino Scattering

Recent results on the Coherent Elastic Neutrino-Nucleus Scattering (CEνNS) on germanium present significant discrepancies among experiments. We perform a combined analysis of the Dresden-II, CONUS+ and COHERENT data, quantifying the impact of quenching factor uncertainties on their CEνNS cross section measurement. No choice of quenching factor can bring these three data sets into mutual agreement, whereas the combination of COHERENT with either Dresden-II or CONUS+ agrees well albeit for very different quenching factors. We further study the quenching factor dependence on the sensitivity of these experiments to a large neutrino magnetic moment, finding that the constraints can vary by up to an order of magnitude. Our work highlights the importance of reducing this uncertainty on quenching factors in order to probe new physics from neutrinos at the low-energy frontier.

neutrino interactions↗

Portable Distortion Free Large Solid Angle Coverage Detector

Neutron single crystal diffractometers require large solid angle coverage for optimum performance. This can be achieved by tiling flat detectors in a cylindrical or spherical geometry, but results in large gaps in detector coverage and parallax distortion. The detector edges exhibit degraded resolution, distortion, and gamma rejection. Since the detector edge regions are a significant fraction of the detector active area, they must be removed from the experimental data set, requiring extra beam time to collect enough analysis data. A spherical detector with a continuous surface would effectively address this issue while eliminating most boundary ‘dead’ areas. Here we report on the development of a novel hemispherical shaped neutron detector using seamlessly tiled readout modules to form the desired shape. The heart of the detector is a specially developed curved neutron scintillator coupled to high resolution silicon photomultiplier (SiPM) Anger cameras via custom made fiber optic tapers (FOTs). The detector has been assembled and initial tests have been conducted at the High Flux Isotope Reactor (HFIR) beamlines at Oak Ridge National Laboratory (ORNL). Here, in this work, we describe details of the scintillator design, fabrication and characterization, evaluation of individual detector modules, the details of the detector design implementation, and evaluation of the assembled detector at ORNL beamlines.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Predicting the Activity and Selectivity of Bimetallic Metal Catalysts for Ethanol Reforming using Machine Learning

Machine learning is ideally suited for the pattern detection in large uniform datasets, but consistent experimental datasets on catalyst studies are often small. Here we demonstrate how a combination of machine learning and first-principles calculations can be used to extract knowledge from a relatively small set of experimental data. The approach is based on combining a complex machine-learning model trained on an extensive computational library of transition-state energies with simple linear regression models of experimental catalytic activities and selectivities from the literature. Using the combined model, we identify the key C–C bond scission reactions involved in ethanol reforming and perform a computational screening for ethanol reforming on monolayer bimetallic catalysts with architectures TM-Pt-Pt(111) and Pt-TM-Pt(111) (TM = 3d transition metals). The model also predicts four promising catalyst compositions for future experimental studies. In conclusion, the approach is not limited to ethanol reforming but is of general use for the interpretation of experimental observations as well as for the computational discovery of novel catalytic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cosmological analysis of the DESI DR1 Lyα 1D power spectrum

We present the cosmological analysis of the one-dimensional Lyman-α flux power spectrum from the first data release of the Dark Energy Spectroscopic Instrument (DESI). We capture the dependence of the signal on cosmology and intergalactic medium physics using an emulator trained on a cosmological suite of hydrodynamical simulations, and we correct its predictions for the impact of astrophysical contaminants and systematics, many of these not considered in previous analyses. We employ this framework to constrain the amplitude and logarithmic slope of the linear matter power spectrum at k ★ = 0.009 km -1 s and redshift z = 3, obtaining Δ 2 ★ = 0.379 ± 0.032 and n ★ = -2.309 ± 0.019 https://github.com/igmhub/cobaya_lya_p1d. The robustness of these constraints is validated through the analysis of mocks and a large number of alternative data analysis variations, with cosmological parameters kept blinded throughout the validation process. We then combine our results with constraints from DESI BAO and temperature, polarization, and lensing measurements from Planck, ACT, and SPT-3G to set constraints on ΛCDM extensions. While our measurements do not significantly tighten the limits on the sum of neutrino masses from the combination of these probes, they sharpen the constraints on the effective number of relativistic species, N eff = 3.02 ± 0.10, the running of the spectral index, α s = 0.0014 ± 0.0041, and the running of the running, β s = -0.0006 ± 0.0048, by a factor of 1.18, 1.27, and 1.90, respectively. We conclude by outlining the improvements needed to fully reach the level of confidence implied by these uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Validity of Machine Learning in the Quantitative Analysis of Complex Scanning Near-Field Optical Microscopy Signals Using Simulated Data

Scattering-type scanning near-field optical microscope (s-SNOM) is a modern technique for subdiffractional optical imaging and spectroscopy. Over the past two decades, tremendous efforts have been devoted to modeling complex tip-sample interactions in s-SNOM, aimed at understanding the electrodynamics of materials at the nanoscale. However, due to complexities in analytical methods and the limited computation power for fully numerical simulations, compromises must be made to facilitate the modeling of tip-sample interaction, such as using quasistatic approximation or unrealistic tip geometries. Here, we apply a variety of widely utilized machine-learning methods, including k nearest neighbor and feedforward neural network etc. to study the phase-resolved spectroscopic near-field response. With only a small set of training data, which is simulated using the finite-dipole model, we demonstrate that the relation between the experimental near-field signal and sample optical constant can be one to one mapped without the need for tip modeling: for a given material with a moderate dielectric function, its complex near-field spectrum can be accurately determined within the mid-IR spectral range, and vice versa. Our preliminary study sets the stage for future exploration using real experimental data. Our method is beneficial for processing the increasing amount of data accumulated across many research groups and especially useful for user facilities such as synchrotron-based national laboratories where a large amount of data is generated on a daily basis.

36 MATERIALS SCIENCE↗

Photovoltaic Data Acquisition (PVDAQ) Public Datasets

The NREL PVDAQ is a large-scale time-series database containing system metadata and performance data from a variety of experimental PV sites and commercial public PV sites. The datasets are used to perform on-going performance and degradation analysis. Some of the sets can exhibit common elements that effect PV performance (e.g. soiling). The dataset consists of a series of files devoted to each of the systems and an associated set of metadata information that explains details about the system hardware and the site geo-location. Some system datasets also include environmental sensors that cover irradiance, temperatures, wind speeds, and precipitation at the site.

Array↗

Binning is Sinning (Supernova Version): The Impact of Self-calibration in Cosmological Analyses with Type Ia Supernovae

Recent cosmological analyses (e.g., JLA, Pantheon) of Type Ia supernovae (SNe Ia) have propagated systematic uncertainties into a covariance matrix and either binned or smoothed the systematic uncertainty vectors in redshift space. We demonstrate that systematic error budgets of these analyses can be improved by a factor of ~ 1.5 × with the use of unbinned and unsmoothed covariance matrices. To understand this, we employ a separate approach that simultaneously fits for cosmological parameters and additional self-calibrating scale parameters that constrain the size of each systematic. We show that the covariance-matrix approach and scale-parameter approach indeed yield equivalent results, implying that in both cases the data can self-calibrate certain systematic uncertainties, but that this ability is hindered when information is binned or smoothed in redshift space. We review the top systematic uncertainties in current analyses and find that the reduction of systematic uncertainties in the unbinned case depends on whether a systematic is solely degenerate with the cosmological model in redshift space or whether it can be described by additional correlations between supernova properties and luminosity. Furthermore, we show that the power of self-calibration increases with the size of the data set, which presents a tremendous opportunity for upcoming analyses of photometrically classified samples, like those of Legacy Survey of Space and Time (LSST) and the Nancy Grace Roman Telescope (NGRST). However, to take advantage of self-calibration in large, photometrically classified samples, we must first address the issue that binning is required in currently used photometric analysis methods.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine-Learned Manifold-Based Models for Large Eddy Simulation of Turbulent Combustion

Reduced-order manifold combustion models are commonly used to lower the cost of reacting Large Eddy Simulations (LES) and may be derived either from physical principles as in flamelet models or from data-driven methods like Principal Component Analysis (PCA). In either case, neural networks are increasingly used as part of these models to provide a nonlinear mapping between a small set of pre-defined variables that parameterize the manifold and outputs of interest, such as reaction rates. In this work, we propose a new manifold-based modeling approach that combines the definition of the manifold-parameterizing variables (linear combinations of species), the nonlinear mapping to the outputs, and closure of filtered quantities for LES into the structure of a single neural network. This allows the process used to train the neural network to simultaneously optimize both the functional form of the model and the identities of the inputs to the model. The new approach can flexibly incorporate thermochemical data from any combustion system; if trained on data from 1D flames it can be interpreted as an optimized flamelet model, but it can also be used to learn models from data from more complex configurations. This work presents a priori evaluations of the new approach in both contexts. Evaluation using data from 1D premixed flames demonstrates the physical interpretability of the manifold variables generated by the new approach. Evaluation using data from direct numerical simulations of turbulent flames shows improved predictions relative to either flamelet or PCA-based models in a more complex configuration.

47 OTHER INSTRUMENTATION↗

EPCAPE-PT-LANL Measurements: Single Particle Soot Photometer

Coastal cities offer a unique environment for studying aerosol-cloud interactions and the effects of urban emissions on cloud properties. As part of the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE), the Partitioning Thrust by Los Alamos National Laboratory (EPCAPE-PT-LANL) was conducted. Our campaign focused on measuring the optical and chemical properties of aerosols and their interactions within marine stratocumulus clouds in La Jolla, California. EPCAPE-PT-LANL enhances the primary goals of EPCAPE through innovative observations of vapor-phase transitions between aerosols and cloud droplets, the impact of black carbon on aerosol-cloud dynamics, and the effects of cloud processing on aerosol optical properties. Instrument: Single Particle Soot Photometer (Droplet Measurements Technology) Data Notes: Contact us if you want additional data products from this instrument. Reported data is the black carbon (rBC) number and mass concentration with diagnostic flags. The SP2 measures incandescence from particles that is induced with a 1064 Nd-YAG laser. Particles that absorb the laser energy and then emit radiation is assumed to contain black carbon. The amount of intensity of radiation is related to the mass of absorbing material in the particle. Here, we calibrate the incandescence intensity to size selected Regal Black (Cabot) nebulized from solution. Data provided for the EPCAPE campaign used the combined broadband high-gain + low-gain channels (BHBL). Each channel has a lower threshold of detection equivalent to 2-s of the respective channel noise. Thresholding has been applied to select between the high-gain and low-gain channels. The detection limit is 80-540 nm (Deq) or 0.36- 55 fg. It is assumed that rBC is the only aerosol type that is in significant concentration in the sampled atmosphere that absorbs the laser energy. Particle data are integrated over 10 second windows to calculate a rBC number and mass concentration. QC/QA: • Diagnostic flags that impacted measured concentration: - Sample, Sheath, and Purge Flow Rates: Despite observing fluctuations in all flows, the rBC detection and mass quantification are generally observed to be stable, although large changes in sample flow rate did impact detection efficiency. A flag was implemented for sample flow deviations >12 cm3/min from the set point averaged over 30 seconds - Laser Power: Detection efficiency decreases with laser power and deviations in laser power also affect mass quantification. Flag for laser power is set for deviations in laser current from the set value >5 mA. - Primary Detector Threshold: Primary thresholding is the signal value which determines whether a “particle” is recorded. Thresholding errors occur when the threshold value is too HIGH and real particles are ignored. • Several periods without data: - 16-18 Nov: Ultra Zero Air generator failure, flows deviated significantly from set points. rBC # conc recovered but questionable. - 20-21 Nov: Offline for calibrations for several hours each day. - 25-27 Nov: Data is missing. - 30 Nov: Power outage 3 Dec: SP2 hard-drive full. - 3 Dec: Power Outage Header: - BHBL_NumbConc[#/cc]: Refers to the number concentration of refractory black carbon (rBC) measured from combined broadband high-gain and low-gain channels, expressed in particles per cubic centimeter. - BHBL_BCMass_Conc[ug/m3]: Refers to the mass concentration of refractory black carbon (rBC) measured from combined broadband high-gain and low-gain channels, expressed in micrograms per cubic meter. - NoData_Flag[bool]: A boolean flag that indicates whether no data was recorded during a measurement. - NoBC_Flag[bool]: A boolean flag indicating whether no black carbon particles were recorded during the measurement. - Laser_Flag[bool]: A boolean flag indicating deviations in laser current during the measurement. - SampleFlow_Flag[bool]: A boolean flag indicating deviations from the set sample flow rate during the measurement. - PrimThresh_Flag[bool]: A boolean flag indicating deviations from the set primary threshold during the measurement, potentially ignoring real particles. - Manual_Flag[bool]: A boolean flag indicating manual intervention or adjustments during the measurement. - CVI_Flag[bool]: A boolean flag indicating whether the Counterflow Virtual Impactor (CVI) was active or inactive during the measurement.

54 ENVIRONMENTAL SCIENCES↗

Disturbance and diversity: Lichen species richness decreases with increasing anthropogenic disturbance

Anthropogenic disturbance is rapidly increasing through habitat degradation, development, and deforestation. Gaps remain in understanding the effects of this disturbance on diverse and ecologically important organisms such as lichens. In North America, studies have focused on epiphytic macrolichens and catastrophic disturbance, largely ignoring microlichens and less severe disturbances. The Intermediate Disturbance Hypothesis predicts these moderate disturbances will lead to higher species richness. Here we examine the effects of anthropogenic disturbance in the form of land management on overall lichen richness, including microlichens, and on the species richness of eight lichen functional groups. The study draws on a comprehensive data set of 872 species, in 208 one-hectare plots throughout the Southern Appalachian Mountains, a global biodiversity hotspot in eastern North America. A habitat quality index based on an established forestry metric was used as a proxy for anthropogenic disturbance and was quantified using a 10-part score including categories such as percent native tree canopy cover and degree of fragmentation. Linear models were used to compare habitat quality scores to overall species richness and to species richness of functional groups. Rather than following the Intermediate Disturbance Hypothesis, all groups uniformly follow a negative linear relationship: as disturbance increases, species richness decreases. This pattern held even for widely variable functional groups such as morphotype. Effective conservation of lichen richness should prioritize the maintenance of existing older, less-disturbed stands within large, contiguously forested areas.

59 BASIC BIOLOGICAL SCIENCES↗

Handing off the outcome of binary neutron star mergers for accurate and long-term postmerger simulations

In this work, we perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (~1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the HandOff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to harm3d, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the HandOff approach.

79 ASTRONOMY AND ASTROPHYSICS↗