SEARCH · Engineering Papers

Results for “Functional data analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,711 records · Page 95

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Electron Tomography and Machine Learning for Understanding the Highly Ordered Structure of Leafhopper Brochosomes

Insects known as leafhoppers (Hemiptera: Cicadellidae) produce hierarchically structured nanoparticles known as brochosomes that are exuded and applied to the insect cuticle, thereby providing camouflage and anti-wetting properties to aid insect survival. Although the physical properties of brochosomes are thought to depend on the leafhopper species, the structure–function relationships governing brochosome behavior are not fully understood. Brochosomes have complex hierarchical structures and morphological heterogeneity across species, due to which a multimodal characterization approach is required to effectively elucidate their nanoscale structure and properties. In this work, we study the structural and mechanical properties of brochosomes using a combination of atomic force microscopy (AFM), electron microscopy (EM), electron tomography, and machine learning (ML)-based quantification of large and complex scanning electron microscopy (SEM) image data sets. This suite of techniques allows for the characterization of internal and external brochosome structures, and ML-based image analysis methods of large data sets reveal correlations in the structure across several leafhopper species. Our results show that brochosomes are relatively rigid hollow spheres with characteristic dimensions and morphologies that depend on leafhopper species. Nanomechanical mapping AFM is used to determine a characteristic compression modulus for brochosomes on the order of 1–3 GPa, which is consistent with crystalline proteins. Altogether, this work provides an improved understanding of the structural and mechanical properties of leafhopper brochosomes using a new set of ML-based image classification tools that can be broadly applied to nanostructured biological materials.

Chemical structure↗

DOE OSTI2023

Search for doubly charged Higgs boson production in multi-lepton final states using 139 fb$^{-1}$ of proton–proton collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for pair production of doubly charged Higgs bosons ($H$ ±± ), each decaying into a pair of prompt, isolated, and highly energetic leptons with the same electric charge, is presented. The search uses a proton–proton collision data sample at a centre-of-mass energy of 13 TeV corresponding to an integrated luminosity of 139 fb -1 recorded by the ATLAS detector during Run 2 of the Large Hadron Collider (LHC). This analysis focuses on same-charge leptonic decays, $H$ ±± → $\ell$ ± $\ell$ '± where $\ell$,$\ell$ ' = $e, μ, τ$, in two-, three-, and four-lepton channels, but only considers final states which include electrons or muons. No evidence of a signal is observed. Corresponding upper limits on the production cross-section of a doubly charged Higgs boson are derived, as a function of its mass $m$($H$ ±± ), at 95% confidence level. Assuming that the branching ratios to each of the possible leptonic final states are equal, $\mathcal{B}$($H$ ±± → $e$ ± $e$ ± ) = $\mathcal{B}$($H$ ±± → $e$ ± $μ$ ± ) = $\mathcal{B}$($H$ ±± → $μ$ ± $μ$ ± ) = $\mathcal{B}$($H$ ±± → $e$ ± $τ$ ± ) = $\mathcal{B}$($H$ ±± → $μ$ ± $τ$ ± ) = $\mathcal{B}$($H$ ±± → $τ$ ± $τ$ ± ) = 1/6, the observed (expected) lower limit on the mass of a doubly charged Higgs boson is 1080 GeV (1065 GeV) within the left right symmetric type-II seesaw model, which is the strongest limit to date produced by the ATLAS Collaboration. Additionally, this paper provides the first direct test of the Zee– Babu neutrino mass model at the LHC, yielding an observed (expected) lower limit of $m$($H$ ±± ) = 900 GeV (880 GeV).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Data for Sugar Accumulation Enhancement in Sorghum Stem is Associated with Reduced Reproductive Sink Strength and Increased Phloem Unloading Activity

Sweet sorghum has emerged as a promising source of bioenergy mainly due to its high biomass and high soluble sugar yield in stems. Studies have shown that loss-of-function Dry locus alleles have been selected during sweet sorghum domestication, and decapitation can further boost sugar accumulation in sweet sorghum, indicating that the potential for improving sugar yields is yet to be fully realized. To maximize sugar accumulation, it is essential to gain a better understanding of the mechanism underlying the massive accumulation of soluble sugars in sweet sorghum stems in addition to the Dry locus. We performed a transcriptomic analysis upon decapitation of near-isogenic lines for mutant (d, juicy stems, and green leaf midrib) and functional (D, dry stems and white leaf midrib) alleles at the Dry locus. Our analysis revealed that decapitation suppressed photosynthesis in leaves, but accelerated starch metabolic processes in stems. SbbHLH093 negatively correlates with sugar levels supported by genotypes (DD vs. dd), treatments (control vs. decapitation), and developmental stages post anthesis (3d vs.10d). D locus gene SbNAC074A and other programmed cell death-related genes were down regulated by decapitation, while sugar transporter-encoding gene SbSWEET1A was induced. Both SbSWEET1A and Invertase 5 were detected in phloem companion cells by RNA in situ assay. Loss of the SbbHLH093 homolog, AtbHLH093, in Arabidopsis led to a sugar accumulation increase. This study provides new insights into sugar accumulation enhancement in bioenergy crops, which can be potentially achieved by reducing reproductive sink strength and enhancing phloem unloading.

Transcriptomics↗

DOE OSTI2026

ICE Calculator 2: Final Report for Phase 1 and 2 of the National Initiative to Update the Interruption Cost Estimate (ICE) Calculator

ICE 2.0 Phase 2 Final Report In 2021, Berkeley Lab and Resource Innovations, Inc. launched the “ICE Calculator 2 Initiative” – a national study to refresh the underlying data and enhance the functionality of the ICE Calculator. The Initiative involves Berkeley Lab contracting with sponsoring utilities to administer identical, updated and comprehensive interruption cost surveys to statistically representative samples of each utility’s customers. Berkeley Lab and Resource Innovations then pool the survey results across the utilities and use them to update the analytical engines that drive the ICE Calculator. The ICE Calculator 2 Initiative is being conducted in phases. Each phase involves the administration of interruption cost surveys to the customers of sponsoring utilities, followed by an update to the ICE Calculator based on analysis of the pooled survey results. This report describes the activities and findings from Phase 1 and 2 of the ICE Calculator 2 Initiative. Phase 1 was sponsored by eight utilities: American Electric Power, Commonwealth Edison, Dominion Energy, Duke Energy, DTE Electric, Exelon, National Grid, and Puget Sound Energy. Phase 2 was sponsored by six utilities: Empire District Electric Company, Evergy Missouri, Pacific Gas & Electric, San Diego Gas & Electric, Southern California Edison, and Union Electric. Phase 1 and 2 involved 15 customer interruption cost survey activities representing a total of 30 electricity distribution service territories. ICE Calculator Version 2.0 and 2.2 Comparison This memorandum describes–at a high-level–the improvements in interruption cost estimates for the version 2.2 of the ICE Calculator (released February 2026) compared to version 2.0 (released in April 2025). Version 2.2 of the ICE Calculator corresponds to Phase 2 of the initiative, while version 2.0 corresponds to Phase 1. The improvements in version 2.2 result from both a significant increase in the number of customer responses that have been collected and the identification of seven additional factors that help estimate customer interruption costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Precision QCD, Hadronic Structure & Forward QCD, Heavy Ions: Report of Energy Frontier Topical Groups 5, 6, 7 submitted to Snowmass 2021

This report was prepared on behalf of three Energy Frontier Topical Groups of the Snowmass 2021 Community Planning Exercise. It summarizes the status and implications of studies of strong interactions in high-energy experiments and QCD theory. We emphasize the rich landscape and broad impact of these studies in the decade ahead. Hadronic interactions play a central role in the high-luminosity Large Hadron Collider (LHC) physics program, and strong synergies exist between the (HL-)LHC and planned or proposed experiments at the U.S. Electron-Ion Collider, CERN forward physics experiments, high-intensity facilities, and future TeV-range lepton and hadron colliders. Prospects for precision determinations of the strong coupling and a variety of nonperturbative distribution and fragmentation functions are examined. We also review the potential of envisioned tests of new dynamical regimes of QCD in high-energy and high-density scattering processes with nucleon, ion, and photon initial states. The important role of the high-energy heavy-ion program in studies of nuclear structure and the nuclear medium, and its connections with QCD involving nucleons are summarized. We address ongoing and future theoretical advancements in multi-loop QCD computations, lattice QCD, jet substructure, and event generators. Cross-cutting connections between experimental measurements, theoretical predictions, large-scale data analysis, and high-performance computing are emphasized.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Comparative analysis of atmospheric parameters from high-resolution spectroscopic sky surveys: APOGEE, GALAH, Gaia -ESO

SDSS-IV APOGEE-2, GALAH, and Gaia-ESO are high-resolution, ground-based, multi-object spectroscopic surveys providing fundamental stellar atmospheric parameters and multiple elemental abundance ratios for hundreds of thousands of stars of the Milky Way. Data from these and other surveys contribute to investigations of the history and evolution of the Galaxy. We undertake a comparison between the most recent data releases of these surveys to investigate the accuracy and precision of derived parameters by placing the abundances on an absolute scale. We also discuss the correlations in parameter and abundance differences as a function of main parameters. Uncovering the variants provides a basis to continue the efforts of future sky surveys. Quality samples from the APOGEE-GALAH (15 537 stars), APOGEE-GES (804 stars), and GALAH-GES (441 stars) overlapping catalogs were collected. We investigated the mean variants between the surveys, and linear trends were also investigated. We compared the slope of correlations and mean differences with the reported uncertainties. The average and scatter of v rad , T eff , log g, [M/H], and v micro , along with numerous species of elemental abundances in the combined catalogs, show that in general there is a good agreement between the surveys. We find large radial velocity scatters ranging from 1.3 km s -1 to 4.4 km s -1 when comparing the three surveys. We observe some weak trends (e.g., in ΔT eff vs. Δlog g for the APOGEE-GES stars) and a clear correlation in the v micro - Δv micro planes in the APOGEE-GALAH common sample. For [α/H], [Ti/H] (APOGEE-GALAH giants), and [Al/H] (APOGEE-GALAH dwarfs) potential strong correlations are discovered as a function of the differences in the main atmospheric parameters, and we also find weak trends for other elements. In general we find good agreement between the three surveys within their respective uncertainties. However, there are certain regimes in which strong variants exist, which we discuss. There are still offsets larger than 0.1 dex in the absolute abundance scales.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

A spectroscopic, structural, and computational study of Ag–oxo interactions in Ag + /UO 2 2+ complexes

For this work, twelve novel Ag + /UO 2 2+ heterometallic complexes have been prepared and characterized via structural, spectroscopic, and computational methods to probe the effects of Ag–oxo interactions on bonding and photophysical properties of the uranyl cation. Structural characterization reveals Ag–oxo interaction distances ranging from 2.475(3) Å to 4.287(4) Å. These interactions were probed using luminescence and Raman spectroscopy which displayed little effect on the luminescence intensity and the energy of the Raman active U═O symmetric stretch peak as compared to previously reported Pb–oxo interactions. Computational efforts via density functional theory-based natural bond orbital analysis revealed that the highest stabilization energy associated with the Ag–oxo interaction had a value of only 11.03 kcal mol –1 and that all other energy values fell at 7.05 kcal mol –1 or below indicating weaker interactions relative to those previously reported for Pb 2+ /UO 2 2+ heterometallic compounds. In contrast, quantum theory of atoms in molecules analysis of bond critical point electron density values indicated higher electron density in Ag–oxo interactions as compared to Pb–oxo interactions which suggests more covalent character with the Ag + . Overall, this data indicates that Ag + has a less significant effect on UO 2 2+ bonding and photophysical properties as compared to other Pb 2+ , likely due to the high polarizability of the cation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Sensitivity of an integrated experiment to uncertainty in the high explosive equations of state

Traditionally, hydrodynamics simulations are performed with a single equation-of-state (EOS) to describe each material. These EOSs typically have a physics-informed functional form with adjustable parameters that are calibrated in order to replicate small-scale data. However, because the calibration data have uncertainty and there are typically inherent degeneracies in fitting the EOS, there are actually multiple EOSs that might be consistent with calibration data. In this work, we perform uncertainty quantification (UQ) for the reactant and product equations of state for the high explosive PBX 9501 to yield an ensemble of EOSs that match the uncertain small-scale calibration data. We then simulate an experiment of an explosively formed penetrator repeatedly with different EOSs to both validate the UQ analysis and determine the effects of EOS uncertainty on the prediction of quantities of interest in the experiment. In general, we find good agreement between the simulation predictions and the experimental measurements, and we identify an EOS variable that contributes most directly to the spread in the predictions as the EOSs are varied.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Genetic switching by the Lac repressor is based on two-state Monod–Wyman–Changeux allostery

High-resolution NMR spectroscopy enabled us to characterize allosteric transitions between various functional states of the dimeric Escherichia coli Lac repressor. In the absence of ligands, the dimer exists in a dynamic equilibrium between DNA-bound and inducer-bound conformations. Binding of either effector shifts this equilibrium toward either bound state. Analysis of the ternary complex between repressor, operator DNA, and inducer shows how adding the inducer results in allosteric changes that disrupt the interdomain contacts between the inducer binding and DNA binding domains and how this in turn leads to destabilization of the hinge helices and release of the Lac repressor from the operator. Based on our data, the allosteric mechanism of the induction process is in full agreement with the well-known Monod–Wyman–Changeux model.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Ultra-Long Distance BOTDA Sensor System Employing Hybrid Amplification and Advanced Noise Reduction Techniques

Brillouin Optical Time Domain Analysis (BOTDA) sensor system play a pivotal role in distributed sensing, which enables precise measurements of strain and temperature across extensive fiber lengths. Nonetheless, challenges emerge as distances grow due to signal attenuation and noise interference resulting in measurement errors. This research offers a comprehensive strategy to extend the sensing range of BOTDA systems beyond 10’s of kilometers while maintaining high spatial resolutions. Such enhanced sensing is realized through the integration of distributed Raman amplification, inline amplification using erbium-doped fiber amplifiers (EDFA), and advanced noise reduction techniques. Leveraging inherent redundancy in measured data as a function of frequency and fiber distance, the non-local means (NLM) filter removes noise while preserving essential physical information. This approach proves particularly advantageous in BOTDA systems, where accurate measurement of Brillouin scattering signals is paramount for long-range sensing, while concurrently safeguarding high spatial resolutions. In summary, this research has shown a holistic exploration of extending BOTDA's distance sensing capabilities up to 150 km with spatial resolutions of 8 meters.

Bhatta, Hari↗

DOE OSTI2026

NOvA joint $ν^{e}$ + $ν^{µ}$ oscillation results in neutrino and antineutrino modes

NOvA is an experiment devoted to studying neutrino oscillations in the NuMI neutrino beam from FNAL (USA). It is a long-baseline experiment consisting of two functionally identical, finely granulated detectors which are separated by 810 km of Earth crust and sited at 14 mrad off the beam axis. By measuring the transition probabilities P(νμ→νe \nu_\mu \rightarrow \nu_e ) and P(νμ→νμ \nu_\mu \rightarrow \nu_\mu ) NOvA is able to extract oscillation parameters: Δm232 \Delta m^2_{32} , mixing angle θ23 \theta_{23} , CP violating phase δCP \delta_{CP} and neutrino mass hierarchy. This analysis will be the first to include both neutrino (9⋅1020 9 \cdot 10^{20} POT) and antineutrino (7⋅1020 7 \cdot 10^{20} POT) data, which helps to resolve degeneracies in the oscillation probability. In this poster, the lastest NOvA oscillation results will be discussed, as well as important intermediate steps in νe \nu_e analysis like signal and background predictions, projected sensitivities in upcoming analyses will be presented.

Back, Ashley [Iowa State U.]↗

DOE OSTI2023

Towards the extraction of the crystal cell parameters from pair distribution function profiles

The approach based on atomic pair distribution function (PDF) has revolutionized structural investigations by X-ray/electron diffraction of nano or quasi-amorphous materials, opening up the possibility of exploring short-range order. However, the ab initio crystal structural solution by the PDF is far from being achieved due to the difficulty in determining the crystallographic properties of the unit cell. A method for estimating the crystal cell parameters directly from a PDF profile is presented, which is composed of two steps: first, the type of crystal cell is inferred using machine-learning approaches applied to the PDF profile; second, the crystal cell parameters are extracted by means of multivariate analysis combined with vector superposition techniques. The procedure has been validated on a large number of PDF profiles calculated from known crystal structures and on a small number of measured PDF profiles. The lattice determination step has been benchmarked by a comprehensive exploration of different classifiers and different input data. The highest performance is obtained using the k -nearest neighbours classifier applied to whole PDF profiles. Descriptors calculated from the PDF profiles by recurrence quantitative analysis produce results that can be interpreted in terms of PDF properties, and the significance of each descriptor in determining the prediction is evaluated. The cell parameter extraction step depends on the cell metric rather than its type. Monometric, dimetric and trimetric cells have top-1 estimates that are correct 40, 20 and 5% of the time, respectively. Promising results were obtained when analysing real nanocrystals, where unit cells close to the true ones are found within the top-1 ranked solution in the case of monometric cells and within the top-6 ranked solutions in the case of dimetric cells, even in the presence of a crystalline impurity with a weight fraction up to 40%.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A comparative analysis of different van der Waals treatments for molecular adsorption on the basal plane of 2H-MoS 2

The binding energy of hydrogen sulfide, ammonia, ethane, ethylene, butadiene, benzene, toluene, pyridine, pyrrole, and thiophene on the basal plane of the semi-conducting 2H-molybdenum sulfide (MoS 2 ) was calculated with the following flavors of Density Functional Theory (DFT): GGA-PW91, PBE-D3, vdW-DF, optPBE, optB86b, optB88, vdW-TS, and BEEF-vdW. The GGA-PW91 binding energies are negligible (<0.07 eV in magnitude) in all cases. The predictions with vdW-DF and PBE-D3 are the closest (error <0.05 eV) to the isosteric heats of adsorption calculated from reported temperature programmed desorption data for thiophene and butadiene. For all dispersion flavors examined here, the magnitude of the dispersion contribution to the binding energy increases linearly with the number of heavy atoms in the adsorbate, with each atom contributing 0.05 eV (BEEF-vdW) – 0.09 eV (optB88-vdW). Further, this implies that the calculated adsorption constants of molecules larger than acridine (i.e., comprising > 14 non-heavy atoms) can vary by more than four orders of magnitude at industrial conditions depending on the chosen method of dispersion correction. Further, dispersion effects fall off rapidly (>0.03 eV/ non-hydrogen atom/Å) as the adsorbate-surface distance increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

A priori analysis on deep learning of subgrid-scale parameterizations for Kraichnan turbulence

In the present study, we investigate different data-driven parameterizations for large eddy simulation of two-dimensional turbulence in the a priori settings. These models utilize resolved flow field variables on the coarser grid to estimate the subgrid-scale stresses. We use data-driven closure models based on localized learning that employs a multilayer feedforward artificial neural network with point-to-point mapping and neighboring stencil data mapping, and convolutional neural network fed by data snapshots of the whole domain. The performance of these data-driven closure models is measured through a probability density function and is compared with the dynamic Smagorinksy model (DSM). The quantitative performance is evaluated using the cross-correlation coefficient between the true and predicted stresses. We analyze different frameworks in terms of the amount of training data, selection of input and output features, their characteristics in modeling with accuracy, and training and deployment computational time. We also demonstrate computational gain that can be achieved using the intelligent eddy viscosity model that learns eddy viscosity computed by the DSM instead of subgrid-scale stresses. Here, we detail the hyperparameters optimization of these models using the grid search algorithm.

42 ENGINEERING↗

DOE OSTI2021

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Furthermore, fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Crustal and uppermost mantle structure near the Gloria Fault, North Atlantic, from ocean bottom seismometer surface wave observations

SUMMARY In this work, we present both 1-D and 3-D shear wave velocity (Vs) models of the oceanic crust and uppermost mantle below the Deep OCean Test ARray area, located ∼ 70 km north of the central section of the Gloria Fault, in the eastern North-Atlantic Ocean. The velocity models are inferred from the dispersion of surface waves recorded on ocean bottom seismometers. Dispersion measurements are obtained from the analysis of ambient seismic noise at short periods (< 14 s) and teleseismic surface waves at long periods (> 14 s) using the two-station method. The 1-D Vs model is inferred from the joint inversion of Rayleigh wave phase velocities and Love wave group and phase velocities. The 3-D tomographic model is obtained by inversion of 2-D Love wave group velocity maps as a function of depth, further constrained by the average of Love wave phase velocities obtained from ambient noise (4–9 s) and the average Rayleigh and Love wave phase velocities calculated from teleseismic data (14–44 s). The 1-D Vs model shows a sediment layer with a low velocity of 1.05 km s−1, similar to previous studies in the region. Below the sediments, we find an oceanic crust with velocities ranging from 3.3 to 4.5 km s−1. The model reaches an unusually high velocity of 4.9 km s−1 in a 20 km thick layer at depths between 16 and 36 km. We interpret this fast velocity layer as indicative of the presence of harzburgite, a residue of enhanced melting that might have been formed by the proximity between the Mid-Atlantic Ridge and the Azores mantle plume. At greater depths the velocity decreases, forming a low-velocity zone that reaches a minimum at ∼ 70 km depth, which we interpret as the maximum depth for the lithosphere–asthenosphere boundary. The 3-D model shows a structure that is mostly horizontally layered, with Vs isocontours at 3.5–4.5 km s−1 highlighting oscillations of the crustal structure with wavelengths of ∼25–30 km. These oscillations may be due to changes in the rate of mantle upwelling and magma supply rate.

Pinzón, Juan I. (ORCID:0000000338078068)↗

DOE OSTI2024

guppy i : a code for reducing the storage requirements of cosmological simulations

ABSTRACT As cosmological simulations have grown in size, the permanent storage requirements of their particle data have also grown. Even modest simulations present a major logistical challenge for the groups which run these boxes and researchers without access to high performance computing facilities often need to restrict their analysis to lower quality data. In this paper, we present guppy, a compression algorithm and code base tailored to reduce the sizes of dark matter-only cosmological simulations by approximately an order of magnitude. guppy is a ‘lossy’ algorithm, meaning that it injects a small amount of controlled and uncorrelated noise into particle properties. We perform extensive tests on the impact that this noise has on the internal structure of dark matter haloes, and identify conservative accuracy limits which ensure that compression has no practical impact on single-snapshot halo properties, profiles, and abundances. We also release functional prototype libraries in C, Python, and Go for reading and creating guppy data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Transcript and metabolite network perturbations in lignin biosynthetic mutants of Arabidopsis

Abstract Lignin, one of the most abundant polymers in plants, is derived from the phenylpropanoid pathway, which also gives rise to an array of metabolites that are essential for plant fitness. Genetic engineering of lignification can cause drastic changes in transcription and metabolite accumulation with or without an accompanying development phenotype. To understand the impact of lignin perturbation, we analyzed transcriptome and metabolite data from the rapidly lignifying stem tissue in 13 selected phenylpropanoid mutants and wild-type Arabidopsis (Arabidopsis thaliana). Our dataset contains 20,974 expressed genes, of which over 26% had altered transcript levels in at least one mutant, and 18 targeted metabolites, all of which displayed altered accumulation in at least one mutant. We found that lignin biosynthesis and phenylalanine supply via the shikimate pathway are tightly co-regulated at the transcriptional level. The hierarchical clustering analysis of differentially expressed genes (DEGs) grouped the 13 mutants into 5 subgroups with similar profiles of mis-regulated genes. Functional analysis of the DEGs in these mutants and correlation between gene expression and metabolite accumulation revealed system-wide effects on transcripts involved in multiple biological processes.

Plant Sciences↗