Search NASA⌕ Search

SEARCH · Search NASA

Results for “United States”

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.

1,328 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Quantum Filtering and Analysis of Multiplicities in Eigenvalue Spectra

Fine-grained spectral properties of quantum Hamiltonians, including both eigenvalues and their multiplicities, provide useful information for characterizing many-body quantum systems as well as for understanding phenomena such as topological order. Extracting such information with small additive error is #BQP-complete in the worst case. In this work, we introduce QFAMES (quantum filtering and analysis of multiplicities in eigenvalue spectra), a quantum algorithm that efficiently identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions, which allows us to bypass worst-case complexity barriers. QFAMES also enables the estimation of observable expectation values within targeted energy clusters, providing a powerful tool for studying quantum phase transitions and other physical properties. We validate the effectiveness of QFAMES through numerical demonstrations, including its applications to characterizing quantum phases in the transverse-field Ising model and estimating the ground-state degeneracy of a topologically ordered phase in the two-dimensional toric code model. We also generalize QFAMES to the setting of mixed initial states. Our approach offers rigorous theoretical guarantees and significant advantages over existing subspace-based quantum spectral analysis methods, particularly in terms of the sample complexity and the ability to resolve degeneracies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Cave-based chronology of valley incision in the glaciated Northern Carpathians

While glaciers leave a mark on mountain landscapes, the importance of glacial erosion in mountains at a global scale remains insufficiently constrained despite over a century of research aimed at quantifying the rates and magnitudes of such erosion. Measuring rates of landscape change across the Quaternary in repeatedly glaciated mountain landscapes requires the ability to accurately partition the relative role of glacial versus fluvial processes in forming relief; however, the modification of landforms by surface erosion hinders quantification of process rates over the last few million years. This study presents 32 paired 26 Al and 10 Be cosmogenic nuclide burial ages from allochthonous, quartz-bearing sediment deposited within a paleophreatic karst system in the northern Carpathians. These data quantify Pliocene-to-present timing of valley incision and karst system reactivation, as well as erosion rates over time. Post-burial, muon production-correction yields 24 fully resolved ages, six minimum ages, and two unresolvable ages. Burial ages are interpreted to document the tempo of valley incision via phreatic passage formation and reveal that ∼80% of current valley relief (∼700 m) had developed by approximately 1.6 Ma. Paleoerosion rates derived from this dataset suggest little change in the processes generating sediment washed into the karst networks since the Pliocene. These results suggest that river incision dominated relief creation in Carpathian valleys that bear the imprint of glaciation. Our findings suggest the importance of investigating other mid-latitude mountainous regions worldwide to determine if they have similar erosion and topographic histories.

Carpathians↗
DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Hybrid Bio-Based Composites: Enabling Cellulose Nanofiber (CNF) Incorporation into Composites via Macroscale Natural Fiber Carriers

Cellulose nanofibers (CNFs) have significant potential in composites as additives to improve mechanical properties, melt rheology, and more. However, agglomeration of CNFs is a key challenge in composite melt processing as obtaining nano-level dispersion of CNFs often requires cost- and energy-intensive processes (e.g., solvent exchange or freeze drying) due to the strong hornification tendencies of CNF. Herein, we avoid these challenges by using a natural fiber carrier method to integrate CNF into thermoplastic composites. Fibers are co-dried to create a hybrid fiber feedstock for compounding in which natural fibers are decorated with dispersed nanofibers. The hybridized fibers result in up to a 24% increase in tensile strength and up to a 35% increase in Young’s modulus compared to composites only containing natural fibers. The lignocellulosic nanofibers are found to outperform their purely cellulosic counterpart, which is theorized to be due to either an increased propensity for fibrillation of the lignocellulosic fibers or the increased hydrophobicity of the fibers due to the presence of lignin. Surface analysis of fiber feedstocks, via streaming potential measurements and dynamic light scattering (DLS), confirmed a significant change in the feedstock hydrophobicity before and after hybridization. While mild additions of CNF (1 wt.% on the macroscale fiber) do not impact the composite melt viscosity, the viscosity is found to increase at higher CNF loadings (5 wt.% on the macroscale fiber), indicating its utility as a rheology modifier. Lastly, use of these materials as novel feedstocks for medium-scale additive manufacturing in high-fidelity part production was demonstrated.

bio-based↗
DOE OSTI2026

Integrating AI Data Centers with the Power Grid

The rapid expansion of artificial intelligence (AI) has triggered an unprecedented surge in electricity demand, with US data center energy use projected to double or triple 2023 levels by 2028. This exponential growth places strain on grid infrastructure, which can hinder timely construction of desired computing capacity. To bridge this supply-demand gap, utilities and AI developers are increasingly turning to demand flexibility, a strategy that incentivizes shifting or reducing power use during peak periods of grid stress. Data centers are uniquely equipped for flexible operations due to their digital workloads, built-in redundancy, and onsite energy assets. This article outlines four primary mechanisms to enable data center flexibility: computational load flexibility (shifting tasks temporally or geographically), flexible use of core facility infrastructure adjustments, energy storage utilization, and onsite electricity generation. To encourage adoption, utilities are deploying new tariff designs, including voluntary interruptible service riders, mandated flexibility requirements, and streamlined interconnection processes for flexible loads. For the highly capitalized and rapidly growing AI industry, the primary motivators for embracing these strategies are expediting facility interconnection, satisfying emerging regulatory mandates, and mitigating community resistance. While demand flexibility cannot substitute the long-term need for new bulk power generation, it serves as an essential, immediate solution for enabling near-term deployment. By transforming data centers from grid stressors into stabilizing assets, flexible operations can ensure reliable grid integration, ease market pressures, and support a resilient power system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Mechanically and Thermally Robust Gel Electrolytes Built from A Charged Double Helical Polymer

Polymer electrolytes have received tremendous interest in the development of solid-state batteries, but often fall short in one or more key properties required for practical applications. Herein, a rigid gel polymer electrolyte prepared by immobilizing a liquid mixture of a lithium salt and poly(ethylene glycol) dimethyl ether with only 8 wt% poly(2,2′-disulfonyl-4,4′-benzidine terephthalamide) (PBDT) is reported. The high charge density and rigid double helical structure of PBDT lead to formation of a nanofibrillar structure that endows this electrolyte with stronger mechanical properties, wider temperature window, and higher battery rate capability compared to all other poly(ethylene oxide) (PEO)-based electrolytes. The ion transport mechanism in this rigid polymer electrolyte is systematically studied using multiple complementary techniques. Li/LiFePO 4 cells show excellent capacity retention over long-term cycling, with thermal cycling reversibility between ambient temperature and elevated temperatures, demonstrating compelling potential for solid-state batteries targeting fast charging at high temperatures and slower discharging at ambient temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Influence of Powder Characteristics and Processing Methods on Creep Performance of Powder Metallurgy Hot Isostatic Pressed SS316

The U.S. nuclear energy expansion goals are driving the demand for manufacturing routes that can rapidly produce large, complex, near net shape components. Powder metallurgy hot isostatic pressing (PM HIP) is an advanced manufacturing technique that can be economically scaled-up, while alleviating the supply chain challenges that forging and casting face in terms of cost and lead time constraints. This makes PM-HIP a viable technology to aid and accelerate large-scale part manufacturing for nuclear applications. However, large-scale qualification and deployment of this technology require a thorough understanding of the influence of powder feedstock quality, powder handling history, hot isostatic pressing (HIP) parameters, and subsequent heat treatment on microstructural evolution and elevated temperature mechanical performance. The present work focusses on 316 austenitic stainless steel (SS316) which is most commonly used in high temperature environments for nuclear applications Results from this study show that PM HIPed SS316 meets ASME tensile requirements at room temperature and at elevated temperature. However, creep performance of PM-HIPed SS316 remains inferior to its wrought counterpart, demonstrating that tensile performance alone is not a reliable metric for long duration high temperature integrity. Further, this report delineates powder derived microstructural features that govern creep damage, with key evidences pointing to “microstructural inheritance” from gas atomized powder feedstocks. Commercial SS316 powders of varying chemical compositions and recycling histories were studies, and the results showed large differences in elemental segregation, oxide surface layers and secondary phase distributions. Multi-scale characterization revealed segregation of chromium, molybdenum, manganese and silicon at the boundaries and the precipitation of manganese-, silicon-, and molybdenum-oxides. During HIP consolidation, these surface oxides transform into decorated prior particle boundaries (PPBs) and grain boundary inclusions that persist through conventional post-HIP solution annealing treatment. The retained oxides in post-HIP microstructures were found to influence grain growth, precipitation behavior, and ultimately creep cavitation and fracture. Such post-HIP heat treatments are therefore limited by a complex trade-off between grain growth, and oxide coarsening which aggravate creep damage by acting as nucleation sites for cavities. The objective of this work is to establish an integrated processing–structure–property framework for PM-HIP 316 stainless steel by investigating the influence of powder feedstock characteristics in pre- and post-HIP processing as well as to understand the significance of post-HIP heat treatment on microstructural evolution and creep properties. The results from this report emphasize the significance of powder feedstock integrity in improving creep performance of PM-HIPed SS316, by highlight the effect of rapid solidification, elemental segregation, oxide formation and powder recycling on microstructural defect inheritance following HIP consolidation. Rather than considering HIP processing, solution annealing, and mechanical performance independently, this report treats powder production, HIP consolidation, post-HIP thermal processing, and creep deformation as interconnected stages within a continuous metallurgical process. The resulting framework will provide a scientific basis for developing feedstock engineering strategies capable of improving long-term reliability of PM-HIP stainless steels and accelerating their qualification for advanced nuclear applications.

Ajjarapu, Pavan [Oak Ridge National Laboratory (OR↗
DOE OSTI2026

Spectrograph stabilization using a single-delay interferometer on the Hale Telescope

We describe a technique for spectrograph stabilization useful when conventional mitigation techniques of vacuum tanks, thermal insulation, and laser frequency comb may be impractical, expensive, heavy, or bulky. This includes spectrographs on airborne platforms or mounted on telescopes where they suffer a changing gravity vector or other drifts. Placing a fixed-delay interferometer in series with a spectrograph forms an externally dispersed interferometer (EDI). This produces a uniform sinusoidal comb multiplying input spectrum, creating (through heterodyning) beats (moiré patterns). In Fourier space for low frequencies up to the comb frequency, the moiré generated signal counter-rotates to ordinary spectra under an unknown disperser wavenumber drift Δx. This generates a large negative feedback signal useful in a conceptual control loop, to converge rapidly to a stable spectrum and yield Δx. A modified EDI data analysis algorithm (“crossfading”) combines frequency-weighted moiré with conventional spectrum to cancel net output spectrum reaction to Δx. Needing only a single-delay, this is a practical improvement over prior crossfading analyses requiring multiple delays. We test crossfading on ThAr data near 4850 cm−1 taken on Hale telescope in an earlier project. In a single pass, we reduce drift 20 times. Using seven iterations, we reduce 0.5 cm−1 (31 km/s Doppler equivalent) drift to 4×10−7 cm−1 (2.5 cm/s). The interferometer delay can wander, because linearity of phase versus wavenumber interpolates science features between bracketing calibrating spectral references. Second, mathematically reversing the heterodyning effect doubles effective spectral resolution without changing disperser slit.

Erskine, David J [Lawrence Livermore National Labo↗
DOE OSTI2026

In situ electric field X-ray total scattering reveals composition-dependent electromechanical strain mechanisms in (1 − x )BiFe 2/8 Ti 3/8 Mg 3/8 O 3 – x PbTiO 3 ceramics

Ferroelectric materials find many applications in energy and aerospace industries. A major challenge for ferroelectric materials is maintaining their properties at elevated temperatures. The (1 − x )BiFe 2/8 Ti 3/8 Mg 3/8 O 3 – x PbTiO 3 (BFTM–xPT) ferroelectric solid solution is a promising candidate for high-temperature applications as it has a high piezoelectric response while also maintaining its ferroelectric phase until a high Curie temperature (T C ) of 650 °C. In this work, we investigate the piezoelectric mechanism in BFTM–xPT, a novel high-T C ferroelectric ceramic. To elucidate the origin of its enhanced piezoelectric performance, relative to other piezoceramics of similar T C , in situ electric field X-ray total scattering experiments were performed. Total scattering combines Bragg scattering (diffraction) and diffuse scattering (pair distribution function), providing insight on various length scales. With information spanning different length scales, extrinsic contributions (i.e., domain wall motion and interphase boundary motion) can be separated from intrinsic contributions (i.e., piezoelectric lattice strain). We show that, for compositions within the morphotropic phase boundary (MPB, x = 0.325), the piezoelectric response is dominated by intrinsic piezoelectric lattice strain, whereas outside the MPB (x = 0.375) the major component of the piezoelectric response is from extrinsic mechanisms including an electric-field-induced phase transition and tetragonal domain wall motion. This article reveals that the origin of piezoelectric response in BFTM–xPT is composition-dependent and shows that high-performance materials may also be located outside the MPB. These findings can help guide material design to optimize properties for specific applications.

Richtik, Brooke N. [University of Calgary, AB (Can↗
DOE OSTI2025

Mechanically Accelerated Depolymerization of Entangled Linear Polymer Melts

Mechanical forces can enhance the chemical depolymerization of synthetic polymers when shear flow accelerates chain scission. To quantify the extent of mechanically-accelerated scission, the effect of simple shear flow (duration and strength) with low Weissenberg and Deborah numbers was investigated by considering the impact of applied work in both simple shear and shear dominated mixed flows. Hydrogenated polyisoprene was chosen as a model linear, entangled system. The conditions (strain amplitude, frequency, and shearing time) necessary to increase chain scission were assessed in the rubbery melt. Shear flow accelerated chain scission at higher temperatures, suggesting an activated process. Isothermal scission versus work curves were superposed by applying shift factors a T,S , whose Arrhenius-like temperature dependence gave an apparent activation energy for chain scission of ~ 110 kJ/mol, which is likely a combination of the activation energy of viscosity and bond energy. This work provides a base for quantifying the impact of shear on depolymerization of polymer melts and highlight the connection between viscous dissipation and scission chemistry.

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

Enzymatic Nylon Deconstruction: Enzyme Discovery, Engineering, and Opportunities

Nylons are widely used synthetic polyamides valued for their strength, versatility, and durability across diverse applications. However, their petrochemical origin and energy-intensive production underscore the need for efficient, circular solutions. Conventional recycling methods remain limited by incomplete recovery, material degradation, and costly sorting requirements. Enzymatic depolymerization offers a selective, low-energy alternative capable of processing mixed waste streams under mild conditions. While significant progress has been achieved for polyesters, enzymatic degradation of polyamides is still at an early stage. The discovery of nylon hydrolases demonstrated the potential of biological systems to evolve catalysts for synthetic polyamides, yet reported depolymerization yields remain low. These limitations reflect both the structural complexity of nylons and the need for improved enzyme discovery and engineering. In conclusion, this review highlights recent advances, key challenges, and future directions for enzymatic nylon recycling, outlining its potential role enabling mixed polymer waste to be used as a green feedstock for remanufacturing.

Amides↗
DOE OSTI2026

Classical-Quantum Algorithm for Solving Stochastic Programs

Stochastic programming provides a rigorous mathematical framework for making decisions under uncertainty in a risk-aware manner. Two-stage stochastic programming is, perhaps, the simplest form of this framework. Here the first-stage variables represent decisions that must be made "here and now" in the face of uncertainty, while the second-stage variables are decisions made after uncertain events. However, the broad adoption of stochastic programming has been hindered by computational challenges caused by the two-stage stochastic programming formulation which requires solving an ensemble of optimization problems. Using quantum amplitude estimation (QAE), quantum computers have shown the theoretic ability to compute expectations with Monte-Carlo methods with quadratically fewer samples than classical methods. In this work, we present a quantum algorithm for computing the expectation term using QAE for given first-stage decisions. Further, we detail methods of computing gradient information from the quantum calculation enabling the application of classical gradient-based optimization techniques. The result is a classical-quantum hybrid method of solving two-stage stochastic programs. These techniques are demonstrated with computational experiments based an engineering optimization problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Spin‐Polarized Interfaces Redirect CO 2 Reduction From CO to Formate

Achieving high product selectivity in electrocatalytic carbon dioxide reduction (CO 2 RR) remains a critical challenge due to competition between multiple proton-coupled electron-transfer pathways on catalyst surfaces. Meanwhile, chirality-induced spin selectivity (CISS), which enables spin-polarized electron transport through chiral interfaces, has recently emerged as a promising strategy to modulate interfacial electrochemical reactions. Although the CISS effect has been shown to enhance selectivity and efficiency in the spin-sensitive oxygen evolution reaction (OER), its role in regulating CO 2 RR pathways and in stabilizing intermediates remains largely unexplored. Here, chiral molecules (R- and S-1,1′-bi-2-naphthyl-2,2′-diyl hydrogen phosphate, BNP) were integrated with SnO 2 to construct chiral-modified catalysts (R-BNP/SnO 2 and S-BNP/SnO 2 ). Compared with bare SnO 2 and racemic BNP-modified SnO 2 (Rac-BNP/SnO 2 ), the chiral catalysts exhibited a pronounced shift in product selectivity from CO toward formate production. Importantly, in-situ attenuated total reflectance surface-enhanced infrared absorption spectroscopy (ATR-SEIRAS) reveals that the chiral interface selectively stabilizes the O-bound *OCHO intermediate associated with the formate pathway and modulates interfacial water structure and hydrogen-bonding dynamics. These findings demonstrate that spin-polarized interfacial electron transfer can regulate CO 2 RR pathway selectivity by modulating the stabilization of key intermediates. More broadly, this work establishes chiral spin-selective interfaces as a new strategy for regulating competitive electrocatalytic reaction pathways.

14 SOLAR ENERGY↗
DOE OSTI2026

Revealing the Hidden Third Dimension of Point Defects in Two-Dimensional MXenes

Point defects govern many important functional properties of two-dimensional (2D) materials. However, resolving the three-dimensional (3D) arrangement of these defects in multi-layer 2D materials remains a fundamental challenge, hindering rational defect engineering. Here, we overcome this limitation using an artificial intelligence-guided electron microscopy workflow to map the 3D topology and clustering of atomic vacancies in Ti3C2TX MXene. Our approach reconstructs the 3D coordinates of vacancies across hundreds of thousands of lattice sites, generating robust statistical insight into their distribution that can be correlated with specific synthesis pathways. This large-scale data enables us to classify a hierarchy of defect structures-from isolated vacancies to nanopores-revealing their preferred formation and interaction mechanisms, as corroborated by molecular dynamics simulations. This work provides a generalizable framework for understanding and ultimately controlling point defects across large volumes, paving the way for the rational design of defect-engineered functional 2D materials.

2D materials↗
DOE OSTI2026

Development of in-situ polymerized intrinsically conductive resin and low-cost carbon pigments offering high conductivity for sensing, EMI shielding and lighting protection

Electrically conductive composites are emerging across diverse industries such as electronic, automotive, aerospace, advanced air mobility, biomedical, infrastructure, defense and security offering static charge dissipation, electromagnetic interference shielding, lighting protection, sensing, dicing, corrosion monitoring, etc. Conductivity enhanced composites provide several advantages compared to conventional metals including weight reduction, corrosion resistance, energy efficient processability, tunable properties and multifunctionality. Polymers are typically insulating in nature and require conducting filler for electron transport. However, dispersion and polymer-filler interphases are critical and often disrupt conducting pathways. Besides, conductive fillers such as graphene, carbon nanotube, MXene, silver nanowire, etc. are expensive, limiting their wide adoption in composite industry. On the other hand, a limited number of intrinsically conductive polymers are available among which polyaniline (PANI) has been widely studied due to its high conductivity, thermal and chemical stability. However, PANI is difficult to process and exhibits weak mechanical properties. In brief, there is a significant demand for electrically conductive polymer formulation with cost-effective conducting fillers that offer processability in scale to expand the market of conductivity enhanced materials.

Kumar, Vipin [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2026

NanoPSD: A software for automatic detection of Nano-Particle Shape Distribution in electron microscopy images

Accurate quantification of the size and morphology of nanoparticles from electron microscopy (EM) images is essential to understand growth mechanisms, surface reactivity, and functional behavior in nanoscale materials. Manual analysis remains slow, subjective, and difficult to reproduce in large datasets. We introduce NanoPSD (Nano-Particle Shape Distribution), an open-source and fully automated framework for quantitative particle detection and morphology analysis from EM images. NanoPSD integrates adaptive contrast enhancement, polarity-agnostic scale-bar detection, Optical Character Recognition (OCR)-based calibration, and classical segmentation via Otsu thresholding with morphological refinement. Particle contours are used to extract geometric descriptors, including equivalent circular diameter, aspect ratio, circularity, and solidity, enabling automated classification into spherical, rod-like, and aggregate morphologies. The framework supports both single-image and batch processing, generating publication-quality visualizations, LaTeX-ready tables, and structured comma-separated values (CSV) datasets. As a demonstration, we applied NanoPSD to plasma-synthesized nanoparticle samples diagnosed via transmission electron microscopy (TEM). The code produced statistically robust size and morphology distributions spanning a few to tens of nanometers with minimal user supervision. The pipeline demonstrates high reproducibility and scalability, processing large image collections with consistent calibration and output formatting. Its modular design enables seamless integration of future deep-learning-based segmentation models, providing a pathway toward intelligent, data-driven electron microscopy analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Heliostat sizing methodology for concentrating solar thermal industrial process heat projects

This study presents a method to obtain a heliostat size that minimizes the levelized cost of heat (LCOH) of a heliostat-based concentrating solar thermal system for applications of solar heating for industrial processes at operating temperatures from 565 to 1550°C. The method extends prior work by embedding a routine for system design that obtains near-optimal subsystem sizes, increasing the fidelity of drive cost functions, and adding an optical performance model to supplement the previously developed cost models, which we update to reflect current pricing trends. An illustrative business case is developed for Daggett, California, targeting specified annual thermal energy outputs of 50 to 400 GWh th . Optical performance is modeled using verified estimates from the literature. A surrogate heliostat cost model, derived from commercial heliostat designs and scaled for production volume, installation, and operations and maintenance costs, is used to develop cost functions. Results show that heliostat size strongly affects the LCOH, producing a characteristic U-shaped trend with a robust near-optimal window of 7-20 m 2 ; the heliostat size producing the lowest project cost in our study grows slightly as the project size increases, and is reduced as the operating temperature increases. The findings in this study are consistent with the general trend of smaller heliostats being deployed at existing projects for high-temperature industrial process heat and reflect the significant reduction in power electronics and other per-heliostat costs. The methodology we propose is general and can be tailored to revised cost curves as the technology continues to evolve.

14 SOLAR ENERGY↗
DOE OSTI2026

EvoDiffMol: evolutionary diffusion framework for 3D molecular design with optimized properties

Designing molecules with specific target properties remains a fundamental challenge in computational chemistry. While existing approaches show promise, most rely on simplified representations like SMILES strings or 2D graphs that lack essential three-dimensional geometric information. We present EvoDiffMol, a computational framework that integrates evolutionary algorithms with three-dimensional diffusion models for property-driven molecular generation. The method operates through adaptive evolutionary optimization, where population-based selection guides the generation process toward desired property landscapes. EvoDiffMol supports both unconstrained molecular design and scaffold-constrained generation that preserves fixed substructures while optimizing complementary regions. Comprehensive evaluation demonstrates exceptional performance, achieving the highest drug-likeness score (0.94) among all compared state-of-the-art methods while maintaining excellent validity, uniqueness, and novelty. Beyond single property optimization, the framework demonstrates flexible multi-property optimization capabilities, simultaneously controlling multiple molecular descriptors including synthetic accessibility, lipophilicity, topological polar surface area, and clinically relevant ADMET properties such as cardiotoxicity (hERG) and intestinal permeability (Caco-2). This adaptability spans from simple descriptors to practical pharmaceutical endpoints without requiring complete model retraining. The framework achieves precise control over target property values, generating molecules with properties closely matching specified targets for both single and multiple descriptors. Scaffold-constrained experiments preserve fixed molecular cores while maintaining effective property optimization. The three-dimensional representation offers advantages in maintaining structural validity during iterative optimization, with potential for geometry-aware applications in materials science and drug discovery.

3D molecular generation↗