Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantitative structure-property relationship”

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 19 records

Prediction of Dielectric Constant in Series of Polymers by Quantitative Structure-Property Relationship (QSPR)

This work is devoted to the investigation of dielectric permittivity which is influenced by electronic, ionic, and dipolar polarization mechanisms, contributing to the material’s capacity to store electrical energy. In this study, an extended dataset of 86 polymers was analyzed, and two quantitative structure–property relationship (QSPR) models were developed to predict dielectric permittivity. From an initial set of 1273 descriptors, the most relevant ones were selected using a genetic algorithm, and machine learning models were built using the Gradient Boosting Regressor (GBR). In contrast to Multiple Linear Regression (MLR)- and Partial Least Squares (PLS)-based models, the gradient boosting models excel in handling nonlinear relationships and multicollinearity, iteratively optimizing decision trees to improve accuracy without overfitting. The developed GBR models showed high R2 coefficients of 0.938 and 0.822, for the training and test sets, respectively. An Accumulated Local Effect (ALE) technique was applied to assess the relationship between the selected descriptors—eight for the GB_A model and six for the GB_B model, and their impact on target property. ALE analysis revealed that descriptors such as TDB09m had a strong positive effect on permittivity, while MLOGP2 showed a negative effect. These results highlight the effectiveness of the GBR approach in predicting the dielectric properties of polymers, offering improved accuracy and interpretability.

Ascencio-Medina, Estefania↗

Elucidation of Structure and Physical Properties of Pyrolytic Sugar Oligomers Derived from Cellulose Depolymerization/Dehydration Reactions: A Density Functional Theory Study

Fast pyrolysis of lignocellulosic materials is a promising research area to produce renewable fuels and chemicals. Dehydration is known to be among the most important reaction families during cellulose pyrolysis; water is the most important product. Together with water, dehydration reactions also form a range of poorly known oligomer species of varying molecular sizes, often collected as part of bio-oil water-soluble (WS) fraction. In this work, we used electronic structure calculations to evaluate the relative thermodynamic stabilities of several oligomer species resulting from up to three consecutive dehydration events from cellulose depolymerization intermediates. A library of the thermodynamically favored candidate molecular structures was compiled. Results revealed that most of the water molecules are eliminated from the non-reducing end, forming thermodynamically more stable conjugated compounds. This is consistent with results reported by other researchers in literature where dehydration reactions occur preferably at the non-reducing ends of oligomers. The physical-chemical properties of the proposed structures were estimated using quantitative structure-property relationships (QSPRs) and quantitative property-property relationships (QPPRs). The anhydro-sugars derived from cellulose are often blamed for coke formation during bio-oil hydrotreatment. Understanding their chemical structure could help to develop rational strategies to mitigate coke formation. Furthermore, the thermo-physical properties reported (boiling point, melting point, Gibb’s free energy of formation, enthalpy of formation, and solubility parameters among others) are also fundamental to conducting first principle engineering calculations to design and analyze new pyrolysis reactors and bio-oil up-grading units.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Artificial Neural Network Models for Octane Number and Octane Sensitivity: A Quantitative Structure Property Relationship Approach to Fuel Design

Octane sensitivity (OS), defined as the research octane number (RON) minus the motor octane number (MON) of a fuel, has gained interest among researchers due to its effect on knocking conditions in internal combustion engines. Compounds with a high OS enable higher efficiencies, especially within advanced compression ignition engines. RON/MON must be experimentally tested to determine OS, requiring time, funding, and specialized equipment. Thus, predictive models trained with existing experimental data and molecular descriptors (via quantitative structure-property relationships (QSPRs)) would allow for the preemptive screening of compounds prior to performing these experiments. Here, the present work proposes two methods for predicting the OS of a given compound: using artificial neural networks (ANNs) trained with QSPR descriptors to predict RON and MON individually to compute OS (derived octane sensitivity (dOS)), and using ANNs trained with QSPR descriptors to directly predict OS. Twenty-five ANNs were trained for both RON and MON and their test sets achieved an overall 6.4% and 5.2% error, respectively. Twenty-five additional ANNs were trained for both dOS and OS; dOS calculations were found to have 15.3% error while predicting OS directly resulted in 9.9% error. A chemical analysis of the top QSPR descriptors for RON/MON and OS is conducted, highlighting desirable structural features for high-performing molecules and offering insight into the inner mathematical workings of ANNs; such chemical interpretations study the interconnections between structural features, descriptors, and fuel performance showing that connectivity, structural diversity, and atomic hybridization consistently drive fuel performance.

09 BIOMASS FUELS↗

Sustainable materials acceleration platform reveals stable and efficient wide-bandgap metal halide perovskite alloys

The vast chemical space of emerging semiconductors, like metal halide perovskites, and their varied requirements for semiconductor applications have rendered trial-and-error environmentally unsustainable. Here, in this work, we demonstrate RoboMapper, a materials acceleration platform (MAP), that achieves 10-fold research acceleration by formulating and palletizing semiconductors on a chip, thereby allowing high-throughput (HT) measurements to generate quantitative structure-property relationships (QSPRs) considerably more efficiently and sustainably. We leverage the RoboMapper to construct QSPR maps for the mixed ion FA 1-y Cs y Pb(I 1-x Br x ) 3 halide perovskite in terms of structure, bandgap, and photostability with respect to its composition. We identify wide-bandgap alloys suitable for perovskite-Si hybrid tandem solar cells exhibiting a pure cubic perovskite phase with favorable defect chemistry while achieving superior stability at the target bandgap of ~1.7 eV. RoboMapper’s palletization strategy reduces environmental impacts of data generation in materials research by more than an order of magnitude, paving the way for sustainable data-driven materials research.

36 MATERIALS SCIENCE↗

Development of efficient solar cells using combination of QSPR and DFT approaches

The awarded research project aimed to design and determine the effective organic dye-sensitizers from the ensemble of chromophores for the dye-sensitized solar cell (DSSC) and fullerene-derivatives (FDs) as acceptor for polymer solar cell (PSC) followed by explore the electron transfer mechanism using the combined quantitative structure-property relationship (QSPR) analysis in conjunction with density functional theory (DFT) and time-dependent DFT (TDDFT)-based calculations. Based on our proposal we have employed different in silico approaches to design photo-efficient organic dye-sensitizers for DSSCs and FDs as acceptor for PSCs with exhaustive screening methodology.

14 SOLAR ENERGY↗

DOE STTR Phase I Final Report Report: Machine-learning Based Prediction of Thermal Limits for Conjugated Organic Materials

The SCANN-DT Phase I STTR project led by NLM Photonics and partnering with the National Renewable Energy Laboratory (NREL) sought to apply machine learning techniques based on graph neural networks (GNNs) towards the prediction of decomposition temperatures (Td) of organic semiconductors, based on prior work on bond dissociation energy (BDE) prediction as implemented in NREL’s ALFABET prediction tool. Using a curated set of experimental decomposition energies, the project examined GNN-based, classical quantitative structure-property relationship (QSPR) based on DFT calculations, and combinations of both methods to predict Td. While the best MAEs in Td achieved were near 40°C, below project targets, the project led to improvements in the ALFABET model, improvements in cloud-based implementations of NWChem software, and a substantial dataset of calculations on medium-sized conjugated organic molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evidence That Sharp Interfaces Suppress Recombination in Thick Organic Solar Cells

Commercialization and scale-up of organic solar cells (OSCs) using industrial solution printing requires maintaining maximum performance at active-layer thicknesses >400 nm – a characteristic still not generally achieved in non-fullerene acceptor OSCs. NT812:PC71BM is a rare system whose performance increases up to these thicknesses due to highly suppressed charge recombination relative to the classic Langevin model. The suppression in this system, however, uniquely depends on device processing, pointing towards the role of nano-morphology. We investigate the morphological origins of this suppressed recombination by combining results from a suite of X-ray techniques. We are surprised to find that while all investigated devices are composed of pure, similarly aggregated nanodomains, Langevin reduction factors can still be tuned from ~2 to >1000. This indicates that pure aggregated phases are insufficient for non-Langevin (reduced) recombination. Instead, we find that large well-ordered conduits and, in particular, sharp interfaces between domains appear to help to keep opposite charges separated and percolation pathways clear for enhanced charge collection in thick active layers. To our knowledge, this is the first quantitative study to isolate donor/acceptor interfacial width correlated to non-Langevin charge recombination. Furthermore, this new structure-property relationship will be key to successful commercialization of printed OSCs at scale.

14 SOLAR ENERGY↗

Atomic-Scale Surface Studies of Bulk Metallic Glasses. Final Report

Bulk metallic glasses (BMGs) are of both scientific and technological interest because the absence of periodic atomic arrangements provides them with unique physical, chemical, and mechanical properties. Their high strength, superior elasticity, and an ability to be easily formed into virtually unlimited shapes with feature sizes from centimeters to Angstroms by blow molding and thermoplastic forming makes them an attractive choice for more and more practical applications and products. Due to their complex internal structure, however, experiments that yield insight into their exact atomic arrangements have been scarce. As a result, glass physics is one of the last remaining unexplored fields of materials science despite the scientific and technological importance of glasses in our daily lives, and the question how to characterize and control matter away from equilibrium, as glasses are, was listed as one of five Grand Challenges in a recent DoE report. The main reason for the slow progress in glass physics is the lack of experimental tools that enable access to atomic-scale structural and behavioral information for disordered materials. Such atomic-scale knowledge is mandatory to establish structure-property relationships that ultimately could allow to custom-design alloys featuring specific desired characteristics. With no such relationships available, theory development in glass remains basic and the few that exist are often untested. The aim of this research was to enable progress in our understanding of BMGs by developing a new approach that will allow a meaningful application of local surface science methods to specially prepared BMG samples to obtain a wealth of quantitative information on their atomic arrangements. Key was the availability of specially prepared samples whose surfaces feature large atomically flat terraces despite being entirely amorphous, which we have produced from a Pt 57.5 Cu 14.7 Ni 5.3 P 22.5 alloy (‘Pt-BMG’) both under ambient conditions as well as in ultrahigh vacuum using a unique setup that has been specially developed within this grant. Our approach starts with the in-situ preparation of oxide crystals that are terminated by large terraces, from which exact mirror images out of BMG will be produced using thermoplastic forming (TPF); for the research within this grant, we have successfully used (001)-oriented SrTiO 3 single crystal surfaces as well as (100)-, (110)-, and (111)-oriented single crystals made from LaAlO 3 . Since the resulting BMG replicas display all features of the original crystal with sub-Angstrom fidelity, thereby mimicking the original crystal’s termination by atomically flat terraces without being crystalline themselves, they are ideally suited for further investigation. The following atomic-scale local studies were then carried within this proposal: (i) high-resolution surface imaging and local spectroscopy using scanning probe microscopy, which showed disordered atom-like features and revealed changes in the gradient of the local surface potential on a 1-2 nm length scale; (ii) characterization of atomic-scale plastic flow with affected volumes as low as 1000 atoms, which showed local hardness near or above the theoretically predicted maximum and, once plastic deformation was initiated, homogeneous flow of the atoms involved; (iii) characterization of surface relaxation processes and the onset of crystallization induced by annealing, which showed that upon heating over the material’s glass transition temperature, the surface rearranges and relaxes towards a more stable, denser packed glass, which increases on-terrace surface roughness, while surface tension smoothens step edges; and (iv) studies that investigate the dependence of the material’s mechanical properties and structure on processing parameters, revealing that relaxed glasses get denser, harder, and more elastic. In combination, this information allows to combine structural models with mechanical properties and preparation history, thereby facilitating the development of preparation-structure-property relationships for metallic glasses. With the availability of such information, bulk metallic glasses can be further optimized to be used in more and more applications in industry.

36 MATERIALS SCIENCE↗

Polarized Resonant X-ray Scattering to Measure Molecular Orientation and Conformation in Organic Nanostructures

The ordering of organic (carbon-based) molecules into nanostructures drives useful materials properties for applications such as organic light emitting diodes, flexible batteries, and polymer nanocarriers. However, there is a current lack of techniques to measure this ordering. X-ray scattering techniques involving molecular resonances have shown promise in identifying molecular types, locations, and orientation. However, resonant soft X-ray scattering (RSoXS) techniques have not been developed even lacking basic optical models. The objectives achieved by this research were the development of methods, models, and analysis techniques that harness these resonant phenomena in X-ray scattering measurements to quantitatively probe ordering of organic molecules within devices and nanostructures. This award resulted in 11 peer reviewed publications with two others soon to be submitted. New physical models for X-ray scattering that accurately describe these complex molecular resonances enable measurements of how molecules assemble, orient and conform at the nanometer scale. These new techniques were applied to both organic electronic devices and aqueous nanoparticles being developed for hydrocarbon sequestration. This new measurement capability will enable characterization of structure-property relationships that lead to disruptive technologies of printable, stretchable, and biocompatible devices based on organic materials.

14 SOLAR ENERGY↗

Bridging the length scales on mechanical property evaluation (Final Report)

The development of small-scale mechanical testing in combination with microstructural investigation is of great interest to the nuclear materials community for both materials development and monitoring applications. Dramatically reducing the sample sizes to reduce radioactivity and obtaining mechanical properties of irradiated samples is truly intriguing. Moreover, such studies promise a range of benefits including cost reduction, fundamental insight in structure-property relationships, increased statistics on less sample material, and reinvestigation of prior irradiated and tested reactor samples while simultaneously enabling the generation of mechanical test data on ion beam irradiated materials with limited penetration depths. Small scale materials testing on sub-sized samples has been studied for several decades, though it has only been after the development of micro-testing based on Focused Ion Beam (FIB) sample manufacturing in that orders of magnitude smaller samples could really be investigated in a quantitative manner. In recent years, small scale mechanical testing techniques at a number of length scales has been developed for both unirradiated and irradiated (ion and neutron) materials. Technological advances made in this field have enabled ex-situ and in-situ transmission electron microscopy (TEM) and scanning electron microscopy (SEM) examination, thus leading to more accurate measurements as well as additional mechanistic information. A recent review of the benefits of these techniques show that these techniques are at a stage to tackle multi-scale ranges of materials investigations and can be utilized to obtain fundamental science-based understanding of nuclear materials. Considering the tremendous advances made, one can see how small-scale mechanical testing techniques combined with modeling can enable true small scale to bulk scale mechanical property correlations. However, for the engineering community to adapt this approach fully, one needs to demonstrate that a) that these techniques can produce results with high fidelity and reproducibility, b) generate engineering stress-strain data that one can utilize to understand bulk behavior, and c) generate new insight into relevant phenomena fostering the true understanding of radiation damage and microstructure in materials for nuclear applications. It is the objective of this proposal to bridge the length scale between macro- and micro- scale mechanical testing of unirradiated and irradiated materials. This involves the development and demonstration of procedures for multi-scale mechanical testing that enable high fidelity reproducibility of data and the generation bulk property data from small-scale mechanical tests. Through this, the proposal aims to enhance the confidence in the obtained data at the smaller length scales and enhance the insight provided from these techniques for bulk scale applications on both unirradiated and irradiated materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Stochastic 3D reconstruction of cracked polycrystalline NMC particles using 2D SEM data

Li-ion battery performance is strongly influenced by the 3D microstructure of its cathode particles. Cracks within these particles develop during calendaring and cycling, reducing connectivity but increasing reactive surface, making their impact on battery performance complex. Understanding these contradictory effects requires a quantitative link between particle morphology and battery performance. However, informative 3D imaging techniques are time-consuming, costly and rarely available, such that analyses often have to rely on 2D image data. This paper presents a novel stereological approach for generating virtual 3D cathode particles exhibiting crack networks that are statistically equivalent to those observed in 2D sections of experimentally measured particles. Consequently, 2D image data suffices for deriving a full 3D characterization of cracked cathodes particles. Such virtually generated 3D particles could serve as geometry input for spatially resolved electro-chemo-mechanical simulations to enhance our understanding of structure-property relationships of cathodes in Li-ion batteries.

36 MATERIALS SCIENCE↗

Comparing Machine Learning and Physics-Based Nanoparticle Geometry Determinations Using Far-Field Spectral Properties

Anisotropic metal nanostructures exhibit polarization-dependent light scattering, a property which has been widely studied and exploited to determine orientations of subwavelength structures using far-field microscopy. Here we explore the use of variational autoencoders (VAEs) to determine the geometries of gold nanorods (NRs) such as in-plane orientation and aspect ratio under linearly polarized dark-field illumination in an optical microscope. We enforce a shared latent space to connect two VAEs trained separately with polarized dark-field scattering spectra and electron microscopy images and achieve image prediction (shape, orientation, and size) of Au NRs using only polarized dark-field scattering spectra. We determine the geometrical parameters of orientational angle and aspect ratio quantitatively via both our dual-VAE and physics-based analysis on the input scattering spectra. We show that orientational angle prediction by dual-VAE performs well with only a small (~300 particle) training set, yielding a mean absolute error (MAE) of 14.4° and a concordance correlation coefficient (CCC) of 0.95. This performance is only marginally worse than the physics-based cos(2?) fitting approach between the scattering intensity and the polarizing angle, which achieves MAE of 8.78° and CCC of 0.99. Aspect ratio determination is also comparable for the dual-VAE and physics-based fitting comparison (MAE of 0.21 vs. 0.23 and CCC of 0.53 vs. 0.68). Here, this dual encoder-decoder architecture effectively exploits the structure-property relationships of plasmonic nanostructures to construct a cross-modal machine learning (ML) approach, providing a pathway to employ ML approaches to address other structure-property relationships in materials science.

Dark-field scattering↗

Structure–property relationships of reduced graphene oxide membranes intercalated with polycyclic aromatics

Graphene oxide (GO) membranes intercalated with various organic moieties have shown excellent potential for a range of water processing applications. However, microstructure–functional property relationships in these structurally disordered membranes are not well understood. We demonstrate a practical methodology for developing such relationships for GO membranes intercalated with molecular species, with polycyclic aromatic toluidine blue O (TBO) as an example functional intercalant. We use solid-state UV–vis absorbance and fluorescence measurements to quantitatively track the arrangements of TBO in a series of TBO-loaded reduced GO (rGO) membranes. This study reveals the evolution of diverse arrangements including TBO monomers, lateral and stacked dimers, and other aggregates, as a function of overall TBO loading. These microstructures are then correlated to changes in overall properties such as interlayer d-spacings, permeate fluxes, and solute rejections. The characterization of these different intercalant microstructures explains non-intuitive flux and rejection trends, which can circumvent flux and solute rejection trade-offs.

02 PETROLEUM↗

Functional group effect of chemically modified microcrystalline methyl cellulose on thermoplastic polyurethane composites

Microcrystalline cellulose (MCC) is a promising bio-based filler for lightweight yet mechanically high-performance eco-polymer composites because of its low density and high biocompatibility. However, intermolecular hydrogen bonding among MCC is stronger than the polymer–filler interactions, which deteriorate the mechanical properties of the composites. Herein, we investigated the effects of functional groups on the mechanical properties of composites by scrutinizing chemically modified microcrystalline methyl cellulose (m-MMC) with three different substitution levels of hydroxyl group to hydroxypropyl (HP) group: no-, low-, and high-level substitution (HP-0, HP-low, and HP-high). The degree of HP substitution of m-MMC was quantitatively measured by CP/MAS NMR analysis. The relatively bulky HP groups interrupted the filler–filler intermolecular interactions and reduced the crystallinity and density of m-MMC, as evident from X-ray diffractometer and pycnometer data, respectively. For scalable production of the composites, the m-MMC were compounded with thermoplastic polyurethane (TPU) by a twin-screw extruder at concentrations between 0.5 and 10 wt%. Despite its low filler concentration, the toughness of m-MMC/TPU composites was remarkably enhanced, up to 28% (229.2 to 294.4 MJ/m 3 ) at 0.5 wt% loading of HP-low, owing to the enhanced polymer–filler interactions. In conclusion, the fundamental understanding on structure–property relationships will provide insights for designing of mechanically robust yet eco-friendly polymer composites.

composites↗

High throughput structure–property relationship for additively manufactured 316L/IN625 alloy mixtures leveraging 2-step Bayesian estimation

While the fabrication of graded materials by directed energy deposition (DED) has led to accelerated materials discovery, the ability to rapidly explore sufficiently large material composition spaces is limited due to the time-intensive nature of conventional materials characterization techniques. The present study investigates the viability of small punch test (SPT) protocols for rapidly evaluating DED-fabricated alloy mixtures of stainless steel 316L (316L) and Inconel 625 (IN625). The SPT protocols evaluated in this study include both the recently established two-step Bayesian estimation framework as well as the empirical relationships established in prior literature. It is shown that these protocols are capable of reliably and quantitatively tracking the changes in the mechanical properties of the alloy mixtures studied. Enhancement of mechanical properties was observed with the addition of IN625 to 316L, which is attributed to the austenite stabilization in the matrix and the formation of fine δ - Ni3Nb precipitates. It is shown that CALPHAD-based Scheil model simulations predicted the formation of different precipitate phases for each composition. The novel protocols presented in this paper open new avenues for high throughput material explorations for additive manufacturing.

2-step Bayesian estimation↗

Exciton engineering of 2D Ruddlesden–Popper perovskites by synergistically tuning the intra and interlayer structures

Designing two-dimensional halide perovskites for high-performance optoelectronic applications requires deep understanding of the structure-property relationship that governs their excitonic behaviors. However, a design framework that considers both intra and interlayer structures modified by the A-site and spacer cations, respectively, has not been developed. Here, we use pressure to synergistically tune the intra and interlayer structures and uncover the structural modulations that result in improved optoelectronic performance. Under applied pressure, (BA) 2 (GA)Pb 2 I 7 exhibits a 72-fold boost of photoluminescence and 10-fold increase of photoconductivity. Based on the observed structural change, we introduce a structural descriptor χ that describes both the intra and interlayer characteristics and establish a general quantitative relationship between χ and photoluminescence quantum yield: smaller χ correlates with minimized trapped excitons and more efficient emission from free excitons. Building on this principle, we design a perovskite (CMA) 2 (FA)Pb 2 I 7 that exhibits a small χ and an impressive photoluminescence quantum yield of 59.3%.

36 MATERIALS SCIENCE↗

Structure–property relations of sodium iron phosphate nuclear waste glasses: Effects of iron redox ratio and glass composition

Iron phosphate glasses, known for their exceptional chemical durability and potential applicability in nuclear waste management, have gained significant attention over the years. The structures of these glasses are complicated by the coexistence of Fe 3+ and Fe 2+ , which plays a crucial role in determining their structures and properties. Here, this work uses molecular dynamics simulations to study the structural changes in Na 2 O–Fe 2 O 3 –P 2 O 5 glasses with varying glass composition and Fe 2+ /Fe 3+ redox ratio. It was found that the redox ratio and modifier contents significantly affected the short-range and medium-range orders in the glasses. Significant changes in the local environments around P 5+ and Fe 3+ were observed, as reflected by the bond distances and coordination numbers. Na + cations are found to preferentially associate with Fe 3+ (rather than Fe 2+ ), whereas Fe 2+ has stronger association with P 5+ than Na + , confirming the structural role of Fe 2+ as a glass modifier. The disruptions in P–O–P linkages upon increasing FeO suggest that FeO causes glass depolymerization. These glasses achieved higher connectivity with increasing Fe 3+ / (Fe 3+ + Fe 2+ ) ratios, conerting phosphorous Q 2 to Q 3 units and iron Q 5 units to Q 4 units. The decrease of nonbridging oxygen fractions with increasing Fe 3+ / (Fe 3+ + Fe 2+ ) ratios, through creating P–O–Fe linkages, is the main reason of enhanced network connectivity. Quantitative structure–property relationship analyses with different structural descriptors were used to correlate with measured properties. The analyses provided valuable insights into structure–property relationships, emphasizing the importance of choosing relevant energy parameters and defining glass network connectivity, particularly in F net descriptors. It was found the Fe–O–P linkage density exhibits strong correlations to measured dissolution rates, supporting the importance of these linkages in improving the chemical durability in iron phosphate glasses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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↗