Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface features”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Evidence for strong electronic correlations in the bulk state of grey arsenic

Here, we investigate the electron band structure of grey arsenic, whose (111) face hosts the topological Shockley state. Interestingly, the bulk band close to the touching point with the surface state exhibits the characteristics of inelastic scattering. Moreover, the band structure analysis reveals linearity in the imaginary part of electron self-energy. These features are analogous to those observed in high-temperature superconductors and marginal Fermi liquid systems, respectively, where strong electronic correlations exist. Our results suggest that correlated many-body states can be connected by non-interacting topological states, providing a viable playground to explore the coupling between topological and correlated states via grey arsenic surface.

ARPES↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2026

Xanthos-Lake Model Source Code

This repository contains the source code for Xanthos-Lake, a lake-modeling extension of the Xanthos framework that introduces a coupled lake component comprising the Xanthos-Lake Snow and Ice Model (xLSIM) and the Xanthos-Lake Water Balance Model (xLWBM). xLSIM is a basin-aware machine-learning model for lake snow, ice, and thermal conditions. It predicts monthly lake ice thickness, snow depth, snow-cover fraction, mixing-layer temperature, and lake ice fraction from meteorological forcing and lake surface-area information. It uses sequence-based deep-learning architectures, including Transformer and hybrid Long Short-Term Memory–Transformer (LSTM–Transformer) models, together with seasonal encoding, multi-lake learning, physical masking, and basin-level cryospheric and non-cryospheric classification. The training workflow uses Ray for scalable execution and includes optional Ray Tune hyperparameter optimization. Model predictions, observations, diagnostics, and feature-importance outputs are written in NetCDF. xLWBM is the water-balance component of the new lake framework. It simulates monthly lake storage, surface area, evaporation, inflow, outflow, and lake–groundwater exchange. It combines physical water-balance equations with calibrated bathymetric relationships, weir-based outlet flow, modified Penman open-water evaporation, groundwater head relaxation, Penman–Monteith snow and ice sublimation, and snow, ice, and thermal conditions supplied by xLSIM. The model calibrates lake parameters against satellite-derived surface-area data, using evaporation-based calibration where surface-area data are unavailable, and supports small, medium, and large lake classes. For large lakes, xLWBM is integrated with the managed-routing workflow so that lake storage and outflow interact directly with downstream river routing and reservoir operations. Together, xLSIM and xLWBM provide Xanthos with a coupled lake-modeling capability. xLSIM supplies the snow, ice, and thermal conditions that affect lake evaporation and snow- and ice-related water exchanges, while xLWBM translates those conditions into dynamic lake storage, surface area, evaporation, and discharge. In return, xLWBM supplies evolving lake surface area to xLSIM. This coupling enables Xanthos to represent lakes as active hydrologic components within basin-scale water-availability and routing simulations.

Machine Learning↗
DOE OSTI2025

A structural analysis of ordered Cs3Sb films grown on single crystal graphene and silicon carbide substrates

Alkali antimonides are well established as high efficiency, low intrinsic emittance photocathodes for accelerators and photon detectors. However, conventionally grown alkali antimonide films are polycrystalline with surface disorder and roughness that can limit achievable beam brightness. Ordering the crystalline structure of alkali antimonides has the potential to deliver higher brightness electron beams by reducing surface disorder and enabling the engineering of material properties at the level of atomic layers. In this report, we demonstrate the growth of ordered Cs3Sb films on single crystal substrates 3C-SiC and graphene-coated 4H-SiC using pulsed laser deposition and conventional thermal evaporation growth techniques. The crystalline structures of the Cs3Sb films were examined using reflection high energy electron diffraction and x-ray diffraction diagnostics, while film thickness and roughness estimates were made using x-ray reflectivity. With these tools, we observed ordered domains in less than 10 nm thick films with quantum efficiencies greater than 1% at 530 nm. Moreover, we identify structural features such as Laue oscillations indicative of highly ordered films. We found that Cs3Sb films grew with flat, fiber-textured surfaces on 3C-SiC and with multiple ordered domains and sub-nanometer surface roughness on graphene-coated 4H-SiC under our growth conditions. We identify the crystallographic orientations of Cs3Sb grown on graphene-coated 4H-SiC substrates and discuss the significance of examining the crystal structure of these films for growing epitaxial heterostructures in future experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Deep Learning Approaches for Predicting the Surface Tension of Ionic Liquids

Ionic liquids (ILs) are a novel class of solvents that have attracted significant attention due to their unique and tunable properties. Among their physiochemical characteristics, surface tension plays a critical role in various industrial applications including electrolytes, heat transfer fluids, and separation processes. However, because of the exploratory nature of IL design and the vast combinatorial space of possible anion–cation pairs, the experimental determination of these properties is often impractical, being both time-consuming and costly. To overcome these challenges, computational approaches are increasingly employed to develop accurate predictive models that can accelerate IL discovery and design. In this study, we present two deep learning (DL) models for predicting the surface tension of ILs across a broad temperature range at a constant pressure. The models use simplified molecular input line entry system, SMILES, representations of ILs to extract molecular features as inputs. Both DL models demonstrate excellent agreement with experimental data, achieving an R 2 value of 0.990 and a root-mean-square error of 0.792 mN/m. In conclusion, these results offer valuable insights for the rapid screening and rational design of ILs with tailored surface tension values.

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

Operando Contactless EFISH Study of the Rate-Determining Step of Light-Driven Water Oxidation on TiO 2 Photoanodes

For many slow solar-fuel-forming reactions, the accumulation of photogenerated minority carriers on the photoelectrode surface leads to light-induced band edge unpinning, affecting the junction properties by decreasing band bending in the semiconductor space charge layer and increasing the driving force of surface reactions in the electric double layer. In this study, we demonstrate a contactless operando electric field-induced second harmonic generation (EFISH) method for measuring the band bending change (δΔΦ SCR L ) on photoelectrodes upon photoexcitation. For ndoped rutile TiO 2 water oxidation photoanodes at pH 7, δΔΦ SCR L increases at more positive potentials or higher illumination power density until it reaches saturation values. We show that under fast mass transport conditions, δΔΦ SCR L is exclusively attributed to the accumulated charged rate-determining species that can be regarded as temporary surface states, and the relationship between the photocurrent and δΔΦ SCR L can be well modeled by assuming that hole trap states function as the reaction center. Kinetic isotope experiments identify proton-coupled electron transfer as the rate-determining step and suggest a possible chemical nature of the key intermediate. We demonstrate that light-induced band edge unpinning is a beneficial feature under high illumination conditions for oxygen evolution reaction on TiO 2 because it maintains the photon-to-current conversion efficiency by enhancing the surface reaction driving force, shedding light on the actual device application.

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

A simple centrifuge cell method for ex situ quantification of electrical conductivity of slurry electrode materials

We present the design, experimental procedure, and experimental evaluation of a system for fast, simple, and ex situ characterization of electrical conductivity of slurry electrode materials. The system uses a custom-designed electrochemical cell compatible with centrifugation in a swing-bucket centrifuge. The cell features cylindrical graphite electrodes that are partially sheathed so as to expose only 2 mm of the electrode surface to the bottom region of the packed particulate pellet. Also presented is a conduction model that provides a shape factor for estimating effective conductivity. We tested aqueous solutions of carbon black (CB), activated carbon (AC), and mixtures thereof. These particles were dispersed in 0.0 and 0.5 M NaCl solutions. Measurements show that the effective conductivity initially increases linearly with pellet mass and then saturates at higher masses. Notably, CB exhibited a fivefold increase in conductivity than AC at equal pellet masses. CB/AC mixtures at a fixed pellet mass were tested with CB mass fractions of 0 to 100%. Interestingly, the mixture conductivity was found to be a non-monotonic function of CB mass fraction, with a maximum conductivity at about 60 % CB mass fraction. At this maximum, the mixture conductivity is approximately 30 % higher than pure CB. NaCl concentration in the slurry solution had no effects on conductivity. These results highlight the interactions between slurry electrode composition and compaction, offering insights for optimizing slurry electrodes. Furthermore, the system and results may also be applicable to evaluation of particulate materials (including slurries) used for Li-ion batteries, capacitive deionization, fuel cells, and flow electrodes.

Capacitive deionization↗
DOE OSTI2024

A Skipper-CCD Light Shield for X-Ray Detection in Space

We present advancements in the design of Skipper-CCD sensors for X-ray detection in environments with high optical backgrounds, such as those expected in space. These packages incorporate a custom-made aluminum shield on the CCD surface that blocks over 99% of visible light while preserving the efficiency for keV X-rays. These features allow us to perform precise X-ray measurements in environments challenged by visible light interference. Furthermore, we briefly discuss the potential implementation of this design concept in frame-transfer CCDs, opening opportunities for broader applications and advancements in imaging technologies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multi-Functional Smart Structures for Smart Vehicles

This report summarizes the development of a new class of recyclable multi-functional composite materials for production of lightweight smart structures and surfaces. Functional high stiffness conductive composites were processed using molding methods that integrated continuous fiber and additively manufactured features. Methods for integration of sensing functionality and controls were also developed to reduce system cost while providing a new capability for structural health monitoring. This new class of composites is applicable to a broad range of vehicle interior, exterior and battery enclosure systems. By way of demonstration, a vehicle instrument panel cross car beam was developed that provided a 38% mass savings compared to steel while maintaining a cost penalty competitive to alternate lightweight material solutions. These technologies were validated for implementation by a uniquely qualified project team comprising a US automotive OEM, Tier 1 and Tier 2 supplier, with key contributions from Oak Ridge National Lab, Purdue University and Michigan State University.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Effect of Solvents on Lignin–Surface Interactions via Molecular Dynamics Simulations

Lignin, an essential building block of lignocellulosic biomass, is a potential abundant source of aromatic monomers for the polymer and chemical industry. Reductive catalytic fractionation (RCF) is one promising process that can produce high yields of phenolic monomers and oligomers from lignin under different catalytic conditions. An important choice in optimizing RCF is the selection of solvent; however, detailed insights into solvent effects on lignin behaviors and interactions remain limited. Here, in this work, we perform all-atom molecular dynamics simulations to study the solvation of lignin, solvent-mediated conformational changes, and the interaction of solvated lignin oligomers with model surfaces. We focus on the behavior of an oligomeric lignin model compound in methanol, ethanol, a binary mixture of ethanol and water, and water at both the RCF reaction temperature (473 K) and room temperature. Analysis of structural features of lignin suggests that these three organic solvent systems favorably solvate lignin, resulting in a more extended conformation suitable for catalytic conversion to valuable chemicals. We further introduce model palladium (Pd) and carbon (C) surfaces to understand how solvent choice impacts adsorption onto a representative catalytic surface and support, and to quantify the competition among the reactant and solvent molecules for the surface. Unbiased simulations suggest strong adsorption of lignin on both Pd and C surfaces at 473 K, with notable solvent-mediated differences in adsorption energies. Additionally, our findings indicate that lignin adsorption is promoted by the entropy change resulting from the displacement of solvent molecules from the surface. This study provides a molecular perspective of adsorption of lignin onto varying surfaces, which is a step towards understanding and optimizing the catalytic conversion of lignin into valuable chemicals.

adsorption↗
DOE OSTI2024

Modifying the Reactivity of Single Pd Sites in a Trimetallic Sn‐Pd‐Ag Surface Alloy: Tuning CO Binding Strength

Abstract Improving control over active‐site reactivity is a grand challenge in catalysis. Single‐atom alloys (SAAs) consisting of a reactive component doped as single atoms into a more inert host metal feature localized and well‐defined active sites, but fine tuning their properties is challenging. Here, a framework is developed for tuning single‐atom site reactivity by alloying in an additional inert metal, which this work terms an alloy‐host SAA. Specifically, this work creates about 5% Pd single‐atom sites in a Pd 33 Ag 67 (111) single crystal surface, and then identifies Sn based on computational screening as a suitable third metal to introduce. Subsequent experimental studies show that introducing Sn indeed modifies the electronic structure and chemical reactivity (measured by CO desorption energies) of the Pd sites. The modifications to both the electronic structure and the CO adsorption energies are in close agreement with the calculations. These results indicate that the use of an alloy host environment to modify the reactivity of single‐atom sites can allow fine‐tuning of catalytic performance and boost resistance against strong‐binding adsorbates such as CO.

Mohrhusen, Lars↗
DOE OSTI2025

WRF-ELM v1.0: a regional climate model to study land–atmosphere interactions over heterogeneous land use regions

Abstract. The Energy Exascale Earth System Model (E3SM) Land Model (ELM) is a state-of-the-art land surface model that simulates the intricate interactions between the terrestrial land surface and other components of the Earth system. Originating from the Community Land Model (CLM) version 4.5, ELM has been under active development, with added new features and functionality, including plant hydraulics, radiation–topography interaction, subsurface multiphase flow, and more explicit land use and management practices. This study integrates ELM v2.1 with the Weather Research and Forecasting (WRF; WRF-ELM) model through a modified Lightweight Infrastructure for Land Atmosphere Coupling (LILAC) framework, enabling affordable high-resolution regional modeling by leveraging ELM's innovative features alongside WRF's diverse atmospheric parameterization options. This framework includes a top-level driver for variable communication between WRF and ELM and Earth System Modeling Framework (ESMF) caps for the WRF atmospheric component and ELM workflow control, encompassing initialization, execution, and finalization. Importantly, this LILAC–ESMF framework demonstrates a more modular approach compared to previous coupling efforts between WRF and land surface models. It maintains the integrity of ELM's source code structure and facilitates the transfer of future developments in ELM to WRF-ELM. To test the ability of the coupled model to capture land–atmosphere interactions over regions with a variety of land uses and land covers, we conducted high-resolution (4 km) WRF-ELM ensemble simulations over the Great Lakes region (GLR) in the summer of 2018 and systematically compared the results against observations, reanalysis data, and WRF-CTSM (WRF coupled with the Community Terrestrial Systems Model). In general, the coupled WRF-ELM model has reasonably captured the spatial distribution of surface state variables and fluxes across the GLR, particularly over the natural vegetation areas. The evaluation results provide a baseline reference for further improvements in ELM in the regional application of high-resolution weather and climate predictions. Our work serves as an example to the model development community for expanding an advanced land surface model's capability to represent fully-coupled land–atmosphere interactions at fine spatial scales. The development and release of WRF-ELM marks a significant advancement for the ELM user community, providing opportunities for fine-scale regional representation, parameter calibration in coupled mode, and examination of new schemes with atmospheric feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Additive Manufacturing of Thermal Energy Storage Composites with Microencapsulated Phase Change Materials Supported in a Multipolymer Matrix

Additive manufacturing (AM) techniques to directly integrate phase change materials (PCMs) are of interest for efficient thermal energy storage (TES) architectures. Complex, high surface-to-volume ratio composites embedded with PCM can improve thermal management with reduced material waste for customizable device fabrication. Reducing feature sizes of TES-integrated heat exchangers using AM can increase heat transfer without thermal conductivity enhancement. Here, composite AM materials containing 60 wt% microencapsulated phase change materials (MEPCM) are fabricated using off-the-shelf printers at common speeds and resolutions. High MEPCM loading in filaments is achieved with powder extrusion using two polymers, thermoplastic-polyurethane (TPU) and polycaprolactone (PCL), that mediate flexibility and rigidity for effective extrusion and printing without filament fracture or buckling. Furthermore, with PCL and TPU at 20 wt% each and 60 wt% MEPCM (P 20 T 20 M 60 ), smooth, form-stable filaments are consistently printed. Powder-based extrusion displays negligible damaging effects on the MEPCM. Printed P 20 T 20 M 60 demonstrates 105 J/g of energy storage with no degradation through 250 thermal cycles, within 5% of the theoretical storage enthalpy. Combining PCL/TPU shows good interfacial adhesion between print layers and produces high surface area objects, like 15% gyroids, and dense, 100% infilled pucks. Prints are also scalable to a 900 cm 3 honeycomb heat exchanger with an estimated 9 Wh energy storage.

25 ENERGY STORAGE↗
DOE OSTI2025

Modeling the Nucleation and Growth of Lead Sulfate Particles on Lead Electrodes

Lead-acid batteries (LABs) play a pivotal role in the energy storage sector with applications spanning from starting-lighting-ignition batteries to grid energy storage. Passivation of lead negative electrodes by PbSO 4 particles is a fundamental mechanism limiting the performance of LAB. In this regard, an electrochemical model is developed that simulates the nucleation and growth (N&G) dynamics of PbSO 4 particles on a flat lead electrode, responsible for its passivation. The model considers the electrochemical reactions between lead electrode and sulfuric acid, N&G of PbSO 4 particles, passivation of the lead surface, and the ternary transport of PbSO 4 (aq), bisulfate, and protons in H 2 SO 4 electrolyte. The model is validated with a dataset of cyclic voltammetry (CV) responses collected at several scan rates and H 2 SO 4 concentrations. The model shows remarkable qualitative and quantitative agreement with the experimental data including CV peak features, discharge capacity, and particle size. The model was employed to explore key N&G quantities, such as supersaturation, nucleation rate, particle count, growth rate, particle size, and surface coverage, and to examine how their interactions influence electrode utilization. Parametric studies were also conducted to evaluate how scan rates and acid concentrations influence the previously mentioned N&G quantities and, subsequently, the utilization of the electrode.

25 ENERGY STORAGE↗
DOE OSTI2026

Electronic correlations and topology in Kondo insulator PuB 6

Utilizing a combination of dynamical mean field theory (DMFT) and density functional theory, it has been theoretically proposed that PuB 6 is a strongly correlated topological insulator characterized by nontrivial 𝐙 2 topological invariants and metallic surface states [X. Deng et al., Phys. Rev. Lett. 111, 176404 (2013)]. Here, we demonstrate through low-temperature magnetotransport measurements and first-principles calculations that PuB 6 exhibits characteristics of a topological Kondo insulating state. These features include a transition in electrical resistivity from high-temperature, thermally activated behavior with a narrow gap at the Fermi level (Δ⁢𝜌 ∼ 20 meV) to a distinctive low-temperature plateau, as well as a surface-to-volume dependence of electrical resistivity at low temperatures. The topological nature of PuB 6 is further supported by the theoretical calculations, which show that GGA + 𝑈 is capable of capturing electronic, topological, and lattice properties of PuB 6 with much lower computational cost than DMFT.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Electrosynthesis of high purity ethylene using high-index facet Cu 2 O nanocrystals electrocatalyst

Electrochemical CO 2 reduction reaction (eCO2RR) to multi-carbon (C 2+ ) products with copper-based catalysts is often limited by poor selectivity. This challenge arises from the concurrent formation of various intermediates, dictated by the atomic arrangement and electronic properties of surface atoms. In this study, we found that copper (I) oxide (Cu 2 O) nanocrystals with 50 facets (50F-NC), predominantly featuring (211) facets that offers high density of under-coordinated sites, demonstrate superior ethylene (C 2 H 4 ) selectivity of 92% ± 2 with an overall current density of 212 mA/cm 2 at -650 mV vs RHE. Furthermore, after one month of storage in a 1 M KOH electrolyte, this catalyst demonstrated a C 2 H 4 Faradaic efficiency of 87% highlighting its stabile structure under strong alkaline environments. Here, operando electrochemical Raman spectroscopy revealed enhanced CO* intermediate coverage on the 50F-NC catalyst, correlating with improved C-C coupling. SEM, TEM, and XPS analyses, along with DFT calculations, suggested that Cu sites on the (211) facet of 50F-NC and those at the Cu/Cu 2 O interface formed in-situ due to the surface reconstruction during the reaction, are likely active sites for effective C-C coupling and sustained high-rate C 2 H 4 production.

Cu nano particle↗
DOE OSTI2025

Hydrologic Model Data for the East Fork Poplar Creek Watershed Simulated with the Advanced Terrestrial Simulator (ATS): Streamflow and Network Expansion–Contraction Dynamics

This dataset supports hydrologic modeling and stream network expansion–contraction analysis for the East Fork Poplar Creek (EFPC) Watershed in Tennessee. It includes a Jupyter notebook for model setup, model configuration files, simulation outputs, and derived products used to evaluate model performance and investigate stream dynamics under varying hydrologic conditions. The dataset was generated using the Watershed Workflow Python package and the Advanced Terrestrial Simulator (ATS), enabling integrated surface–subsurface hydrologic simulations using a stream-aligned mesh. Outputs include high-resolution time series of streamflow, active network length, water table depth, and related hydrologic variables. Also included are spatially explicit stream persistency indices and classifications of reaches as perennial or non-perennial. These data facilitate reproducibility and support further research on stream intermittency and variability in network extent.The model data archive is organized in following directories:1) model_setup_inputsContains the Watershed Workflow Jupyter notebooks (accessed through any open source code editor), selected input datasets, and resulting ATS input files, including XML files (access through any open source code editor), computational mesh (.exo files can be viewed using Paraview), and meteorological forcing files (.h5 files can be accessed through h5py python package and HDFView open source software). 2) model_outputsIncludes ATS simulation outputs relevant to this study. Time series of spatially integrated or averaged variables (e.g., streamflow, water table depth) are provided as CSV files. Select spatial fields (e.g., ponded depth and water table depth) are saved as pickled Python objects to reduce file size, and can be accessed through pickle package in Python. Key geometry objects from Watershed Workflow—such as the surface mesh and river tree—are also included to support analysis of streamflow persistency and expansion–contraction dynamics. These files can also be accessed through Watershed Workflow Python package.3) model_evaluationProvides observed streamflow time series and field survey-based flow regime classifications used to evaluate model performance. Jupyter notebooks for processing ATS outputs and comparing model predictions with observations to build confidence in the model prior to scientific analysis are also included.4) Q_L_relationshipsContains workflows for generating time series of discharge, active network length, and related hydrologic variables used in the stream network expansion–contraction analysis. Includes routines for delineating baseflow-dominated periods. For each catchment, notebooks and processed data (as pickled DataFrames accessed through Pandas Python package) are provided. 5) figure_scriptsProvides the Jupyter notebooks used to generate the figures presented in the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Sub-melt nanosecond pulsed-laser induced densification and strain-field relaxation in single-crystal diamond

Dislocations and polishing-induced defect networks in synthetic diamond introduce local strain fields and broaden Raman features, limiting performance in optical, thermal, and electronic applications. Laser annealing is emerging as a promising approach to repair surface and near-surface defects in diamond without entering the melt regime, yet surface densification, defect-state modification, and associated structural changes have not been well quantified. In this work, we show that sub-melt nanosecond pulsed-laser annealing (PLA) induces near-surface densification and defect-mediated strain relaxation in single-crystal Chemical Vapor Deposition (CVD) diamond. Single- and two-pulse PLA were applied, and structural evolution was quantified using co-registered ISO 25,178 white-light interferometry, depth-resolved Raman spectroscopy, and cross-sectional STEM with geometric phase analysis (GPA). Across a 5 × 6 grid (n = 30), responsive regions exhibit large reductions in local slope (Sdq 45–65%), developed area (Sdr 60–90%), height spread (Sp, Sz 30–65%), void volume (Vv 57–60%), and roughness amplitude (Sa, Sq 48–57%), consistent with densification of ∼4–6.5 nm. Raman profiling shows narrowing of the diamond line and improved spectral uniformity to depths of ∼2–3 μm. Given that the Raman probing depth significantly exceeds the densified layer thickness, this response is interpreted as consistent with long-range strain-field redistribution originating from the near-surface region. STEM-GPA strain maps further support this interpretation, showing smoother strain fields, suppressed hotspots, and redistribution of localized strain concentrations following PLA. These results are consistent with defect-mediated strain relaxation and densification-driven modification of the near-surface energy state. The approach provides a scalable pathway for improving near-surface structural quality in diamond relevant to electronic, photonic, and quantum applications.

Areal surface metrology (ISO 25,178)↗