Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complex Geometry”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

3D printed polymer gaskets for custom-form batteries

Additive manufacturing (AM) processes, like 3D printing, help to facilitate complex and customizable battery geometries which can provide design freedom and enhance volumetric energy density within electronic devices. AM materials must have the thermal and mechanical properties that enable printability, and when used in batteries, AM materials must also be chemically and electrochemically compatible with the battery chemistry. The compatibility between AM materials and the battery is of particular importance for the cell packaging materials which must be inert and are often overlooked. This study systematically studies AM-compatible polymeric materials for use as gaskets in lithium-ion cells. The materials investigated include three thermoplastics suitable for material extrusion printing: polylactic acid (PLA), polycarbonate, and polypropylene/polyethylene copolymer (PPPEC); and two photoresins suitable for vat photopolymerization (VPP) printing: an acrylate-based photoresin and a polyethylene glycol diacrylate photoresin. The AM gasket materials were tested in comparison to a conventional commercial polypropylene gasket. Mechanical testing (swell measurements and material stiffness) and electrochemical testing (linear sweep voltammetry and galvanostatic cycling of full cells) demonstrated that PLA and the VPP polymers were the least compatible with the lithium-ion battery chemistry, despite their prevalent use in studies of AM batteries, and that PPPEC was the most compatible.

3D printing↗
DOE OSTI2025

Photocurrent Nanoscopy of Quantum Hall Bulk

Understanding nanoscale electronic and thermal transport of two-dimensional (2D) electron systems in the quantum Hall regime, particularly in the bulk insulating state, poses considerable challenges. One of the primary difficulties arises from the presence of chiral edge channels, whose transport behavior obscures the investigation of the insulating bulk. Using near-field optical and photocurrent nanoscopy, we probe real-space variations of the optical and thermal dynamics of graphene in the quantum Hall regime without relying on complex sample or electrode geometries. Near the charge-neutrality point, we detect strong optical and photothermal signals from resonant inter-Landau-level (LL) magnetoexciton excitations between the zeroth and ± 1 st LLs, which gradually weaken with increasing doping due to Pauli blocking. Interestingly, at higher doping levels and full-integer LL fillings, photothermal signals reappear across the entire sample over an approximately 10 − μ m scale, indicating unexpectedly long cooling lengths and nonlocal photothermal heating through the insulating bulk. This observation suggests that thermal conductivity persists for the localized states even as electronic transport is suppressed—a clear violation of the Wiedemann-Franz law. Our experiments provide novel insights into nanoscale thermal and electronic transport in incompressible 2D gases, highlighting the roles of magnetoexcitons and chiral edge states in the thermo-optoelectric dynamics of the Dirac quantum Hall state. Published by the American Physical Society 2025

36 MATERIALS SCIENCE↗
DOE OSTI2024

UHT-CAMANCHE: Ultra-High Temperature Ceramic Additively Manufactured Compact Heat Exchangers

The conceptual basis for this project is the convergence of advanced ultra-high temperature ceramic materials and additive manufacturing technologies to produce compact ceramic heat exchangers with complex internal flow path geometries. Task areas were broadly divided into materials and manufacturing development, heat exchanger design, component testing, and techno-economic analysis. Technical challenges included the design and commissioning of new test facilities, improving feature resolution and deposition rate of ceramic additive manufacturing techniques, establishing process-structure-property relationships in additively manufactured ultra-high temperature ceramics, and assessing high temperature materials compatibility in CO 2 environments. The primary candidate material evaluated in this work is a composite comprising zirconium diboride (ZrB2) with 30 vol. % silicon carbide (SiC) which was selected based on its desirable combination of high temperature mechanical properties, high thermal and electrical conductivities, and oxidation resistance. High solids loaded ZrB2-SiC pastes suitable for extrusion-based additive manufacturing were developed for the first time as part of this work. Materials compatibility studies indicate this material oxidizes in CO 2 to form a protective borosilicate scale which transforms to pure silica above 1000°C. Parts made by additive manufacturing displayed enlarged grain sizes produced by pressureless sintering as compared to hot-press sintering. Increases in microstructural coarseness have outsized effect on oxidation performance up to 1400°C due to incomplete oxidation of coarse large diameter SiC particles resulting in lower amounts of silica that apparently inhibit protective scale formation. Additive manufacturing as a forming technique did not appear to significantly affect thermal conductivity, hardness, or elastic modulus, though flexural strength was reduced by half or more as compared to traditionally hot-pressed materials. This effect was attributed to the presence of strength-limiting flaws (ca. 40 microns in size) originating from extrudate inhomogeneities that could potentially be eliminated with further process improvements. Attempts to attain economies of scale for production of multi-kilowatt scale heat exchangers by ceramic additive manufacturing proved difficult. Lack of automation and a modest extrusion rate while retaining fine feature resolution made the overall process labor intensive and limited experimental throughput. A number of full-scale components were taken through post-process heat treatments including drying, binder burnout, and sintering; however, none survived without significant flaws or cracks. Therefore, no operational data from a newly installed heat exchanger test loop were able to be obtained during the performance period. Continued research and development is recommended to improve economic feasibility of ceramic additive manufacturing by standardizing the use of advanced sensors, artificial intelligence, and automation tools to reduce associated labor costs and accelerate production rates. The materials and manufacturing techniques demonstrated in this work are likely to find applications in defense and energy applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Statistically Resolved Planetary Boundary Layer Height Diurnal Variability Using Spaceborne Lidar Data

The Planetary Boundary Layer Height (PBLH) significantly impacts weather, climate, and air quality. Understanding the global diurnal variation of the PBLH is particularly challenging due to the necessity of extensive observations and suitable retrieval algorithms that can adapt to diverse thermodynamic and dynamic conditions. This study utilized data from the Cloud-Aerosol Transport System (CATS) to analyze the diurnal variation of PBLH in both continental and marine regions. By leveraging CATS data and a modified version of the Different Thermo-Dynamics Stability (DTDS) algorithm, along with machine learning denoising, the study determined the diurnal variation of the PBLH in continental mid-latitude and marine regions. The CATS DTDS-PBLH closely matches ground-based lidar and radiosonde measurements at the continental sites, with correlation coefficients above 0.6 and well-aligned diurnal variability, although slightly overestimated at nighttime. In contrast, PBLH at the marine site was consistently overestimated due to the viewing geometry of CATS and complex cloud structures. The study emphasizes the importance of integrating meteorological data with lidar signals for accurate and robust PBLH estimations, which are essential for effective boundary layer assessment from satellite observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

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

Advanced manufacturing techniques, such as additive manufacturing (AM), that can directly integrate phase change materials (PCMs) have garnered interest in recent years due to their potential for development of highly efficient thermal energy storage architectures. Complex, high surface area geometries embedded with PCMs that are only feasible with AM can improve thermal management with reduced material waste. Our work focuses on developing composite filaments with microencapsulated phase change materials (MEPCM) bound within a single or dual polymer matrix that can be processed through standard filament extruders and additively manufactured using off-the-shelf 3D printers. Polymer powders, rather than polymer pellets, were key to homogenously mixed filaments achieving high MEPCM loadings with no deterioration in thermal energy storage (TES) capability during extrusion. Composite filaments contain upwards of 60 wt% MEPCM and were printed without loss in feature resolution, print speed, or layer adhesion. Storage enthalpies of printed composites range from 100 - 130 kJ/kg, which were within 5% of the theoretical enthalpy based on weight fraction of MEPCM and maintained enthalpies within 1% over 500 thermal cycles. We can reliably manufacture low density, high surface area structures like 15% gyroid infill, along with dense, compact pucks at a 100% concentric infill. Prints were also scalable to a 900 cm3 honeycomb infill heat exchanger model that has an estimated energy storage capacity of 9 Wh.

3D printing↗
DOE OSTI2024

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

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

Fluid Mechanics Challenges in Direct-Ink-Writing Additive Manufacturing

Direct-ink writing (DIW) has rapidly become a versatile 3D fabrication method due to its ability to deposit a wide range of complex fluids into customizable 3D geometries. This review highlights key fundamental fluid mechanics and soft matter challenges across the different stages of the DIW printing process. The rheology of fluids and suspensions governs the flow behavior through narrow nozzles, posing questions about extrudability, confined flow dynamics, and clogging mechanisms. Downstream, the formation and deposition of extruded filaments involve extensional flows and potential instabilities, while postdeposition dynamics introduces complexities related to yield stress and structural stability. These stages are inherently interdependent, as optimizing material composition without considering filament stability risks compromising the final structure. As DIW applications expand through advanced ink formulations, developing fundamental fluid mechanics frameworks is essential to replace trial-and-error approaches with predictive design methodologies to enable more precise control and improved reliability of the printing process.

3D printing↗
DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Improved manufacturability of vacuum thermoforming molds by segmentation

Thin polymer shell components used in consumer goods, the automotive industry, and marine applications are often produced by thermoforming. Depending on the quantity of parts, molding processes can increase throughput significantly when compared to processes such as milling or additive manufacturing. One disadvantage to mold-based manufacturing methods, however, is that molds must be produced to form the required geometry. Mold production is often complex and requires a large capital investment. For large-scale molding processes, mold production is sometimes complicated by the requirement for large working volume machine tools which are capable of machining the cavities to the appropriate tolerances for the selected molding process. Here, this paper focuses on a segmented design for a large draw ratio (>3:1 molded surface area to sheet surface area) thermoforming mold and a cost-benefit analysis of segmentation versus monolithic design. Discussion of the thermoforming tests are included. The primary goals were to improve manufacturability and reduce the overall cost of the mold.

Manufacturability↗
DOE OSTI2025

Building confidence in models for complex barrier systems for radionuclides

The modeling and simulation of the Cement-clay Interaction-Diffusion field (CI-D) experiment at the Mont Terri site in Switzerland presented here demonstrates that it is possible to capture the multiscale physical and chemical features of natural and engineered barrier systems for radionuclides. The simulations are successfully carried out with the newly developed CrunchODiTi high-performance computing software that accounts for multiple continua, including a continuum representing the electrical double layer (EDL) developed along negatively charged clay particles in clay rock. The simulation also accounts for both the complex three-dimensional (3D) geometry, expected as the norm in a geological waste repository, and the anisotropy of the geological formation. In addition, the high resolution of the model makes it possible to include "skin effects" developed at the interface between highly reactive materials, in this case between the high pH cement and the circumneutral but electrostatic Opalinus Clay. The successful history matching with the field experiment demonstrates that the distinct geochemical and physical properties of the cement and the Opalinus Clay in the CI-D experiment can be accounted for. Such analyses are essential for developing a defensible safety case for the underground storage of radioactive waste.

Sarsenbayev, Dauren↗
DOE OSTI2025

Fast and Invertible Simplicial Approximation of Magnetic‐Following Interpolation for Visualizing Fusion Plasma Simulation Data

We introduce a fast and invertible approximation for fusion plasma simulation data represented as 2D planar meshes with connectivities approximating magnetic field lines along the toroidal dimension in deformed 3D toroidal spaces. Scientific variables (e.g., density and temperature) in these fusion data are interpolated following a complex magnetic-field-line-following scheme in the toroidal space represented by a cylindrical coordinate system. This deformation in the 3D space poses challenges for root-finding and interpolation. To this end, we propose a novel paradigm for visualizing and analyzing such data based on a newly developed algorithm for constructing a 3D simplicial mesh within the deformed 3D space. Our algorithm generates a tetrahedral mesh that connects the 2D meshes using tetrahedra while adhering to the constraints on node connectivities imposed by the magnetic field-line scheme. Specifically, we first divide the space into smaller partitions to reduce complexity based on the input geometries and constraints on connectivities. Then, we independently search for a feasible tetrahedralization of each partition, considering nonconvexity. We demonstrate our method with two X-Point Gyrokinetic Code (XGC) simulation datasets on the International Thermonuclear Experimental Reactor (ITER) and Wendelstein 7-X (W7-X), and use an ocean simulation dataset to substantiate broader applicability of our method. An open source implementation of our algorithm is available at https://github.com/rcrcarissa/DeformedSpaceTet.

Ren, Congrong [The Ohio State Univ., Columbus, OH ↗
DOE OSTI2024

Accessing bands with extended quantum metric in kagome Cs 2 Ni 3 S 4 through soft chemical processing

Flat bands that do not merely arise from weak interactions can produce exotic physical properties, such as superconductivity or correlated many-body effects. The quantum metric can differentiate whether flat bands will result in correlated physics or are merely dangling bonds. A potential avenue for achieving correlated flat bands involves leveraging geometrical constraints within specific lattice structures, such as the kagome lattice; however, materials are often more complex. In these cases, quantum geometry becomes a powerful indicator of the nature of bands with small dispersions. We present a simple, soft-chemical processing route to access a flat band with an extended quantum metric below the Fermi level. By oxidizing Ni-kagome material Cs 2 Ni 3 S 4 to CsNi 3 S 4 , we see a two orders of magnitude drop in the room temperature resistance. However, CsNi 3 S 4 is still insulating, with no evidence of a phase transition. Using experimental data, density functional theory calculations, and symmetry analysis, our results suggest the emergence of a correlated insulating state of unknown origin.

Science & Technology - Other Topics↗
DOE OSTI2024

FY2024 Q4 L2 Milestone 8237 for Opacity-on-NIF

This document addresses parts 2 and 3 of the original milestone request. In the following writeup, the development of OpSpecTR at LLNL by LLNL and NNSS personnel will be detailed first, followed by a write-up of the development of Vme-resolved simulaVons of the opacity-on-NIF experiment. Capsule backlight simulaVons, hohlraum simulaVons, and sample simulaVons from CASSIO are combined to produce the line of sight from the capsule to an effecVve OpSpec posiVon. Spectra are generated by post-processing these simulaVons from either Spect3D, a commercial software by Prism ComputaVonal Sciences, or FESTR (Finite Element Spectroscopic Transport of RadiaVon), a LANL code. While the simulaVons do not at this stage proceed late enough to include the worst of the backgrounds generated from the hohlraum experimentally and each post-processing simulaVon includes only one ray for each spectrometer channel, these efforts represent proof of concept for a new capability to complete Vme-resolved and Vmegated simulaVons for the complete line of sight of the complex opacity-on-NIF geometry. Future work will be discussed at the end of the report for both efforts.

42 ENGINEERING↗
DOE OSTI2025

PyARC Status Report: New Integrations and Upgrades to the Fast Reactor Analysis Workflow Management Tool

PyARC was initially developed as an open source tool to support fast reactor analyses using the Argonne Reactor Computation (ARC) code suite as a part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Workbench initiative in FY17. The goal of this initiative is to provide a common user interface for model generation, real-time validation, execution, output processing, and visualization for all integrated codes. This is accomplished through the reliance on tools available in the Workbench framework and runtime environment. While initially developed to support the ARC codes, PyARC was extended in FY22 to wrap other NEAMS and non-ARC codes, including Griffin and OpenMC, in the supported other neutronics workflows, and support users in the adoption of NEAMS-supported high fidelity analysis codes. Most recently, NUBOW-3D, a recently adopted ARC code, was integrated to support reactor bowing calculations as well. Integration of these codes into the NEAMS Workbench directly benefits the advanced reactor modeling community by: • Providing a set of controlled, maintained, documented and validated scripts to generate inputs, which promotes best practices, reduces the learning curve, and facilitates project collaboration. • Improving the user experience: the Workbench interface provides assistance for building an input through auto-completion, real-time validation, document navigation, and geometry and results visualization. • Automating complex calculations and workflows for reactor analysis. • Helping users transition to using high-fidelity NEAMS codes along-side the ARC codes. In FY22, a progress report was published that described the state of each of the tools integrated into PyARC. Since then, there have been many enhancements and upgrades to the existing integrations as well as entirely new code integrations as well. This report details all new integrations and major developments in PyARC since the version 2.0.0 release highlighted in the FY22 report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Scalable Fabrication of a Fibrous Amine-functionalized Matrix (FAM) Sorbent for Critical Mineral Recovery

We report a novel flat sheet Fibrous Amine-functionalized Matrix (FAM) sorbent platform designed for efficient and selective capture of CM from dilute solutions. The FAM sorbent features crosslinked amine microfilms coated onto/within a glass fiber matrix, providing fast mass transfer and excellent mechanical stability. Systematic batch and flow-through tests with FAM revealed rapid metal uptake kinetics and high capacity for representative species, achieving ~90 mg/g of Gallium, ~100 mg/g of Cobalt, and ~90 mg/g for Neodymium. Moreover, multiple eluents, including mineral acids and complexing agents, enabled highly effective desorption of adsorbed metals, demonstrating the feasibility of regenerating FAM sorbents. Importantly, tests with authentic coal ash leachate demonstrated strong selectivity toward U.S. Department of Energy (DOE)-listed CM and rare earth elements over abundant base cations, confirming the robustness of FAM in realistic complex solutions. The flat sheet geometry was amenable to scaling into durable spiral wound modules, highlighting the potential for future regeneration and reuse. This work establishes FAM sorbents as a promising platform for the recovery of CM from wastewaters, advancing both resource sustainability and environmental stewardship.

critical mineral recovery↗
DOE OSTI2025

Coordination Chemistry of Solvated Metal Ions in Soft Donor Solvents

The structures of hexaammine solvated indium(III) and thallium(III) ions in liquid ammonia solution are determined by EXAFS. Both complexes have regular octahedral coordination geometry with mean In-N and Tl-N bond distances of 2.23(1) and 2.29(2) Å, respectively. Ammine solvated thallium(III) in liquid ammonia is characterized with 205Tl NMR measurements. Solvents such as liquid ammonia, N,N-dimethylthioformamide (DMTF), trialkyl and triphenyl phosphite and phosphine are strong electron pair donors and thereby able to form bonds with a large covalent contribution with strong electron pair acceptors. A survey of reported structures of ammine, DMTF, trialkyl and triphenyl phosphite and phosphine solvated metal ions in the solid state and solution is presented. The M-N and M-S bond distances in ammine and DMTF solvated metal ions are compared with the M-O bond distance in the corresponding metal ion hydrates, expected to form mainly electrostatic interactions with metal ions. The d10 metal ions have high ability to form bonds with a high degree of covalency with increasing ability down the group and with decreasing charge of the metal ion. The difference in M-N and M-O bond distances between ammine solvated and hydrated metal ions with the same coordination geometry decreases significantly with the increasing ability of the metal ion to form bonds with a large covalent contribution. This difference correlates well with the covalent bonding index, γM2*r.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Homoleptic 1,2-benzenedithiolate complexes of thorium and uranium

Reaction of 4 equiv. of [Li(TMEDA)] 2 [1,2-S 2 C 6 H 4 ] with [ThCl 4 (DME) 2 ] or [UCl 4 (THF) 3 ] in THF results in formation of [Li(THF) 2 ] 4 [An(1,2-S 2 C 6 H 4 ) 4 ] (An = Th, 1; An = U, 2), whereas reaction of 4 equiv. of [Li(TMEDA)] 2 [1,2-S 2 C 6 H 4 ] with UCl 4 in Et 2 O results in formation of [Li(TMEDA)] 4 [U(1,2-S 2 C 6 H 4 ) 4 ] (3). Complexes 1–3 represent the first reported benzenedithiolate complexes of the actinides. Here, they were characterized by NMR spectroscopy and X-ray crystallography. In the solid state, complexes 1–3 exhibit triangular dodecahedral geometries about their actinide centers. Additionally, their Li + cations are bound by two sulfur atoms of adjacent [1,2-S 2 C 6 H 4 ] 2− ligands, in addition to two solvent donor atoms. In solution, complexes 2 and 3 exhibit spectral data consistent with S 4 symmetry (and non-exchanging Li + sites), whereas complex 1 exhibits spectral properties consistent with labile Li + cations.

X-ray↗
DOE OSTI2024

Ligand Many-Body Expansion as a General Approach for Accelerating Transition Metal Complex Discovery

Methods that accelerate the evaluation of molecular properties are essential for chemical discovery. While some degree of ligand additivity has been established for transition metal complexes, it is underutilized in asymmetric complexes, such as the square pyramidal coordination geometries highly relevant to catalysis. To develop predictive methods beyond simple additivity, we apply a many-body expansion to octahedral and square pyramidal complexes and introduce a correction based on adjacent ligands (i.e., the cis interaction model). We first test the cis interaction model on adiabatic spin-splitting energies of octahedral Fe(II) complexes, predicting DFT-calculated values of unseen binary complexes to within an average of 1.4 kcal/mol. Uncertainty analysis reveals the optimal basis, comprising the homoleptic and mer symmetric complexes. We next show that the cis model (i.e., the cis interaction model solved for the optimal basis) infers both DFT- and CCSD(T)-calculated model catalytic reaction energies to within 1 kcal/mol on average. The cis model predicts low-symmetry complexes with reaction energies outside the range of binary complex reaction energies. We observe that trans interactions are unnecessary for most monodentate systems but can be important for some combinations of ligands, such as complexes containing a mixture of bidentate and monodentate ligands. Lastly, we demonstrate that the cis model may be combined with Δ-learning to predict CCSD(T) reaction energies from exhaustively calculated DFT reaction energies and the same fraction of CCSD(T) reaction energies needed for the cis model, achieving around 30% of the error from using the CCSD(T) reaction energies in the cis model alone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗