Search NASA⌕ Search

SEARCH · Search NASA

Results for “differentiable programming”

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

Refining the GENEVA method for Higgs boson production via gluon fusion

We describe a number of improvements to the Geneva method for matching NNLO calculations to parton shower programs. In particular, we detail changes to the resummed calculation used in the matching procedure, including disentangling the cross section dependence on factorisation and beam scales, and an improved treatment of timelike logarithms. We also discuss modifications in the implementation of the splitting functions which serve to make the resummed calculation differential in the higher multiplicity phase space. These changes improve the stability of the numerical cancellation of the nonsingular term at small values of the resolution parameter. As a case study, we consider the gluon-initiated Higgs boson production process gg → H. We validate the NNLO accuracy of our predictions against independent calculations, and compare our showered and hadronised results with recent data taken at the ATLAS and CMS experiments in the diphoton decay channel, finding good agreement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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↗
DOE OSTI2017

Sandia QIS Capabilities.

Sandia National Laboratories has developed a broad set of capabilities in quantum information science (QIS), including elements of quantum computing, quantum communications, and quantum sensing. The Sandia QIS program is built atop unique DOE investments at the laboratories, including the MESA microelectronics fabrication facility, the Center for Integrated Nanotechnologies (CINT) facilities (joint with LANL), the Ion Beam Laboratory, and ASC High Performance Computing (HPC) facilities. Sandia has invested $75 M of LDRD funding over 12 years to develop unique, differentiating capabilities that leverage these DOE infrastructure investments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Loss of Vacuum Condition for IER 479 [Slides]

Conservative thermal transient calculations and simulations show that a loss of vacuum condition would lead to a small (~6 °C) temperature differential across the Jemima Plates. Jemima plates would experience very low thermal stresses, well below the temperature differences of concern (>40 °C).

42 ENGINEERING↗
DOE OSTI2017

Effect of compensatory immigration on the genetic structure of coyotes

Despite efforts to reduce their impacts on livestock and native ungulates within the southeastern United States, coyotes (Canis latrans) are able to recover from control programs. It is unknown how coyotes compensate for high mortality following trapping, so there is great interest to identify methods that can provide insight into coyote response to intensive trapping. To investigate if population genetic tools can decipher how coyotes recover from intensive trapping, we combined an empirical test of how genetic differentiation, diversity, and familial structure changed following trapping on the Savannah River Site (SRS), South Carolina with spatially explicit genetic simulations.The empirical dataset found that pre- and post-trapping periods had similar genetic diversities and were not genetically differentiated as expected by either compensatory reproduction or immigration from a single genetic source. The post-trapping coyote populations exhibited weaker signatures of philopatry with little evidence for increased dispersal distances of young coyotes, which suggests immigration caused a decrease in familial structure. Our simulations indicated that spatial autocorrelation coefficients and observed heterozygosities change as immigration increases while population differentiation, allelic richness, and displacement distances do not. Collectively, our results suggest that coyotes recover from intensive trapping via both reproduction and immigration, which likely makes preventing compensation difficult. Furthermore, monitoring post-trapping populations may offer more insight into maximizing the effectiveness of control efforts, and based on our simulations, population genetics can provide critical information about the amount of compensatory immigration following trapping.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Functional importance of cardiac enhancer-associated noncoding RNAs in heart development and disease

We report here that the key information processing units within gene regulatory networks are enhancers. Enhancer activity is associated with the production of tissue-specific noncoding RNAs, yet the existence of such transcripts during cardiac development has not been established. Using an integrated genomic approach, we demonstrate that fetal cardiac enhancers generate long noncoding RNAs (lncRNAs) during cardiac differentiation and morphogenesis. Enhancer expression correlates with the emergence of active enhancer chromatin states, the initiation of RNA polymerase II at enhancer loci and expression of target genes. Orthologous human sequences are also transcribed in fetal human hearts and cardiac progenitor cells. Through a systematic bioinformatic analysis, we identified and characterized, for the first time, a catalog of lncRNAs that are expressed during embryonic stem cell differentiation into cardiomyocytes and associated with active cardiac enhancer sequences. RNA-sequencing demonstrates that many of these transcripts are polyadenylated, multi-exonic long noncoding RNAs. Moreover, knockdown of two enhancer-associated lncRNAs resulted in the specific downregulation of their predicted target genes. Interestingly, the reactivation of the fetal gene program, a hallmark of the stress response in the adult heart, is accompanied by increased expression of fetal cardiac enhancer transcripts. Altogether, these findings demonstrate that the activity of cardiac enhancers and expression of their target genes are associated with the production of enhancer-derived lncRNAs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Kinetic Model of HoxEFU reduction by NADH [SWR-26-087]

This repository is used to release code generated for manuscripts on the Photosynthetic Energy Transduction core program. This code simulates the reduction of HoxEFU by NADH. The electro transfer rate constants for the simulation are specified in the .csv files. The two .csv files correspond tot he two models described in Dawson et al. Cell. Rep. Phys. Sci. 2026. The code utilizes a chemical master equation, a set of differential equations, defining the time evolution of the oxidation and reduction kinetics of NAD+, NADH, a FMN flavin, and a set of iron sulfur clusters. The kinetics of HoxEFU reduction by NADH are evaluated by numerical integration of the chemical master equation using a variable-time-step Runge-Kutta algorithm.

Dahl, Peter [National Laboratory of the Rockies (N↗
DOE OSTI2011

Precise predictions for (nonstandard) W{gamma}+ jet production

We report on a detailed investigation of the next-to-leading-order QCD corrections to W{gamma}+jet production at the Tevatron and the LHC using a fully flexible parton-level Monte Carlo program. We include the full leptonic decay of the W, taking into account all off-shell and finite width effects, as well as nonstandard WW{gamma} couplings. We find particularly sizable corrections for the currently allowed parameter range of anomalous couplings imposed by LEP data. In total, the next-to-leading-order differential distributions reveal a substantial phase-space dependence of the corrections, leaving considerable sensitivity to anomalous couplings beyond scale uncertainty at large momentum transfers in the anomalous vertex.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2012

Functionality and operation of fluoroscopic automatic brightness control/automatic dose rate control logic in modern cardiovascular and interventional angiography systems: A Report of Task Group 125 Radiography/Fluoroscopy Subcommittee, Imaging Physics Committee, Science Council

Task Group 125 (TG 125) was charged with investigating the functionality of fluoroscopic automatic dose rate and image quality control logic in modern angiographic systems, paying specific attention to the spectral shaping filters and variations in the selected radiologic imaging parameters. The task group was also charged with describing the operational aspects of the imaging equipment for the purpose of assisting the clinical medical physicist with clinical set-up and performance evaluation. Although there are clear distinctions between the fluoroscopic operation of an angiographic system and its acquisition modes (digital cine, digital angiography, digital subtraction angiography, etc.), the scope of this work was limited to the fluoroscopic operation of the systems studied. The use of spectral shaping filters in cardiovascular and interventional angiography equipment has been shown to reduce patient dose. If the imaging control algorithm were programmed to work in conjunction with the selected spectral filter, and if the generator parameters were optimized for the selected filter, then image quality could also be improved. Although assessment of image quality was not included as part of this report, it was recognized that for fluoroscopic imaging the parameters that influence radiation output, differential absorption, and patient dose are also the same parameters that influence image quality. Therefore, this report will utilize the terminology ''automatic dose rate and image quality'' (ADRIQ) when describing the control logic in modern interventional angiographic systems and, where relevant, will describe the influence of controlled parameters on the subsequent image quality. A total of 22 angiography units were investigated by the task group and of these one each was chosen as representative of the equipment manufactured by GE Healthcare, Philips Medical Systems, Shimadzu Medical USA, and Siemens Medical Systems. All equipment, for which measurement data were included in this report, was manufactured within the three year period from 2006 to 2008. Using polymethylmethacrylate (PMMA) plastic to simulate patient attenuation, each angiographic imaging system was evaluated by recording the following parameters: tube potential in units of kilovolts peak (kVp), tube current in units of milliamperes (mA), pulse width (PW) in units of milliseconds (ms), spectral filtration setting, and patient air kerma rate (PAKR) as a function of the attenuator thickness. Data were graphically plotted to reveal the manner in which the ADRIQ control logic responded to changes in object attenuation. There were similarities in the manner in which the ADRIQ control logic operated that allowed the four chosen devices to be divided into two groups, with two of the systems in each group. There were also unique approaches to the ADRIQ control logic that were associated with some of the systems, and these are described in the report. The evaluation revealed relevant information about the testing procedure and also about the manner in which different manufacturers approach the utilization of spectral filtration, pulsed fluoroscopy, and maximum PAKR limitation. This information should be particularly valuable to the clinical medical physicist charged with acceptance testing and performance evaluation of modern angiographic systems.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Thermal Effects on Far-Field Distributed Acoustic Strain-Rate Sensors

Summary Fiber-optic cables cemented outside of the casing of an unconventional well measure crosswell strain changes during fracturing of neighboring wells with low-frequency distributed acoustic sensing (LF-DAS). As a hydraulic fracture intersects an observation well instrumented with fiber-optic cables, the fracture fluid injected at ambient temperatures can cool a section of the sensing fiber. Often, LF-DAS and distributed temperature sensing (DTS) cables are run in tandem, enabling the detection of such cooling events. The increasing use of LF-DAS for characterizing unconventional hydraulic fracture completions demands an investigation of the effects of temperature on the measured strain response by LF-DAS. Researchers have demonstrated that LF-DAS can be used to extract the temporal derivative of temperature for use as a differential-temperature-gradient sensor. However, differential-temperature-gradient sensing is predicated on the ability to filter strain components out of the optical signal. In this work, beginning with an equation for optical phase shift of LF-DAS signals, a model relating strain, temperature, and optical phase shift is explicitly developed. The formula provides insights into the relative strength of strain and temperature effects on the phase shift. The uncertainty in the strain-rate measurements due to thermal effects is estimated. The relationship can also be used to quantify uncertainties in differential-temperature-gradient sensors due to strain perturbations. Additionally, a workflow is presented to simulate the LF-DAS response accounting for both strain and temperature effects. Hydraulic fracture geometries are generated with a 3D fracture simulator for a multistage unconventional completion. The fracture width distributions are imported by a displacement discontinuity method (DDM) program to compute the strain rates along an observation well. An analytic model is used to approximate the temperature in the fracture. Using the derived formulae for optical phase shift, the model outputs are then used to compute the LF-DAS response at a fiber-optic cable, enabling the generation of waterfall plots including both strain and thermal effects. The model results suggest that before, during, and immediately following a fracture intersecting a well instrumented with fiber, the strain on the fiber drives the LF-DAS signal. However, at later times, as completion fluid cools the observation well, the temperature component of the LF-DAS signal can be equal to or exceed the strain component. The modeled results are compared to a published field case in an attempt to enhance the interpretation of LF-DAS waterfall plots. Finally, we propose a sensing configuration to identify the events when “wet fractures” (fractures with fluids) intersect the observation well.

Engineering↗
DOE OSTI2025

Drought adaptation index (DAI) based on BLUP as a selection approach for drought-resilient switchgrass germplasm

This study introduces a Drought Adaptation Index (DAI), derived from Best Linear Unbiased Prediction (BLUP), as a method to assess drought resilience in switchgrass (Panicum virgatum L.). A panel of 404 genotypes was evaluated under drought-stressed (CV) and well-watered (UC) conditions over four consecutive years (2019–2022). BLUP-estimated biomass yields were used to calculate the DAI, which enabled classification of genotypes into four adaptation groups: very well-adapted, well-adapted, adapted, and unadapted. The DAI was compared with conventional drought tolerance indices, including the Stress Susceptibility Index (SSI), Stress Tolerance Index (STI), Geometric Mean Productivity (GMP), and Yield Stability Index (YSI). Correlation analyses demonstrated strong agreement between DAI and these indices, supporting its validity and consistency. Biplot analyses using the Genotype plus Genotype-by-Environment Interaction (GGE) and Additive Main Effects and Multiplicative Interaction (AMMI) models revealed significant genotype-by-environment interactions (GEI) and identified J222.A, J463.A, and J295.A. A as high-performing genotypes, with J222.A exhibiting greater yield stability across treatments and years. Additionally, DAI isoline curves provided a graphical representation of differential genotype performance under drought and control conditions. These visualizations aided in distinguishing genotypes with stable and superior biomass yield across contrasting environments. Overall, the BLUP-based DAI is a robust and practical selection tool that improves the accuracy of identifying drought-resilient, high-yielding switchgrass genotypes. Its integration into breeding programs offers a comprehensive framework for improving biomass productivity and stress adaptation under variable climatic conditions. The application of DAI supports the development of climate-resilient cultivars and contributes to sustainable bioenergy and forage production systems.

BLUP↗
DOE OSTI2012

AVESTAR Center for Operational Excellence of Electricity Generation Plants

To address industry challenges in attaining operational excellence for electricity generation plants, the U.S. Department of Energy’s (DOE) National Energy Technology Laboratory (NETL) has launched a world-class facility for Advanced Virtual Energy Simulation Training and Research (AVESTARTM). This presentation will highlight the AVESTARTM Center simulators, facilities, and comprehensive training, education, and research programs focused on the operation and control of high-efficiency, near-zero-emission electricity generation plants. The AVESTAR Center brings together state-of-the-art, real-time, high-fidelity dynamic simulators with full-scope operator training systems (OTSs) and 3D virtual immersive training systems (ITSs) into an integrated energy plant and control room environment. AVESTAR’s initial offering combines--for the first time--a “gasification with CO2 capture” process simulator with a “combined-cycle” power simulator together in a single OTS/ITS solution for an integrated gasification combined cycle (IGCC) power plant with carbon dioxide (CO2) capture. IGCC systems are an attractive technology option for power generation, especially when capturing and storing CO2 is necessary to satisfy emission targets. The AVESTAR training program offers a variety of courses that merge classroom learning, simulator-based OTS learning in a control-room operations environment, and immersive learning in the interactive 3D virtual plant environment or ITS. All of the courses introduce trainees to base-load plant operation, control, startups, and shutdowns. Advanced courses require participants to become familiar with coordinated control, fuel switching, power-demand load shedding, and load following, as well as to problem solve equipment and process malfunctions. Designed to ensure work force development, training is offered for control room and plant field operators, as well as engineers and managers. Such comprehensive simulator-based instruction allows for realistic training without compromising worker, equipment, and environmental safety. It also better prepares operators and engineers to manage the plant closer to economic constraints while minimizing or avoiding the impact of any potentially harmful, wasteful, or inefficient events. The AVESTAR Center is also used to augment graduate and undergraduate engineering education in the areas of process simulation, dynamics, control, and safety. Students and researchers gain hands-on simulator-based training experience and learn how the commercial-scale power plants respond dynamically to changes in manipulated inputs, such as coal feed flow rate and power demand. Students also analyze how the regulatory control system impacts power plant performance and stability. In addition, students practice start-up, shutdown, and malfunction scenarios. The 3D virtual ITSs are used for plant familiarization, walk-through, equipment animations, and safety scenarios. To further leverage the AVESTAR facilities and simulators, NETL and its university partners are pursuing an innovative and collaborative R&D program. In the area of process control, AVESTAR researchers are developing enhanced strategies for regulatory control and coordinated plant-wide control, including gasifier and gas turbine lead, as well as advanced process control using model predictive control (MPC) techniques. Other AVESTAR R&D focus areas include high-fidelity equipment modeling using partial differential equations, dynamic reduced order modeling, optimal sensor placement, 3D virtual plant simulation, and modern grid. NETL and its partners plan to continue building the AVESTAR portfolio of dynamic simulators, immersive training systems, and advanced research capabilities to satisfy industry’s growing need for training and experience with the operation and control of clean energy plants. Future dynamic simulators under development include natural gas combined cycle (NGCC) and supercritical pulverized coal (SCPC) plants with post-combustion CO2 capture. These dynamic simulators are targeted for use in establishing a Virtual Carbon Capture Center (VCCC), similar in concept to the DOE’s National Carbon Capture Center for slipstream testing. The VCCC will enable developers of CO2 capture technologies to integrate, test, and optimize the operation of their dynamic capture models within the context of baseline power plant dynamic models. The objective is to provide hands-on, simulator-based “learn-by-operating” test platforms to accelerate the scale-up and deployment of CO2 capture technologies. Future AVESTAR plans also include pursuing R&D on the dynamics, operation, and control of integrated electricity generation and storage systems for the modern grid era. Special emphasis will be given to combining load-following energy plants with renewable and distributed generating supplies and fast-ramping energy storage systems to provide near constant baseload power.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Low Cost Glass-Ceramic Matrix Composite Heat Exchanger

As part of ARPA-E’s High Intensity Thermal Exchange through Materials and Manufacturing Processes (HITEMMP) program, this project sought to develop novel heat exchanger (HX) capabilities to enable efficient and power dense power generation cycles. This class of HX comes under the category of ceramic/composite materials with the higher temperature goal in the program of ≥1100 °C inlet temperature operation. The enabling capability of this effort is the use of glass-ceramic matrix composite (GCMC) material which provides the high temperature durability of a ceramic, the flaw tolerance of a composite, a significantly faster and lower cost manufacturing process than conventional matrix CMCs and very low porosity levels < 0.5%. For thin-walled HX structures and the need to minimize leakage, the low porosity differentiator is particularly important. RTRC has prior experience with this material system and in the current project advanced the component design and manufacturing methods into new territory to produce features required for effective heat exchange under high pressures. In this approach, silicon carbide fiber is fabricated into a fiber preform using various textile processes. Graphite tooling is used both during the build-up of the fiber preform (interior tooling) and after the fiber preform has been completed (exterior tooling). This tooling assembly is heated to high temperature in an environment that has been evacuated and backfilled with inert gas. A reservoir of specialty glass is present and once the desired temperature has been reached to achieve the desired glass viscosity, an actuator distributes the glass throughout the fiber preform using passageways which are part of the tooling design in a process known as glass transfer molding. After the tooling has been removed, the composite is heat treated to convert the amorphous glass to a crystalline ceramic, providing improved properties. The project was divided into three phases focusing on the following: 1) 10 kW HX design and coupon-level tube sheet fabrication, 2) 10 kW HX fabrication, 3) 50 kW HX fabrication. During Budget Period 1 (BP1), additional risks were encountered and the need for additional funds was agreed upon by ARPA-E program leadership. Due to a variety of factors, the contract modification required nominally 18 months to execute at which time the HITEMMP program was effectively concluding. Because of this and the time that would be required to perform BP2 tasks, it was decided to conclude the project at the end of BP1. During the design of the 10 kW HX, manufacturing constraints were learned and incorporated, leading to a revised configuration for the fiber preform and HX. Heat exchange and pressure drop predictions also played a role in modifying the original design concept to be a higher aspect ratio shell-and-tube HX, simplifying the manufacturing process and improving the heat exchanger performance. Good gravimetric and volumetric thermal power densities of 11.2 kW/kg and 10,200 kW/m3 for the entire HX were projected that involved thermo-structural Finite Element Analysis to determine the structural mass needed for the high operation pressures of 250 bar cold inlet and 80 bar hot inlet. Fiber preforms using textile processes were produced for multiple headered tube sheets. Additional challenges were encountered during the glass transfer molding step for which solutions were identified, but programmatics did not allow them to be implemented in BP1. While complete HX test articles were not fabricated, the benefits of this GCMC material for a variety of high temperature applications remain.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2013

Absolute properties of the eclipsing binary VV CORVI

We have obtained red-wavelength spectroscopy and Johnson B and V differential photoelectric photometry of the eclipsing binary VV Crv = HR 4821. The system is the secondary of the common proper motion double star ADS 8627, which has a separation of 5.''2. VV Crv has an orbital period of 3.144536 days and a low but non-zero eccentricity of 0.085. With the Wilson-Devinney program we have determined a simultaneous solution of our spectroscopic and photometric observations. Those orbital elements produce masses of M {sub 1} = 1.978 ± 0.010 M {sub ☉} and M {sub 2} = 1.513 ± 0.008 M {sub ☉}, and radii of R {sub 1} = 3.375 ± 0.010 R {sub ☉} and R {sub 2} = 1.650 ± 0.008 R {sub ☉} for the primary and secondary, respectively. The effective temperatures of the two components are 6500 K (fixed) and 6638 K, so the star we call the primary is the more massive but cooler and larger component. A comparison with evolutionary tracks indicates that the components are metal rich with [Fe/H] = 0.3, and the system has an age of 1.2 Gyr. The primary is near the end of its main-sequence lifetime and is rotating significantly faster than its pseudosynchronous velocity. The secondary is still well ensconced on the main sequence and is rotating more slowly than its pseudosynchronous rate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Evaluation of genetic diversity, agronomic traits, and anthracnose resistance in the NPGS Sudan Sorghum Core collection

The United States Department of Agriculture (USDA) National Plant Germplasm System (NPGS) sorghum core collection contains 3011 accessions randomly selected from 77 countries. Genomic and phenotypic characterization of this core collection is necessary to encourage and facilitate its utilization in breeding programs and to improve conservation efforts. In this study, we examined the genome sequences of 318 accessions belonging to the NPGS Sudan sorghum core set, and characterized their agronomic traits and anthracnose resistance response. We identified 183,144 single nucleotide polymorphisms (SNPs) located within or in proximity of 25,124 annotated genes using the genotyping-by-sequencing (GBS) approach. The core collection was genetically highly diverse, with an average pairwise genetic distance of 0.76 among accessions. Population structure and cluster analysis revealed five ancestral populations within the Sudan core set, with moderate to high level of genetic differentiation. In total, 171 accessions (54%) were assigned to one of these populations, which covered 96% of the total genomic variation. Genome scan based on Tajima’s D values revealed two populations under balancing selection. Phenotypic analysis showed differences in agronomic traits among the populations, suggesting that these populations belong to different ecogeographical regions. A total of 55 accessions were resistant to anthracnose; these accessions could represent multiple resistance sources. Genome-wide association study based on fixed and random model Circulating Probability (farmCPU) identified genomic regions associated with plant height, flowering time, panicle length and diameter, and anthracnose resistance response. Integrated analysis of the Sudan core set and sorghum association panel indicated that a large portion of the genetic variation in the Sudan core set might be present in breeding programs but remains unexploited within some clusters of accessions. The NPGS Sudan core collection comprises genetically and phenotypically diverse germplasm with multiple anthracnose resistance sources. Population genomic analysis could be used to improve screening efforts and identify the most valuable germplasm for breeding programs. The new GBS data set generated in this study represents a novel genomic resource for plant breeders interested in mining the genetic diversity of the NPGS sorghum collection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

High spatial resolution mass spectrometry imaging reveals the genetically programmed, developmental modification of the distribution of thylakoid membrane lipids among individual cells of maize leaf

Metabolism in plants is compartmentalized among different tissues, cells and subcellular organelles. Mass spectrometry imaging (MSI) with matrix-assisted laser desorption ionization (MALDI) has recently advanced to allow for the visualization of metabolites at single-cell resolution. Here we applied 5- and 10 μm high spatial resolution MALDI-MSI to the asymmetric Kranz anatomy of Zea mays (maize) leaves to study the differential localization of two major anionic lipids in thylakoid membranes, sulfoquinovosyldiacylglycerols (SQDG) and phosphatidylglycerols (PG). The quantification and localization of SQDG and PG molecular species, among mesophyll (M) and bundle sheath (BS) cells, are compared across the leaf developmental gradient from four maize genotypes (the inbreds B73 and Mo17, and the reciprocal hybrids B73 × Mo17 and Mo17 × B73). SQDG species are uniformly distributed in both photosynthetic cell types, regardless of leaf development or genotype; however, PG shows photosynthetic cell-specific differential localization depending on the genotype and the fatty acyl chain constituent. Overall, 16:1-containing PGs primarily contribute to the thylakoid membranes of M cells, whereas BS chloroplasts are mostly composed of 16:0-containing PGs. Furthermore, PG 32:0 shows genotype-specific differences in cellular distribution, with preferential localization in BS cells for B73, but more uniform distribution between BS and M cells in Mo17. Maternal inheritance is exhibited within the hybrids, such that the localization of PG 32:0 in B73 × Mo17 is similar to the distribution in the B73 parental inbred, whereas that of Mo17 × B73 resembles the Mo17 parent. As a result, this study demonstrates the power of MALDI-MSI to reveal unprecedented insights on metabolic outcomes in multicellular organisms at single-cell resolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Simulation and fitting of complex reaction network TPR: The key is the objective function

In this research, a method has been developed for finding improved fits during simulation and fitting of data from complex reaction network temperature programmed reactions (CRN-TPR). It was found that simulation and fitting of CRN-TPR presents additional challenges relative to simulation and fitting of simpler TPR systems. The method used here can enable checking the plausibility of proposed chemical mechanisms and kinetic models. The most important finding was that when choosing an objective function, use of an objective function that is based on integrated production provides more utility in finding improved fits when compared to an objective function based on the rate of production. The response surface produced by using the integrated production is monotonic, suppresses effects from experimental noise, requires fewer points to capture the response behavior, and can be simulated numerically with smaller errors. For CRN-TPR, there is increased importance (relative to simple reaction network TPR) in resolving of peaks prior to fitting, as well as from weighting of experimental data points. Using an implicit ordinary differential equation solver was found to be inadequate for simulating CRN-TPR. Lastly, the method employed here was capable of attaining improved fits in simulation and fitting of CRN-TPR when starting with a postulated mechanism and physically realistic initial guesses for the kinetic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗