Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural 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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Phase transitions of correlations in black hole geometries

We study the holographic realization of optimized correlation measures—measures of quantum correlation that generalize elementary entropic formulas—in two-dimensional thermal states dual to spacetimes with a black hole horizon. We consider the symmetric bipartite optimized correlation measures: the entanglement of purification, Q-correlation, R-correlation, and squashed entanglement, as well as the mutual information, a nonoptimized correlation measure, and identify the bulk surface configurations realizing their geometric duals over the parameter space of boundary region sizes and the black hole radius. This parameter space is divided into phases associated with given topologies for these bulk surface configurations, and first-order phase transitions occur as a new topology of bulk surfaces becomes preferred. The distinct phases can be associated with different degrees of correlation between the boundary regions and the thermal environment. The Q-correlation has the richest behavior, with a structure of nested optimizations leading to two topologically distinct bulk surface configurations being equally valid as geometric duals at generic points in the phase diagram.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

G 0 W 0 Ionization Potentials of First-Row Transition Metal Aqua Ions

Here, we report computations of the vertical ionization potentials within the GW approximation of the near-complete series of first-row transition metal (V–Cu) aqua ions in their most common oxidation states, i.e., V 3+ , Cr 3+ , Cr 2+ , Mn 2+ , Fe 3+ , Fe 2+ , Co 2+ , Ni 2+ , and Cu 2+ . The d-orbital occupancy of these systems spans a broad range from d 2 to d 9 . All of the structures were first optimized at the density functional theory level using a large cluster of explicit water molecules that are embedded in a continuum solvation model. Vertical ionization potentials were computed with the one-shot G 0 W 0 approach on a range of transition metal ion clusters (6, 18, 40, and 60 explicit water molecules), wherein the convergence with respect to the basis set size was evaluated using the systems with 40 water molecules. We assess the results using three different density functional approximations as starting points for the vertical ionization potential calculations, namely, G 0 W 0 @PBE, G 0 W 0 @PBE0, and G 0 W 0 @r 2 SCAN. While the predicted ground-state structures are similar to all three exchange-correlation functionals, the vertical ionization potentials were in closer agreement with experiment when using the G 0 W 0 @PBE0 and G 0 W 0 @r 2 SCAN approaches, with the r 2 SCAN-based calculations being significantly less expensive. Computed bond distances and vertical ionization potentials for all structures are in good agreement with available experimental data.

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

Schottky Anomalies in the Low-Temperature Specific Heat of YBa{sub 2}Cu{sub 3}O{sub y} HTSC

The evolution of the Schottky anomalies in the temperature dependence of the specific heat (in the temperature range T = 2–10 K and the magnetic field range H = 0–9 T) of polycrystalline YBa{sub 2}Cu{sub 3}O{sub y} HTSCs is studied when the oxygen content in coarse-crystalline samples changes and when the degree of structural disordering in optimally doped fine-crystalline samples of the same compound changes. The number of paramagnetic centers with uncompensated spin moments, which generate the Schottky anomalies, and the splitting of the lower energy levels are estimated. A correlation between these parameters and both the number of oxygen vacancies in chain planes and the method of oxygen ordering in chains is found. The experimental data are compared with the results of studying the magnetic susceptibility at T > T{sub c} in the same samples. Based on this comparison, we can conclude that the nature of the Schottky centers is related to the formation of paramagnetic centers at breaks in the Cu1–O4 copper–oxygen chains in the basal planes of the crystal structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Silica accelerates the selective hydrogenation of CO 2 to methanol on cobalt catalysts

The reaction pathways on supported catalysts can be tuned by optimizing the catalyst structures, which helps the development of efficient catalysts. Such design is particularly desired for CO 2 hydrogenation, which is characterized by complex pathways and multiple products. Here, we report an investigation of supported cobalt, which is known for its hydrocarbon production and ability to turn into a selective catalyst for methanol synthesis in CO 2 hydrogenation which exhibits good activity and stability. The crucial technique is to use the silica, acting as a support and ligand, to modify the cobalt species via Co-O-SiOn linkages, which favor the reactivity of spectroscopically identified *CH 3 O intermediates, that more readily undergo hydrogenation to methanol than the C-O dissociation associated with hydrocarbon formation. Cobalt catalysts in this class offer appealing opportunities for optimizing selectivity in CO 2 hydrogenation and producing high-grade methanol. By identifying this function of silica, we provide support for rationally controlling these reaction pathways.

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

Search for Catalysts by Inverse Design: Artificial Intelligence, Mountain Climbers, and Alchemists

In silico catalyst design is a grand challenge of chemistry. Traditional computational approaches have been limited by the need to compute properties for an intractably large number of possible catalysts. Recently, inverse design methods have emerged, starting from a desired property and optimizing a corresponding chemical structure. Techniques used for exploring chemical space include gradient-based optimization, alchemical transformations, and machine learning. Though the application of these methods to catalysis is in its early stages, further development will allow for robust computational catalyst design. This review provides an overview of the evolution of inverse design approaches and their relevance to catalysis. In conclusion, the strengths and limitations of existing techniques are highlighted, and suggestions for future research are provided.

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

Structure and dynamics of small polyimide oligomers with silicon as a function of aging

The effect of UV curing and shearing on the structure and behavior of a polyimide (PI) binder as it disperses silicon particles in a battery electrode slurry was investigated. PI dispersant effectiveness increases with UV curing time, which controls the overall binder molecular weight. Here, the shear force during electrode casting causes higher molecular weight PI to agglomerate, resulting in battery anodes with poorly dispersed Si particles that do not cycle well. It is hypothesized that when PI binder is added above a critical amount, it conformally coats the silicon particles and greatly impedes Li ion transport. There is an “interzonal region” for binder loading where it disperses silicon well and provides a coverage that facilitates Li transport through the anode material and into the silicon particles. These results have implications in ensuring reproducible electrode manufacturing and increasing cell performance by optimizing the PI structure and coordination with the silicon precursor.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On the scaling of wind turbine rotors

This paper formulates laws for scaling wind turbine rotors. Although the analysis is general, the article primarily focuses on the subscaling problem, i.e., on the design of a smaller-sized model that mimics a full-scale machine. The present study considers both the steady-state and transient response cases, including the effects of aerodynamic, elastic, inertial, and gravitational forces. The analysis reveals the changes to physical characteristics induced by a generic change of scale, indicates which characteristics can be matched faithfully by a subscaled model, and states the conditions that must be fulfilled for desired matchings to hold. Based on the scaling laws formulated here, the article continues by considering the problem of designing scaled rotors that match desired indicators of a full-scale reference. To better illustrate the challenges implicit in scaling and the necessary tradeoffs and approximations, two different approaches are contrasted. The first consists in a straightforward geometric zooming. An analysis of the consequences of zooming reveals that, although apparently simple, this method is often not applicable in practice, because of physical and manufacturing limitations. This motivates the formulation of scaling as a constrained optimal aerodynamic and structural matching problem of wide applicability. Practical illustrations are given considering the scaling of a large reference 10 MW wind turbine of about 180 m in diameter down to three different sizes of 54, 27, and 2.8 m. Results indicate that, with the proper choices, even models characterized by very significant scaling factors can accurately match several key performance indicators. Additionally, when an exact match is not possible, relevant trends can at least be captured.

17 WIND ENERGY↗
DOE OSTI2024

Dosage optimization for reducing tumor burden using a phenotype-structured population model with a drug-resistance continuum

Abstract Drug resistance is a significant obstacle to effective cancer treatment. To gain insights into how drug resistance develops, we adopted a concept called fitness landscape and employed a phenotype-structured population model by fitting to a set of experimental data on a drug used for ovarian cancer, olaparib. Our modeling approach allowed us to understand how a drug affects the fitness landscape and track the evolution of a population of cancer cells structured with a spectrum of drug resistance. We also incorporated pharmacokinetic (PK) modeling to identify the optimal dosages of the drug that could lead to long-term tumor reduction. We derived a formula that indicates that maximizing variation in plasma drug concentration over a dosing interval could be important in reducing drug resistance. Our findings suggest that it may be possible to achieve better treatment outcomes with a drug dose lower than the levels recommended by the drug label. Acknowledging the current limitations of our work, we believe that our approach, which combines modeling of both PK and drug resistance evolution, could contribute to a new direction for better designing drug treatment regimens to improve cancer treatment.

Life Sciences & Biomedicine - Other Topics↗
DOE OSTI2018

Geometries of edge and mixed dislocations in bcc Fe from first principles calculations

We use density functional theory (DFT) to compute the core structures of $a_0[100](010)$ edge, $a_0[100](011)$ edge, $a_0/2[\bar{1}\bar{1}1](1\bar{1}0)$ edge, and $a_0/2[111](1\bar{1}0)$ $71^{\circ}$ mixed dislocations in body-centered cubic (bcc) Fe. The calculations are performed using flexible boundary conditions (FBC), which effectively allow the dislocations to relax as isolated defects by coupling the DFT core to an infinite harmonic lattice through the lattice Green function (LGF). We use the LGFs of the dislocated geometries in contrast to most previous FBC-based dislocation calculations that use the LGF of the bulk crystal. The dislocation LGFs account for changes in the topology of the crystal in the core as well as local strain throughout the crystal lattice. A simple bulk-like approximation for the force constants in a dislocated geometry leads to dislocation LGFs that optimize the core structures of the $a_0[100](010)$ edge, $a_0[100](011)$ edge, and $a_0/2[111](1\bar{1}0)$ $71^{\circ}$ mixed dislocations. This approximation fails for the $a_0/2[\bar{1}\bar{1}1](1\bar{1}0)$ dislocation however, so in this case we derive the LGF from more accurate force constants computed using a Gaussian approximation potential. The standard deviations of the dislocation Nye tensor distributions quantify the widths of the dislocation cores. The relaxed cores are compact, and the local magnetic moments on the Fe atoms closely follow the volumetric strain distributions in the cores. We also compute the core structures of these dislocations using eight different classical interatomic potentials, and quantify symmetry differences between the cores using the Fourier coefficients of their Nye tensor distributions. Most of the core structures computed using the classical potentials agree well with the DFT results. Furthermore the DFT core geometries provide benchmarking for classical potential studies of work-hardening, as well as substitutional and interstitial sites for computing solute-dislocation interactions that serve as inputs for mesoscale models of solute strengthening and solute diffusion near dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystal structure of toceranib, C 22 H 25 FN 4 O 2

The crystal structure of toceranib has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Toceranib crystallizes in space group P2 1 /c (#14) with a = 10.6899(6), b = 24.5134(4), c = 7.8747(4) Å, β = 107.7737(13)°, V = 1965.04(3) Å 3 , and Z = 4. The crystal structure consists of stacks of approximately planar molecules, with N–H∙∙∙O hydrogen bonds between the layers. The commercial reagent sample was a mixture of two or more phases with toceranib being the dominant phase. The difference between the Rietveld-refined and DFT-optimized structures is larger than usual. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗
DOE OSTI2018

A first principles study of energetics and electronic structural responses of uranium-based coordination polymers to Np incorporation

Abstract Recently developed coordination polymers (CPs) and metal organic frameworks (MOFs) may find applications in areas such as catalysis, hydrogen storage, and heavy metal immobilization. Research on the potential application of actinide-based CPs (An-CP/MOFs) is not as advanced as transition metal-based MOFs. In order to modify their structures necessary for optimizing thermodynamic and electronic properties, here, we described how a specific topology of a particular actinide-based CP or MOF responds to the incorporation of other actinides considering their diverse coordination chemistry associated with the multiple valence states and charge-balancing mechanisms. In this study, we apply a recently developed DFT-based method to determine the relative stability of transuranium incorporated CPs in comparison to their uranium counterpart considering both solid and aqueous state sources and sinks to understand the mechanism and energetics of charge-balanced Np 5+ incorporation into three uranium-based CPs. The calculated Np 5+ +H + incorporation energies for these CPs range from 0.33 to 0.52 eV, depending on the organic linker, when using the solid oxide Np source Np 2 O 5 and U sink UO 3 . Incorporation energies of these CPs using aqueous sources and sinks increase to 2.85–3.14 eV. The thermodynamic and structural analysis in this study aides in determining, why certain MOF topologies and ligands are selective for some actinides and not for others. This means that once this method is extended across a variety of CPs with their respective linker molecules and different actinides, it can be used to identify certain CPs with certain organic ligands being specific for certain actinides. This information can be used to construct CPs for actinide separation. This is the first determination of the electronic structure (band structure, density of states) of these uranium- and transuranium-based CPs which may eventually lead to design CPs with certain optical or catalytic properties. While the reduction of the DFT-determined-bandgap goes from 3.1 eV to 2.4 eV when going from CP1 to CP3, showing the influence of the linker, Np 6+ incorporation reduces the bandgap for CP1 and CP3, while increasing it for CP2. The coupled substitution of U 6+ →Np 5+ +H + reduces the bandgap significantly, but only for CP3.

Chemistry↗
DOE OSTI2019

Building a Consistent and Reproducible Database for Adsorption Evaluation in Covalent–Organic Frameworks

We present a workflow that traces the path from the bulk structure of a crystalline material to assessing its performance in carbon capture from coal's postcombustion flue gases. This workflow is applied to a database of 324 covalent-organic frameworks (COFs) reported in the literature, to characterize their CO 2 adsorption properties using the following steps: (1) optimization of the crystal structure (atomic positions and unit cell) using density functional theory, (2) fitting atomic point charges based on the electron density, (3) characterizing the pore geometry of the structures before and after optimization, (4) computing carbon dioxide and nitrogen isotherms using grand canonical Monte Carlo simulations with an empirical interaction potential, and finally, (5) assessing the CO 2 parasitic energy via process modeling. The full workflow has been encoded in the Automated Interactive Infrastructure and Database for Computational Science (AiiDA). Both the workflow and the automatically generated provenance graph of our calculations are made available on the Materials Cloud, allowing peers to inspect every input parameter and result along the workflow, download structures and files at intermediate stages, and start their research right from where this work has left off. In particular, our set of CURATED (Clean, Uniform, and Refined with Automatic Tracking from Experimental Database) COFs, having optimized geometry and high-quality DFT-derived point charges, are available for further investigations of gas adsorption properties. We plan to update the database as new COFs are being reported.

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

Durable and High-Performance SOECs Based on Proton Conductors for Hydrogen Production

Proton-conducting solid oxide electrolysis cells (P-SOECs) are a promising technology for cost-effective and efficient production of green hydrogen. Breakthroughs in materials development, optimization of cell structure, and achievement of high performance and durability are essential to significantly increase the commercial competitiveness of these technologies. The main objective of this project is to gain scientific knowledge for the rational design, fabrication, and demonstration of a robust, highly efficient, and low-cost SOEC technology based on a proton-conducting electrolyte membrane for hydrogen production. We focused on better understanding the degradation mechanisms of proton-conducting electrolytes, air electrodes, and catalyst materials under electrolysis mode to develop an effective strategy for rationalizing new materials that are vital for enhancing cell performance and durability. The scope includes enhancing the performance and durability of the electrolyte and electrode materials under realistic operating conditions, developing highly active and robust catalysts to minimize electrode losses while improving tolerance to contaminant poisoning, revealing the mechanism of enhanced activity and stability of the catalyst, and understanding the underlying degradation mechanisms. In addition, various characterization techniques were employed to gain a fundamental understanding of the materials’ behavior and their impact on cell performance, providing vital information to guide materials discovery and cell design. After defect chemistry engineering, the optimized donor and acceptor co-doped electrolytes BaMo/W 0.03 Ce 0.71 Yb 0.26 O 3-δ (BM/W03) showed substantially improved chemical stability against high concentrations of CO 2 and H 2 O compared to the state-of-the-art electrolyte (BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3-δ , BZCYYb1711) while maintaining comparable ionic conductivity and ionic transference number. To bypass the inherent trade-off between conductivity and chemical stability, we fabricated a bi-layer electrolyte composed of BZCYYb1711 coated with a highly-stable thin layer of BaHf 0.83 Yb 0.17 O 3-δ (BHYb). This bi-layer electrolyte displayed excellent chemical stability against high concentration CO 2 ; there was no detectable formation of BaCO 3 after exposure to 97% CO 2 (with 3% H 2 O) at 500 °C for 1000 hours and the rate of degradation in resistance was about 0.4% per 1,000 hours (kh). In contrast, the same BZCYYb1711 electrolyte without a BHYb coating degraded significantly under the same testing conditions; the degradation rate was increased to 5.1%/kh. In addition, a triple conducting air electrode Ba 0.9 Pr 0.1 Hf 0.1 Y0.1Co 0.8 O 3-δ (BPHYC) was developed by heavily doping transition metal ions into a proton-conducting material. This air electrode material, composed of 3 distinct phases, exhibits superior electrocatalytic activity due to the synergistic effect from the three component phases. Moreover, an active and durable catalyst, La 2 Ni 0.5 Fe 0.5 O 4+δ (LNF), was developed, showing excellent catalytic activity and contaminant tolerance, with a degradation rate of only 0.49%/kh when exposed to high concentrations of steam and Cr. Finally, single cells were constructed from the best electrolytes, electrodes, and catalyst coatings developed in this project. These cells demonstrated superior high current density at a given cell voltage, high roundtrip efficiency, and remarkable durability (up to 1000 hours of operation).

08 HYDROGEN↗
DOE OSTI2023

Extend an innovative HPC-Compatible Multiple Temporal-spatial Resolution Concurrent Finite Element Modeling Approach to Guide Laser Powder Bed Fusion Additive

Laser power bed fusing (PBF) additive manufacturing is a key enabling technology to manufacture highly complex and integrated automotive structures. However, the geometric complexity of PBF-AM technique also leads to highly non-uniform heating and cooling rate in the manufactured part, which may cause flaw formation and produce excessive and nonuniform residual stresses, which increase quality uncertainties and manufacture issues, leading to increases in cost and energy consumption in the form of rejected parts. In this research project, we developed an innovative Multi-Spatial-Temporal-Resolution Finite Element (MUST-FE) method and completed the corresponding high performance computation (HPC) platform-based in-house code, which enables high accuracy prediction of temperature and residual stress fields for component-scale PBF-AM manufacture in efficient computation time. The MUST-FE model is calibrated and validated with a “2D pad” AlSi10Mg experiments by matching the melt pool shape and dimension, and with a “XY-cross” AlSi10Mg experiment by matching the thermal distortion and residual stress. The innovative multi-resolution and concurrent modeling approach adopted in this code ensures accuracy and computational efficiency, which will enable energy-efficient and high-yield, low-cost manufacturing of optimized, qualifiable automotive structures and contribute towards reaching technical targets outlined in AMO’s Program Plan to develop additive manufacturing systems that deliver consistently reliable parts with predictable properties.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Elastic Properties of Lithium and Sodium Amides

This paper presents the optimization of the structural parameters of tetragonal and orthorhombic crystals α-LiNH{sub 2} and α-NaNH{sub 2} and computation of elastic constants and flexibility coefficients using the all-electron implementation of the Gaussian and augmented plane wave density functional method. The obtained elastic constants satisfy all conditions of mechanical stability of LiNH{sub 2} and α-NaNH{sub 2} crystals. From the viewpoint of elastic properties, α-LiNH{sub 2} crystal is isotropic, while α-NaNH{sub 2} crystal is rather anisotropic. Anisotropy under compression in α-NaNH{sub 2} crystal is higher than that under shear. The former is brittle, the latter is more plastic.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Synthesis and structural characterization of the new Zintl phases Ba 3 Cd 2 P 4 and Ba 2 Cd 2 P 3 . Rare example of small gap semiconducting behavior with negative thermopower within the range 300 K-700 K

The new Zintl phases Ba 3 Cd 2 P 4 and Ba 2 Cd 2 P 3 have been synthesized using Pb flux, which allowed for the growth of 4-5 mm large crystals. The structures were determined utilizing single-crystal X-ray diffraction methods. Both compounds crystalize in the monoclinic crystal system (space group C2/m (No. 12)) and their structures are closely related. The structure of Ba 3 Cd 2 P 4 can be seen as being comprised of divalent Ba atoms and conjoined CdP 4 tetrahedra in the form of [Cd 2 P 4 ] 6- layers. Within the layers, homoatomic P–P bonds are present, which if cleaved, leave two infinite [CdP 3 ] 7- chains running along the crystallographic b-axis. The other structure, that of Ba 2 Cd 2 P 3 , can be rationalized as also having divalent Ba atoms and conjoined CdP4 tetrahedra in the form of [Cd 2 P 3 ] 6- layers. These layers, again, can be visualized as chains that run down the crystallographic b-axis, which are further connected by P-P dimers. Electronic band structure calculations show that each structure has an optimal number of valence electrons, and therefore conform to the Zintl-Klemm concept. Accordingly, the two compounds can be considered small band gap semiconductors, with band gaps of ca. 0.1 eV and 0.6 eV for Ba3Cd2P4 and Ba 2 Cd 2 P 3 , respectively. Electrical resistivity measurements show that Ba3Cd2P4 displays a large resistivity value at room temperature and an experimental band gap of ca. 0.05 eV, which fits reasonably well with the theoretical predictions. Thermopower measurements show that throughout the temperature range 300 K-700 K, Ba 3 Cd 2 P 4 displays a negative Seebeck coefficient. Here, the extremum value of -84 μV is reached at 630 K, suggestive of an n-type semiconductor, a rarity among Zintl phases.

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

First considerations on the supporting structures of FCC-ee booster and collider in the arc regions

In 2022, the FCC Feasibility Study management mandated a working group to analyse the best configuration of the FCC-ee tunnel in the arc regions, in view of the construction of a mock-up of the arc half-cell. One of the main and most challenging goals of the study, named FCC-ee Arc Half-Cell Mock-up Project, was to perform a preliminary investigation on the principles of supporting the Short-Straight Sections and dipoles of the half-cells, both for the booster and for the collider machines. This is an important input needed for the choice of the best configuration of the relative placement of the booster with respect to the collider. The structural stiffness, mass and stability of the supporting structures must be optimized to minimize the vibrations transmitted/transferred to the magnetic system of the accelerators by elements such as pumps, water cooling system, beam thermomechanical stresses, powering elements, etc. To perform the study, tools such as CAD software, FEM and analytical techniques were employed. This paper summarizes the preliminary design concepts and the results of the simulations performed.

Accelerator Subsystems and Technologies↗
DOE OSTI2023

Controlling Molecular Structure and Spin with Multiconfigurational Quantum Chemistry (Final Technical Report)

For many first-row transition metal complexes, structure-property relationships can be obtained from high-level molecular geometry optimizations and subsequent electronic structure studies. The project funded under this award utilized newly implemented fully internally contracted (FIC) nuclear gradients for extended multi-state (XMS) complete active space second-order multireference perturbation theory (CASPT2) to explore geometry changes in first-row transition metal coordination complexes. At the start of the project, only one full geometry optimization using FIC-CASPT2 analytical gradients had been reported (J. Chem. Theory Comput., 2016, 12 (8), 3781), and many open-questions regarding the performance and achievable accuracy in applying such computations to larger complexes persisted. Key deliverables in the report include the implementation of the numerical Hessian and subsequent vibrational analysis in the BAGEL program package. This includes both full Hessian vibrational analysis (FHVA) and a partial Hessian vibrational analysis (PHVA). The latter of which allows us to reduce the significant cost of computing the vibrational frequencies in a large molecule by focusing on the modes of highest interest. The gradient code also proved useful in the development of an approach to improve the Hubbard-U correction by removing bias in predicting spin-splitting in Fe(II) complexes via plane-wave density functional theory (DFT). Finally, we showed for three families of complexes (spin-crossover (SCO) complexes, metallocorroles complexes, and those with metal-metal bonds) that CASPT2 can result in good molecular geometries and established best practices in undertaking this work. We are now using this approach in collaborative projects with experimental groups, which would not have been possible at the start of this project. A common theme also arose through this work that static and dynamic correlation must be recovered in a balanced way to yield quantitative results. By systematically varying how both effects were included, we were able to make key chemical insights. This work supported the training of two postdoctoral scholars, one graduate student, and one undergraduate student in applying high-level wavefunction based methods to challenging systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗