Search NASA⌕ Search

SEARCH · Search NASA

Results for “Double diffusion maps”

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.

285 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Li 21 Ge 8 P 3 S 34 : New Lithium Superionic Conductor with Unprecedented Structural Type

Abstract Lithium superionic conductors are pivotal for enabling all‐solid‐state batteries, which aim to replace liquid electrolytes and enhance safety. Herein, we report the discovery of an unprecedented lithium superionic conductor, Li 21 Ge 8 P 3 S 34 , featuring a novel structural type and a new composition in the Li–Ge–P–S system. This material exhibits high lithium ionic conductivity of approximately 1.0 mS cm −1 at 303 K with a low activation energy of 0.20(1) eV. It's unique crystal structure was elucidated using three‐dimensional electron diffraction (3D ED) and further refined through combined powder X‐ray and neutron diffraction analyses. The structure consists of alternating two‐dimensional slabs: one of corner‐sharing GeS 4 tetrahedra and the other of isolated PS 4 tetrahedra, enabling efficient lithium‐ion transport through a tetrahedrally interconnected network of 1D, 2D, and 3D diffusion pathways. This distinctive structural motif provides a novel design strategy for next‐generation solid electrolytes, broadening the structural landscape of lithium superionic conductors. With further advancements in compositional tuning and interfacial engineering, Li 21 Ge 8 P 3 S 34 could contribute to the development of high‐performance all‐solid‐state batteries.

Chemistry↗
DOE OSTI2026

Unlocking superplasticity in medium and high-entropy alloys

Superplasticity, the capacity of materials to sustain extraordinary tensile elongations at elevated temperatures, underpins a range of advanced metal-forming technologies. Conventionally, it is achieved in fine-grained alloys where deformation is dominated by grain-boundary sliding, accommodated by diffusion and dislocation activity. The advent of medium- and high-entropy alloys (M/HEAs), with their high chemical complexity and unconventional phase stability, offers new pathways to superplastic behavior beyond traditional alloy systems. Although investigated only recently, several M/HEAs already exhibit elongations that rival or exceed those of classical superplastic materials, particularly when ultrafine or metastable microstructures are engineered. Here, we review progress in understanding superplastic deformation in M/HEAs, emphasizing the interplay among composition, initial microstructure, thermomechanical processing, and microstructural evolution during high-temperature deformation. We discuss approaches to generating the fine-grained structures necessary for grain-boundary sliding, including severe plastic deformation and tailored heat treatments. We further highlight dynamic phenomena such as phase transformations, evolving grain-boundary chemistry, and deformation-induced grain refinement that can enhance plasticity in these systems. These mechanisms often shift the balance of deformation processes, enabling large elongations even outside classical criteria. Finally, we outline key challenges for application, including cost, scalability, recyclability, and microstructural stability.

klenam, Desmond [University of the Witwatersrand, ↗
DOE OSTI2026

High-resolution modeling of indoor radon exposure with uncertainty quantification in Utah

Indoor radon accounts for 37% of population-level exposure to ionizing radiation in the United States. However, radon metrics are typically reported at coarse spatial scales, potentially obscuring meaningful local variation. We developed a high-resolution modeling framework to estimate indoor radon concentrations across Utah while explicitly quantifying predictive uncertainty. A total of 19,497 residential radon measurements collected between 2006 and 2017 were combined with environmental and housing characteristics and analyzed using a geospatial neural network that accommodates spatial dependence and nonlinear associations. Predictions were generated on a uniform hexagonal grid at 0.73 km2 resolution (H3 level 8). Out-of-sample predictions aggregated to the H3 level 8 grid showed good agreement with observed concentrations (Pearson r=0.64), while household-level predictions exhibited more moderate agreement (r=0.45). The model produced well-calibrated uncertainty estimates, with 24.1% of held-out observations exceeding the predicted 75th-percentile threshold. Maps of predicted radon concentrations and the probability of exceeding the U.S. EPA action level of 148 Bq/m3 (4 pCi/L) revealed substantial fine-scale spatial heterogeneity that was not apparent in conventional coarse-resolution summaries, with greater local variability observed in densely monitored urban counties than in sparsely sampled regions. High-resolution radon models that explicitly quantify uncertainty provide a useful framework for characterizing the spatial distribution of indoor radon and identifying areas of elevated exceedance risk. These findings highlight the value of fine-scale monitoring data and uncertainty-aware modeling approaches for radon exposure assessment, environmental risk characterization, and radon-related health research.

Wu, Yunhan [ORNL] (ORCID:0000000178842994)↗
DOE OSTI2026

Automated scanning probe microscopy of combinatorial ferroelectric libraries: Gaussian-process-guided exploration and noise-aware experiment planning

Combinatorial materials libraries provide an efficient route for mapping composition–property relationships, but their broader impact depends on rapid, quantitative, and functionally relevant characterization. Scanning Probe Microscopy (SPM), including piezoresponse force microscopy (PFM), offers significant potential for quantitative, functionally relevant combi-library readouts. Here, we implement a fully automated SPM workflow for ferroelectric combinatorial libraries and benchmark Gaussian-process-based Bayesian optimization strategies for autonomous experiment planning. The workflow integrates automated probe motion, contact optimization, imaging, and dual amplitude resonance tracking-PFM spectroscopy, and uses scalarized spectroscopic observables to guide subsequent measurements. Stage motion, probe engagement, in-contact tuning, imaging, spectroscopy, and the choice of the next measurement location all proceed without human input. We demonstrate the approach on Sm-doped BiFeO 3 and Zn x Mg 1−x O libraries. By comparing vanilla Bayesian optimization with a measured-noise variant, we show that explicit treatment of local reproducibility can improve modeling of composition-dependent response when the measured variance is physically meaningful, but can also reduce robustness when variability is dominated by outliers or topographic artifacts. Furthermore, these results establish automated SPM as a bridge between combinatorial synthesis and quantitative functional characterization.

Liu, Yu [University of Tennessee, Knoxville, TN (U↗
DOE OSTI2023

Atomistic‐Level Effects of Noncovalent Interactions and Crystalline Packing for Organic Material Structural Integrity upon Exposure to Gamma Radiation

Developing an atomistic understanding of ionizing radiation induced changes to organic materials is necessary for intentional design of greener and more sustainable materials for radiation shielding and detection. Cocrystals are promising for these purposes, but a detailed understanding of how the specific intermolecular interactions within the lattice upon exposure to radiation affect the structural stability of the organic crystalline material is unknown. This study evaluates atomistic-level effects of γ radiation on both single- and multicomponent organic crystalline materials and how specific noncovalent interactions and packing within the crystalline lattice enhance structural stability. Dose studies were performed on all crystalline systems and evaluated via experimental and computational methods. Changes in crystallinity were evaluated by p-XRD and free radical formation was analyzed via EPR spectroscopy. Type of intermolecular interactions and packing within the crystal lattice was delineated and related to the specific free radical species formed and the structural integrity of each material. Periodic DFT and HOMO-LUMO surface mapping calculations provided atomistic-level identifications of the most probable sites for the radicals formed upon exposure to γ radiation and relate intermolecular interactions and molecular packing within the crystalline lattice to experimental results.

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

Decoupling thermal and irradiation effects on grain boundary segregation

Radiation-induced segregation (RIS) is most often measured by peak solute concentration at a boundary. However, this may give an incomplete picture of segregation quantity and phenomena. Radiation-induced and thermal segregation at grain boundaries was investigated in Fe-9.6 at.% Cr after 9 MeV Fe 3+ ion irradiation at 400 °C. The experimental results were compared to kinetic Monte Carlo (KMC) simulations. The study revealed that Cr enrichment (peak segregation) at the grain boundaries was comparable in both the irradiated and non-irradiated conditions, although irradiation resulted in broader segregation profiles in both experiments and simulations, indicating an overall increase in grain boundary segregation due to irradiation. This broadening is attributed to back diffusion into the grain interior. While it is an established phenomenon, this study offers a quantitative evaluation using experimental data and KMC modeling. Further, these results emphasize the importance of analyzing the entire segregation profile and decoupling the thermal and irradiation contributions to solute segregation.

Grain boundary segregation↗
DOE OSTI2026

The impact of chemistry on anion migration in bixbyite-structured lanthanide oxides

This manuscript describes atomistic calculations of oxygen vacancy and interstitial migration in bixbyite structured lanthanide oxides. We examine two types of compounds, one in which only one type of lanthanide cation is present and a second class in which two lanthanides are present in a 3:1 ratio as dictated by the symmetry of the bixbyite crystal structure. Using temperature accelerated dynamics and the nudged elastic band method, we quantify the role of chemistry on the energy barriers for the most important pathways for both vacancy and interstitial migration. We then analyze the impact of these variations on the overall diffusivity of each defect, quantifying the contribution of each pathway using the theory of kinosons. We find that vacancy mobility can vary by as much as three orders of magnitude through changes in chemistry at 500 K. Changes in interstitial mobility are more modest but can still vary by an order of magnitude. This points to the ability to tune the mass transport characteristics of these compounds through appropriate choices in chemistry. We have also included supplementary information containing the atomic structures of the relevant pathways.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Circumventing data imbalance in magnetic ground state data for magnetic moment predictions

Abstract Magnetic materials play a crucial role in the transition to more sustainable forms of energy and electric vehicles. There is an anticipated shortage in magnetic materials in the future, and as a result there is an urgent need to discover and design new magnetic materials. Computational magnetic material design using density functional theory is daunting because of the challenge in identifying magnetic ground states from a combinatorially large set of possibilities. Machine learning offers a path forward by enabling efficient surrogate models that can more readily enumerate these states, but there is a dearth of training data available, and what is available tends to be imbalanced with too much non-magnetic data. In this work we show that the discrete and previously tackled data imbalance that exists at the level of the magnetic ordering leads to an imbalanced continuous distribution with many zeros when the data is unraveled at the atomic magnetic moment level, which subsequently leads to models with low accuracy for magnetic properties. We mitigate this by using a two-part model framework. Our scheme is able to classify atoms into magnetic and non-magnetic with an F1 score and Matthew’s correlation coefficient (MCC) of ~91% and then to provide an implicit embedding representation that maps directly onto the magnitude of the magnetic moment with a mean absolute error of 0.1 μ B . Beyond screening for new magnetic materials, we demonstrate an additional practical use case of our scheme: the provision of good initial guesses for magnetic moments in first-principles electronic relaxations. Such initialization is shown to lead to faster convergence to configurations that lie closer to the ground state.

Computer Science↗
DOE OSTI2025

Quantum graph models for transport in filamentary switching

The formation of metallic nanofilaments bridging two electrodes across an insulator is a mechanism for resistive switching. Examples of such phenomena include atomic synapses, which constitute a distinct class of memristive devices the behavior of which is closely tied to the properties of the filament. Until recently, experimental investigation of the low-temperature regime and quantum transport effects has been limited. However, with growing interest in understanding the true impacts of the filament on device conductance, comprehending quantum effects has become crucial for quantum neuromorphic hardware. Here, we discuss quantum transport resulting from filamentary switching in a narrow region where the continuous approximation of the contact is not valid, and only a few atoms are involved. In this scenario, the filament can be represented by a graph depicting the adjacency of atoms and the overlap between atomic orbitals. Using the theory of quantum graphs with locally diffusive node scattering, we calculate the scattering amplitude of charge carriers on this graph and explore the interplay between filamentary formation and quantum transport effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2025

Airports as Energy Nodes Activity Summary

Advanced aircraft concepts that use non-traditional aviation energy storage methods such as batteries or cryogenic hydrogen are in development and expected to enter regular service at airports worldwide within the next decade. The energy needs for these aircraft may quickly overwhelm the existing energy infrastructure at airports, particularly at smaller and more remote facilities. Without energy upgrades, these airports will not be able to host these advanced vehicles, but without the advanced vehicle traffic, these airports will not have the rationale or funding to build up their energy infrastructure. The Airports as Energy Nodes (ÆNodes) activity, a collaboration between the National Aeronautics and Space Administration (NASA) and the National Renewable Energy Laboratory (NREL), was executed to understand and model the energy needs that advanced aircraft concepts may levy on these smaller airports, determine cost-effective approaches to enhance the airport energy infrastructure, and demonstrate the enhanced resilience of these energy infrastructure upgrades to the airport and surrounding community via “digital twin” simulation at relevant energy and dynamic time scales. The ÆNodes team also investigated future reference aircraft designs and materials to enable cryogenic hydrogen storage for aircraft. The ÆNodes team conducted analysis at two U.S. airport partner sites — Winchester Regional Airport in Winchester, Virginia, and Tweed/New Haven Airport in New Haven, Connecticut. The goal of this partnership was to develop data and reference infrastructure designs that could accommodate advanced aircraft in the future at these airports while also enhancing the resiliency of the energy supply to the surrounding airport community, which could be used to capture funding to enable the infrastructure upgrades. Over the course of the study, a method was developed to estimate air traffic requiring advanced energy services over the course of a year using a mix of historical data and companion studies on advanced aircraft transportation networks. The study has concluded at NASA but continues at NREL, who will develop a final report discussing the energy infrastructure upgrades and digital twin results. Preliminary results indicate that unrestricted adoption of advanced battery-electric aircraft may double traffic at these airports and increase peak daily power usage by an order of magnitude, while increase electricity energy needs by a factor of two to four. The infrastructure upgrades necessary to accommodate these increased energy needs could be used to provide enhanced energy services to the airport community to offset the cost and increase the utility of the upgrades, which will be described in the NREL final report.

Airports↗
NASA NTRSDate not supplied

SERENE: Saturn Enceladus Return Explorer with Nuclear Electric Propulsion

A ‘quick’ Enceladus sample return mission concept was developed based on the scientist recommendations at the recent ‘Accelerating Space Science with Nuclear Technology Workshop’. The Nuclear Electric Propulsion spacecraft assumed a follow-on 40 kWe nuclear reactor using the demonstrated 40 kWe Fission Surface Power system, expected in the early 2030s. The NEP vehicle also utilized a set of NEXT-C ion thrusters as well as planned Artemis commercial launchers. By launching the 40 kW NEP vehicle on a Starship and adding the propellants of 15 tankers, the 27t probe could be sent on a direct trajectory to Saturn (no Earth or Jupiter flybys) where NEP was used for Saturn capture, spiral down, spiral up and return to the Earth. A small lander obtained the surface Enceladus sample. Using the 40 kW NEP provided a round-trip time of only 16.5 years. A second option was more attractive from a science perspective whereby the NEP vehicle would deliver a large, 6t chemical lander to low Enceladus orbit where it would grab and return a sample to Earth (similar to the recent Orbilander design but in reverse). After deploying the lander, the NEP vehicle would stay in Saturn space performing a moon tour by orbiting four more moons and mapping the large moon of Titan. This option took slightly longer (18.5 yrs) due to the chemical return leg limitations. Both options demonstrated the agility, payload capability, and sample return goals the workshop recommended. An all-chemical option with two stages was roughly analyzed but took 21.5 yrs and required a Jupiter gravity assist.

Nuclear Electric Propulsion↗
DOE OSTI2024

Thermal stability of 3D interface Cu/Nb nanolaminates

Nanocrystalline alloys are promising structural materials yet lack thermal stability in many cases. Recent work shows that interface structure has an outsize effect on the thermal behavior of nanostructured alloys. Here, this work focuses on the role of controlled heterophase interface structure in the thermal evolution of model Cu/Nb nanolaminates. We introduce 3D interfaces containing nanoscale heterogeneities in all spatial dimensions between Cu and Nb, forming 3D Cu/Nb. TEM, nanoindentation, and DSC are used in tandem to establish thermal stability and to identify shifts in microstructure as a function of static annealing temperature. 3D interfaces are shown to survive annealing to 300 °C for 1 hr., while 3D Cu/Nb microstructure evolves to form low-density and voided regions correlating to the onset of layer pinch-off between 500 and 600 °C annealing temperatures. A diffusivity- and vacancy energetics-based mechanism is developed to explain void formation driven by 3D interface degradation at elevated temperature.

3D interfaces↗
DOE OSTI2025

Active Fluids Form System-Spanning Filamentary Networks

Recent experimental realizations of liquid-liquid phase separation of active liquid crystals have offered an insight into the interaction between phase separation, ubiquitous in soft matter and biology, and chaotic active flows. Here, in this Letter, we use continuum theory to examine phase separation of an active liquid crystal and a passive fluid and report two new results. First, we provide an analytical derivation of the activity-induced suppression of the phase boundary of the coexistence region—a result first reported in simulations and experiments. We show that the shift in the critical point is a result of the balance between self-stirring active flows and phase-separating diffusive fluxes. Second, we show that this same balance is responsible for dramatically changing the morphology of the phase separated state, resulting in the emergence of a new mixed active phase consisting of a dynamical filamentous active network that invades the entire system area, trapping droplets of passive material. This structure exists even for very low volume fractions of active material. Our work provides an important step towards the goal of understanding how to use activity as a new handle for sculpting interfaces.

active nematics↗
DOE OSTI2025

Chemical Bond Covalency in Superionic Halide Solid‐State Electrolytes

Abstract Halide solid‐state electrolytes (SSEs) are promising superionic conductors with high oxidative stability and ionic conductivity, making them attractive for all‐solid‐state lithium‐ion batteries. However, most studies have focused on ion‐stacking structures, overlooking the role of bond characteristics in ionic transport. Here, we investigate bond dynamics and the superionic transition (SIT) in bromide electrolyte, Li 3 InBr 6 , using synchrotron X‐ray techniques and ab initio molecular dynamics (AIMD) simulations. We demonstrate that the SIT in halide SSEs is driven by a thermally induced transition in bonding character (ionic to covalent) rather than a change in crystal phase. AIMD simulations further reveal enhanced Li⁺ diffusion and collective anion motion at elevated temperatures. Expanding our study to Li 3 LnBr 6 (Ln = Gd, Tb, Ho, Tm, and Lu), we confirm the widespread occurrence of SIT in this material class, with Li 3 GdBr 6 exhibiting the highest ionic conductivity (5.2 mS cm −1 at 298 K). More importantly, the ionic‐covalent transition is highly tunable through electrolyte modifications, such as cation/anion substitution and synthesis methods. Our findings provide a new perspective on ionic transport, highlighting the critical role of chemical bond characteristics in halide SSEs.

Chemistry↗
DOE OSTI2025

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems↗
DOE OSTI2025

Internally Catalyzed Hydrogen Atom Transfer (I-CHAT)—A New Class of Reactions in Combustion Chemistry

The current paradigm of low-T combustion and autoignition of hydrocarbons is based on the sequential two-step oxygenation of fuel radicals. The key chain-branching occurs when the second oxygenation adduct (OOQOOH) is isomerized releasing an OH radical and a key ketohydroperoxide (KHP) intermediate. The subsequent homolytic dissociation of relatively weak O–O bonds in KHP generates two more radicals in the oxidation chain leading to ignition. Based on the recently introduced intramolecular “catalytic hydrogen atom transfer” mechanism (J. Phys. Chem. 2024, 128, 2169), abbreviated here as I-CHAT, we have identified a novel unimolecular decomposition channel for KHPs to form their classical isomers—enol hydroperoxides (EHP). The uncertainty in the contribution of enols is typically due to the high computed barriers for conventional (“direct”) keto–enol tautomerization. Remarkably, the I-CHAT dramatically reduces such barriers. The novel mechanism can be regarded as an intramolecular version of the intermolecular relay transfer of H-atoms mediated by an external molecule following the general classification of such processes (Catal. Rev.-Sci. Eng. 2014, 56, 403). Here, we present a detailed mechanistic and kinetic analysis of the I-CHAT-facilitated pathways applied to n-hexane, n-heptane, and n-pentane models as prototype molecules for gasoline, diesel, and hybrid rocket fuels. We particularly examined the formation kinetics and subsequent dissociation of the γ-enol-hydroperoxide isomer of the most abundant pentane-derived isomer γ-C5-KHP observed experimentally. To gain molecular-level insight into the I-CHAT catalysis, we have also explored the role of the internal catalyst moieties using truncated models. All applied models demonstrated a significant reduction in the isomerization barriers, primarily due to the decreased ring strain in transition states. In addition, the longer-range and sequential H-migration processes were also identified and illustrated via a combined double keto–enol conversion of heptane-2,6-diketo-4-hydroperoxide as a potential chain-branching model. To assess the possible impact of the I-CHAT channels on global fuel combustion characteristics, we performed a detailed kinetic analysis of the isomerization and decomposition of γ-C5-KHP comparing I-CHAT with key alternative reactions—direct dissociation and Korcek channels. Calculated rate parameters were implemented into a modified version of the n-pentane kinetic model developed earlier using RMG automated model generation tools (ACS Omega, 2023, 8, 4908). Simulations of ignition delay times revealed the significant effect of the new pathways, suggesting an important role of the I-CHAT pathways in the low-T combustion of large alkanes.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Field-tunable BKT and quantum phase transitions in spin-$\frac{1}{2}$ triangular lattice antiferromagnet

Quantum magnetism is one of the most active fields for exploring exotic phases and phase transitions. The recently synthesized Na 2 BaCo(PO 4 ) 2 (NBCP) is an ideal material incarnation of the spin-$\frac{1}{2}$ easy-axis triangular lattice antiferromagnet (TLAF). Experimental evidence shows that NBCP hosts the spin supersolid state with a giant magnetocaloric effect. Theory further predicts that magnetic fields can drive NBCP through Berezinskii-Kosterlitz-Thouless (BKT) and other richer quantum phase transitions. However, detecting these transitions is challenging, as they onset at ultralow temperatures near 60 mK and require high magnetization sensitivity. Using a newly developed gradient force magnetometer in a dilution refrigerator, we mapped the magnetic susceptibility phase diagram down to 30 mK. Our results provide a more comprehensive and accurate understanding of BKT melting of spin supersolidity and several field-tunable quantum phase transitions, which establish NBCP as a model platform for frustrated magnetism and highlight potential applications of its giant magnetocaloric effects.

BKT transition↗
DOE OSTI2023

Overcoming the Entropy Penalty of Direct Air Capture for Efficient Gigatonne Removal of Carbon Dioxide

Atmospheric carbon poses an existential threat to civilization via global climate change. Hundreds of gigatonnes of carbon dioxide must be removed from earth’s atmosphere in the next three decades, necessitating a low-cost, energy-efficient process to extract low concentrations of carbon dioxide for conversion to a stable material permanently stored for thousands of years. In this work, the challenge of removing gigatonnes of CO 2 is described via the scale of effort and the thermodynamics of collecting and reducing this diffuse chemical, the accumulation of which imparts a substantial entropy penalty on any atmospheric carbon capture process. The methods of CO 2 reduction combined with upstream direct air capture (DAC) including absorption, membrane separation, and adsorption are compared with biomass torrefaction and permanent burial (BTB). A Monte Carlo model assesses the mass, energy, and economics of the full process of biomass torrefaction from biomass collection and transport to stable carbon burial to determine that 95% of scenarios could remove carbon for less than $200 per CO 2 -tonne-equivalent. Torrefied carbon is further discussed for its long-term stability and availability at the scale required to substantially mitigate the threat of climate change.

biomass↗