Search NASA⌕ Search

SEARCH · Search NASA

Results for “Transportation Energy data book”

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.

979 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Expanding the horizon of bio-naphtha beyond gasoline blend: property characterization and conversion opportunity assessment through technoeconomic and life-cycle analyses

Bio-naphtha, a common by-product of biorefineries, is expected to experience substantial growth in supply due to increasing demands for renewable diesel and synthetic aviation fuel (SAF). However, demand for bio-naphtha itself as a gasoline blendstock is limited because of the electrification of light-duty vehicles. This work investigated valorization opportunities for bio-naphtha from catalytic fast pyrolysis, hydrothermal liquefaction, Fischer–Tropsch synthesis, and hydrotreated esters and fatty acids pathways. These opportunities include producing polymer-grade olefin via steam cracking, SAF via steam cracking followed by olefin oligomerization, and renewable aromatics benzene, toluene, and xylene (BTX) and hydrogen via catalytic reforming. Process models were developed in Aspen Plus V14 and Aspen HYSYS V14 to calculate the mass and energy balances for each conversion step. Technoeconomic assessment and life-cycle analysis were conducted to evaluate the minimum fuel/product selling price, conversion cost, and life-cycle CO2 equivalent (CO2e) emission reduction. Technoeconomic assessment results suggest a minimum fuel/product selling price as low as $1.9/kg of olefins, $6.30/gal of gasoline-equivalent SAF, and $1.2/kg of BTX without any incentives. For all pathways, these prices are dominated by bio-naphtha feedstock costs, which account for at least 76% of the total cost. Compared with petroleum baselines, bio-naphtha-derived SAF, olefins, and BTX can achieve significant CO2e emission reductions from the use of renewable carbon resources. The results of life-cycle analysis and subsequent technoeconomic assessment, incorporating carbon credits, indicate the economic viability of using bio-naphtha for polymer-grade olefin and BTX production, with product costs comparable to market prices.

Xu, Shuang↗
DOE OSTI2025

Nuclear Reactor Heat Extraction for Synthetic Fuel Plants

This report looks at the viability and best approach when producing synthetic fuels using heat and power supplied by advanced nuclear energy systems. The report looks at a low temperature integration pathway with four types of advanced reactors: a pressurized water reactor (PWR), an advanced light water reactor (A-LWR), a sodium fast reactor (SFR), and a high-temperature gas-cooled reactor (HTGR). A failure modes and effects analysis (FMEA) of the coupling system between the nuclear plant and the synthetic fuel systems is performed with the goal of identifying the reliability of such a thermal delivery system. Furthermore, a high temperature pathway is investigated for synfuel coupling to determine if this is more efficient and more cost effective as a coupling approach.

10 - SYNTHETIC FUELS↗
DOE OSTI2023

Electron Microscopy Studies of Soft Nanomaterials

This review highlights recent efforts on applying electron microscopy (EM) to soft (including biological) nanomaterials. We will show how developments of both the hardware and software of EM have enabled new insights into the formation, assembly, and functioning (e.g., energy conversion and storage, phonon/photon modulation) of these materials by providing shape, size, phase, structural, and chemical information at the nanometer or higher spatial resolution. Specifically, we first discuss standard real-space two-dimensional imaging and analytical techniques which are offered conveniently by microscopes without special holders or advanced beam technology. The discussion is then extended to recent advancements, including visualizing three-dimensional morphology of soft nanomaterials using electron tomography and its variations, identifying local structure and strain by electron diffraction, and recording motions and transformation by in situ EM. On these advancements, we cover state-of-the-art technologies designed for overcoming the technical barriers for EM to characterize soft materials as well as representative application examples. Here, the even more recent integration of machine learning and its impacts on EM are also discussed in detail. With our perspectives of future opportunities offered at the end, we expect this review to inspire and stimulate more efforts in developing and utilizing EM-based characterization methods for soft nanomaterials at the atomic to nanometer length scales in academic research and industrial applications.

Imaging↗
DOE OSTI2025

Manganese‐Based Spinel Cathodes: A Promising Frontier for Solid‐State Lithium‐Ion Batteries

Recently, all-solid-state lithium-ion batteries (ASSLIBs), which exhibit improved safety and enhanced energy density compared to conventional commercialized lithium-ion batteries (LIBs), thereby have garnered extensive research interest. Among the promising cathode candidates, Mn-based spinel cathodes LiMn 2 O 4 (LMO) and LiNi 0.5 Mn 1.5 O 4 (LNMO), with the unique characteristics of low cost, structural stability, and 3D Li-ion diffusion channels, have demonstrated excellent performance in LIBs and presented great potential in ASSLIBs applications. However, several challenges, including structural degradations, poor interfacial contact, large interfacial resistance, and Mn-dissolution/diffusion during the electrochemical cycling, hinder their practical applications and commercialization in the ASSLIBs. Particularly, the high-voltage LNMO cathodes suffer from the challenge of electrochemical incompatibility with most of the solid-state electrolytes (SSEs). Herein, the spinel structure, the electrochemical behavior, and the structural degradation of the LMO/LNMO are explored. The characteristics and recent progress of the mitigating strategies to the challenges of various SSEs, including polymer-, oxide-, composite-, sulfide-, halide-, and LiPON-based SSEs, are introduced when paired with LMO/LNMO. Finally, the directions for future research to advance Mn-based spinel cathodes and fulfill the requirements of the next-generation ASSLIBs are also discussed.

Dou, Yu [Concordia University, Montreal, QC (Canad↗
DOE OSTI2024

Cation valency in water-in-salt electrolytes alters the short- and long-range structure of the electrical double layer

Highly concentrated aqueous electrolytes (termed water-in-salt electrolytes, WiSEs) at solid-liquid interfaces are ubiquitous in myriad applications including biological signaling, electrosynthesis, and energy storage. This interface, known as the electrical double layer (EDL), has a different structure in WiSEs than in dilute electrolytes. Here, we investigate how divalent salts [zinc bis(trifluoromethylsulfonyl)imide, Zn(TFSI) 2 ], as well as mixtures of mono- and divalent salts [lithium bis(trifluoromethylsulfonyl)imide (LiTFSI) mixed with Zn(TFSI) 2 ], affect the short- and long-range structure of the EDL under confinement using a multimodal combination of scattering, spectroscopy, and surface forces measurements. Raman spectroscopy of bulk electrolytes suggests that the cation is closely associated with the anion regardless of valency. Wide-angle X-ray scattering reveals that all bulk electrolytes form ion clusters; however, the clusters are suppressed with increasing concentration of the divalent ion. To probe the EDL under confinement, we use a Surface Forces Apparatus and demonstrate that the thickness of the adsorbed layer of ions at the interface grows with increasing divalent ion concentration. Multiple interfacial layers form following this adlayer; their thicknesses appear dependent on anion size, rather than cation. Importantly, all electrolytes exhibit very long electrostatic decay lengths that are insensitive to valency. It is likely that in the WiSE regime, electrostatic screening is mediated by the formation of ion clusters rather than individual well-solvated ions. This work contributes to understanding the structure and charge-neutralization mechanism in this class of electrolytes and the interfacial behavior of mixed-electrolyte systems encountered in electrochemistry and biology.

Science & Technology - Other Topics↗
DOE OSTI2026

Atomic faulting drives exceptional toughness in low thermal expansion chromium alloys

Endowing functional properties with mechanical responses in traditional metals has been a frontier topic, akin to transforming base metal into gold. Chromium and its alloys, with their functional deficiencies and limited ductility, serve as typical examples. Herein, we report a Cr96Fe4Ge1.3B1 alloy that unifies low thermal expansion (LTE, αl = 1.79 × 10-6 K-1, 200 − 315 K) with exceptional toughness (240.2 J·cm-3). The enhancement in mechanical responses is primarily attributed to layered Cr2B intermetallic precipitates, which ameliorate interfacial cohesion and simultaneously refine the grain structure. The weakened interlayer interactions within the Cr-B layers facilitate the nucleation and movement of numerous tiny stacking faults in precipitates, efficiently alleviating strain energy and resulting in marked work-hardening ability. Additionally, antiferromagnetic fluctuations in the BCC matrix contribute to the unique LTE behavior. This paves the way for the design of high-performance alloys featuring layered-symmetry precipitates.

Yu, Chengyi [University of Science and Technology ↗
DOE OSTI2022

Efficient electrocatalytic valorization of chlorinated organic water pollutant to ethylene

Electrochemistry can provide an efficient and sustainable way to treat environmental waters polluted by chlorinated organic compounds. However, the electrochemical valorization of 1,2-dichloroethane (DCA) is currently challenged by the lack of a catalyst that can selectively convert DCA in aqueous solutions into ethylene. Here, in this work, we report a catalyst comprising cobalt phthalocyanine molecules assembled on multiwalled carbon nanotubes that can electrochemically decompose aqueous DCA with high current and energy efficiencies. Ethylene is produced at high rates with unprecedented ~100% Faradaic efficiency across wide electrode potential and reactant concentration ranges. Kinetic studies and density functional theory calculations reveal that the rate-determining step is the first C–Cl bond breaking, which does not involve protons—a key mechanistic feature that enables cobalt phthalocyanine/carbon nanotube to efficiently catalyse DCA dechlorination and suppress the hydrogen evolution reaction. The nanotubular structure of the catalyst enables us to shape it into a flow-through electrified membrane, which we have used to demonstrate >95% DCA removal from simulated water samples with environmentally relevant DCA and electrolyte concentrations.

Choi, Chungseok [Yale Univ., New Haven, CT (United↗
DOE OSTI2026

Pinning ångström-size solid ionic channels for rare-earth element separation

High-purity rare-earth elements are essential for modern technologies, yet current solvent extraction processes are energy-intensive and environmentally harmful because of inadequate selectivity and ligand toxicity. Although combining size exclusion and binding affinity can improve lanthanide separation, the role of long-range confinement remains underexplored. Here we report lanthanide separation in aqueous systems using extremely confined manganese oxide solid ionic channels with optimized layer spacing. Different lanthanides induce distinct solid-state phase transformations in manganese oxide, creating a strong driving force for separation. Two lanthanide groups, differing by ~1.4 Å in spacing, were identified and confirmed to be stable by density functional theory. The narrower confinement of heavier Group II lanthanides improves cross-group separation by increasing the dehydration barrier for lighter Group I lanthanides without inducing strong binding. Here, we further developed a strategy to pin the confinement dimensions and enhance same-group separation, increasing enrichment factors for La–Nd and La–Pr pairs from 1.6 ± 0.1 and 1.5 ± 0.1 to 5.4 ± 0.1 and 4.2 ± 0.1, respectively.

Chemical engineering↗
DOE OSTI2026

Terahertz‐Nanoscale Visualization of the Microscopic Spin‐Charge Architecture of Colossal Magnetoresistive Switching

Resolving sub-10 nm spin switching and the associated terahertz (THz) electrodynamics during the colossal magnetoresistance (CMR) transition is a definitive frontier in reaching the fundamental spatial, temporal, and energy-dissipation limits of spin-electronics. Yet, simultaneous control of high magnetic field, cryogenic environment, and nanometer resolution has remained an elusive benchmark for THz nanoscopy, leaving the local THz dynamics of these transitions largely unexplored. Here, we overcome these limitations by utilizing a custom-built cryogenic magneto-THz scattering-type scanning near-field optical microscopy (cm-THz-sSNOM) to resolve the near-field THz spectroscopic evolution of the magnetic field-driven CMR transition in a manganite single crystal. Our measurements provide a nanoscale visualization of the THz conductivity, capturing the moment that magnetic-field-induced spin switching triggers the transition from an antiferromagnetic insulator to a ferromagnetic metal. An ellipsoidal near-field model reveals a multi-scale transition initiated by 1–2 nm isolated spin-flip sites at low magnetic fields, which coalesce into ∼15 nm conducting regions as the threshold field is approached. These results provide an nano-THz view of CMR switching, establishing an analysis framework for mapping spin–charge–lattice–orbit–coupled dynamics at spatial scales that transcend the nominal sSNOM resolution.

Haeuser, Samuel [Iowa State University, Ames, IA (↗
DOE OSTI2026

Observed Land Surface Influence on Atmospheric Heat and Moisture Profiles During Interstorms

Land-atmospheric (L-A) feedbacks have historically been studied using models whose structure and parameterizations influence outcomes and insights. The representation of L-A feedbacks based on observations alone remains an ongoing challenge for understanding boundary layer development and precipitation. To address this gap, we use ground-based passive remote sensing and in-situ observations to present an analysis of the atmosphere during 103 interstorm soil moisture drydown events spanning nine warm seasons (2016–2024) in the U.S. Southern Great Plains region. By separating events based on local L-A coupling signals and characterizing the profiles of atmospheric heat and moisture to surface energy flux behavior, we investigate the physical mechanisms linking land surface processes to boundary layer development. We find that during interstorm drydowns, the atmospheric column follows a consistent pattern: moisture increases within the boundary layer, peaks near its top, and declines rapidly above, while warming occurs through the depth. Drydowns that shift toward evaporation produce stronger and deeper thermodynamic responses than cases dominated by sensible heating, which are weaker and shallower. Additionally, moisture is accumulated faster within the boundary layer during shorter drydowns, with longer drydowns representing slower, moisture-limited growth. Drydowns with wetter initial soil moisture will sustain stronger moistening within and above the boundary layer, accelerating buoyancy growth and convective potential toward the next storm. These results provide observational evidence linking surface flux evolution to boundary layer thermodynamics and offer a process-level benchmark for evaluating coupled L-A representations in models and demonstrating the influence of soil moisture on short-term weather forecasting skill.

Zhang, M. S. [Massachusetts Inst. of Technology (M↗
DOE OSTI2023

AutoSourceID-Classifier: Star-galaxy classification using a convolutional neural network with spatial information

Aims.Traditional star-galaxy classification techniques often rely on feature estimation from catalogs, a process susceptible to introducing inaccuracies, thereby potentially jeopardizing the classification’s reliability. Certain galaxies, especially those not manifesting as extended sources, can be misclassified when their shape parameters and flux solely drive the inference. We aim to create a robust and accurate classification network for identifying stars and galaxies directly from astronomical images. Methods.The AutoSourceID-Classifier (ASID-C) algorithm developed for this work uses 32x32 pixel single filter band source cutouts generated by the previously developed AutoSourceID-Light (ASID-L) code. By leveraging convolutional neural networks (CNN) and additional information about the source position within the full-field image, ASID-C aims to accurately classify all stars and galaxies within a survey. Subsequently, we employed a modified Platt scaling calibration for the output of the CNN, ensuring that the derived probabilities were effectively calibrated, delivering precise and reliable results. Results.We show that ASID-C, trained on MeerLICHT telescope images and using the Dark Energy Camera Legacy Survey (DECaLS) morphological classification, is a robust classifier and outperforms similar codes such as SourceExtractor. To facilitate a rigorous comparison, we also trained an eXtreme Gradient Boosting (XGBoost) model on tabular features extracted by SourceExtractor. While this XGBoost model approaches ASID-C in performance metrics, it does not offer the computational efficiency and reduced error propagation inherent in ASID-C’s direct image-based classification approach. ASID-C excels in low signal-to-noise ratio and crowded scenarios, potentially aiding in transient host identification and advancing deep-sky astronomy.

Astronomy & Astrophysics↗
DOE OSTI2026

The 10 September 2025 M w 4.1 Earthquake in Northeastern Utah, United States: An Archetypal Continental Mantle Event

The 10 September 2025 M w 4.1 earthquake in northeastern Utah, United States, had a focal depth 68 km beneath sea level, which is ∼20–25 km greater than estimates of local crustal thickness, making it a rare example of a continental mantle earthquake (CME). The focal depth is well resolved from arrival-time inversion (nearest station ∼13 km away) and moment tensor inversion of regional waveforms. Similar to other CMEs in the Intermountain West, there were no obvious aftershocks or foreshocks, and the waveforms were enriched in high-frequency energy. Spectral modeling gives a stress drop of ∼80 MPa and a radiation efficiency of ∼0.08, albeit with large uncertainties. The high stress drop and low radiation efficiency are consistent with a dissipative source process such as thermal runaway. Also similar to previous Intermountain West CMEs, the event occurred along the boundary of the Archean Wyoming craton, where pressure–temperature conditions favor ductile deformation. We hypothesize that edge-driven or regional-scale mantle convection produces increased strain rates near the craton boundary that make either conventional brittle failure or thermal runaway feasible at relatively high pressure–temperature conditions. High conductivity inferred around the edge of the craton may suggest that fluids also contribute to CME occurrence.

Koper, Keith D. [Univ. of Utah, Salt Lake City, UT↗
DOE OSTI2026

Insulating ground state and 2−𝑘 magnetic structure of candidate Weyl hydrogen-atom K2⁢Mn3⁢(AsO4)3

The ideal Weyl "hydrogen-atom" semimetal exhibits only a single pair of Weyl nodes and no other trivial states at the Fermi energy. Such a material would be a panacea in the study of Weyl quasi-particles, allowing direct unambiguous observation of their topological properties. The alluaudite-like K2⁢Mn3⁢(AsO4)3 compound was recently proposed as such a material. Here, we use comprehensive experimental work and first-principle calculations to assess this prediction. We find K2⁢Mn3⁢(AsO4)3 crystallizes in the 𝐶⁢2/𝑐 symmetry with a quasi-one-dimensional Mn sublattice, growing as small needle-like crystals. Bulk property measurements reveal magnetic transitions at ≈8 and ≈4 K, which neutron scattering experiments show correspond to two distinct magnetic orders, first a partially ordered ferrimagnetic 𝐤𝟏=(0,0,0) structure at 8 K and a second transition of 𝐤𝟐=(1,0,0) at 4 K to a fully ordered state. Below the second transition, both ordering vectors are necessary to describe the complex magnetic structure with modulated spin magnitudes. Both of the best-fit magnetic structures in this work are found to break the symmetry necessary for the generation of Weyl nodes, though one of the magnetic structures allowed by 𝐤𝟏 does preserve this symmetry. However, the crystals are optically transparent and ellipsometry measurements reveal a large band gap, undermining expectations of semimetallic behavior. Density functional theory calculations predict an insulating antiferromagnetic ground state, in contrast to previous reports, and suggest potential frustration on the magnetic sublattice. Given the wide tunability of the alluaudite structure, we consider ways to push the system closer to a semimetallic state.

Taddei, Keith M [Argonne National Laboratory]↗
DOE OSTI2026

Reaction Optimization for Enzymatic Deconstruction of Industrially Relevant Nylon Composites

Plastics such as polyamides (PAs) possess unique physicochemical properties that make them indispensable in modern society. However, their energy‐intensive production and challenging end‐of‐life management highlight the urgent need for efficient recycling or remanufacturing solutions. Enzymatic depolymerization offers a promising route toward circular recycling, yet remains constrained by limited enzyme characterization, lack of validation under industrially relevant conditions and substrates, and overall performance. Here, we optimized the reaction conditions for three recently discovered nylon‐degrading enzymes. One of them, Nyl12, achieved product titers with PA6 and PA66 that exceed previously reported values, without enzyme engineering or substrate pretreatment. We further demonstrated the scalability of the process and its application to complex PA‐based materials used in microelectronic components. Analysis of substrate features, including surface area and particle size, revealed key parameters governing enzymatic activity and provided a framework for future pretreatment and process optimization efforts. In combination, these efforts provide a new benchmark for enzymatic nylon recycling.

nylon↗
DOE OSTI2024

Localizing the non-thermal X-ray emission of PSR J2229+6114 from its multi-wavelength pulse profiles

Context.Pulsars are detected over the whole electromagnetic spectrum, from radio wavelengths up to very high energies, in the GeV-TeV range. While the radio emission site for young pulsars is well constrained to occur at altitudes about several percent of the light-cylinder radius andγ-ray emission is believed to be produced in the striped wind, outside the light cylinder, the non-thermal X-ray production site remains unknown. Aims.The aim of this Letter is to localize the non-thermal X-ray emission region based on a multi-wavelength pulse profile fitting for PSR J2229+6114, which stands as a particularly good candidate given its high X-ray brightness. Methods.Based on the geometry deduced from the joint radio andγ-ray pulse profiles, we fixed the magnetic axis inclination angle and the line-of-sight (LoS) inclination angle. However, we left the region of X-ray emission unlocalized, setting it somewhere between the surface and the light cylinder. We localized this region and its extension by fitting the X-ray pulse profile as observed by the NICER,NuSTAR, and RXTE telescopes in the ranges of 2–7 keV, 3–10 keV, and 9.4–22.4 keV, respectively. Results.We constrained the non-thermal X-ray emission to arise from altitudes between 0.2 r L and 0.55 r L wherer L is the light-cylinder radius. The magnetic obliquity is approximatelyα ≈ 45° −50° and the LoS inclination angle isζ ≈ 32° −48°. Conclusions.This Letter is among the first works to offer a tight constraint on the location of non-thermal X-ray emission from pulsars. We plan to apply this procedure to several other promising candidates to confirm this new result.

Astronomy & Astrophysics↗
DOE OSTI2026

Fabrication of freeform potassium dihydrogen phosphate crystals by belt-on-wheel polishing for spatially tailored polarization control in high-power lasers

Large-aperture (>30 cm) optical components that provide spatially tailored control of the properties of a laser beam, such as phase, polarization, and amplitude, have been the focus of much research and development in recent years. Such optics can improve energy throughput and target implosion efficiency in high-power laser systems conducting inertial confinement fusion research. Here, a method is demonstrated for fabricating freeform surface topographies into the birefringent crystalline material, potassium dihydrogen phosphate (KDP). Belt-on-wheel polishing with an oil-based fluid, developed to mitigate the deliquescent properties of KDP, is used to deterministically polish a freeform topography into KDP, enabling the fabrication of a wave plate with spatially arbitrary retardance. Testing of the laser-induced damage threshold of belt-on-wheel polished KDP was conducted at a wavelength of 351 nm and a pulse duration of 1 ns. The results revealed that the polished surfaces are highly resistant to laser-induced damage, making them suitable for large-aperture, high-power laser applications.

Urban, Nathaniel D. [Univ. of Rochester, NY (Unite↗
DOE OSTI2023

A Self‐Healing System for Polydicyclopentadiene Thermosets

Self-healing offers promise for addressing structural failures, increasing lifespan, and improving durability in polymeric materials. Implementing self-healing in thermoset polymers faces significant manufacturing challenges, especially due to the elevated temperature requirements of thermoset processing. To introduce self-healing into structural thermosets, the self-healing system must be thermally stable and compatible with the thermoset chemistry. This article demonstrates a self-healing microcapsule-based system stable to frontal polymerization (FP), a rapid and energy-efficient manufacturing process with a self-propagating exothermic reaction (≈200 °C). A thermally latent Grubbs-type complex bearing two N-heterocyclic carbene ligands addresses limitations in conventional G2-based self-healing approaches. Under FP's elevated temperatures, the catalyst remains dormant until activated by a Cu(I) co-reagent, ensuring efficient polymerization of the dicyclopentadiene (DCPD) upon damage to the polyDCPD matrix. The two-part microcapsule system consists of one capsule containing the thermally latent Grubbs-type catalyst dissolved in the solvent, and another capsule containing a Cu(I) coagent blended with liquid DCPD monomer. Using the same chemistry for both matrix fabrication and healing results in strong interfaces as demonstrated by lap-shear tests. In an optimized system, the self-healing system restores the mechanical properties of the tough polyDCPD thermoset. Self-healing efficiencies greater than 90% via tapered double cantilever beam tests are observed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Growth of deuterium supersaturated surface layer with increasing ion flux and fluence in plasma-exposed tungsten

Deuterium supersaturated surface layer (DSSL) in tungsten, a few nm thick layer exhibiting extremely high deuterium content (> 5%), has been studied as a function of deuterium plasma ion flux and fluence. Tungsten samples were exposed at 400 K and deuterium ion energy of ∼ 60 eV. The deuterium ion flux spanned over an order of magnitude, being 7.3 × 10 20 , 4.2 × 10 21 , and 3.8 × 10 22 D/m 2 . The fluence ranged from 3.4 × 1024 to 3.4 × 1025 D/m2. The samples and their deuterium content were analyzed by nuclear reaction analysis (NRA), scanning transmission electron microscopy (STEM), and thermal desorption spectroscopy (TDS). The largest thickness of DSSL was found to be around 10.5 nm and was observed in the case of the highest exposure flux and fluence, whereas the layer was only about 3.8 nm thick in the case of the lowest value of the two parameters. The thickness of the DSSL was found to monotonically increase with both deuterium flux and fluence. Finally, similar to the thickness, the estimated D concentration seems to follow the same trend, increasing from the lowest value of around 3 at.% to the highest value of around 5 at.%.

Deuterium↗