Search NASA⌕ Search

SEARCH · Search NASA

Results for “root development”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Fast and Invertible Simplicial Approximation of Magnetic‐Following Interpolation for Visualizing Fusion Plasma Simulation Data

We introduce a fast and invertible approximation for fusion plasma simulation data represented as 2D planar meshes with connectivities approximating magnetic field lines along the toroidal dimension in deformed 3D toroidal spaces. Scientific variables (e.g., density and temperature) in these fusion data are interpolated following a complex magnetic-field-line-following scheme in the toroidal space represented by a cylindrical coordinate system. This deformation in the 3D space poses challenges for root-finding and interpolation. To this end, we propose a novel paradigm for visualizing and analyzing such data based on a newly developed algorithm for constructing a 3D simplicial mesh within the deformed 3D space. Our algorithm generates a tetrahedral mesh that connects the 2D meshes using tetrahedra while adhering to the constraints on node connectivities imposed by the magnetic field-line scheme. Specifically, we first divide the space into smaller partitions to reduce complexity based on the input geometries and constraints on connectivities. Then, we independently search for a feasible tetrahedralization of each partition, considering nonconvexity. We demonstrate our method with two X-Point Gyrokinetic Code (XGC) simulation datasets on the International Thermonuclear Experimental Reactor (ITER) and Wendelstein 7-X (W7-X), and use an ocean simulation dataset to substantiate broader applicability of our method. An open source implementation of our algorithm is available at https://github.com/rcrcarissa/DeformedSpaceTet.

Ren, Congrong [The Ohio State Univ., Columbus, OH ↗
DOE OSTI2022

Application of Orthogonal Defect Classification for Software Reliability Analysis

The modernization of existing and new nuclear power plants with digital instrumentation and control systems (DI&C) is a recent and highly trending topic. However, there lacks strong consensus on best-estimate reliability methodologies by both the United States (U.S.) Nuclear Regulatory Commission (NRC) and the industry. This has resulted in hesitation for further modernization projects until a more unified methodology is realized. In this work, we develop an approach called Orthogonal-defect Classification for Assessing Software Reliability (ORCAS) to quantify probabilities of various software failure modes in a DI&C system. The method utilizes accepted industry methodologies for software quality assurance that are also verified by experimental or mathematical formulations. In essence, the approach combines a semantic failure classification model with a reliability growth model to predict (and quantify) the potential failure modes of a DI&C software system. The semantic classification model is used to address the question: How do latent defects in software contribute to different software failure root causes? The use of reliability growth models is then used to address the question: Given the connection between latent defects and software failure root causes, how can we quantify the reliability of the software? A case study was conducted on a representative I&C platform (ChibiOS) running a smart sensor acquisition software developed by Virginia Commonwealth University (VCU). The testing and evidence collection guidance in ORCAS was applied, and defects were uncovered in the software. Qualitative evidence, such as condition coverage, was used to gauge the completeness and trustworthiness of the assessment while quantitative evidence was used to determine the software failure probabilities. The reliability of the software was then estimated and compared to existing operational data of the sensor device. It is demonstrated that by using ORCAS, a semantic reasoning framework can be developed to justify if the software is reliable (or unreliable) while still leveraging the strength of the existing methods.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Application of Orthogonal Defect Classification for Software Reliability Analysis

The modernization of existing and new nuclear power plants with digital instrumentation and control systems (DI&C) is a recent and highly trending topic. However, there lacks strong consensus on best-estimate reliability methodologies by both the United States (U.S.) Nuclear Regulatory Commission (NRC) and the industry. This has resulted in hesitation for further modernization projects until a more unified methodology is realized. In this work, we develop an approach called Orthogonal-defect Classification for Assessing Software Reliability (ORCAS) to quantify probabilities of various software failure modes in a DI&C system. The method utilizes accepted industry methodologies for software quality assurance that are also verified by experimental or mathematical formulations. In essence, the approach combines a semantic failure classification model with a reliability growth model to predict (and quantify) the potential failure modes of a DI&C software system. The semantic classification model is used to address the question: How do latent defects in software contribute to different software failure root causes? The use of reliability growth models is then used to address the question: Given the connection between latent defects and software failure root causes, how can we quantify the reliability of the software? A case study was conducted on a representative I&C platform (ChibiOS) running a smart sensor acquisition software developed by Virginia Commonwealth University (VCU). The testing and evidence collection guidance in ORCAS was applied, and defects were uncovered in the software. Qualitative evidence, such as condition coverage, was used to gauge the completeness and trustworthiness of the assessment while quantitative evidence was used to determine the software failure probabilities. The reliability of the software was then estimated and compared to existing operational data of the sensor device. It is demonstrated that by using ORCAS, a semantic reasoning framework can be developed to justify software reliability (or unreliability) while still leveraging the strength of the existing methods.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Coupling High-Throughput Experiments and Regression Algorithms to Optimize PGM-Free ORR Electrocatalyst Synthesis

Over the past decades, significant improvement has been achieved in the performance of platinum group metal-free (PGM-free) materials as an alternative to Pt-based electrocatalysts for oxygen reduction reaction (ORR). However, further progress in ORR activity requires evaluation of precursors and synthesis approaches. In response to this challenge, we generated a first of its kind experimental data set of 36 samples using high-throughput synthesis and activity measurements. Several control parameters (e.g., Fe precursor identity, the precursor content, and pyrolysis temperature) were varied. We then developed several state-of-the-art machine learning (ML) based regression models to predict ORR activity, dependent on selected synthesis variables. Through an iterative algorithm, higher prediction accuracy (smaller root-mean-square error) was achieved. We identified that gradient boosting regression (GBR) and support vector regression (SVR), among several methods, work best for this data set. Aided by our ML-based surrogate models, we decided to alter catalyst synthesis conditions, which resulted in a 36% increase in measured ORR activity in comparison to the maximum ORR mass activity value of 21.9 A/g catalyst in the original data set. Overall, this combined experiment and machine learning approach represents a promising path forward toward developing highly efficient next-generation ORR electrocatalysts and, more generally, functional materials.

25 ENERGY STORAGE↗
DOE OSTI2021

Dynamic P-Q Capability and Abnormal Operation Analysis of a Wind Turbine with Doubly-Fed Induction Generator

After several accidents occurred in wind farms in the United States and around the world, the subsynchronous oscillation (SSO) issues in grid-connected wind farms have gained serious attention. Particularly, these issues have caused significant challenges to wind turbines with a doubly fed induction generator (DFIG) because it is connected to the grid via both its stator and rotor paths. Traditionally, a P-Q capability chart is utilized to assure the safe operation boundary for a synchronous generator. But the energy conversion characteristics of a DFIG wind turbine are completely different. A critical factor to affect the reliable operation of a DFIG wind turbine is the rated current and pulse-width modulation (PWM) saturation constraints of its power converters. These constraint conditions can be affected by the wind turbine rotating speed and grid conditions. However, a detailed study of DFIG P-Q capability from these perspectives has not been conducted, which has hindered adequate understanding of many abnormal wind turbine operations reported in the literature and the development of advanced control technologies to overcome the challenges. The proposed study in this paper considers vector control implementation to DFIG power electronic converters in the dq reference frame, and the models and algorithms developed for the P-Q capability study have addressed specific DFIG power converter constraints that are different from those of a traditional synchronous generator. The paper especially focuses on exploring the dynamic natures of DFIG P-Q capability under uncertain and variable conditions to explore the root causes of many abnormal operations of DFIG wind turbines reported in the literature. The proposed study is validated through an electromagnetic transient simulation model of a grid-connected DFIG wind turbine. The proposed study has the potential to lead to the development of new DFIG control technologies that can help overcome the challenges for many abnormal operations of DFIG wind turbines.

17 WIND ENERGY↗
DOE OSTI2024

Transgenic Sugarcane–Oilcane: An Alternative Feedstock for the Production of Drop-in Fuel and Value-Added Bioproducts

The utilization of plants and other agricultural produce can partly offset petroleum dependency for energy requirements and can potentially provide sustainable solutions to global environmental problems. To this end, synthetic biology has shown great potential in developing transgenic bioenergy grasses such as sugarcane, sorghum, miscanthus, and energy cane that hyperaccumulate energy-rich lipid molecules in their vegetative tissues, such as leaves, stems, and roots. These perennial high-biomass transgenic C4 grasses are not targeted to grow on prime agricultural land and can be dedicatedly used to produce biofuels and other value-added bioproducts. Recently, sugarcane has been metabolically engineered to sequestrate carbon from juice towards biosynthesis of lipid molecules in the vegetative tissues. Transgenic sugarcane is referred to as “oilcane.” Transgenic sugarcane–oilcane has improved energy density due to an elevated lipid content in the vegetative tissues. Transgenic bioenergy crops cater to both cellulosic sugars and vegetative lipids, hence resulting in higher biofuel yield (biodiesel/renewable diesel and bioethanol) per unit area of cultivable land. In this book chapter, the genetic engineering of sugarcane and its bioprocessing are discussed to illustrate its development and use as an alternative feedstock for the production of biofuels (bio-jet fuel, biodiesel, and bioethanol) and value-added bioproducts.

Maitra, Shraddha↗
DOE OSTI2022

Essential amino acids in the Plant-Conserved and Class-Specific Regions of cellulose synthases

The Plant-Conserved Region (P-CR) and the Class-Specific Region (CSR) are two plant-unique sequences in the catalytic core of cellulose synthases (CESAs) for which specific functions have not been established. Here, we used site-directed mutagenesis to replace amino acids and motifs within these sequences predicted to be essential for assembly and function of CESAs. We developed an in vivo method to determine the ability of mutated CesA1 transgenes to complement an Arabidopsis (Arabidopsis thaliana) temperature-sensitive root-swelling1 (rsw1) mutant. Replacement of a Cys residue in the CSR, which blocks dimerization in vitro, rendered the AtCesA1 transgene unable to complement the rsw1 mutation. Examination of the CSR sequences from 33 diverse angiosperm species showed domains of high-sequence conservation in a class-specific manner but with variation in the degrees of disorder, indicating a nonredundant role of the CSR structures in different CESA isoform classes. The Cys residue essential for dimerization was not always located in domains of intrinsic disorder. Expression of AtCesA1 transgene constructs, in which Pro417 and Arg453 were substituted for Ala or Lys in the coiled-coil of the P-CR, were also unable to complement the rsw1 mutation. Despite an expected role for Arg457 in trimerization of CESA proteins, AtCesA1 transgenes with Arg457Ala mutations were able to fully restore the wild-type phenotype in rsw1. Our data support that Cys662 within the CSR and Pro417 and Arg453 within the P-CR of Arabidopsis CESA1 are essential residues for functional synthase complex formation, but our data do not support a specific role for Arg457 in trimerization in native CESA complexes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Variation in forest root image annotation by experts, novices, and AI

Abstract Background The manual study of root dynamics using images requires huge investments of time and resources and is prone to previously poorly quantified annotator bias. Artificial intelligence (AI) image-processing tools have been successful in overcoming limitations of manual annotation in homogeneous soils, but their efficiency and accuracy is yet to be widely tested on less homogenous, non-agricultural soil profiles, e.g., that of forests, from which data on root dynamics are key to understanding the carbon cycle. Here, we quantify variance in root length measured by human annotators with varying experience levels. We evaluate the application of a convolutional neural network (CNN) model, trained on a software accessible to researchers without a machine learning background, on a heterogeneous minirhizotron image dataset taken in a multispecies, mature, deciduous temperate forest. Results Less experienced annotators consistently identified more root length than experienced annotators. Root length annotation also varied between experienced annotators. The CNN root length results were neither precise nor accurate, taking ~ 10% of the time but significantly overestimating root length compared to expert manual annotation ( p = 0.01). The CNN net root length change results were closer to manual ( p = 0.08) but there remained substantial variation. Conclusions Manual root length annotation is contingent on the individual annotator. The only accessible CNN model cannot yet produce root data of sufficient accuracy and precision for ecological applications when applied to a complex, heterogeneous forest image dataset. A continuing evaluation and development of accessible CNNs for natural ecosystems is required.

Handy, Grace↗
DOE OSTI2022

Generalized Bootstrap AMG and AIR-AMG for coupled PDE systems with a focus on spacetime discretizatoins

The Pennsylvania State University (“Subcontractor”) has worked on the design of new algebraic, parallel, multilevel methods that obtain the full space and time solution of systems of PDEs. In particular, the PI and his collaborators explored semi-intrusive approaches based on algebraic multigrid (AMG). The focus of the research has been on the development of these techniques for the Euler equations in 1d and 2d. The overall research focused on the development of adaptive AIR (approximate ideal restriction) AMG solvers for these problems. The PI also explored the use of smoothed aggregation and root-node energy-based AMG solvers for these problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A robust approach to Gaussian process implementation

Abstract. Gaussian process (GP) regression is a flexible modeling technique used to predict outputs and to capture uncertainty in the predictions. However, the GP regression process becomes computationally intensive when the training spatial dataset has a large number of observations. To address this challenge, we introduce a scalable GP algorithm, termed MuyGPs, which incorporates nearest-neighbor and leave-one-out cross-validation during training. This approach enables the evaluation of large spatial datasets with state-of-the-art accuracy and speed in certain spatial problems. Despite these advantages, conventional quadratic loss functions used in the MuyGPs optimization, such as root mean squared error (RMSE), are highly influenced by outliers. We explore the behavior of MuyGPs in cases involving outlying observations and, subsequently, develop a robust approach to handle and mitigate their impact. Specifically, we introduce a novel leave-one-out loss function based on the pseudo-Huber function (LOOPH) that effectively accounts for outliers in large spatial datasets within the MuyGPs framework. Our simulation study shows that the LOOPH loss method maintains accuracy despite outlying observations, establishing MuyGPs as a powerful tool for mitigating unusual observation impacts in the large data regime. In the analysis of US ozone data, MuyGPs provides accurate predictions and uncertainty quantification, demonstrating its utility in managing data anomalies. Through these efforts, we advance the understanding of GP regression in spatial contexts.

Mukangango, Juliette↗
DOE OSTI2022

Active rotor coning for a 25 MW downwind offshore wind turbine

A two-bladed downwind turbine system was upscaled from 13.2 MW to 25 MW by redesigning aerodynamics, structures, and controls. In particular, three 25-MW rotors were developed, and the final version is a fully redesigned model of the original rotor. Despite their radically large sizes, it was found that these 25-MW turbine rotors satisfy this limited set of structural design drivers at the rated condition and that larger blade lengths are possible with conewise load-alignment. In addition, flapwise morphing (varying the cone angle with a wind-speed schedule) was investigated to minimize mean and fluctuating blade root bending loads using steady inflow proxies for the maximum and lifetime damage equivalent load moments. Compared to the fixed coned rotor case, morphing can provide an Annual Energy Production (AEP) increase of 6%, and the maximum blade root flapwise bending moment increases 21% (still under the constraint, i.e., 10% of the ultimate moments) as a trade-off. The resulting series of 25-MW rotors can be a valuable baseline for further development and assessment of ultra-large-scale wind turbines.

17 WIND ENERGY↗
DOE OSTI2022

Revenue prediction for integrated renewable energy and energy storage system using machine learning techniques

Revenue estimation for integrated renewable energy and energy storage systems is important to support plant owners or operators’ decisions in battery sizing selection that leads to maximized financial performances. A common approach to optimizing revenues of a hybrid hydro and energy storage system is using mixed-integer linear programming (MILP). Although MILP models can provide accurate production cost estimations, they are typically very computationally expensive. To provide a fast yet accurate first-step information to hydropower plant owners or operators who consider integrating energy storage systems, we propose an innovative approach to predicting optimal revenues of an integrated energy generation and storage system. In this study, we examined the performance of two prediction techniques: Generalized Additive Models (GAMs) and machine learning (ML) models developed based on artificial neural networks (ANN). Predictive equations and models are generated based on optimized solutions from a market participation optimization model, the Conventional Hydropower Energy and Environmental Resource System (CHEERS) model. The two predicting techniques reduce the computational time to evaluate annual revenue for one set of battery configurations from 3 h to 1 to 4 min per run while also being implementable with significantly less data. The model validation prediction errors of developed GAMs and ML models are generally below 5%; for model testing predictions, the ML models consistently outperform the regression equations in terms of root mean square errors. This new approach allows plant owners, operators, or potential investors to quickly access multiple battery configurations under different energy generation and market scenarios. This new revenue prediction method will therefore help reduce the barriers, and thereby promoting the deployment of battery hybridization with existing renewable energy sources.

13 HYDRO ENERGY↗
DOE OSTI2024

A novel conditional generative model for efficient ensemble forecasts of state variables in large-scale geological carbon storage

Integrating monitoring data to efficiently update reservoir pressure and CO 2 plume distribution forecasts presents a significant challenge in geological carbon storage (GCS) applications. Inverse modeling techniques are commonly used to fuse observational data and refine reservoir model parameters, thereby improving state variable forecasts. However, these techniques often rely on linear or Gaussian assumptions, which can limit their effectiveness in accurately predicting state variables. Moreover, simulating large-scale three-dimensional (3D) GCS problems is computationally expensive, making iterative runs in inverse problems prohibitive. To address these challenges, we propose a conditional generative model utilizing the score-based diffusion method for real-time 3D pressure and saturation field distribution predictions. Our approach involves solving the score function with a mini-batch-based Monte Carlo estimator to generate labeled data. This data is subsequently employed to train a fully connected neural network, enabling it to learn the conditional sample generator within a supervised learning framework. This method enables the rapid generation of a large ensemble of predictions, facilitating comprehensive uncertainty quantification of state variables. Here we applied our method to forecast the dynamic 3D distributions of pressure and saturation fields over a 30-year injection period. The statistical assessment with low root mean square error (RMSE) values demonstrates that our method can accurately predict the spatiotemporal distributions of both pressure and saturation fields. Moreover, the developed conditional generative model shows high computational efficiency by generating 100 ensemble forecasts of 3D state variables in less than 10 min. The consistency between ensemble averages and ground truth values further illustrates the model’s capability to capture state variable dynamics during the CO 2 plume injection process. Notably, the ground truth values fall within the ensemble forecasts, indicating that our uncertainty quantification effectively captures variability and potential noise in the observations. Thus, the developed conditional generative model proves to be a more efficient, accurate, and practical tool for GCS applications, facilitating timely risk analysis and informed decision-making.

58 GEOSCIENCES↗
DOE OSTI2023

Morphologically-Directed Raman Spectroscopy as an Analytical Method for Subvisible Particle Characterization in Therapeutic Protein Product Quality

Abstract Subvisible particles (SVPs) are a critical quality attribute of injectable therapeutic proteins (TPs) that needs to be controlled due to potential risks associated with drug product quality. The current compendial methods routinely used to analyze SVPs for lot release provide information on particle size and count. However, chemical identification of individual particles is also important to address root-cause analysis. Herein, we introduce Morphologically-Directed Raman Spectroscopy (MDRS) for SVP characterization of TPs. The following particles were used for method development: (1) polystyrene microspheres, a traditional standard used in industry; (2) photolithographic (SU-8); and (3) ethylene tetrafluoroethylene (ETFE) particles, candidate reference materials developed by NIST. In our study, MDRS rendered high-resolution images for the ETFE particles (> 90%) ranging from 19 to 100 μm in size, covering most of SVP range, and generated comparable morphology data to flow imaging microscopy. Our method was applied to characterize particles formed in stressed TPs and was able to chemically identify individual particles using Raman spectroscopy. MDRS was able to compare morphology and transparency properties of proteinaceous particles with reference materials. The data suggests MDRS may complement the current TPs SVP analysis system and product quality characterization workflow throughout development and commercial lifecycle.

Science & Technology - Other Topics↗
DOE OSTI2020

Advances in catalytic production processes of biomass-derived vinyl monomers

Plastic industries currently source the majority of monomers from crude oil substances. Although we have witnessed a significant interest of companies in the utilization of sustainable feedstock materials for bio-based compound synthesis in the past decade, the transition to photosynthetic or chemosynthetic plant-based production in a circular carbon economy is largely slow due to complex biomass processing, costs and related reaction factors. The upgrade of the separated platform chemicals with a deeper supporting understanding of processes, models and a root–cause analysis about the underlying distribution challenges of engineered transformation mechanisms, catalyzed conversion pathways and selectivity would be beneficial to advance applied scientific development, bio-refining and manufacturing output amount. This review provides a summary, assessment and perspective for three important polymer-forming vinyl molecules, i.e. acrylic acid, methacrylic acid and styrene. These provide a backbone to produce acrylates, polystyrene, resins, rubbers, protective surface coatings, adhesives, textiles and other obtained copolymers. A succinct analytical overview on the thermo-catalytic intermediate routes for property-wise drop-in alternatives is presented. Sugars, acrolein, allyl alcohol, ethylbenzene, glycerol, 3-hydroxypropionic acid, isobutene, itaconic acid and lactic acid are considered as the main starting reactants. Catalysts span mixed metal oxides, silicates, native or impregnated zeolite frameworks (HBEA, HZSM, and MFI) and hetero-poly acids as well as homogeneous base hydroxides or platinum group metals, supported on carbon, alumina and sulfates. Finally, the article concludes with a brief state sum-up of the results, topics and opportunities for systematic future research or scaling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Improved node culture methods for rapid vegetative propagation of switchgrass (Panicum virgatum L.)

Abstract Background Switchgrass ( Panicum virgatum L.) is an important bioenergy and forage crop. The outcrossing nature of switchgrass makes it infeasible to maintain a genotype through sexual propagation. Current asexual propagation protocols in switchgrass have various limitations. An easy and highly-efficient vegetative propagation method is needed to propagate large natural collections of switchgrass genotypes for genome-wide association studies (GWAS). Results Micropropagation by node culture was found to be a rapid method for vegetative propagation of switchgrass. Bacterial and fungal contamination during node culture is a major cause for cultural failure. Adding the biocide, Plant Preservative Mixture (PPM, 0.2%), and the fungicide, Benomyl (5 mg/l), in the incubation solution after surface sterilization and in the culture medium significantly decreased bacterial and fungal contamination. In addition, “shoot trimming” before subculture had a positive effect on shoot multiplication for most genotypes tested. Using the optimized node culture procedure, we successfully propagated 330 genotypes from a switchgrass GWAS panel in three separate experiments. Large variations in shoot induction efficiency and shoot growth were observed among genotypes. Separately, we developed an in planta node culture method by stimulating the growth of aerial axillary buds into shoots directly on the parent plants, through which rooted plants can be generated within 6 weeks. By circumventing the tissue culture step and avoiding application of exterior hormones, the in planta node culture method is labor- and cost-efficient, easy to master, and has a high success rate. Plants generated by the in planta node culture method are similar to seedlings and can be used directly for various experiments. Conclusions In this study, we optimized a switchgrass node culture protocol by minimizing bacterial and fungal contamination and increasing shoot multiplication. With this improved protocol, we successfully propagated three quarters of the genotypes in a diverse switchgrass GWAS panel. Furthermore, we established a novel and high-throughput in planta node culture method. Together, these methods provide better options for researchers to accelerate vegetative propagation of switchgrass.

Wang, Yongqin↗
DOE OSTI2020

High-frequency and Vertical Variations of Soil Moisture and their Impact on Ecosystem Dynamics

The vertical distribution of subsurface moisture is a major control on ecosystem dynamics and climate. The Mediterranean climate of California, with wet winters and dry summers, provides a testbed for understanding how transpiration is sustained through extended dry periods. The project developed and tested a new subsurface hydrology model that includes a parameterization of stochastic hydraulic conductivity to capture fracture flow, a diagnostic model of rooting profile that depends on the subsurface hydrology, and a representation of hydraulic redistribution by plant roots. A sensitivity study shows that it is plants with rooting depths >8 m that could access the subsurface moisture. The insight is applied to understand differential tree mortality during the 2012-2016 drought in California. The analysis shows that greater mortality trees on plutonic rocks than on sedimentary rocks. The study demonstrates, for the first time, the importance of lithology for the distribution and health of ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Simulating low-energy neutrino interactions with MARLEY

Monte Carlo event generators are a critical tool for the interpretation of data obtained by neutrino experiments. Several modern event generators are available which are well-suited to the GeV energy scale used in studies of accelerator neutrinos. However, theoretical modeling differences make their immediate application to lower energies difficult. In this paper, I present a new event generator, MARLEY, which is designed to better address the simulation needs of the low-energy (tens of MeV and below) neutrino community. The code is written in C++14 with an optional interface to the popular ROOT data analysis framework. The current release of MARLEY (version 1.2.0) emphasizes simulations of the reaction 40 Ar ( ν e , e − ) 40 K ⁎ but is extensible to other channels with suitable user input. This paper provides detailed documentation of MARLEY's implementation and usage, including guidance on how generated events may be analyzed and how MARLEY may be interfaced with external codes such as Geant4. Further information about MARLEY is available on the official website at http://www.marleygen.org. Program title:MARLEY 1.2.0 CPC Library link to program files:https://doi.org/10.17632/4v7zxnc8j3.1 Developer's respository link:http://github.com/MARLEY-MC/marley Code Ocean capsule:https://codeocean.com/capsule/9868179 Licensing provisions: GNU General Public License 3.0 Programming language: C++14 External routines/libraries used: GNU Scientific Library [1,2] (required), ROOT [3,4] (optional) Nature of problem: Simulation of neutrino-nucleus scattering events at energies of tens-of-MeV and below Solution method: Initial two-to-two scattering kinematics are sampled using the allowed approximation differential cross section and tables of precomputed nuclear matrix elements. Subsequent de-excitations of the remnant nucleus are simulated using a Monte Carlo implementation of the Hauser-Feshbach statistical model and tabulated γ-ray decay schemes for discrete nuclear levels. Additional comments including restrictions and unusual features: Input data are provided with the code that are suitable for producing simulations of the charged-current reaction 40 Ar ( ν e , e − ) 40 K ⁎ , coherent elastic neutrino-nucleus scattering on spin-zero target nuclei, and neutrino-electron elastic scattering on any atomic target. Preparation of new reaction input files (whose format is documented in Appendix B) would enable other reaction channels and nuclear targets to be handled by the existing code framework. Although there is no maximum neutrino energy enforced by the code itself, realistic neutrino-nucleus scattering events may be generated up to roughly 50 MeV. Above this energy, the effects of forbidden nuclear transitions, which are neglected in the current treatment of the cross sections (see section 2.1), become increasingly important.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗