Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GaP”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Overcoming barriers to improved decision-making for battery deployment in the clean energy transition

Decarbonization plans depend on the rapid, large-scale deployment of batteries to sufficiently decarbonize the electricity system and on-road transport. This can take many forms, shaped by technology, materials, and supply chain selection, which will have local and global environmental and social impacts. Current knowledge gaps limit the ability of decision-makers to make choices in facilitating battery deployment that minimizes or avoids unintended environmental and social consequences. These gaps include a lack of harmonized, accessible, and up-to-date data on manufacturing and supply chains and shortcomings within sustainability and social impact assessment methods, resulting in uncertainty that limits incorporation of research into policy making. These gaps can lead to unintended detrimental effects of large-scale battery deployment. To support decarbonization goals while minimizing negative environmental and social impacts, we elucidate current barriers to tracking how decision-making for large-scale battery deployment translates to environmental and social impacts and recommend steps to overcome them.

25 ENERGY STORAGE↗
DOE OSTI2025

Unveiling the nature of Ga-based chalcogenides for electrical switching selectors

Three-dimensional phase-change memory with stackable crossbar architecture is a promising technology to meet the urgent demands for high-density storage and rapid information processing in the era of explosive data growth. The performance depends strongly on the properties of ovonic threshold switching (OTS) selectors, which control the on/off states of memory units. Amorphous GaS serves as an outstanding OTS material, distinguished by its sizable mobility gap and high crystallization temperature, while the underlying mechanism continues to be inadequately comprehended. Here, in this work, we systematically studied the structural and electronic properties of amorphous Ga-X (X = S/Se/Te) using first-principles calculations. The results show that Ga atoms adopt tetrahedral motifs, while S/Se/Te atoms predominantly exhibit the structure of a distorted triangular pyramid. This structural arrangement is ascribed to the substantial dative bonds formed by the lone-pair electrons of the anions and the vacant sp3 orbitals around Ga atoms. Large mobility gaps (e.g., GaS: 2.43 eV, GaSe: 1.76 eV, GaTe: 1.26 eV) and distinct mid-gap states (e.g., ∼0.66 eV above valence band tail) ensure that these three chalcogenide glasses can be switched on under an external electric field while effectively suppressing leakage current without a bias, and the defect electronic states originate from short, robust Ga-Ga bonds due to the formation of distorted chain-like local structures. Our research elucidates the mechanisms of amorphous Ga-X as OTS materials, enriching the spectrum of electrical switching selectors by incorporating III-VI chalcogenides. This inclusion offers novel opportunities for the refinement and optimization of high-density integrated memory systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Correlation of Early-Stage Growth Process Conditions with Dislocation Evolution in MOCVD-Based GaP/Si Heteroepitaxy

To identify the complex relationships between early-stage growth processes and the resultant defect microstructure in GaP/Si heteroepitaxy, a holistic study of several key metal-organic chemical vapor deposition (MOCVD) parameters was conducted, focusing on Si surface preparation and GaP atomic layer epitaxy (ALE) based nucleation processes. Here, crystalline defects related to the lattice mismatch and/or interfacial heterovalency, namely misfit dislocations (MD), threading dislocations (TD), and stacking fault pyramids (SFP), were quantitatively characterized via electron channeling contrast imaging (ECCI) and correlated against the different process variations. Choice of Si surface preparation method between the two examined (dilute SiH4 annealing versus Si2H6 based homoepitaxy) had little impact on resultant GaP film morphology and defect content, whereas differing GaP ALE nucleation conditions produced much more substantial changes. In particular, the initial precursor species (tert-butylphosphine versus triethygallium) and ALE cycle purge times both yielded significant influence over threading dislocation densities (TDD) in thin (100 nm), post-critical thickness GaP/Si films, with TDD spanning two orders of magnitude, from 6.7×10 7 cm -2 to 7.1×10 5 cm -2 , depending on the specific process conditions employed. SFP densities were also found to follow a similar trend, ranging from 2.0×10 7 cm -2 to 1.8×10 5 cm -2 , but with no apparent causal relationship between SFP density and TDD. To help explain the dramatic differences observed, detailed, large-area MD network characterization was used to provide statistically-relevant quantitative analyses of the critical dislocation dynamics (introduction rates and glide velocities) associated with the different process variants. These extracted values are then correlated against the ALE process variants to provide insight into the potential mechanistic roles of the different growth processes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

New insights on the effects of chemistry and temperature on α’ precipitation during aging of FeCrAl alloys

Here, the effect of Al on α' precipitation in FeCrAl alloys was studied through thermal aging of several model binary FeCr and FeCrAl alloys with Cr content of 13 wt.%, 17 wt.%, and 25 wt.% with and without 5.5 wt.% Al. Aging was performed up to 1,000h at varying temperatures from 400-500°C. At 400°C no age hardening was observed due to slow kinetics at this temperature. For both the 17Cr and 25Cr alloys, the addition of Al shows a lowering of the miscibility gap, consistent with other reports in literature. Interestingly, however, for the 25Cr alloys the addition of Al in the FeCrAl ternary alloy accelerated α' precipitation below the miscibility gap. Such enhanced precipitation has also been predicted by our atomistic kinetic Monte Carlo (AKMC) simulations. While previous literature has often focused on Al suppressing α', here we show that while Al can lower the miscibility gap in the Fe-Cr ferritic system thermodynamically, it may also enhance the kinetics of precipitation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nanochannel electrodes facilitating interfacial transport for PEM water electrolysis

Proton-exchange membrane water electrolyzers (PEMWEs) are a promising technology for green hydrogen production; however, interfacial transport behaviors are poorly understood, hindering device performance and longevity. Here, we first utilized finite-gap electrolyzer to demonstrate the possibility of proton transfer through water in PEMWEs. The measured high-frequency resistances (HFRs) exhibit a linear trend with increasing gap distance, where extrapolation shows a lower value compared with HFRs in regular zero-gap electrolyzers, indicating that ohmic resistance could be further reduced. We introduce nanochannels to facilitate mass transport, as evidenced by both liquid-fed and vapor-fed electrolysis. Nanochannel electrodes achieve a voltage reduction of 190 mV at 9 A·cm –2 compared with the Ir-PTEs without nanochannels. Furthermore, nanochannel electrodes show negligible degradation through 100,000 accelerated-stress tests and over 2,000 h of operation at 1.8 A·cm –2 with a decay rate of 11.66 μV·h –1 . These results provide new insights into localized transport dynamics for PEMWEs and highlight the significance of interfacial engineering for electrochemical devices.

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

Influence of grain size on α' Cr precipitation in an isothermally aged Fe-21Cr-5Al alloy

Cr-rich α' precipitation during aging typically leads to hardening and accordingly embrittlement of FeCrAl alloys, which needs to be suppressed. The influence of grain size on α' precipitation was studied by aging coarse-grained (CG), ultra-fine grained (UFG), and nanocrystalline (NC) ferritic Kanthal-D [KD; Fe-21Cr-5Al (wt.%) alloy] at 450, 500 and 550 °C for 500 h. After aging at 450 and 500 °C, less hardening was observed in the UFG KD than in CG KD. Atom probe tomography indicated a lower number density and larger sized intragranular α' in the UFG versus the CG alloy. The smaller grain size and higher defect (vacancy and dislocation) density in the UFG KD facilitated diffusion and accordingly enhanced precipitation kinetics, leading to coarsening of precipitates, as well as saturation of precipitation at lower temperatures, as compared to those in CG KD. No hardening occurred in UFG and CG KD after aging at 550 °C, indicating that the miscibility gap is between 500 and 550 °C. NC KD exhibited softening after aging owing to grain growth. α' precipitation occurred in NC KD aged at 450 °C but not at 500 °C, indicating that miscibility gap is between 450 and 500 °C. Thus, the significantly smaller grain size in NC KD decreased the miscibility gap, as compared to that in CG and UFG KD. Finally, this is attributed to the absorption of vacancies by migrating grain boundaries during aging, suppressing α' nucleation and enhancing Cr solubility.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Efficacy of vibro-impact energy harvesting absorbers on controlling dynamical systems under vortex-induced vibrations and base excitation

Here this study investigates improving the efficacy of an energy harvesting absorber's ability to control a structure under vortex-induced vibrations, base excitation, and a combination of the two by including mechanical amplitude stoppers. The nonlinear reduced-order model is developed through modifying trilinear spring models to represent the impact forces, a modified van der Pol oscillator to represent the forcing due to the vortex-induced vibrations and using the Euler-Lagrange principle to express the equations of motion. It is seen that a soft stopper stiffness and a 5mm gap performs the most effectively of increasing the power generated from the absorber while still greatly reducing the primary structure's amplitude. By changing the stopper's location towards the middle of the energy harvesting absorber, the large effects of the impact forces are reduced and improves the efficacy of medium and hard stopper stiffnesses to generate near the amount of power the soft stopper does, while greatly improving the control of the primary structure. When the system is under combined loadings, the large oscillations of the synchronization region cause the effective configuration to be that of a 27.5 mm gap with soft stiffnesses. The results shows that medium stiffness stoppers with small gaps generate large aperiodic regions due to the high impact force. When the oscillations are close to the stoppers, the beating phenomenon is observed and is not overpowered by the vibro-impact force.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Passive heat removal in horizontally oriented micro-HTGRs

There are novel, horizontally oriented high-temperature gas-cooled micro-reactor (HTGR) designs under consideration because of their portability, simplicity, and reliable operation. These features of the micro-HTGR design can meet the fission battery attributes, such as economic, standardize, installed, unattended and reliable. Although there have been HTGR related safety studies, micro-HTGRs when horizontally oriented are expected to exhibit significantly different thermal physics. During the unavailability of helium circulation, the internal reactor core is designed to cool by block-to-block conduction and radiation, and the reactor vessel surface is cooled by the ambient air. This scenario is anticipated during the transport of the micro-HTGR in a shipping container. The conduction and radiation between the prismatic micro-HTGR blocks in the core can be influenced by variances in the thermal contacts. This work investigated the conduction within a simulated horizontal HTGR core. An experimental setup was used to validate a numerical model. The experimental setup consisted of a hexagonal assembly with scaled prismatic blocks placed within a high-temperature vacuum environment. The gaps between the blocks were well controlled and monitored. The experimental setup was designed in such a way that the temperature variation in the axial direction was minimal, such that the experiment could be observed as a 2D (r, ) heat transfer problem. The pressurized tube was equipped with Infra-Red (IR) transparent windows, and an IR camera was used to capture the spatio-temporal temperature evolution in the hexagonal block assembly during the cooldown experiments. The experimental scenario was computationally modeled with a finite element analysis (FEA) program. Once validated, the computational model was used to investigate the impact of gap conductance on overall decay heat removal. Using a conservative estimate for gap conductance value (100 W/m 2 -K) between the prismatic blocks, there is a negligible increase in temperature observed during decay heat generation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Heat transfer coefficients of moving particle beds from flow-dependent thermal conductivity and near-wall resistance

Accurate determination of heat transfer coefficients for flowing packed particle beds is essential to the design of particle heat exchangers and other thermal and thermochemical equipment. While such dense granular flows mostly fall into the well-known plug-flow regime, the discrete nature of granular materials alters the thermal transport processes in both the near-wall and bulk regions of flowing particle beds from their stationary counterparts. As a result, heat transfer correlations based on the stationary particle bed thermal conductivity could be inadequate for flowing particles in a heat exchanger. Most earlier works have achieved a reasonable agreement with experiments by treating granular heat transfer media as a plug-flow continuum with a near-wall thermal resistance in series. However, the thermal conductivity values of the continuum were often obtained from measurements on stationary beds owing to the difficulty of flowing bed measurements. In this work, it was found that the properties of a stationary bed are highly sensitive to the method of particle packing and there is a decrease in the particle bed thermal conductivity and increase in the near-wall thermal resistance, measured as an effective air gap thickness, on the onset of particle flow. These variations in thermal conductivity of stationary and flowing particle beds can lead to errors in heat transfer coefficient calculations. Therefore, the heat transfer coefficients for granular flows were calculated using experimentally determined flowing particle bed thermal conductivity and near-wall air gap for ceramic particles – CARBO CP 40/100 (mean diameter = 275 µm), HSP 40/70 (404 µm) and HSP 16/30 (956 µm); at velocities of 5–15 mm·s –1 ; and temperatures of 300–650 °C. The thermal conductivity and air gap values for CP 40/100 and HSP 40/70 were further used to calculate heat transfer coefficients across different particle bed temperatures and velocities for different parallel-plate heat exchanger dimensions. Furthermore, these calculations, which show good agreement with measured HTC values reported in literature, can be used as a guide for heat exchanger designs. Graphical abstract

14 SOLAR ENERGY↗
DOE OSTI2023

Modeling and design of III-V heterojunction solar cells for enhanced performance

Heterojunctions can increase the efficiency of solar cell devices relative to homojunctions, but there is a large parameter space with significant tradeoffs that must be considered. Here, we present an experimental and computational study of III-V heterojunction solar cells and show how the emitter doping, emitter band gap, and heteroband offsets impact device efficiency. Efficiency is maximized by pushing the junction depletion region into the wider band gap material while minimizing the effects of heteroband offsets through optimized choice of emitter band gap, emitter electron affinity, and/or emitter doping density. We use these results to guide optimization of devices grown by halide vapor phase epitaxy, achieving 27% efficiency in a GaAs/GaInPAs heterojunction device. We also show that heterojunctions yield proportionally larger efficiency improvements in lower-quality materials. Although the modeling was developed and validated using III-V materials, the results are theoretically applicable to materials systems outside III-Vs.

14 SOLAR ENERGY↗
DOE OSTI2026

Lattice-Mismatched van der Waals Epitaxy and Photoluminescence of Two-Dimensional GaxIn1-xSe Alloys on Si(111)

GaxIn1-xSe (GIS) alloys are two-dimensional (2D) layered materials with band gaps and lattice parameters of interest for many energy and electronic applications. They can be fabricated using van der Waals epitaxy, which is an emerging technique that offers unprecedented opportunities for 2D optoelectronic devices and epitaxy processes. This work has demonstrated van der Waals epitaxy of GIS alloys for the first time. Films with x compositions of 0, 0.062, 0.164, 0.680, 0.894, and 1 and tunable lattice constants were grown on Si(111) substrates, and characterized by X-ray diffraction pole figure and transmission electron microscope analysis. In spite of the lattice mismatches (InSe is 4.1% too large and GaSe is 2.8% too small), these alloys grow epitaxially, with Si(111) || GIS(001) and Si[1-10] || GIS[100] orientation. Photoluminescence was used to measure tunable band gaps in the absorber-relevant 1.3-2.0 eV range as a function of x composition and showed GIS did not degrade after capping with Se and prolonged storage. Therefore, GIS alloys exhibit a technologically advantageous combination of tunable band gap and photoluminescence with relaxed lattice parameter and rotational registry with the substrate.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ligand-Dependent Tuning of Interband and Intersubband Transitions of Colloidal CdSe Nanoplatelets

Although surface ligands of colloidal nanocrystals are known to adjust the absolute energy levels of valence and conduction bands of semiconductor nanocrystals, they typically have only minor influence on the band gap or effective masses. This changes in nanoplatelets. Ligand exchange of CdSe colloidal nanoplatelets induces large (up to 300 meV) bathochromic shifts of both interband and intersubband transitions. Here, three families of ligands-halides, thiolates, and phosphonates-are used to tune interband transitions, reflecting electron and hole confinement, across visible wavelengths and intersubband transitions, reflecting electron confinement, across the near-infrared spectral window. Furthermore, careful examination shows that delocalization from expansion of the nanoplatelet short axis, which was reported previously, cannot alone explain observed red shifts. Instead, comparison of intersubband, interband, and hole energy levels shows that ligand head group chemistry confers specific, idiosyncratic adjustments of the contribution of conduction and valence bands to the observed bathochromic shifts. Phosphonate ligands show the largest band gap reductions but the smallest red shift of intersubband transition energies; halide-exchanged samples displayed smaller reductions in band gap but large red shifts of intersubband transitions; thiolates fall in between. A related specificity is observed in hole states, which implicates ligand-responsive valence band curvature as an additional contribution driving optical changes. For nanoplatelets, surface ligand chemistry offers not only a tool to adjust the absolute energy level of conduction and valence bands but also an alternative route to preferential electron or hole band engineering that is normally achieved with inorganic shells.

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

Local Structural Disorder in Metavanadates MV 2 O 6 (M = Zn and Cu) Synthesized by the Deep Eutectic Solvent Route: Photoactive Oxides with Oxygen Vacancies

Metavanadates MV 2 O 6–δ (M = Zn and Cu) are synthesized by using a deep eutectic solvent (DES), a mixture of hydrogen bond donor and acceptor, as a reaction medium. Dissolution of stable binary metal oxides in a DES followed by a heat treatment yields phase-pure vanadates. According to in situ powder X-ray diffraction, ternary phases (α-Zn 2 V 2 O 7 and β-Cu 2–x V 2 O 7 , x ~ 0.27) are intermediates in the reaction pathway taken. Identifying a polymorphic phase transformation temperature for CuV 2 O 6 as well as the narrow temperature range between formation and decomposition for both metavanadates allows for tackling the challenge of the synthesis of these materials. The oxygen vacancy introduced by the DES route is accompanied by the formation of reduced V 4+ and Cu + in the oxide matrix, based on X-ray photoelectron spectroscopy. These oxygen vacancies modify the vibrational modes in the corresponding Raman spectra and are also responsible for broad optical absorptions in the 1.8–1.1 eV range. The optical band gaps of the materials are found at 1.8 eV (CuV 2 O 6 ) and 2.2 eV (ZnV 2 O 6 ), approximately 0.1–0.3 eV below the values reported in the literature. The reduced band gaps and sub-band gap photon absorption are key features of the oxygen-deficient metavanadates. Surface photovoltage spectroscopy reveals that all the synthesized vanadates are n-type materials with electrons as the majority charge carriers, and photoelectrochemical measurements confirm photoanodic currents for methanol oxidation. Furthermore, these results provide insight into the synthesis and structure–property relationship of the metavanadates for their potential use as photoanodes.

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

FeSi 4 P 4 and CoSi 3 P 3 : Hidden Gems of Ternary Tetrel Pnictides with Outstanding Nonlinear Optical Properties

Metal silicon phosphides have shown promise as nonlinear optical materials. To be practically useful and cheap, earth-abundant 3d transition metals are preferred over their scarcer and more expensive 4d and 5d counterparts. Here, we developed a synthetic method to produce polycrystalline bulk powders and millimeter-sized single crystals of ternary compounds FeSi 4 P 4 and CoSi 3 P 3 . Both studied compounds have noncentrosymmetric and chiral crystal structures with ordered Si/P arrangements as was confirmed by single-crystal X-ray diffraction and solid-state NMR. Despite the presence of the transition metal, FeSi 4 P 4 and CoSi 3 P 3 are semiconductors with direct band gaps of 1.3 and 1.6 eV, respectively, indicating low-spin d 6 electronic configuration for octahedral Fe 2+ and Co 3+ . Relative to reported sulfide materials, FeSi 4 P 4 and CoSi 3 P 3 small band gap semiconductors demonstrate an outstanding combination of second-harmonic generation (SHG) activity and laser damage threshold (LDT). Both studied materials are phase-matchable with a 2.09 μm laser and not only exhibit 2.5–3.0 times stronger SHG signal than that of the state-of-the-art AgGaS 2 standard but also demonstrate an LDT response of 2.3–2.5 times higher than that of AgGaS 2 (at 1.09 μm laser with a pulse width of 10 ns)-which is unprecedented for small band gap semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Understanding and Tuning the Electronic Properties of Prussian Blue Analogues

Prussian blue analogues (PBAs) have attracted increasing interest owing to their potential applications in various fields such as energy storage and conversion, neuromorphic computing, and magnetic switching. With a general formula of A x M N [M C (CN) 6 ], they feature an open framework that provides abundant channels for diffusion of alkali metal ions A and allows flexible compositional control of transition metal ions M N and M C . The oxidation states of transition metal ions can be tuned by adjusting the amount (x) of alkali ions A. Here, we carried out density functional theory calculations combined with experimental measurements to investigate the effects of transition metal ions, alkali ions, and oxidation states on the electronic properties of PBAs. Our calculations found that the band gaps of PBAs can be tuned from close to 0 eV to more than 4 eV. Experimentally, we introduced the synthesis/characterization of five previously unreported PBAs (M N = Ru, Os; M C = Fe, Ru, and Os) to complete the nine stable M N :M C transition metal combinations in group VIII of the periodic table. The optically measured intervalence charge transfer excitation energies of group VIII PBAs are consistent with calculated band gaps. They demonstrate wide band gap tunability by adjusting transition metals and oxidation states, enabling semiconductor-to-metal transitions for memristor applications and enhancing electronic conductivity for battery applications. In conclusion, this work provides a computational/experimental database of electronic properties versus structural compositions for PBAs.

Electrical conductivity↗
DOE OSTI2022

Cubic Crystal Structure Formation and Optical Properties within the Ag–B II –M IV –X (B II = Sr, Pb; M IV = Si, Ge, Sn; X = S, Se) Family of Semiconductors

Quaternary chalcogenide semiconductors are promising materials for energy conversion and nonlinear optical applications, with properties tunable primarily by varying the elemental composition and crystal structure. Here, we first analyze the connections among several cubic crystal structure types, as well as the orthorhombic Ag 2 PbGeS 4 -type structure, reported for select members within the Ag–B II –M IV –X (B II = Sr, Pb; M IV = Si, Ge, Sn; X = S, Se) compositional space. Focusing on the Ag–Pb–Si–S and Ag–Sr–Sn–S systems, we show that one structure type, with the formulas Ag 2 Pb 3 Si 2 S 8 and Ag 2 Sr 3 Sn 2 S 8 , is favored. We have prepared powder and single-crystal samples of Ag 2 Pb 3 Si 2 S 8 and Ag 2 Sr 3 Sn 2 S 8 , showing that each takes on the noncentrosymmetric cubic space group I$\bar{43}$d and is isostructural to the previously reported compound Ag 2 Sr 3 Ge 2 Se 8 . Through hybrid density functional theory calculations, these cubic compounds are demonstrated to be (quasi-)direct band gap semiconductors with high densities of states at the band maxima. Furthermore, the band-gap energies are measured by reflectance spectroscopy as 1.95(3) and 2.66(4) eV for Ag 2 Pb 3 Si 2 S 8 and Ag 2 Sr 3 Sn 2 S 8 , respectively. We further measure the optical properties and show the electronic band structures of three other isostructural A I –B II –M IV –X-type materials, i.e., Ag 2 Sr 3 Si 2 S 8 , Ag 2 Sr 3 Ge 2 S 8 , and Ag 2 Sr 3 Ge 2 Se 8 , showing that the band gaps can be predictably tuned by element substitution. Detailed visual analyses of the different structures and of their relationships with other members of the Ag–B II –M IV –X compositional family provide a basis for a broader understanding of the structure formation and optoelectronic properties within the quaternary chalcogenide semiconductor family.

14 SOLAR ENERGY↗
DOE OSTI2022

The Effect of Geometry, Spin, and Orbital Optimization in Achieving Accurate, Correlated Results for Iron–Sulfur Cubanes

Iron–sulfur clusters comprise an important functional motif in the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from correlated, large active space calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. Here, to provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of the [Fe(II/III) 4 S 4 (SMe) 4 ] -2 clusters using two complementary, correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (<1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the singlet and the spin state 2S + 1 = 19, corresponding to a maximal number of Fe-d electrons being unpaired and of parallel spin, is strongly geometry dependent, changing by a factor of 3 upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap compared to typical mean-field orbital approaches. Our results show the need for performing large active space calculations to unveil the challenging electronic structure of these complex catalytic centers and should serve as accurate starting points for fully correlated treatments upon inclusion of dynamical correlation outside the active space.

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

Quantum Vibronic Effects on the Electronic Properties of Molecular Crystals

We present a study of molecular crystals, focused on the effect of nuclear quantum motion and anharmonicity on their electronic properties. We consider a system composed of relatively rigid molecules, a diamondoid crystal, and one composed of floppier molecules, NAI-DMAC, a thermally activated delayed fluorescence compound. We compute fundamental electronic gaps at the density functional theory (DFT) level of theory, with the Perdew–Burke–Erzenhof (PBE) and strongly constrained and approximately normed (SCAN) functionals, by coupling first-principles molecular dynamics with a nuclear quantum thermostat. We find a sizable zero-point renormalization (ZPR) of the band gaps, which is much larger in the case of diamondoids (0.6 eV) than for NAI-DMAC (0.22 eV). We show that the frozen phonon (FP) approximation, which neglects intermolecular anharmonic effects, leads to a large error (~50%) in the calculation of the band gap ZPR. Instead, when using a stochastic method, we obtain results in good agreement with those of our quantum simulations for the diamondoid crystal. However, the agreement is worse for NAI-DMAC where intramolecular anharmonicities contribute to the ZPR. Our results highlight the importance of accurately including nuclear and anharmonic quantum effects to predict the electronic properties of molecular crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗