Search NASA⌕ Search

SEARCH · Search NASA

Results for “transform limit”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Limit Cycle Oscillating Detonation Wave Behavior Analysis Within a Rotating Detonation Engine

Limit cycle oscillation (LCO) detonation wave behaviors are presented and analyzed for test times exceeding 20 s in a water-cooled rotating detonation engine (RDE). LCO detonation waves exhibit cyclic acceleration and deceleration, resulting in oscillating wave spacing at unique process conditions. In previous RDE studies, similar behaviors have been studied as microsecond-scale instabilities leading to ascending or descending modal transitions. In the current work, however, LCO waves are considered a persistent wave mode, occupying unique portions of the operational envelope adjacent to those of their equally spaced counterparts. These occurrences of LCO waves are repeatable, enduring behaviors. A method to generate shifted contour surfaces specifically intended to extract and analyze wave spacing variation through time, termed limit cycle oscillation visualization (LCOV) surfaces, is presented. LCOV surfaces transform data into the reference frame of a primary traveling wave and are used to analyze quasi-steady, short-timescale, and transitional LCO modes. Results are leveraged to understand the relationship between fill height, wave strength, local wave acceleration, and subsequent LCO wave spacing for individual wave sets. In conclusion, quasi-steady LCO waves display wave spacing oscillations between equal spacing values associated with ±1 wave across runs exceeding 18 s.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Aircraft-derived particle fluxes distinguish entrainment zone and decoupled layer nucleation in marine boundary layers

The vertical distribution of freshly nucleated aerosol particles in the marine boundary layer remains poorly constrained, limiting our ability to represent new particle formation in climate models. Here we characterize 3–10 nm particle events, termed small particle events (SPEs), by deriving their vertical turbulent fluxes from aircraft measurements during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign. To overcome stationarity limitations of traditional eddy covariance methods, we applied continuous wavelet transform analysis to data collected during June–July 2017 and January–February 2018 flights over the Azores. Our flux-based analysis revealed two distinct SPE scenarios with different vertical structures and spatial extents. The first featured nucleation in the entrainment zone, where free tropospheric air entrains into the boundary layer. The second showed nucleation in the decoupled layer, a stratified region between the well-mixed surface layer and cloud-topped upper boundary layer. In both cases, convergence of air masses from different layers diluted preexisting aerosol surface area to very low levels, creating conditions favorable for nucleation and generating strong downward particle fluxes. SPEs occurred in 15 % of flights, challenging prevailing theoretical expectations that new particle formation should rarely occur in marine boundary layers due to high condensation and coagulation sink capacity of sea spray aerosols. Aircraft-derived particle fluxes provide first observational constraints on the vertical location and source strength of likely nucleation regions in the remote marine boundary layer, improving aerosol source representations in climate models and reducing uncertainties in aerosol-cloud interactions.

Rasheeda Satheesh, Ajmal [North Carolina State Uni↗
DOE OSTI2020

Dynamic re-equilibration controlled multi-step transformations in (Mn, Fe) 2 (P, Si) alloys

In (Mn, Fe) 2 (P, Si) alloys crystallographic first-order phase transformations enable strong coupling of magnetic and entropic properties, potentially leading to high-efficiency energy generation and refrigeration applications. Although hysteresis losses that limit these applications can be reduced through careful control of alloy composition, compositional tuning can also unfavorably influence transformation temperatures and magnetocaloric coupling strength. Hence, exploration of additional processing variables enabling independent control of transformation properties is crucial. In this work, we investigate the role of thermal history as an additional processing variable, exploiting thermally-activated mechanisms to control properties of non-diffusive transformations in (Mn, Fe) 2 (P, Si) alloys. In so doing, we report an unusual transformation-splitting phenomenon following annealing at intermediate times, where a single well-defined magneto-structural transformation evolves towards a multi-step transformation with individual steps occurring at multiple distinct temperatures. Here, on longer annealing at the same temperatures, single-step transformation behavior is recovered. Through additional magnetic and crystallographic characterization, we show that the thermal history-controlled multi-step behavior results from sluggish thermally-activated diffusion. The two-step transformation corresponds to non-equilibrium bimodal composition distributions in the transforming phase, and these develop through a dynamic re-equilibration process as the alloy passes relatively slowly between different thermal equilibria. Together, these results suggest that thermal history primarily controls the transformation properties of (Mn, Fe) 2 (P, Si) alloys indirectly through the composition of one or more transforming phases. Additional investigations are needed to develop thermal history processing for decoupling hysteresis control from other transformation properties.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

High-Throughput Microfluidic Electroporation (HTME): A Scalable, 384-Well Platform for Multiplexed Cell Engineering

Electroporation-mediated gene delivery is a cornerstone of synthetic biology, offering several advantages over other methods: higher efficiencies, broader applicability, and simpler sample preparation. Yet, electroporation protocols are often challenging to integrate into highly multiplexed workflows, owing to limitations in their scalability and tunability. These challenges ultimately increase the time and cost per transformation. As a result, rapidly screening genetic libraries, exploring combinatorial designs, or optimizing electroporation parameters requires extensive iterations, consuming large quantities of expensive custom-made DNA and cell lines or primary cells. To address these limitations, we have developed a High-Throughput Microfluidic Electroporation (HTME) platform that includes a 384-well electroporation plate (E-Plate) and control electronics capable of rapidly electroporating all wells in under a minute with individual control of each well. Fabricated using scalable and cost-effective printed-circuit-board (PCB) technology, the E-Plate significantly reduces consumable costs and reagent consumption by operating on nano to microliter volumes. Furthermore, individually addressable wells facilitate rapid exploration of large sets of experimental conditions to optimize electroporation for different cell types and plasmid concentrations/types. Use of the standard 384-well footprint makes the platform easily integrable into automated workflows, thereby enabling end-to-end automation. We demonstrate transformation of E. coli with pUC19 to validate the HTME's core functionality, achieving at least a single colony forming unit in more than 99% of wells and confirming the platform's ability to rapidly perform hundreds of electroporations with customizable conditions. This work highlights the HTME's potential to significantly accelerate synthetic biology Design-Build-Test-Learn (DBTL) cycles by mitigating the transformation/transfection bottleneck.

Gaillard, William R↗
DOE OSTI2023

Development of Scalable Reactive Transport Framework for PFAS

Research on PFAS chemicals is extensive and covers all aspects, including analytical quantification, determination of properties, toxicology, ex situ treatment, and in situ remediation. PFAS chemicals are very stable and thus persistent/recalcitrant in the environment. Although there are many unknowns about PFAS chemicals, degradation pathways, reaction rates, etc., many different sorption, oxidation, reduction, biological, and innovative treatment approaches are being developed. Sorption with activated carbon is currently the only fully available in situ treatment technology for PFAS-impacted groundwater. Given the wide array of PFAS chemicals and transformation products, remediation may need multi-step treatment trains to fully address the PFAS contamination. The work here provides kinetic reaction modules that represent an initial cut at functionality representing PFAS migration and reaction in groundwater aquifer flow and transport models. One reaction kinetics module provides a method to model kinetically limited adsorption using a mass transfer model. The second reaction module represents biological transformation of 8:2 FTOH and daughter species, illustrating how a complex reaction pathway network can be represented. Both reaction modules allow for spatially variable parameter values so that a variety of remediation approaches (e.g., a PRB or volumetric treatment or variations in geochemical conditions) can be simulated. The intent with these PFAS reaction modules is to provide tools for practitioners to aid in the selection, design, and assessment of potential in situ PFAS remediation strategies. It is anticipated that, as PFAS remediation technologies and scientific understanding advances, these modules would be refined or replaced to match new knowledge.

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

Machine learning enables reconstruction of past fire regimes from charcoal-derived fire intensity and fuel composition

Background Fire is a foundational ecological process that shapes ecosystem structure, diversity, and resilience. Quantifying paleofire regime attributes such as frequency, severity, and intensity is essential for understanding the historical range of variability in fire behavior and its ecological effects. While frequency and severity are often reconstructed in paleofire studies, quantitative reconstructions of fire intensity remain limited. Recent work has shown that maximum pyrolysis temperature—a proxy for fire intensity—and plant species type can be inferred from charcoal using transmission Fourier-transform infrared (FTIR) spectroscopy. However, the sample preparation for transmission FTIR is destructive and time-consuming, limiting application and reuse of materials for other analyses. We evaluated reflectance FTIR spectroscopy as a non-destructive alternative for reconstructing combustion temperature and plant species from laboratory-generated charcoal. We also examined the influence of contrasting airflow environments (ambient air versus nitrogen-rich) on pyrolysis temperature and plant species reconstruction prediction accuracies and compared predictive performance between a novel, neural network–based deep learning model with the traditional modern analogue technique (MAT) using k-nearest neighbor functions. As proof of concept, we apply our enhanced methodology to ancient charcoal to demonstrate applicability at improving long-term fire regime reconstructions and the ability to link paleofire records with contemporary fire ecology. Results Our analysis shows that transmission and reflectance FTIR spectra yield comparable spectral profiles. However, sample preparation for reflectance FTIR is minimal and non-destructive, unlike transmission FTIR which is destructive. We demonstrate that oxygen environments improved reconstruction accuracy relative to nitrogen-rich conditions. Finally, our deep learning neural network (DL) achieved testing accuracies of 98.7% for temperature and 96.2% for species identification, outperforming MAT’s k-NN approach (89.8% and 65.9%, respectively). A Shapley importance analysis identified 5 key spectral regions that greatly influenced the model’s temperature or species categorization. When applied to ancient charcoal, our results show historic fires from the most recent past primarily burned at low intensities (400–500 °C), reflective of natural fire regimes in ponderosa pine forests. Our results corroborate charcoal morphology data that suggests all ancient charcoal originated from burned woody plant types. Conclusions By combining reflectance FTIR spectroscopy with a deep learning approach, we provide the first accuracies high enough to confidently identify both species and temperature from laboratory-produced charcoal, improving quantitative reconstructions of fire intensity and fuel composition from paleofire records. This opens a wide range of research into the link between fire and larger drivers (i.e., climate or human) and greater ecological understanding of fire regimes beyond that of burn scars or recent observations. These methodological improvements have direct relevance for fire management by improving interpretation of historical fire behavior, informing fuel–fire relationships, and providing a scalable analytical framework applicable to both long-term ecological studies and contemporary fire science.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Optimal Spectral Resolution for Infrared Studies of Solids and Liquids

Due to a legacy originating in the limited capability of early computers, the spectroscopic resolution used in Fourier transform infrared spectroscopy and other systems has largely been implemented using only powers of two for more than 50 years. In this study, we investigate debunking the spectroscopic lore of, e.g., using only 2, 4, 8, or 16 cm −1 resolution and determine the optimal resolution in terms of both (i) a desired signal-to-noise ratio and (ii) efficient use of acquisition time. The study is facilitated by the availability of solids and liquids reference spectral data recorded at 2.0 cm −1 resolution and is based on an examination in the 4000–400 cm −1 range of 61 liquids and 70 solids spectra, with a total analysis of 4237 peaks, each of which was also examined for being singlet/multiplet in nature. Of the 1765 liquid bands examined, only 27 had widths <5 cm −1 . Of the 2472 solid bands examined, only 39 peaks have widths <5 cm −1 . For both the liquid and solid bands, a skewed distribution of peak widths was observed: For liquids, the mean peak width was 24.7 cm −1 but the median peak width was 13.7 cm −1 , and, similarly, for solids, the mean peak width was 22.2 cm −1 but the median peak width was 11.2 cm −1 . While recognizing other studies may differ in scope and limiting the analysis to only room temperature data, we have found that a resolution to resolve 95% of all bands is 5.7 cm −1 for liquids and 5.3 cm −1 for solids; such a resolution would capture the native linewidth (not accounting for instrumental broadening) for 95% of all the solids and liquid bands, respectively. After decades of measuring liquids and solids at 4, 8, or 16 cm −1 resolution, we suggest that, when accounting only for intrinsic linewidths, an optimized resolution of 6.0 cm −1 will capture 91% of all condensed-phase bands, i.e., broadening of only 9% of the narrowest of bands, but yielding a large gain in signal-to-noise with minimal loss of specificity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Exploration of Electrochemical Reactions at Organic–Inorganic Halide Perovskite Interfaces via Machine Learning in In Situ Time-of-Flight Secondary Ion Mass Spectrometry

The instability of hybrid organic–inorganic perovskite (HOIP) devices is one of the significant challenges preventing commercialization. Exploring these phenomena is severely limited by the complexity of the intrinsic electrochemistry of HOIPs, the presence of multiple volatile and mobile ionic species, and the possible role of environmentally induced reactions at surfaces and triple-phase junctions. Here, in situ studies of the electrochemistry of methylammonium lead bromide perovskite with the Au electrode interface are reported via light- and voltage-dependent time-of-flight secondary ion mass spectrometry (ToF-SIMS) imaging of lateral perovskite heterostructures. While ToF-SIMS allows for the visualization of the chemical composition along the surface and its evolution with light and electrical bias, the interpretation of the multidimensional data obtained is often limited due to strong correlations between chemical signatures and the need to track multiple peaks at once. Here, a machine learning workflow combining the Hough transform and non-negative matrix factorization and no-negative tensor decomposition is developed to avoid this limitation and extract salient features of associated chemical changes and to separate the light- and voltage-dependent dynamics. Finally, combining these in situ characterizations and the machine learning workflow provides comprehensive information on the chemical nature of moving species, ion accumulation, and interfacial electrochemical reactions in HOIP devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

scPlantAnnotate: an accurate and robust transformer-based model for plant cell type annotation

Accurate cell type annotation remains a major bottleneck in plant single-cell RNA sequencing (scRNA-seq), where existing tools are often adapted from animal studies and perform sub-optimally on plant data. The lack of plant-specific computational frameworks limits the construction of plant cell atlases and downstream biological discovery. We develop and evaluate scPlantAnnotate, a Transformer-based reference annotation framework tailored for plant scRNA-seq data, and benchmark it against state-of-the-art deep learning and conventional methods across multiple plant species. Species-specific scPlantAnnotate models were trained using curated datasets from Arabidopsis thaliana, Zea mays, Oryza sativa, and Glycine max. We compared scPlantAnnotate with leading baselines under both standard random-split evaluation and a more stringent leave-one-dataset-out setting, which tests robustness to completely unseen datasets and tissue types. scPlantAnnotate consistently outperforms existing approaches across all four species under random-split evaluation. In the leave-one-dataset-out setting for A. thaliana, where performance drops markedly for all methods due to strong batch effects and dataset heterogeneity, scPlantAnnotate nonetheless achieves the highest Accuracy, Macro-F1, Balanced Accuracy, and Macro-AUROC on average and ranks first on most held-out datasets. These results demonstrate improved robustness to dataset shifts, a critical yet underexplored challenge in plant scRNA-seq analysis. A freely accessible web server enables users to annotate their own datasets using pretrained models. scPlantAnnotate provides a plant-specific, Transformer-based framework for single-cell annotation that delivers state-of-the-art performance and enhanced robustness to unseen datasets. By addressing limitations of existing tools and enabling scalable reference-based annotation, scPlantAnnotate supports the development of comprehensive plant cell atlases and facilitates broader use of single-cell genomics in plant biology.

Bioinformatics↗
DOE OSTI2025

Compressing Vision Transformers in Geospatial Transfer Learning with Manifold-Constrained Optimization

Deploying geospatial foundation models on resource-constrained edge devices demands compact architectures that maintain high downstream performance. However, their large parameter counts and the accuracy loss often induced by compression limit practical adoption.In this work, we leverage manifold-constrained optimization framework DLRT to compress large vision transformer–based geospatial foundation models during transfer learning. By enforcing structured low-dimensional parameterizations aligned with downstream objectives, this approach achieves strong compression while preserving task-specific accuracy. We show that the method outperforms of-the-shelf low-rank methods as LoRA. Experiments on diverse geospatial benchmarks confirm substantial parameter reduction with minimal accuracy loss, enabling high-performing, on-device geospatial models.

Snyder, Thomas [Yale University]↗
DOE OSTI2020

Vulnerable window of yield strength for swelling-driven fracture of phase-transforming battery materials

Abstract Despite numerous experimental and theoretical investigations of the mechanical behavior of high-capacity Si and Ge Li-ion battery anodes, our basic understanding of swelling-driven fracture in these materials remains limited. Existing theoretical studies have provided insights into elasto-plastic deformations caused by large volume change phase transformations, but have not modeled fracture explicitly beyond Griffith’s criterion. Here, we use a multi-physics phase-field approach to model self-consistently anisotropic phase transformation, elasto-plastic deformation, and crack initiation and propagation during lithiation of Si nanopillars. Our computational results reveal that fracture occurs within a “vulnerable window” inside the two-dimensional parameter space of yield strength and fracture energy and highlight the importance of taking into account the surface localization of plastic deformation to accurately predict the magnitude of tensile stresses at the onset of fracture. They further demonstrate how the increased robustness of hollow nanopillars can be understood as a direct effect of anode geometry on the size of this vulnerable window. Those insights provide an improved theoretical basis for designing next-generation mechanically stable phase-transforming battery materials undergoing large volume changes.

25 ENERGY STORAGE↗
DOE OSTI2024

High-Entropy Metal-Organic Frameworks (HEMOFs): A New Frontier in Materials Design for CO 2 Utilization

High-entropy materials (HEMs) emerged as promising candidates for a diverse array of chemical transformations, including CO 2 utilization. However, traditional HEMs catalysts are nonporous, limiting their activity to surface sites. Designing HEMs with intrinsic porosity can open the door toward enhanced reactivity while maintaining the many benefits of high configurational entropy. Here, in this study, a synergistic experimental, analytical, and theoretical approach to design the first high-entropy metal-organic frameworks (HEMOFs) derived from polynuclear metal clusters is implemented, a novel class of porous HEMs that is highly active for CO 2 fixation under mild conditions and short reaction times, outperforming existing heterogeneous catalysts. HEMOFs with up to 15 distinct metals are synthesized (the highest number of metals ever incorporated into a single MOF) and, for the first time, homogenous metal mixing within individual clusters is directly observed via high-resolution scanning transmission electron microscopy. Importantly, density functional theory studies provide unprecedented insight into the electronic structures of HEMOFs, demonstrating that the density of states in heterometallic clusters is highly sensitive to metal composition. This work dramatically advances HEMOF materials design, paving the way for further exploration of HEMs and offers new avenues for the development of multifunctional materials with tailored properties for a wide range of applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Celestial operator product expansions and w 1+∞ symmetry for all spins

The operator product expansion of massless celestial primary operators of arbitrary spin is investigated. Poincaré symmetry is found to imply a set of recursion relations on the operator product expansion coefficients of the leading singular terms at tree-level in a holomorphic limit. The symmetry constraints are solved by an Euler beta function with arguments that depend simply on the right-moving conformal weights of the operators in the product. These symmetry-derived coefficients are shown not only to match precisely those arising from momentum-space tree-level collinear limits, but also to obey an infinite number of additional symmetry transformations that respect the algebra of w 1+∞ . In tree-level minimally-coupled gravitational theories, celestial currents are constructed from light transforms of conformally soft gravitons and found to generate the action of w 1+∞ on arbitrary massless celestial primaries. Results include operator product expansion coefficients for fermions as well as those arising from higher-derivative non-minimal couplings of gluons and gravitons.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

High performance Nb/TiNi nanocomposites produced by packaged accumulative roll bonding

Taking full advantage of intrinsic high strength of nano-reinforcements has been proven difficult in a composite until the concept of strain matching was experimental verified in a nanowire Nb/TiNi in-situ composite, which was obtained via solidifying eutectic Nb/TiNi and subsequent severe wire drawing. However, the volume fraction of the nano-reinforcement was dictated by eutectic reaction, which severely limited the development of high-performance composites. Here, the martensitic transformation nanocomposites were successfully acquired with increased volume fraction and changeable morphology of the nano-reinforcement through hot packaged ARB (PARB) and subsequent wire drawing. The nano-reinforcement with lamellar configuration in the composite could exhibit a large lattice strain of similar to 3.8%, which demonstrated the applicability and effectiveness of the strain matching principle in nano-lamellar composite. The obtained composites show a yield strength over 2.0 GPa and a quasi-linear elastic strain as larger as similar to 6.5% with adjustable apparent elastic modulus from 30 to 50 GPa.

42 ENGINEERING↗
DOE OSTI2025

Deep Learning with Reflection High-Energy Electron Diffraction Images to Predict Cation Ratio in Sr 2 x Ti 2(1– x ) O 3 Thin Films

Machine learning (ML) with in-situ diagnostics offers a transformative approach to accelerate, understand, and control thin film synthesis by uncovering relationships between synthesis conditions and material properties. In this study, we demonstrate the application of deep learning to predict the stoichiometry of Sr 2x Ti 2(1–x) O 3 thin films using reflection high-energy electron diffraction images acquired during pulsed laser deposition. A gated convolutional neural network trained for regression of the Sr atomic fraction achieved accurate predictions with a small dataset of 31 samples. Explainable AI techniques revealed a previously unknown correlation between diffraction streak features and cation stoichiometry in Sr 2x Ti 2(1–x) O 3 thin films. Here, our results demonstrate how ML can be used to transform a ubiquitous in-situ diagnostic tool, that is usually limited to qualitative assessments, into a quantitative surrogate measurement of continuously valued thin film properties. Such methods are critically needed to enable real-time control, autonomous workflows, and accelerate traditional synthesis approaches.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

A multi-omics approach to lignocellulolytic enzyme discovery reveals a new ligninase activity from Parascedosporium putredinis NO1

Lignocellulose, the structural component of plant cells, is a major agricultural byproduct and the most abundant terrestrial source of biopolymers on Earth. The complex and insoluble nature of lignocellulose limits its conversion into value-added commodities, and currently, efficient transformation requires expensive pretreatments and high loadings of enzymes. In this work, we report on a fungus from the Parascedosporium genus, isolated from a wheat-straw composting community, that secretes a large and diverse array of carbohydrate-active enzymes (CAZymes) when grown on lignocellulosic substrates. We describe an oxidase activity that cleaves the major β-ether units in lignin, thereby releasing the flavonoid tricin from monocot lignin and enhancing the digestion of lignocellulose by polysaccharidase mixtures. We show that the enzyme, which holds potential for the biorefining industry, is widely distributed among lignocellulose-degrading fungi from the Sordariomycetes phylum.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Autoregressive long-horizon prediction of plasma edge dynamics *

Accurate modeling of scrape-off layer (SOL) and divertor-edge dynamics is vital for designing plasma-facing components in fusion devices. High-fidelity edge fluid/neutral codes such as SOLPS-ITER capture SOL physics with high accuracy, but their computational cost limits broad parameter scans and long transient studies. We present transformer-based, autoregressive surrogates for efficient prediction of 2D, time-dependent plasma edge state fields. Trained on SOLPS-ITER spatiotemporal data for the KSTAR tokamak, the surrogates forecast electron temperature, electron density, and radiated power over extended horizons. We evaluate model variants trained with increasing autoregressive horizons (1–100 steps) on short- and long-horizon prediction tasks. Longer-horizon training systematically improves rollout stability and mitigates error accumulation, enabling stable predictions over hundreds to thousands of steps and reproducing key dynamical features such as the motion of high-radiation regions. Measured end-to-end wall-clock times show the surrogate is orders of magnitude faster than SOLPS-ITER, enabling rapid parameter exploration. Prediction accuracy degrades when the surrogate enters physical regimes not represented in the training dataset, motivating future work on data enrichment and physics-informed constraints. Overall, this approach provides a fast, accurate surrogate for computationally intensive plasma edge simulations, supporting rapid scenario exploration, control-oriented studies, and progress toward real-time applications in fusion devices.

autoregressive deep learning↗
DOE OSTI2021

encore : an O ( N g2) estimator for galaxy N -point correlation functions

ABSTRACT We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal {O}(N_\mathrm{g}^2)$ for Ng particles, or $\mathcal {O}(N_\mathrm{FFT}\log N_\mathrm{FFT})$ when using a Fast Fourier Transform with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in ${\sim}100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

79 ASTRONOMY AND ASTROPHYSICS↗