Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Simulation”

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 1,621 records · Page 90

Plant roots stimulate the decomposition of complex, but not simple, soil carbon

Roots release carbon into soil and can alleviate energy limitation of microbial organic matter decomposition. We know little about the effects of roots on microbial decomposition of different organic matter substrates, despite the importance for soil carbon stocks and turnover. Through implementing root-microbe interactions, the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model was previously shown to represent dynamics of total soil carbon in temperate forest field experiments. However, the model permits alternative hypotheses concerning microbial-substrate affinity. We investigated how root inputs affect decomposition of soil organic carbon (SOC) with variable decomposability. We simulated SOC stocks in CORPSE and compared microbial degradation of two substrates types with varying root-microbe interactions under two alternative hypotheses that varied in microbial-substrate affinity. We compared our modeled hypotheses to a forest field experiment where we quantified decomposition of isotopically-labeled starch and leaf tissues in soils with manipulated root access to microbes. Here, we tested the hypothesis that decomposition of leaves would be more sensitive to root inputs than decomposition of starch, corresponding to the alternative model hypothesis. In the field study, leaf decomposition increased with root density while starch decomposition was unchanged by root density. Microbial biomass and enzyme activity consistently increased with root inputs in CORPSE and the field study. Our field experiment supported the CORPSE simulations with high microbial-substrate affinity. Roots stimulated microbial growth and enzyme production, which increased degradation of more complex substrates such as leaf tissues. Substrates that were easily decomposed, such as starch, may already be degrading at a maximum rate in the absence of rhizosphere influence because their decomposition rate was unchanged by root inputs. We found that the degree to which roots stimulate microbial decomposition depends on the substrate being decomposed, and that root-microbe interactions influenced SOC stocks in both our model and field experiment. Environmental changes that alter root-microbe interactions could, therefore, alter soil C stocks and biogeochemical cycling, and models of these interactions should incorporate differential influence of rhizosphere inputs on different substrates.

59 BASIC BIOLOGICAL SCIENCES↗

Geothermal Fault Zone and Fluid Imaging through Joint Airborne ZTEM and Ground MT Data Inversion Analysis

This project has aimed to achieve detailed electrical resistivity resolution at geothermal reservoir scales by combining airborne natural electromagnetic (EM) field surveying (ZTEM) with ground magnetotelluric (MT) measurements to approximate an airborne MT geophysical method. MT alone is relatively expensive and may have permitting challenges in sensitive areas. Airborne ZTEM field data contains only the magnetic field, requires a background assumption, and has been limited to relatively high frequencies, thus suffering uniqueness problems. Based on proto-type 2D simulations, ZTEM ambiguities may be reduced through formal incorporation with possibly sparse ground MT soundings, which we pursued in full 3D for this project. The methodology was tested at the high-temperature Roosevelt Hot Springs geothermal system, Utah, which was considered advantageous given the near total exposure of crystalline reservoir rocks across the project area. ZTEM and ground MT survey data were acquired in 2017, subcontracted to outside parties with which we have worked in the past. These included 80 remote-referenced tensor MT soundings over the Mineral Mountains and adjacent Roosevelt Hot Spring producing geothermal system. These MT stations abut later coverage of a similar number of MT stations taken for the Utah FORGE project providing excellent total data aperture to re-solve structure beneath both project areas better than either set alone. The airborne ZTEM survey covered 704 line kilometers in E-W flight lines with a 250 m line spacing. Although this survey was timed during a maintenance-related shutdown of power production at the Roosevelt Hot Springs, other noise sources difficult to identify but including two high-voltage state-scale transmission lines compromised the ZTEM survey badly leading to unusable responses. Thus, with DOE management concurrence, the project proceeded to emphasize inversion and interpretation of the joint SubTER-FORGE MT data sets with regard to the Roosevelt Hot Springs reservoir recharge and to deep heat sources for both it and the Utah FORGE EGS project area. We also investigated the joint ZTEM-MT sampling concept with data sets from the Eleven Mile Canyon prospect area donated by the U.S. Navy (A. Sabin, PoC). Inversion of the SubTER-FORGE MT data using the HexMT 3D finite element algorithm reveals a large, low-resistivity anomaly extending sub-vertically through the depth range of the crust beneath the western Mineral Mountains. The steep conductive zone connects in the lower crust to a more tabular conductor characteristic of much of the Great Basin that generally is ascribed to current mafic magmatic underplating, hybridization and fluid release. The location of the resolved anomaly relative to the recent (0.5-0.8 Ma) eruptive centers of the Mineral Mountains implicates it as remnants of the magma body which fed these centers. This structure appears to be currently feeding heat and fluids upward into the Roosevelt Hot Springs hydrothermal system, as well as heat laterally to the FORGE project area. Separate and joint inversion models were carried out for the donated Eleven Mile Canyon MT-ZTEM data set to demonstrate concept. ZTEM only inversion showed two main alteration zones in the western portion of the project area known from geological mapping. Joint inversion including an E-W profile of MT soundings sharpened these features considerably. It also resolved in much greater detail the graben related normal faulting structure of the central project area which lies at depths exceeding the sensitivity of ZTEM alone. The sparse number of MT da-ta relative to the ZTEM required upweighting the former by a factor of several, but an exact procedure awaits future research. Our final impression is that sparse MT data can improve resolution of the subsurface over that of ZTEM alone. However, well sampled MT data are to be preferred and offer the simplicity of interpreting just one data type, and possess the superior resolution capability coming with the electric field everywhere, and from their high bandwidth.

15 GEOTHERMAL ENERGY↗

New technology for an ancient fish: A lamprey life cycle modeling tool with an R Shiny application

Lampreys (Petromyzontiformes) are an ancient group of fishes with complex life histories. We created a life cycle model that includes an R Shiny interactive web application interface to simulate abundance by life stage. This will allow scientists and managers to connect available demographic information in a framework that can be applied to questions regarding lamprey biology and conservation. We used Pacific lamprey ( Entosphenus tridentatus ) as a case study to highlight the utility of this model. We applied a global sensitivity analysis to explore the importance of individual life stage parameters to overall population size, and to better understand the implications of existing gaps in knowledge. We also provided example analyses of selected management scenarios (dam passage, fish translocations, and hatchery additions) influencing Pacific lamprey in fresh water. These applications illustrate how the model can be applied to inform conservation efforts. This tool will provide new capabilities for users to explore their own questions about lamprey biology and conservation. Simulations can hone hypotheses and predictions, which can then be empirically tested in the real world.

Gomes, Dylan G. E. (ORCID:0000000226423728)↗

Process design and techno-economic analysis for bio-based graphite and liquid hydrocarbons production from lignocellulosic biomass

The worldwide demand for graphite, as the main anode material for Li-ion batteries, is expected to double by 2028 since it supports the use of electricity, including transient renewable sources, for energy storage, sustainable mobility, and automation. However, the dependence on non-renewable and external resources jeopardizes the world supply chain. This study explores the technical and economic performance of transforming lignocellulosic biomass into biographite and fuel-grade hydrocarbons through pyrolysis bio-oil upgrading. According to simulation results, the total power demand for the biorefinery reached 10,784 kWh per tonne of biographite, of which 36 % can be supplied by the heat integration network and power plant. Sensitivity and risk analyses were conducted to evaluate the economics, with process yields identified as the most relevant indicators to the minimum selling price (MSP). The analysis revealed a promising cost-competitive range for biographite MSP against fossil-based graphite (medium quality synthetic graphite Chinese market price ~$\$$4.2/kg). Case D, which includes biofuels as a byproduct, presents the best metrics, reaching a MSP of $\$$3.3/kg of anode-grade biographite with a profit margin of 27 %. While including biofuels in the product slate provides the best economic performance, the uncertainty associated with the big capital investment makes its risk 13 % higher to attain an IRR >20 % than the case in which biographite is produced as a standalone product. Overall, this study demonstrates that integrated biorefineries can produce a cost-competitive bio-based anode material for Li-ion batteries.

09 BIOMASS FUELS↗

Multilevel peel-off patterning of a prototype semitransparent organic photovoltaic module

Semitransparent organic photovoltaics (ST-OPVs) with applications to power generating windows have shown substantial increases in power conversion efficiency (PCE) and average photopic transmission (APT) at the laboratory scale. Here, the demonstration of similarly efficient, large-scale ST-OPV modules with geometric fill factors (GFF) approaching 100%, however, remains a challenge. Here, we employ a multilevel peel-off patterning method that can achieve micron-scale resolution without exposing chemically sensitive organic materials to solvents. Eight, 4 cm × 0.4 cm cells are connected in series to realize a prototype ST-OPV module with GFF = 95.8%, with PCE = 7.3 ± 0.2% under simulated AM 1.5G illumination at 1 sun intensity, APT = 41.8 ± 1.4%, and a light utilization efficiency of LUE = 3.1 ± 0.1%. A neutral color STOPV module is also demonstrated with 1.7 ± 0.1% and International Commission on Illumination (CIE) LAB coordinates of (L*, a*, b*) = (53.7, -1.9, -3.9).

14 SOLAR ENERGY↗

The Work of Mechanical Degradation in Elongating Polymer Melts

Molecular dynamics simulations are used to study the mechanical degradation of well-entangled polymer melts during uniaxial extensional flow. Simulations measure the transient rise in extensional stresses and relate them to the molecular alignment and scission of chain backbones. Intermolecular entanglements couple chain scission in space and time, making degradation sensitive to deformation history and strain rate in ways not displayed by dilute polymer solutions. The rate of chain scission is nonmonotonic and peaks at strains corresponding to the maximum extensibility of entanglement segments but prior to the full extension of chain backbones. We measure a specific work per scission event w* and decompose it into separate contributions associated with chain alignment, chemical bond breaking, and scission-induced plasticity. We find chain scission in melts requires activating plastic dissipation that is multiple orders of magnitude larger than the chemical work required to break a covalent backbone bond. Our findings underscore the critical need to consider bulk polymer mechanics and rheology in designing efficient mechanical degradation and mechanochemical processes.

36 MATERIALS SCIENCE↗

Understanding the Effect of Lead Iodide Excess on the Performance of Methylammonium Lead Iodide Perovskite Solar Cells

The presence of unreacted lead iodide in organic-inorganic lead halide perovskite solar cells is widely correlated with an increase in power conversion efficiency. We investigate the mechanism for this increase by identifying the role of surfaces and interfaces present between methylammonium lead iodide perovskite films and excess lead iodide. We show how type I and II band alignments arising under different conditions result in either passivation of surface defects or hole injection. Through first-principles simulations of solid-solid interfaces, we find that lead iodide captures holes from methylammonium lead iodide and modulates the formation of defects in the perovskite, affecting recombination. Using surface-sensitive optical spectroscopy techniques, such as transient reflectance and time-resolved photoluminescence, we show how excess lead iodide affects the diffusion and surface recombination velocity of charge carriers in methylammonium lead iodide films. Our coupled experimental and theoretical results elucidate the role of excess lead iodide in perovskite solar cells.

14 SOLAR ENERGY↗

Probing lithium mobility at a solid electrolyte surface

Abstract Solid-state electrolytes overcome many challenges of present-day lithium ion batteries, such as safety hazards and dendrite formation 1,2 . However, detailed understanding of the involved lithium dynamics is missing due to a lack of in operando measurements with chemical and interfacial specificity. Here we investigate a prototypical solid-state electrolyte using linear and nonlinear extreme-ultraviolet spectroscopies. Leveraging the surface sensitivity of extreme-ultraviolet-second-harmonic-generation spectroscopy, we obtained a direct spectral signature of surface lithium ions, showing a distinct blueshift relative to bulk absorption spectra. First-principles simulations attributed the shift to transitions from the lithium 1 s state to hybridized Li- s /Ti- d orbitals at the surface. Our calculations further suggest a reduction in lithium interfacial mobility due to suppressed low-frequency rattling modes, which is the fundamental origin of the large interfacial resistance in this material. Our findings pave the way for new optimization strategies to develop these electrochemical devices via interfacial engineering of lithium ions.

25 ENERGY STORAGE↗

Charge transfer as a mechanism for chlorophyll fluorescence concentration quenching

Highly concentrated solutions of chlorophyll display rapid fluorescence quenching. The same devastating energy loss is not seen in photosynthetic light-harvesting antenna complexes, despite the need for chromophores to be in close proximity to facilitate energy transfer. A promising, though unconfirmed mechanism for the observed quenching is energy transfer from an excited chlorophyll monomer to a closely associated chlorophyll pair that subsequently undergoes rapid nonradiative decay to the ground state via a short-lived intermediate charge-transfer state. In this work, we make use of newly emerging fast methods in quantum chemistry to assess the feasibility of this proposed mechanism. We calculate rate constants for the initial charge separation, based on Marcus free-energy surfaces extracted from molecular dynamics simulations of solvated chlorophyll pairs, demonstrating that this pathway will compete with fluorescence (i.e., drive quenching) at experimentally measured quenching concentrations. We show that the rate of charge separation is highly sensitive to interchlorophyll distance and the relative orientations of chromophores within a quenching pair. We discuss possible solvent effects on the rate of charge separation (and consequently the degree of quenching), using the light-harvesting complex II (LH2) protein from rps. acidophila as a specific example of how this process might be controlled in a protein environment. Crucially, we reveal that the LH2 antenna protein prevents quenching, even at the high chlorophyll concentrations required for efficient energy transfer, by restricting the range of orientations that neighboring chlorophyll pairs can adopt.

59 BASIC BIOLOGICAL SCIENCES↗

NuSTAR tests of sterile-neutrino dark matter: New Galactic bulge observations and combined impact

In this work, we analyze two dedicated NuSTAR observations with exposure ~190 ks located ~10° from the Galactic plane, one above and the other below, to search for x-ray lines from the radiative decay of sterile-neutrino dark matter. These fields were chosen to minimize astrophysical x-ray backgrounds while remaining near the densest region of the dark matter halo. We find no evidence of anomalous x-ray lines in the energy range 5–20 keV, corresponding to sterile neutrino masses 10–40 keV. Interpreted in the context of sterile neutrinos produced via neutrino mixing, these observations provide the leading constraints in the mass range 10–12 keV, improving upon previous constraints in this range by a factor ~2. We also compare our results to Monte Carlo simulations, showing that the fluctuations in our derived limit are not dominated by systematic effects. An updated model of the instrumental background, which is currently under development, will improve NuSTAR’s sensitivity to anomalous x-ray lines, particularly for energies 3–5 keV.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterizing the expected behavior of non-Poissonian template fitting

We have performed a systematic study of the statistical behavior of non-Poissonian template fitting (NPTF), a method designed to analyze and characterize unresolved point sources in general counts datasets. Here, in this paper, we focus on the properties and characteristics of the Fermi-LAT gamma-ray data set. In particular, we have simulated and analyzed gamma-ray sky maps under varying conditions of exposure, angular resolution, pixel size, energy window, event selection, and source brightness. We describe how these conditions affect the sensitivity of NPTF to the presence of point sources, for inner-galaxy studies of point sources within the Galactic Center excess, and for the simplified case of isotropic emission. We do not find opportunities for major gains in sensitivity from varying these choices, within the range available with current Fermi-LAT data. We provide an analytic estimate of the NPTF sensitivity to point sources for the case of isotropic emission and perfect angular resolution, and find good agreement with our numerical results for that case.

79 ASTRONOMY AND ASTROPHYSICS↗

Electron-positron, parton-parton, and photon-photon production of τ -lepton pairs: Anomalous magnetic and electric dipole moments spin effects

Anomalous contributions to the electric and magnetic dipole moments of the τ lepton from new physics scenarios have brought renewed interest in the development of new charge-parity violating signatures in τ -pair production at Belle II energies, and also at higher energies of the Large Hadron Collider and the Future Circular Collider. In this paper, we discuss the effects of spin correlations, including transverse degrees of freedom, in the τ -pair production and decay. These studies include calculating analytical formulas, obtaining numerical results, and building semirealistic observables sensitive to the transverse spin correlations induced by the dipole moments of the τ lepton. The effects of such anomalous contributions to the dipole moments are introduced on top of precision simulations of e − e + → τ − τ + , q q ¯ → τ − τ + and γ γ → τ − τ + processes, involving multibody final states. The τ decays are simulated along with radiative corrections, in particular electroweak box contributions of W W and Z Z exchanges are taken into account. Respective extensions of the Standard Model amplitudes and the reweighting algorithms are implemented into the Monte Carlo, which is used to simulate τ -pair production in e − e + collisions, and the aupinner program, which is used to reweight events with τ pair produced in p p collisions. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Improving neutrino energy estimation of charged-current interaction events with recurrent neural networks in MicroBooNE

We present a deep learning-based method for estimating the neutrino energy of charged-current neutrino-argon interactions. We employ a recurrent neural network (RNN) architecture for neutrino energy estimation in the MicroBooNE experiment, utilizing liquid argon time projection chamber (LArTPC) detector technology. Traditional energy estimation approaches in LArTPCs, which largely rely on reconstructing and summing visible energies, often experience sizable biases and resolution smearing because of the complex nature of neutrino interactions and the detector response. The estimation of neutrino energy can be improved after considering the kinematics information of reconstructed final-state particles. Utilizing kinematic information of reconstructed particles, the deep learning-based approach shows improved resolution and reduced bias for the muon neutrino Monte Carlo simulation sample compared to the traditional approach. In order to address the common concern about the effectiveness of this method on experimental data, the RNN-based energy estimator is further examined and validated with dedicated data-simulation consistency tests using MicroBooNE data. We also assess its potential impact on a neutrino oscillation study after accounting for all statistical and systematic uncertainties and show that it enhances physics sensitivity. This method has good potential to improve the performance of other physics analyses. Published by the American Physical Society 2024

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Robust Anthropogenic Signal Identified in the Seasonal Cycle of Tropospheric Temperature

Previous work identified an anthropogenic fingerprint pattern in T AC (x, t), the amplitude of the seasonal cycle of mid- to upper-tropospheric temperature (TMT), but did not explicitly consider whether fingerprint identification in satellite T AC (x, t) data could have been influenced by real-world multidecadal internal variability (MIV). Here we address this question here using large ensembles (LEs) performed with five climate models. LEs provide many different sequences of internal variability noise superimposed on an underlying forced signal. Despite differences in historical external forcings, climate sensitivity, and MIV properties of the five models, their T AC (x, t) fingerprints are similar and statistically identifiable in 239 of the 240 LE realizations of historical climate change. Comparing simulated and observed variability spectra reveals that consistent fingerprint identification is unlikely to be biased by model underestimates of observed MIV. Even in the presence of large (factor of 3–4) intermodel and inter-realization differences in the amplitude of MIV, the anthropogenic fingerprints of seasonal cycle changes are robustly identifiable in models and satellite data. This is primarily due to the fact that the distinctive, global-scale fingerprint patterns are spatially dissimilar to the smaller-scale patterns of internal T AC (x, t) variability associated with the Atlantic multidecadal oscillation and El Niño–Southern Oscillation. The robustness of the seasonal cycle detection and attribution results shown here, taken together with the evidence from idealized aquaplanet simulations, suggest that basic physical processes are dictating a common pattern of forced T AC (x, t) changes in observations and in the five LEs. The key processes involved include GHG-induced expansion of the tropics, lapse-rate changes, land surface drying, and sea ice decrease.

54 ENVIRONMENTAL SCIENCES↗

Modeling Molten Salt Spreading and Heat Transfer using MELTSPREAD – Model Development Updates

This report summarizes spreading and heat transfer results calculated using MELTSPREAD for a scenario of molten FLiNaK spilling onto a flat stainless steel substrate that was inspired by the “maximum credible accident” scenario of the Molten Salt Reactor Experiment (MSRE). MELTSPREAD was developed at Argonne to model the one-dimensional flowing and freezing of molten corium and was applied to model molten salts for the first time during FY21. The molten salt MELTSPREAD model that was developed during FY21 was updated as part of this work using a corrected heat of fusion value of FLiNaK and to accommodate a larger salt spill volume. The model was run with and without the inclusion of contributions from decay heat and a sensitivity analysis of initial spill conditions was performed to determine the importance of those factors on model outcome. The results provide insight into the expected spreading and heat transfer behavior of simulated and irradiated fuel salt that has been spilled and indicate aspects of the model to be revised in future versions of the MELTSPREAD model for applications to molten salt.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

High-Resolution Wind Resource Data Set of the Greater Puerto Rico Region

In February 2022, the U.S. Department of Energy and six national laboratories launched the Puerto Rico Grid Resilience and Transitions to 100% Renewable Energy Study (PR100). PR100 aims to provide a comprehensive analysis of possible pathways for Puerto Rico's energy future, with a goal of 100% renewable energy by 2050. As a part of the renewable energy potential assessment in this project, we developed 20 years (2001-2020) of data using a numerical weather prediction (NWP) model for onshore and offshore wind resource assessment for the Puerto Rico. The research steps in developing the long-term wind resource data sets based on the NWP model were: 1. Model wind resource based on the Weather Research and Forecasting (WRF) model. 2. Develop WRF model configurations for Puerto Rico. 3. Test WRF with 11 different physics parameterizations for planetary boundary layer (PBL). 4. Assess WRF output from the different PBL schemes against observations. 5. Select a final model configuration which can produce the modeled wind speed with sufficient accuracy. 6. Produce 20 years wind resource data sets for Puerto Rico region. In the first stage of our framework for developing wind resource data, we developed a WRF model configuration using two nested domains (9 km and 3 km) to cover Puerto Rico and U.S. Virgin Islands and downscale the ERA5 reanalysis data (0.25 degrees x 0.25 degrees; hourly interval) to a 3-km domain. For the second stage, we implemented one-year simulations focused on using 11 different PBL physics parameterizations to find a combination of WRF physics parameterizations that could provide accurately modeled wind speed for Puerto Rico. We also analyzed the sensitivity of the modeled wind speed to PBL schemes for onshore and offshore locations. The WRF output resulting from the 11 WRF experiments using different PBL parameterizations were evaluated against observations obtained from the National Data Buoy Center (NDBC) as well as at hub height for a location for which measurements were available. A final model setup selected through the validation with observational data was used to produce 20 years of data with 3-km spatial and 5-minute temporal resolution. The WRF model output was post-processed to include wind profiles and basic atmospheric variables in a format that can be easily used for downstream modeling. The 20 years of wind resource data will be made available through NREL and support the estimation of wind energy development costs for the PR100 study.

17 WIND ENERGY↗

Road to PROSPECT-II

The Precision Reactor Oscillation and SPECTrum (PROSPECT) experiment is based in a segmented liquid scintillator antineutrino detector situated approximately 7 meters from the highly enriched High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory. Its main goal is to investigate short-baseline antineutrino oscillations. The first phase of data collection, known as PROSPECT-I, was held from 2018 to 2019 and was used for several high-precision analyses, including multiple measurements of the $^{235}U$ antineutrino spectrum and searches for $eV$-scale sterile antineutrino oscillations. The collaboration is now preparing for the second phase, PROSPECT-II, which features an upgraded detector design. This advancement will enhance sensitivity and statistical power, allowing for a broader range of analyses beyond those achieved in PROSPECT-I. As we transition into this new phase, new questions have arisen concerning background simulation and its potential differences from those conducted during the initial phase of the experiment. Moreover, it is essential to ascertain, through simulations, the positive effects that an improved detector could have on the study of oscillations. This information is crucial for justifying the proposed enhancements, and I will cover all of this in the talk.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Data and Code for Atomic Scale Etching of Diamond: Insights from Molecular Dynamics Simulations

This work investigates the effects of argon ions, hydrogen atoms, and hydrogen ions on the diamond (100) surface using classical molecular dynamics simulations. The purpose of this investigation was to asses plasma processing techniques for applications in quantum device manufacturing. The simulations suggest that combining argon ion smoothing with selective, near threshold energy H removal of amorphous C could be an effective strategy for diamond surface engineering, leading to more reliable and sensitive diamond color center devices. Results were found to differ significantly with interatomic potential, and an analysis of these differences was also carried out. Included in this repository are LAMMPS source files, input scripts, and plotting scripts required to reproduce the data. Also included are the output data required to make all the plots included in the associated publication.

Brenner↗