Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Process modelling”

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 1,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Investigating microstructure evolution in block copolymer membranes

Block copolymer self-assembly in conjunction with nonsolvent-induced phase separation (SNIPS) has been increasingly leveraged to fabricate integral-asymmetric membranes. The large number of formulation and processing parameters associated with SNIPS, however, has prevented the reliable construction of high performance membranes. In this study, we apply dynamical self-consistent field theory to model the SNIPS process and investigate the effect of various parameters on the membrane morphology: solvent selectivity, nonsolvent selectivity, initial film composition, and glass transition composition. We examine how solvent selectivity and concentration of polymers in the film impact the structure of micelles that connect to form the membrane matrix. In particular, we find that preserving the order in the surface layer and forming a connection between the supporting and surface layer are nontrivial and sensitive to each parameter studied. Furthermore, the effect of each parameter is discussed, and suggestions are made for successfully fabricating viable block copolymer membranes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling Noise: Paths toward AI-Enabled Stochastic Earth System Models and Parameterizations

Some of the key challenges in Earth system prediction arise from an uncertain representation of unpredictable natural variability, across time and space scales. This variability, hereafter referred to as noise, represents the quantity against which the strength of a signal of interest is measured to assess predictability. Underestimating or overestimating this noise in Earth system models can lead to issues ranging between overconfidence in an erroneous prediction and a lack thereof in an otherwise accurate one. Noise is part of the physical processes in the Earth system and can amplify or damp a signal of interest in a complex manner that is not systematically characterized. Specifically, while there has been significant progress in our understanding of multiscale interactions among known modes of variability (potentially predictable signals), how these are impacted by the noise is not clear. Recent effort toward stochastic parametrization schemes that provide a representation of noise due to uncertain sub-grid processes in climate models (Berner et al., 2017 and references therein) has shown promise not only in reducing biases and improving probabilistic prediction but also in improving the representation of natural variability in the model toward what is observed. However, today’s standard schemes for stochastic parametrization are still simplistic as they do not take the state of the atmosphere into account (other than by the use of multiplicative noise for model tendencies). Furthermore they require hand-tuning and trial-and-error testing for the magnitude of the added noise, which is difficult in a chaotic system that is as complex as the Earth system.

58 GEOSCIENCES↗
DOE OSTI2022

Finding Electronic Structure Machine Learning Surrogates without Training

A myriad of phenomena in materials science and chemistry rely on quantum-level simulations of the electronic structure in matter. While moving to larger length and time scales has been a pressing issue for decades, such large-scale electronic structure calculations are still challenging despite modern software approaches and advances in high-performance computing. The silver lining in this regard is the use of machine learning to accelerate electronic structure calculations – this line of research has recently gained growing attention. The grand challenge therein is finding a suitable machine-learning model during a process called hyperparameter optimization. This, however, causes a massive computational overhead in addition to that of data generation. We accelerate the construction of machine-learning surrogate models by roughly two orders of magnitude by circumventing excessive training during the hyperparameter optimization phase. We demonstrate our workflow for Kohn-Sham density functional theory, the most popular computational method in materials science and chemistry.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A 1D Model for Nucleation of Ice From Aerosol Particles: An Application to a Mixed‐Phase Arctic Stratus Cloud Layer

Abstract Mixed‐phase clouds (MPCs) have been identified as significant contributors to uncertainties in climate projections, attributable to model representation of processes controlling the formation and loss of supercooled water droplets and ice particles from the atmosphere. Arctic MPCs are commonly widespread and long‐lived, with sustained ice crystal formation processes that challenge current understanding. This study examines the ice‐nucleating particle (INP) reservoir dynamics governing immersion‐mode heterogeneous freezing in an observed case of Arctic MPCs using a simplified 1D aerosol‐cloud model. The model setup includes prescribed dynamical forcings and thermodynamic profiles, and represents INPs as multicomponent and polydisperse particle size distributions. Diagnostic and prognostic approaches to immersion freezing parameterization are compared, including time‐independent (singular) number‐ and surface area‐based descriptions and a time‐dependent description following classical nucleation theory (CNT). The choice of freezing parameterization defines the size of the INP reservoir. The CNT‐based description yields an orders of magnitude larger INP reservoir than the singular parameterizations, which is the dominant factor for sustained ice crystal formation. The efficiency of the freezing process and cloud cooling are of secondary importance. A diagnostic treatment neglecting INP loss is only accurate when the INP reservoir size is large and INP depletion weak. Since a larger INP reservoir sustains ice crystal formation substantially longer, and ice water path scales with ice crystal concentrations for the conditions considered, resolving the source of differences in INP reservoir dynamics due to model implementation is a high priority for advancing climate model physics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Meshfree simulation and prediction of recrystallized grain size in friction stir processed 316L stainless steel

Friction stir processing (FSP) is a promising solid-phase microstructural modification technique that can repair and enhance damaged stainless steel surfaces exposed to harsh environments. The quality of the repaired material is closely correlated to the recrystallized grain size in the stir zone (SZ), which is influenced by the thermomechanical conditions dictated by FSP process parameters. Thus, establishing a reliable relationship between these parameters and recrystallized grain size in the SZ is crucial for optimizing repair quality. However, existing experimental approaches often rely on indirect temperatures measured far from the SZ, along with rough strain rate estimations, which are imprecise and time-consuming. Meanwhile, existing mesh-based modeling methods usually face numerical challenges when dealing with the large material deformations inherent in FSP. Here, to address these issues, this study introduces a meshfree process model for FSP based on the smoothed particle hydrodynamics (SPH) method, aimed at predicting process conditions under different parameters. The model is validated using experimental data from 11 combinations of tool traverse and rotation speeds on 316 L stainless steel. Correlations between process parameters, material flow, temperature, strain, strain rate, and recrystallized grain size are revealed through SPH simulations and electron backscatter diffraction (EBSD) imaging. The results show that in situ SZ temperatures range from 1071 to 1322°C, which exceed the tool temperature by over 300°C. Furthermore, SZ temperature, strain rate, and grain size increase monotonically with higher tool temperature and faster traverse speed. A relationship is then established between the model-predicted Zener-Hollomon parameter and the recrystallized grain size based on EBSD data, expressed as ln(d) = -0.364 ln(Z) + 14.673. Finally, this relationship exhibits satisfactory accuracy with errors of less than 26.9% in predicting grain sizes at various SZ locations, which offers valuable insights for optimizing FSP repair processes for 316 L stainless steel.

316L stainless steel↗
DOE OSTI2024

Feasibility of Formulating Ecosystem Biogeochemical Models From Established Physical Rules

Abstract To improve the predictive capability of ecosystem biogeochemical models (EBMs), we discuss the feasibility of formulating biogeochemical processes using physical rules that have underpinned the many successes in computational physics and chemistry. We argue that the currently popular empirically based approaches, such as multiplicative empirical response functions and the law of the minimum, will not lead to EBM formulations that can be continuously refined to incorporate improved mechanistic understanding and empirical observations of biogeochemical processes. Instead, we propose that EBM parameterizations, as a lossy data compression problem, can be better formulated using established physical rules widely used in computational physics and chemistry, and different biogeochemical processes can be more robustly integrated within a reactive‐transport framework. Through several examples, we demonstrate how mathematical representations derived from physical rules can improve understanding of relevant biogeochemical processes and enable more effective communication between modelers, observationalists, and experimentalists regarding essential questions, such as what measurements are needed to meaningfully inform models and how can models generate new process‐level hypotheses to test in empirical studies. Finally, while empirical models with more parameters are often less robust, physical rules‐based models can be more robust and show lower predictive equifinality, stemming from their enhanced consistency in representations of processes, interactions and spatial scaling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modeling Novel Aqueous Particle and Cloud Chemistry Processes of Biomass Burning Phenols and Their Potential to Form Secondary Organic Aerosols

Phenols emitted from biomass burning contribute significantly to secondary organic aerosol (SOA) formation through the partitioning of semivolatile products formed from gas-phase chemistry and multiphase chemistry in aerosol liquid water and clouds. The aqueous-phase SOA (aqSOA) formed via hydroxyl radical (˙OH), singlet molecular oxygen ( 1 O 2 *), and triplet excited states of organic compounds ( 3 C*), which oxidize dissolved phenols in the aqueous phase, might play a significant role in the evolution of organic aerosol (OA). However, a quantitative and predictive understanding of aqSOA has been challenging. Here, we develop a stand-alone box model to investigate the formation of SOA from gas-phase ˙OH chemistry and aqSOA formed by the dissolution of phenols followed by their aqueous-phase reactions with ˙OH, 1 O 2 *, and 3 C* in cloud droplets and aerosol liquid water. Here, we investigate four phenolic compounds, i.e., phenol, guaiacol, syringol, and guaiacyl acetone (GA), which represent some of the key potential sources of aqSOA from biomass burning in clouds. For the same initial precursor organic gas that dissolves in aerosol/cloud liquid water and subsequently reacts with aqueous phase oxidants, we predict that the aqSOA formation potential (defined as aqSOA formed per unit dissolved organic gas concentration) of these phenols is higher than that of isoprene-epoxydiol (IEPOX), a well-known aqSOA precursor. Cloud droplets can dissolve a broader range of soluble phenols compared to aqueous aerosols, since the liquid water contents of aerosols are orders of magnitude smaller than cloud droplets. Further, our simulations suggest that highly soluble and reactive multifunctional phenols like GA would predominantly undergo cloud chemistry within cloud layers, while gas-phase chemistry is likely to be more important for less soluble phenols. But in the absence of clouds, the condensation of low-volatility products from gas-phase oxidation followed by their reversible partitioning to organic aerosols dominates SOA formation, while the SOA formed through aqueous aerosol chemistry increases with relative humidity (RH), approaching 40% of the sum of gas and aqueous aerosol chemistry at 95% RH for GA. Our model developments of biomass-burning phenols and their aqueous chemistry can be readily implemented in regional and global atmospheric chemistry models to investigate the aqueous aerosol and cloud chemistry of biomass-burning organic gases in the atmosphere.

09 BIOMASS FUELS↗
DOE OSTI2024

A comparative analysis of residual stresses from friction stir processing of aluminum cast 380 and wrought 7075 alloy sheets: experimental characterization and modeling

Residual stresses are often overlooked in friction stir processing (FSP), but their significant impact on fatigue performance necessitates their consideration in optimizing processing parameters. The first step in this effort is understanding how process conditions influence residual stress distributions, especially across different alloys. This study focuses on determining and explaining the through-thickness residual stress variations and the effect of process temperature on the residual stress magnitude in wrought AA7075 and cast AA380.0 alloys. Additionally, for AA380.0, the impact of a second FSP pass was investigated. To achieve this, hole-drilling electronic speckle pattern interferometry (ESPI) and the thermal pseudo-mechanical (TPM) model within finite element analysis were employed to study the 3D distributions of in-plane residual stresses in processed samples under various conditions. A key finding was the varying impact of process temperatures on residual stress magnitudes. Higher process temperatures reduced stresses in AA380.0 but increased them in AA7075. Additionally, the through-thickness stress distributions differed between the two alloys. Further analysis revealed that yield stresses are crucial in explaining these phenomena and the effects of additional FSP passes. Further, this fundamental understanding will be vital in guiding the efforts to mitigate residual stresses and assess their impact on the performance of FSP aluminum alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Sea–Level Rise Impacts on Tidal Marshes and Estuarine Biogeochemical Processes

We used a numerical model to investigate the effects of sea-level rise (SLR) on the biogeochemical processes in the York River Estuary with extensive tidal marshes. The fully-coupled hydrodynamic-water quality-marsh model accounts for the spatial and temporal variations of physical-biogeochemical interactions between the tidal marshes and surrounding waters. This study focuses on an SLR scenario where the vertical accretion of tidal marshes keeps pace with the rising sea levels. Results show that SLR amplifies the tidal range and prolongs flooding duration, which results in enhanced porewater exchanges of materials between the tidal marshes and the surrounding waters. The increased availability of shallow-water habitats and enhanced light utilization in the shallow areas under SLR promote phytoplankton production (PP) in the shallow-water regions of the York River. Consequently, the organic carbon in the open water is fueled by the contributions from shallow waters and the enhanced export of organic carbon from the marshes. The change in the dissolved oxygen (DO) budget in the York River Estuary is attributed to changes in water column respiration, net metabolism of the benthic layer, reaeration, PP, and increased stratification under SLR. The net DO flux out of the York River increases at the York River mouth. Diel DO variation, especially in the marshes in the upper estuary, promotes phosphorus release from the sediment. The changes in dissolved nitrogen under SLR are relatively minimal.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Multi-Fidelity Gaussian Process Regression Method for Probabilistic Wind Farm Power Curve Estimation

Accurate estimation of the power curve for wind turbines or wind farms is crucial to ensure their efficient operation and management. However, conventional methods for power curve estimation rely either on expensive and infrequent measurements or on low-quality numerical simulations. Moreover, the majority of previous studies on power curve estimation for wind turbines or wind farms focused on deterministic estimation, which provides a point estimate of the relationship between wind speed and power generation. Nevertheless, the deterministic approach fails to consider the inherent uncertainty associated with wind energy production resulting from varying turbine characteristics. This can lead to inaccurate power generation estimation and suboptimal decisions regarding energy management. In this paper, a kernel density estimation (KDE) based Multi-Fidelity Gaussian Process Regression (MFGPR) model is proposed to fuse theoretical power curve data and the ground true measurements to create a mapping of wind speed and wind power. By conducting a case study on an actual wind farm in China, the efficacy of the proposed MFGPR model was demonstrated in characterizing the variability of wind power. The probabilistic MFGPR model was also able to generate confidence intervals that encompassed the measured power, thereby improving the accuracy and confidence in wind power estimation or wind resource assessment. Overall, the proposed MFGPR model offers a reliable approach to integrate high-fidelity ground measurements and theoretical power curve data, resulting in precise wind resource assessment and power estimation.

Gaussian process regression↗
DOE OSTI2022

Hydro-mechanical behavior of heated bentonite buffer for geologic disposal of high-level radioactive waste: A bench-scale X-ray computed tomography investigation

The behavior of heated bentonite buffer is critical for the security and long-term performance of a geological repository for high-level radioactive waste (HLW). While laboratory column experiments have been conducted to investigate compacted bentonite and coupled THMC (thermal-hydro-mechanical and chemical) processes for a moderate temperature range of up to 100 °C, data for a higher temperature range are limited. Understanding bentonite behavior and coupled THMC processes under higher temperatures (e.g., up to 200 °C) could allow for a more economic repository design and would expand the data and knowledge base for more reliable modeling. In this study, a bench-scale experiment was conducted in a compacted bentonite column experiencing both heating up to 200 °C in the center and hydration from a sand-clay boundary surrounding the column. During the experiment run for 1.5 years, frequent X-ray computed tomography (CT) scanning of bentonite provided insights into the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. After the experiment, a comprehensive post-dismantling characterization of bentonite samples was conducted. Results showed that the bentonite hydration was axi-symmetrical despite the initial heterogeneity due to packing, confirming the ability of bentonite to seal fast flow/transport paths. Further, compared to a non-heated control experiment, the heated column showed greater CT density variations along the radial distance, indicating that homogenization of bentonite might be more difficult if a temperature gradient is maintained in the repository. Precipitation of an anhydrite layer occurred in the inner hot zone, pointing to potential concerns about salt precipitation causing canister corrosion. Ultimately, the experiments provided a high-resolution window into the strongly dynamic and coupled behavior of bentonite exposed to heating, hydration and swelling, which will be valuable for improving modeling of coupled processes, especially for the early state of a HLW repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Optimisation of the Kaplan hydropower system via PID 2 and digital twin

Here, this paper proposes a proportional–integral-double–derivative (PID 2 ) optimisation method for the Kaplan hydropower system by building a digital twin. The study first uses one multilayer perceptron (MLP) to model the hydroturbine dynamic and then adopts three connected MLPs to model the generator dynamic, both in an open-loop fashion. Inspired by stochastic distribution control (SDC) theory, we regard the training of the turbine's neural network model as a process control problem, and we propose minimising entropy loss to update the network parameters. The next step is to build the digital twin by connecting the neural network models with a PID 2 controller and a lead-lag exciter and run the whole model in a closed-loop fashion. After that, a binary search approach is applied to optimise the PID 2 parameters based on the obtained digital twin model. The simulation results show that the proposed method can reduce the mean square tracking error by more than 90%. Furthermore, the method is extended to jointly optimise the PID 2 controller and excitation system gains through multiobjective optimisation, leveraging Pareto frontier analysis to balance active power and voltage tracking performance. Simulation results confirm the effectiveness of the proposed method, achieving a 83.46% reduction in relative mean square error of active power, a 47.13% reduction in terminal voltage tracking error, and an 82.78% improvement in the overall scalarized objective.

Hydropower system↗
DOE OSTI2023

IM3 GO WEST Parameter Search Dataset

GO WEST is an open-source power grid modeling framework for U.S. Western Interconnection, which allows users to tailor the model depending on their research study and science questions. It is developed to address weather and water dynamics, and associated vulnerabilities in this bulk power system. It covers 28 balancing authorities (BA) and 12 states in U.S. Western Interconnection. GO WEST allows users to select different number of nodes and come up with a simplified network by utilizing 10,000 nodal topology of U.S. Western Interconnection created by Texas A&M University. Users can try and select different number of nodes, mathematical formulations (linear programming vs. mixed-integer linear programming), transmission line limit scaling factors, and hurdle rate scaling factors. GO WEST offers a unit commitment and economic dispatch (UC/ED) module to simulate grid operations on an hourly scale. In this sense, users can calibrate and validate their model versions by comparing model outputs to historical datasets. Therefore, GO WEST can help researchers to strike a balance between model fidelity (i.e. accuracy) and computational complexity (i.e. runtime). This dataset includes model inputs and outputs from 600 model versions for each 2019, 2020, and 2021. The folder naming convention is as follows: Exp{Number of Nodes}_{Mathematical Formulation}_{Transmission Line Limit Scaling Factor in MW}_{Hurdle Rate Scaling Factor in %}_{Year}. Linear programming is designated with "simple" label whereas mixed-integer linear programming is designated with "coal" label. For example, "Exp100_simple_1000_50_2019" folder contains inputs and outputs from 2019 model version with 100 nodes, linear programming, +1000 MW transmission line limit scaling factor, and +50% hurdle rate scaling factor. GO WEST GitHub repository hosts all raw datasets, processing scripts, and model scripts. Please refer to the README file for a detailed description of the included files.

Economics↗
DOE OSTI2025

HYBRD (High Resolution HYBrid Regional Downscaling) Model: Input data and Code

The HYBRD (HYBrid Regional Downscaling) model is a high-resolution urban land downscaling model that can be used to downscale intermediate urban land use and land cover (LULC) products into a high-resolution (30-meters). HYBRD uses a sequential hybrid process, combining statistical models with cellular-automata-based spatial algorithms. This repository contains all the necessary model code and inputs needed to successfully run HYBRD for Los Angeles, California. The repo also contains example outputs of each model step, except the final simulated raster outputs. Examples of simulated raster outputs for multiple scenarios for Los Angeles are available at DOI: 10.57931/2575233. Please refer to Related Works below.

Land↗
DOE OSTI2022

Deep Learning of Dark Energy Spectroscopic Instrument Mock Spectra to Find Damped Lyα Systems

We have updated and applied a convolutional neural network (CNN) machine-learning model to discover and characterize damped Ly α systems (DLAs) based on Dark Energy Spectroscopic Instrument (DESI) mock spectra. We have optimized the training process and constructed a CNN model that yields a DLA classification accuracy above 99% for spectra that have signal-to-noise ratios (S/N) above 5 per pixel. The classification accuracy is the rate of correct classifications. This accuracy remains above 97% for lower S/N ≈1 spectra. This CNN model provides estimations for redshift and H i column density with standard deviations of 0.002 and 0.17 dex for spectra with S/N above 3 pixel -1 . Also, this DLA finder is able to identify overlapping DLAs and sub-DLAs. Further, the impact of different DLA catalogs on the measurement of baryon acoustic oscillations (BAO) is investigated. The cosmological fitting parameter result for BAO has less than 0.61% difference compared to analysis of the mock results with perfect knowledge of DLAs. This difference is lower than the statistical error for the first year estimated from the mock spectra: above 1.7%. We also compared the performances of the CNN and Gaussian Process (GP) models. Our improved CNN model has moderately 14% higher purity and 7% higher completeness than an older version of the GP code, for S/N > 3. Both codes provide good DLA redshift estimates, but the GP produces a better column density estimate by 24% less standard deviation. A credible DLA catalog for the DESI main survey can be provided by combining these two algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

DUET - Distribution of Understory using Elliptical Transport: A mechanistic model of leaf litter and herbaceous spatial distribution based on tree canopy structure

Heterogeneity in surface fuels produced by overstory trees and understory vegetation is a major driver of fire behavior and ecosystem dynamics. Previous attempts at predicting tree leaf and needle litter accumulation over time have been constrained in scope to probabilistic models that consider a limited number of key factors influencing tree litter dispersal patterns and decomposition processes. Here we present a mechanistic model for estimating variation in surface fuels called the Distribution of Understory using Elliptical Transport (DUET). DUET uses a pre-generated voxelated canopy array and distributes the leaf litter from the trees in an elliptical shape based on tree species characteristics, wind data, and location-specific features. DUET then calculates grass growth between trees and around the litter patterns and sets moisture levels based on relative humidity. In this study, we use DUET to simulate spatial distributions of leaf litter and illustrate model sensitivity to several parameters, providing inputs for testing fuel spatial heterogeneity impacts on fire behavior using process-based models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Utilizing waste heat in wastewater treatment plants for water desalination: Modeling and Multi-Objective optimization of a Multi-Effect desalination system using Decision Tree Regression and Pelican optimization algorithm

This paper examines the feasibility of using waste heat from wastewater treatment plants (WWTPs) for water desalination. A model was developed to utilize waste heat from the gensets at As Samra WWTP in Jordan, using real data and TRNSYS® software to calculate available waste heat. The desalination process was then modeled with ASPEN PLUS® software, focusing on multi-effect desalination (MED). Both series and parallel configurations for the MED system were compared. The study investigated the effects of system feeding flow rate, feeding pressure, and heat input on productivity, performance ratio, and recovery ratio. The study also introduces a novel optimization technique combining machine learning and modern optimization algorithms to maximize system productivity and performance. Initially, a decision tree regression (DTR) model is developed to establish relationships between key independent variables (flow rate, feed pressure, and heat input) and dependent variables (productivity, performance ratio, and recovery ratio). The Pelican Optimization Algorithm (POA) is then used to identify the optimal values of the independent variables for maximum productivity and performance. The results show that using a series configuration yields a system productivity of 3984.2 kg/hr, a performance ratio of 3.78, and a recovery ratio of 0.991 at a feed flow rate of 4000 kg/hr, feed pressure of 3 bars, and heat input of 719 kW. Optimal productivity (4421 kg/hr), performance ratio (3.81), and recovery ratio (0.851) are achieved at a feed flow rate of 5166 kg/hr, feed pressure of 3.2 bars, and heat input of 794 kW. In conclusion, the techno-economic assessment indicates a levelized cost of water of 1.63 USD/m 3 for parallel configurations and 1.65 USD/m 3 for series configurations, with a payback period of less than two years.

42 ENGINEERING↗
DOE OSTI2024

Geometry-complete diffusion for 3D molecule generation and optimization

Abstract Generative deep learning methods have recently been proposed for generating 3D molecules using equivariant graph neural networks (GNNs) within a denoising diffusion framework. However, such methods are unable to learn important geometric properties of 3D molecules, as they adopt molecule-agnostic and non-geometric GNNs as their 3D graph denoising networks, which notably hinders their ability to generate valid large 3D molecules. In this work, we address these gaps by introducing the Geometry-Complete Diffusion Model (GCDM) for 3D molecule generation, which outperforms existing 3D molecular diffusion models by significant margins across conditional and unconditional settings for the QM9 dataset and the larger GEOM-Drugs dataset, respectively. Importantly, we demonstrate that GCDM’s generative denoising process enables the model to generate a significant proportion of valid and energetically-stable large molecules at the scale of GEOM-Drugs, whereas previous methods fail to do so with the features they learn. Additionally, we show that extensions of GCDM can not only effectively design 3D molecules for specific protein pockets but can be repurposed to consistently optimize the geometry and chemical composition of existing 3D molecules for molecular stability and property specificity, demonstrating new versatility of molecular diffusion models. Code and data are freely available on GitHub .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗