Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geometry optimization”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The FMO-DFTB Method

Although the fragment molecular orbital (FMO) method enables electronic structure calculations with near-linear scaling behavior with respect to system size, the computational cost of ab initio methods typically employed in conjunction with FMO is still prohibitive for routine calculations of very large systems or long timescale molecular dynamics simulations. We, therefore, combined the FMO and density-functional tight-binding (DFTB) method, which is one of the emerging semi-empirical quantum chemical methods, and have demonstrated that FMO-DFTB is capable of performing geometry optimizations for systems containing up to one million atoms using limited computational resources. In this chapter, we will review the basics of the DFTB method first before introducing FMO-DFTB, focusing on the relationship with density functional theory and other FMO methodologies. We also demonstrate the latest scalings of FMO-DFTB2 and DFTB3 using three-dimensional water clusters, showing that the most favorable scaling is O(N1.16). Applications of FMO-DFTB to various systems are briefly summarized, and an outlook to future applications is provided.

Nishimoto, Yoshio↗
DOE OSTI2022

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. So, we describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Agentic framework for programmatic crystal structure generation using a fine-tuned worker–supervisor large language model

Platinum group metals (PGMs) underpin many catalytic technologies but face severe supply constraints, motivating the search for alternative materials and computational methods to accelerate discovery. While atomistic simulation tools such as Pymatgen and ASE have streamlined structure manipulation, they require detailed inputs, limiting accessibility for experimentalists and slowing early-stage exploration. Here, in this study, we present an AI-driven agentic framework that orchestrates worker–supervisor large language models (LLMs). The worker translates natural-language prompts of varying abstraction into valid crystallographic structures using a compact LLM fine-tuned with low-rank adaptation on a curated text–code–CIF dataset, emphasizing energy-efficient training. Benchmarking against the baseline CodeGen-350M-mono model shows that fine-tuning reduces hallucination rates from 100% to as low as 5% and improves structural match accuracy to up to 82% for fully specified inputs. Accuracy declines with decreasing prompt detail but remains nontrivial even when only stoichiometry and space group are provided, underscoring the LLM’s capacity for crystallographic inference. The supervisor Claude LLM evaluates the outputs and triggers iterative refinement through the worker’s built-in structure manipulation capabilities (e.g., supercell scaling, strain, vacancy, and substitution operations). We further demonstrate use cases for technologically relevant catalysts, including IrO 2 , pyrochlore Pb 2 Ir 2 O 7 , Ni 2 FeO 4 , and Ni 3 Mo, where the framework generates physically consistent structures that can be refined via geometry optimization. This work introduces a low-energy, language-driven pathway for integrating human and machine intelligence in materials design, paving the way for AI-assisted synthesis planning and high-throughput screening of complex oxides.

AI agent↗
DOE OSTI2022

A three-dimensional pore-scale model for redox flow battery electrode design analysis

A three-dimensional (3-D) pore-scale model has been developed to construct multiscale fibrous electrodes for redox flow batteries (RFB). New designs, such as biporous electrodes modify electrode structures by creating secondary pores on single carbon fiber to reduce internal battery resistance. Existing pore-scale models only resolve electrodes to the single carbon fiber scale and cannot incorporate recent multiscale electrode designs into numerical models. Our new model aims to bridge this gap and provide a tool to rapidly screen new electrode configurations. Two multiscale electrodes, laser-perforated and biporous electrodes, were investigated at varying operation conditions with the proposed 3-D pore-scale models. The laser-perforated electrode exhibits a reduced pressure drop, but follows the same permeability correlation compared to the corresponding pristine electrode. For the biporous electrode, the added specific surface area and faster reaction kinetics from the secondary pores are the most influential factors leading to improved battery efficiency. However, operating the biporous electrode in limiting current density conditions should be avoided due to the decreased mass transfer efficiency and a more significant voltage loss. Finally, we believe that our 3-D pore-scale model can accelerate the flow battery electrode design process and provide new insights into electrode geometry optimizations.

25 ENERGY STORAGE↗
DOE OSTI2023

Development of an accelerator-based neutron source to prototype Mo-99 production, part I: A liquid LBE windowless target

In this article, Molybdenum-99 (Mo-99)’s decay product, technetium-99 (Tc-99 m), is one of the most critical isotopes for medical diagnostics. To provide U.S. domestic supply of Mo-99 without using high-enriched uranium (HEU), a subcritical uranium target assembly (UTA) is irradiated by an accelerator-based neutron source to create Mo-99 through fission. This study discusses the development of the accelerator-based neutron source. The high-energy electrons from the accelerator irradiate a liquid lead-bismuth eutectic (LBE) target to produce neutrons. Part I of this work focuses on numerical and experimental analysis towards the development of a liquid LBE windowless target. Unlike the existing windowless targets in literature, the current design creates a vertical free surface for a beam to irradiate. First, a hydrodynamic analysis of the LBE windowless target is performed. Simplified analytical calculations are assisted by 2D computational fluid dynamics (CFD) simulations to design the target, with the focus on eliminating recirculation zones and avoiding cavitation. With the optimized geometry, the experimental study is performed to investigate the flow hydrodynamics using liquid LBE. The experiments (1) compare pressure drop in the system to correlation predictions; (2) visualize the free surface liquid LBE flow from the beam view; (3) validate the LBE flow profile using temperature sensitive paint from the side view; and (4) validate the liquid LBE film thickness using gamma densitometer measurements. Second, the power handling capability of the designed windowless target is investigated. The divider plate in the current design is susceptible to overheating due to the thin LBE film in front. As LBE erosion and corrosion is likely to occur at an LBE velocity of 2.0 m/s and temperature above 500 °C, a power limit of 10 kW of beam power was established to prevent this corrosion from occurring, which is calculated by a Nusselt number correlation. The divider plate surface temperature at 10 kW agrees well with the 3D CFD simulation results. Part I demonstrates the fundamental physics in liquid LBE windowless target design and associated testing. A companion paper, Part II will demonstrate how to couple this windowless target into the Mo-99 production system, including an accelerator system operating under an ultra-high vacuum and the UTA cooled by water at room temperature.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Photon counting with intensified charge coupled device (ICCD) – I. In-depth detector characterization

While the adsorption properties of transition metal catalysts have been widely studied, leading to the discovery of various scaling relations, descriptors of catalytic activity, and well-established computational models, a similar understanding of semiconductor catalysts has not yet been achieved. In this work, we present a high-throughput density functional theory investigation into the adsorption properties of 5 oxides of interest to the photocatalytic CO 2 reduction reaction: TiO 2 (rutile and anatase), SrTiO 3 , NaTaO 3 , and CeO 2 . Using a systematic approach, we exhaustively identify unique surfaces and construct adsorption structures to undergo geometry optimizations. We then perform a data-driven analysis, which reveals the presence of weak adsorption energy scaling relations, the propensity of adsorbates of interest to interact with oxygen surface sites, and the importance of slab deformation upon adsorption. Our findings are presented in the context of experimental observations and in comparison to previously studied classes of catalysts, such as pure metals and tellurium-containing semiconductors, and reinforce the need for a comprehensive approach to the study of site-specific surface phenomena on semiconductors.

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

Thermomechanical analysis and modeling of a high-temperature light trapping planar cavity solar receiver

The development of durable particle-based high-temperature solar receivers is critical for advancing concentrating solar-thermal (CST) technologies to enable high-efficiency power generation and industrial process heat. Here, this study presents a computational framework to evaluate the thermomechanical performance of a proposed enclosed light-trapping planar cavity receiver designed for particle-based thermal energy systems. The receiver incorporates absorptive cavities and fluidized particle-bed channels to enhance heat capture and reduce thermal losses. Finite element analysis (FEA) is employed to assess stress, strain, and creep-fatigue behavior under concentrated solar flux using realistic thermal boundary conditions derived from coupled system models and experimental assembly parameters. The analysis investigates the influence of particle-to-wall heat transfer coefficients (HTC) ranging from 800 to 1800 W/m 2 .K on the thermomechanical response of six candidate high-temperature alloys: Alloy 740H, Alloy 282, Alloy 617, 316H, Alloy 230, and 800H. Results show that increasing HTC reduces thermal gradients, leading to lower stresses and strains and extended minimum predicted creep life. While all materials satisfy fatigue life requirements under the investigated conditions, significant differences in creep resistance are observed. Alloy 740H consistently exhibits the longest minimum predicted creep life and the most favorable durability margins, followed by Alloy 282, with the remaining materials showing reduced creep resistance under identical loading. The reported creep lives are conservative lower-bound estimates intended for comparative material evaluation. This framework highlights the critical roles of material selection and geometry optimization in improving mechanical durability and reliability of solar-thermal receivers, forming a foundation for future experimental validation and design optimization.

14 SOLAR ENERGY↗
DOE OSTI2025

Systematic computational study of oxide adsorption properties for applications in photocatalytic CO 2 reduction

While the adsorption properties of transition metal catalysts have been widely studied, leading to the discovery of various scaling relations, descriptors of catalytic activity, and well-established computational models, a similar understanding of semiconductor catalysts has not yet been achieved. In this work, we present a high-throughput density functional theory investigation into the adsorption properties of 5 oxides of interest to the photocatalytic CO 2 reduction reaction: TiO 2 (rutile and anatase), SrTiO 3 , NaTaO 3 , and CeO 2 . Using a systematic approach, we exhaustively identify unique surfaces and construct adsorption structures to undergo geometry optimizations. We then perform a data-driven analysis, which reveals the presence of weak adsorption energy scaling relations, the propensity of adsorbates of interest to interact with oxygen surface sites, and the importance of slab deformation upon adsorption. Our findings are presented in the context of experimental observations and in comparison to previously studied classes of catalysts, such as pure metals and tellurium-containing semiconductors, and reinforce the need for a comprehensive approach to the study of site-specific surface phenomena on semiconductors.

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

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modeling of transverse stimulated Raman scattering in KDP/DKDP in large-aperture plates suitable for polarization control

Transverse stimulated Raman scattering (TSRS) in potassium dihydrogen phosphate (KDP) and deuterated potassium dihydrogen phosphate (DKDP) plates for large-aperture, inertial confinement fusion (ICF)-class laser systems is a well-recognized limitation giving rise to parasitic energy conversion and laser-induced damage. The onset of TSRS is manifested in plates exposed to the ultraviolet section of the beam. TSRS amplification is a coherent process that grows exponentially and is distributed nonuniformly in the crystal and at the crystal surfaces. To understand the growth and spatial distribution of TSRS energy in various configurations, a modeling approach has been developed to simulate the operational conditions relevant to ICF-class laser systems. Specific aspects explored in this work include (i) the behavior of TSRS in large-aperture crystal plates suitable for third-harmonic generation and use as wave plates for polarization control in current-generation ICF-class laser system configurations; (ii) methods, and their limitations, of TSRS suppression and (iii) optimal geometries to guide future designs.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Exploring the impacts of conformer selection methods on ion mobility collision cross section predictions

The prediction of structure dependent molecular properties, such as collision cross sections as measured using ion mobility spectrometry, are crucially dependent on the selection of the correct population of molecular conformers. Here, we report an in-depth evaluation of multiple conformation selection techniques, including simple averaging, Boltzmann weighting, lowest energy selection, low energy threshold reductions, and similarity reduction. Generating 50,000 conformers each for 18 molecules, we used the In Silico Chemical Library Engine (ISiCLE) to calculate the collision cross sections for the entire dataset. First, we employed Monte Carlo simulations to understand the variability between conformer structures as generated using simulated annealing. Then we employed Monte Carlo simulations to the aforementioned conformer selection techniques applied on the simulated molecular property—the ion mobility collision cross section. Based on our analyses, we found Boltzmann weighting to be a good tradeoff between precision and theoretical accuracy. Combining multiple techniques revealed that energy thresholds and root-mean-squared deviation-based similarity reductions can save considerable computational expense while maintaining property prediction accuracy. Molecular dynamic conformer generation tools like AMBER can continue to generate new lowest energy conformers even after tens of thousands of generations, decreasing precision between runs. This reduced precision can be ameliorated and theoretical accuracy increased by running density functional theory geometry optimization on carefully selected conformers.

Nielson, Felicity F.↗
DOE OSTI2021

Computational Spectroscopy of the Cr–Cr Bond in Coordination Complexes

In this work, we report the accurate computational vibrational analysis of the Cr–Cr bond in dichromium complexes using second-order multireference complete active space methods (CASPT2), allowing direct comparison with experimental spectroscopic data both to facilitate interpreting the low-energy region of the spectra and to provide insights into the nature of the bonds themselves. Recent technological development by the authors has realized such computation for the first time. Accurate simulation of the vibrational structure of these compounds has been hampered by their notorious multiconfigurational electronic structure that yields bond distances that do not correlate with bond order. Some measured Cr–Cr vibrational stretching modes, ν(Cr 2 ), have suggested weaker bonding, even for so-called ultrashort Cr–Cr bonds, while others are in line with the bond distance. Here, we optimize geometries and compute ν(Cr 2 ) with CASPT2 for three well-characterized complexes, Cr 2 (O 2 CCH 3 ) 4 (H 2 O) 2 , Cr 2 (mhp) 4 , and Cr 2 (dmp) 4 . We obtain CASPT2 harmonic ν(Cr 2 ) modes in good agreement with experiment at 282 cm –1 for Cr 2 (mhp) 4 and 353 cm –1 for Cr 2 (dmp) 4 , compute 50 Cr and 54 Cr isotope shifts, and demonstrate that the use of the so-called IPEA shift leads to improved Cr–Cr distances. Additionally, normal mode sampling was used to estimate anharmonicity along ν(Cr 2 ), leading to an anharmonic mode of 272 cm –1 for Cr 2 (mhp) 4 and 333 cm –1 for Cr 2 (dmp) 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling Molecules under Pressure with Gaussian Potentials

The computational modeling of molecules under high pressure is a growing research area that augments experimental high-pressure chemistry. Here, a new electronic structure method for modeling atoms and molecules under pressure, Gaussians On Surface Tesserae Simulate HYdrostatic Pressure (GOSTSHYP) approach, is introduced. In this method, a set of Gaussian potentials is distributed evenly on the van der Waals surface of the investigated chemical system, leading to a compression of the electron density and the atomic scaffold. Since no parameters other than pressure need to be specified, GOSTSHYP allows straightforward geometry optimizations and ab initio molecular dynamics simulations of chemical systems under pressure for nonexpert users. In this work, calculated energies, bond lengths, and dipole moments under pressure fall within the range of established computational methods for high-pressure chemistry. A Diels-Alder reaction and the cyclotrimerization of acetylene showcase the ability of GOSTSHYP to model pressure-induced chemical reactions. The connection to mechanochemistry is pointed out.

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

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An Accurate, Affordable Density Functional Tight-Binding Model for Excited State Hydrocarbon Polymer Molecular Dynamics

We have developed a density functional tight-binding model for hydrocarbon excited state dynamics by referencing to high-level electronic structure theory and incorporating a many-body repulsive energy. We then validate our model against n-octane geometry optimizations, bond dissociation scans, and vibrational frequencies. Our model is approximately 10 3 times more efficient than hybrid time-dependent density functional theory calculations with comparable accuracy. In conclusion, our efforts enable longer timescale excited state simulations of photochemistry and scattering of incident radiation.

Chemical calculations↗
DOE OSTI2026

Analytic Nuclear Gradients Including Oriented External Electric Fields in a Molecule-Fixed Frame

Electric-field-assisted chemistry has attracted much attention in recent years, particularly in the context of oriented external electric fields for controlling molecular structure and reactivity. Such fields have been explored in a wide range of applications, including switching materials, nanoparticles, controllable catalysts, medicines, and clinical therapies. However, the determination of fixed fields in the laboratory frame becomes ineffective for flexible molecules, as conformational changes can significantly alter the relative orientation between the applied field and molecular structure. In this work, we propose two molecular reference frames─the principal axis frame and the local reference frame─to define oriented electric fields within the molecular framework. These coordinate systems powerfully eliminate ambiguities in the relative orientation between the applied field and the molecule. Analytic nuclear gradients in the presence of external electric fields are derived and implemented, with an initial application to field-dependent geometry optimizations of cis - and trans -formanilide. Analysis of the resulting field-induced equilibrium structures reveals distinct structural responses, validating the accuracy and robustness of the proposed formalism. The analytic gradient framework enables systematic investigations of molecular properties and reactivity under arbitrarily oriented electric fields, opening new opportunities for computational modeling and rational design in electric-field-controlled chemistry.

electric fields↗
DOE OSTI2020

Electron Affinity and Electronic Structure of Hexafluoroacetone (HFA) Revealed by Photodetaching the [HFA] •– Radical Anion

A great deal of effort has been focused on developing a metal-free catalytic system for epoxidation of unreactive alkenes. Fluoroketones are thought as remarkably promising catalysts for epoxidation reactions. The combination of fluorinated alcohols and catalytic amounts of hexafluoroacetone (HFA) gives a versatile and effective medium for epoxidation of various olefins with hydrogen peroxide. However, the fundamental physicochemical properties of HFA remained largely unclear, although they were very important to understand the related interactions. Here, in this study, we performed a joint study on the electron affinity and electronic structure of HFA employing negative ion photoelectron (NIPE) spectroscopy and quantum chemistry calculations. Two distinct bands with complicated vibrational progressions were observed in the 193 nm NIPE spectrum. The adiabatic/vertical detachment energies (ADE/VDE) were derived to be 1.42/2.06 and 4.43/4.86 eV for the ground singlet state and excited triplet state, respectively. Using the optimized geometries and vibrational frequencies of the anion and the neutral, the Franck–Condon factors were calculated for electron detachments to produce HFA in its lowest singlet and triplet states. Good agreements are obtained hereby for both bands between the experimental and calculated NIPE spectra, when taking into account combination vibrational excitations, unequivocally revealing that HFA possesses a singlet ground state with a giant singlet–triplet energy difference (ΔE ST ). The electron affinity (EA) and ΔE ST of HFA were therefore determined to be EA = 1.42 ± 0.02 eV and ΔE ST = −3.01 eV.

Wang, Lei↗
DOE OSTI2023

Excited State Dynamics of a Conformationally Fluxional Copper Coordination Complex

The conversion of solar energy into chemical fuel represents a capstone goal of the 21st century and has the potential to supply terawatts of power in a globally distributed manner. However, the disparate time scales of photodriven charge separation (~fs) and steps in chemical reactions (~μs) represent an inherent bottleneck in solar-to-fuels technology. Here to address this discrepancy, we are developing earth-abundant coordination complexes that undergo light-induced conformational rearrangements such that charge separation (CS) is hastened, while charge recombination (CR) is slowed. To these ends, we report the preparation and characterization of a new series of conformationally fluxional copper coordination complexes that contain a twisted intramolecular charge transfer (TICT) fluorophore as part of their ligand scaffold. Structural and spectroscopic characterization of the Cu(I) and Cu(II) complexes formed with these ligands in their ground states establish oxidation state-dependent conformational dynamicity, while time-resolved emission and transient absorption spectroscopies define the photopysical parameters of photo-induced excited states. Building on initial reports with a related set of molecules, the improved ligand design presented here greatly simplifies the observed photophysics, effectively shutting down unwanted ligand-centered excited states previously observed. Finally, time-dependent density functional theory (TDDFT) analyses reveal an unusual metal-to-TICT electronic transition only reported once before, and though the formation of a CS state is not observed directly through experiments, TDDFT geometry optimizations in the excited states support the formation of transient Cu(II) CS species, lending credence to the potential success of our approach. These studies establish a clear model for the excited state dynamics at play in proof-of-concept systems and clarify key design parameters for future optimizations toward achieving long-lived CS via photoinduced conformational gating.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗