SEARCH · Engineering Papers

Results for “Functional coupling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,459 records · Page 81

All sources NASA NTRS DOE OSTI

NASA NTRS2004

LSPRAY-II: A Lagrangian Spray Module

LSPRAY-II is a Lagrangian spray solver developed for application with parallel computing and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo Probability Density Function (PDF) solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type for the gas flow grid representation. It is mainly designed to predict the flow, thermal and transport properties of a rapidly vaporizing spray because of its importance in aerospace application. The manual provides the user with an understanding of various models involved in the spray formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers. With the development of LSPRAY-II, we have advanced the state-of-the-art in spray computations in several important ways.

Raju, M. S.↗

NASA NTRS2008

LSPRAY-III: A Lagrangian Spray Module

LSPRAY-III is a Lagrangian spray solver developed for application with parallel computing and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo Probability Density Function (PDF) solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type for the gas flow grid representation. It is mainly designed to predict the flow, thermal and transport properties of a rapidly vaporizing spray because of its importance in aerospace application. The manual provides the user with an understanding of various models involved in the spray formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers. With the development of LSPRAY-III, we have advanced the state-of-the-art in spray computations in several important ways.

Raju, M. S.↗

DOE OSTI2022

Photon detector response function methodology using MCNP and shift hybrid radiation transport code for wide-area contamination assay applications

Here, radiation transport modeling using the Monte Carlo N-Particle (MCNP) radiation transport code and Monte Carlo code, Shift, were employed to model detector responses for a variety of wide-area photon contamination scenarios. In this study, 2" × 2" and 3" × 3" cylindrical NaI(Tl) scintillation detector configurations at source detector-distances of 0.5 cm, 1 cm, 2.54 cm, 10 cm, and 30 cm were modeled. Media of soil, concrete, and steel were evaluated for contamination depths ranging from surface to a depth of an infinite thickness in each medium for photon energies ranging from 20 keV to 3 MeV, which correspond to the energies that current detectors can discern. Monoenergetic photon surface contamination detector responses for each of the media, source–detector distances, and detectors were estimated using MCNP v6.2. Shift was harnessed for improved variance reduction of particle transport in highly attenuating media to obtain average cell fluxes in the two MCNP NaI(Tl) scintillation detector configurations. Average cell flux values in Shift were coupled with detector responses from MCNP to convert average cell flux in a void to energy distribution of pulses in the NaI(Tl) scintillation detector crystal of interest. An optimized detector response function methodology was developed by coupling the MCNP radiation transport method with the Consistent Adjoint Driven Importance Sampling (CADIS) hybrid radiation transport method built into Shift to significantly decrease the runtime of thousands of MCNP pulse height simulations. The methodology may be utilized to quickly and accurately facilitate the assessment of a broad range of wide-area environmental contamination assay and decommissioning cleanup applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Ultrafast infrared transient absorption spectroscopy of gas-phase Ni(CO) 4 photodissociation at 261 nm

Here, we employ ultrafast mid-infrared transient absorption spectroscopy to probe the rapid loss of carbonyl ligands from gas-phase nickel tetracarbonyl following ultraviolet photoexcitation at 261 nm. Here, nickel tetracarbonyl undergoes prompt dissociation to produce nickel tricarbonyl in a singlet excited state; this electronically excited tricarbonyl loses another CO group over tens of picoseconds. Our results also suggest the presence of a parallel, concerted dissociation mechanism to produce nickel dicarbonyl in a triplet excited state, which likely dissociates to nickel monocarbonyl. Mechanisms for the formation of these photoproducts in multiple electronic excited states are theoretically predicted with one-dimensional cuts through the potential energy surfaces and computation of spin–orbit coupling constants using equation of motion coupled cluster methods (EOM-CC) and coupled cluster theory with single and double excitations (CCSD). Bond dissociation energies are calculated with CCSD, and anharmonic frequencies of ground and excited state species are computed using density functional theory (DFT) and time-dependent density functional theory (TD-DFT).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Coupling Subsurface and Above-Surface Models for Optimizing the Design of Borefields and District Heating and Cooling Systems

Accurate dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange for seasonal energy storage. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes without flowing groundwater. While detailed simulation tools for subsurface heat and mass transfer exist, these fall short in simulating above-surface energy systems. To support the design and operation of such systems, the study developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the subsurface. For the first, it uses the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, it uses the TOUGH simulator. The TOUGH family of codes can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. The study validated the coupled modeling approach by comparing the simulation results with one from the g-function based ground response model. It then looked into effects when the water table and the regional groundwater flow are considered in the ground, from the perspective of heat exchange between borehole and ground, and the electrical consumption of the district heating and cooling systems. To access the simulation models, please find the links in the submission: -- For coupled approach validation: see model Buildings.Fluid.Geothermal.Borefields.Examples.BorefieldsWithTough and Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from the "Modelica Building Library" resource, branch issue1495_tough_interface, commit a2667c0. -- For the study of the effect of water table: see model Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO, commit 760de49. -- For the study of the effect of regional groundwater flow: see Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO_3D, commit c2a2d2a. The coupling interface script "GrounResponse.py" can be found from the above links in the folder Buildings/Resources/Python-Sources. Also, the needed files for TOUGH simulation are in the folder Buildings/Resources/Python-Sources/ToughFiles that can be accessed through the above links. A brief description of these files is given below; detailed specifications for the first three files may be found in the TOUGH3 Users Guide (Jung et al., 2018) https://tough.lbl.gov/documentation/tough-manuals/. (1) INCON - initial conditions for each grid block (2) INFILE - main input file with material properties and control parameters (3) MESH - description of the computational grid (4) readsave - Modelica/TOUGH interface program: read the final output of TOUGH simulation after TOUGH time step and prepare for transfer to Modelica for next Modelica time step (5) readsave.inp - input parameters for program readsave (6) writeincon - Modelica/TOUGH interface program: write the output of Modelica after Modelica time step and prepare for transfer to TOUGH as initial conditions for the next TOUGH step (7) writeincon.inp - input parameters for program writeincon

15 GEOTHERMAL ENERGY↗

NASA NTRS1984

Bolometers - Ultimate sensitivity, optimization, and amplifier coupling

Theoretical expressions for Johnson noise and thermal noise in bolometers are considered, and optimization with respect to thermal conductivity and bias power is performed. Numerical approximations are given for the ultimate NEP of bolometers as a function of material parameters and compared with photon noise including photon correlations. A resonating capacitor is shown to improve the coupling to an amplifier, so that the amplifier need not limit performance even for very low temperature bolometers.

Mather, J. C.↗

NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗

NASA NTRS1977

The molecular composition of dense interstellar clouds

Presented in this paper is an ab initio chemical model for dense interstellar clouds that incorporates 598 grain surface reactions, with small grains providing the reaction area. Gas-phase molecules are depleted through collisions with grains. The abundances of 372 chemical species are calculated as a function of time and are found to be of sufficient magnitude to explain most observations. Peak abundances are achieved on time scales of the order of 100,000 to 1 million years, depending on cloud density and kinetic temperature. The reaction rates for ion-molecule chemistry are approximately the same, indicating that surface and gas-phase chemistry may be coupled in certain regions. The composition of grain mantles is shown to be a function of grain radius. In certain grain-size ranges, large molecules containing two or more heavy atoms are more predominant than lighter 'ices' - H2O, NH3, and CH4. It is possible that absorption due to these large molecules in the mantle may contribute to the observed 3-micron band in astronomical spectra.

Allen, M.↗

DOE OSTI2023

Trapping-induced quantum beats in a van-der-Waals heterostructure microcavity observed by two-dimensional micro-spectroscopy

Spatial confinement has been frequently engineered to control the flow and relaxation dynamics of exciton polaritons. While widely investigated in GaAs microcavities, exciton-polariton coupling between discretized polariton modes arising from spatially confined 2D crystals been has been less exhaustively studied. Here, we use coherent 2D photoluminescence-detected micro-spectroscopy to detect oscillating 2D peaks exclusively from a spatial trap in a microcavity with an embedded van-der-Waals heterostructure at room temperature. We observe a wide variation of oscillatory phases as a function of spectral position within the 2D spectrum, which suggests the existence of a coupling between the discretized polariton modes. The latter is accompanied by the generation of coherent phonons.

Materials Science↗

DOE OSTI2022

Development of a Control Co-Design Modeling Tool for Marine Hydrokinetic Turbines: Preprint

This report describes the ongoing and planned development of the software package CT-Opt (Current/Tidal Optimization), a control co-design modeling tool for marine hydrokinetic turbines. The commercialization of these turbines has faced significant challenges due to the complex, multidisciplinary nature of their design and the extreme environmental conditions of their operation. This project aims to create a modeling tool that will enable the efficient design of robust, cost-competitive hydrokinetic turbine systems. Rather than using traditional optimization methods, CT-Opt combines multiple models across a range of fidelities to enable coupled optimization of the system design and system controller via a control co-design approach. With this method, the parameters that affect system performance are considered more comprehensively at every stage of the design process. The lowest-fidelity, frequency-domain model called by CT-Opt is RAFT (Response Amplitudes of Floating Turbines), which was originally developed by the National Renewable Energy Laboratory (NREL) to model response amplitudes of floating offshore wind turbines. The highest-fidelity, time-domain model is OpenFAST, which was developed by NREL for land-based and offshore wind turbines. As part of the CT-Opt project, new functionalities will be added to RAFT and OpenFAST to enable the accurate simulation of fixed and floating marine hydrokinetic turbines. In addition to expanding the capabilities of RAFT and OpenFAST, new midfidelity models will be developed. These models will be based on RAFT and OpenFAST and will consist of linearized, state-space models derived from the fully coupled, nonlinear OpenFAST equations and derivative function surrogate models that approximate the nonlinear system behavior. Each model will be coupled with controllers to allow control co-design methods to be applied both within models and across fidelity levels, enabling efficient system optimization.

control co-design↗

NASA NTRS2004

Improved Ionospheric Electrodynamic Models and Application to Calculating Joule Heating Rates

Improved techniques have been developed for empirical modeling of the high-latitude electric potentials and magnetic field aligned currents (FAC) as a function of the solar wind parameters. The FAC model is constructed using scalar magnetic Euler potentials, and functions as a twin to the electric potential model. The improved models have more accurate field values as well as more accurate boundary locations. Non-linear saturation effects in the solar wind-magnetosphere coupling are also better reproduced. The models are constructed using a hybrid technique, which has spherical harmonic functions only within a small area at the pole. At lower latitudes the potentials are constructed from multiple Fourier series functions of longitude, at discrete latitudinal steps. It is shown that the two models can be used together in order to calculate the total Poynting flux and Joule heating in the ionosphere. An additional model of the ionospheric conductivity is not required in order to obtain the ionospheric currents and Joule heating, as the conductivity variations as a function of the solar inclination are implicitly contained within the FAC model's data. The models outputs are shown for various input conditions, as well as compared with satellite measurements. The calculations of the total Joule heating are compared with results obtained by the inversion of ground-based magnetometer measurements. Like their predecessors, these empirical models should continue to be a useful research and forecast tools.

Weimer, D. R.↗

DOE OSTI2024

Impact Analysis of High-Altitude Electromagnetic Pulse Coupling Effects on Power Grid Protection Relays

Protection relays are important equipment used for protection, control, and metering functions in the power grid. These relays are used to protect critical and difficult-to-replace equipment, like generators, transformers, and capacitor banks. Once the protection devices are disturbed or damaged, a high risk of power generation interruption occurs. Therefore, it is important and necessary to study the relay’s immunity to electromagnetic pulse (EMP) events. As a preliminary step toward empirical experimentation on actual equipment, this manuscript outlines an economical and efficient methodology for evaluating the impact of an EMP. An impedance measurement strategy was employed to model the equipment, setting the stage for subsequent immunity analyses. These analyses included the pulse current injection (PCI) method, which utilized an injecting probe to introduce the transient, and frequency domain electromagnetic (FEKO) simulation, which integrated electromagnetic coupling effects into the transient simulation. The impedance measurement and simulation results in this paper provide a reliable basis for gauging equipment performance in the face of HEMP threats. The results obtained using the PCI and FEKO simulations demonstrated the performance of different port responses under a high-altitude EMP, indicating the requirement for some protection to ensure the reliable operation of relays.

42 ENGINEERING↗

NASA NTRS2015

Assessing and Promoting Functional Resilience in Flight Crews During Exploration Missions

NASA plans to send humans to Mars in about 20 years. The NASA Human Research Program supports research to mitigate the major risks to human health and performance on extended missions. However, there will undoubtedly be unforeseen events on any mission of this nature - thus mitigation of known risks alone is not sufficient to ensure optimal crew health and performance. Research should be directed not only to mitigating known risks, but also to providing crews with the tools to assess and enhance resilience, as a group and individually. We can draw on ideas from complexity theory and network theory to assess crew and individual resilience. The entire crew or the individual crewmember can be viewed as a complex system that is composed of subsystems (individual crewmembers or physiological subsystems), and the interactions between subsystems are of crucial importance for overall health and performance. An understanding of the structure of the interactions can provide important information even in the absence of complete information on the component subsystems. This is critical in human spaceflight, since insufficient flight opportunities exist to elucidate the details of each subsystem. Enabled by recent advances in noninvasive measurement of physiological and behavioral parameters, subsystem monitoring can be implemented within a mission and also during preflight training to establish baseline values and ranges. Coupled with appropriate mathematical modeling, this can provide real-time assessment of health and function, and detect early indications of imminent breakdown. Since the interconnected web of physiological systems (and crewmembers) can be interpreted as a network in mathematical terms, we can draw on recent work that relates the structure of such networks to their resilience (ability to self-organize in the face of perturbation). There are many parameters and interactions to choose from. Normal variability is an established characteristic of a healthy physiological response. Healthy coupling has been investigated less extensively, but there are cases in which too tight or too loose coupling can be problematic. This might be in inter-individual behaviors, such as sleep cycles, coordination of work and meal times, and coupled motions during communication. Less apparent are couplings of physiological systems, nevertheless examples abound of coupled systems which might be monitored: cardio-respiratory rhythms; circadian rhythms, body temperature, and sleep; stress markers and cognition, sleep, and performance; profiles of biochemical markers related to immune function and nutritional status; sensorimotor aspects such as motion sickness, ataxia, reaction time, and manual control. Tools for resilience are then the means to measure and analyze these parameters, incorporate them into appropriate models of normal variability and interconnectedness, and recognize when parameters or their couplings are outside of normal limits. What to do when a problem is identified depends on its nature. Changes can be made to crew procedures, work pacing, interpersonal interactions, sleep cycles, meal timing and content, as guided by the model. Use and continued development of these methods could not only provide tools for resilience, but also meaningful autonomous work for the crew on an extended flight.

Shelhamer, Mark↗

DOE OSTI2020

Effective gravitational couplings of four-dimensional $ \mathcal{N} $ = 2 supersymmetric gauge theories

The low energy effective couplings of a four-dimensional N = 2 supersymmetric gauge theory to topological invariants of the background gravitational field are described by two functions A and B. These two functions play an important role in the study of topologically twisted four-dimensional N = 2 supersymmetric gauge theories and in the computation of central charges of N = 2 superconformal theories. In this paper, we compute A and B from the partition function in the Ω-background for SU(2) gauge theories. Our results not only confirm the predicted expressions of the effective gravitational couplings, but also give the previously undetermined overall multiplicative factors. We also analyze A and B for the SU(N) super-Yang-Mills theory, and confirm all the previous predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2006

Tunable Filter Made From Three Coupled WGM Resonators

A tunable third-order band-pass optical filter has been constructed as an assembly of three coupled, tunable, whispering-gallery-mode resonators similar to the one described in Whispering-Gallery-Mode Tunable Narrow-Band-Pass Filter (NPO-30896), NASA Tech Briefs, Vol. 28, No. 4 (April 2004), page 5a. This filter offers a combination of four characteristics that are desirable for potential applications in photonics: (1) wide real-time tunability accompanied by a high-order filter function, (2) narrowness of the passband, (3) relatively low loss between input and output coupling optical fibers, and (4) a sparse spectrum. In contrast, prior tunable band-pass optical filters have exhibited, at most, two of these four characteristics. As described in several prior NASA Tech Briefs articles, a whispering-gallery-mode (WGM) resonator is a spheroidal, disklike, or toroidal body made of a highly transparent material. It is so named because it is designed to exploit whispering-gallery electromagnetic modes, which are waveguide modes that propagate circumferentially and are concentrated in a narrow toroidal region centered on the equatorial plane and located near the outermost edge. Figure 1 depicts the optical layout of the present filter comprising an assembly of three coupled, tunable WGM resonators. Each WGM resonator is made from a disk of Z-cut LiNbO3 of 3.3-mm diameter and 50-m thickness. The perimeter of the disk is polished and rounded to a radius of curvature of 40 microns. The free spectral range of each WGM resonator is about 13.3 GHz. Gold coats on the flat faces of the disk serve as electrodes for exploiting the electro-optical effect in LiNbO3 for tuning. There is no metal coat on the rounded perimeter region, where the whispering-gallery modes propagate. Light is coupled from an input optical fiber into the whispering-gallery-modes of the first WGM resonator by means of a diamond prism. Another diamond prism is used to couple light from the whispering-gallery-modes of the third WGM resonator to an output optical fiber. The filter operates at a nominal wavelength of 1,550 nm and can be tuned over a frequency range of plus or minus 12 GHz by applying a potential in the range of plus or minus 150 V to the electrodes. The insertion loss (the loss between the input and output coupling optical fibers) was found to be repeatable at 6 dB. The resonance quality factor (Q) of the main sequence of resonator modes was found to be 5 x 10(exp 6), which corresponds to a bandwidth of 30 MHz. The filter can be shifted from one operating frequency to another within a tuning time less than or equal to 30 micro seconds. The transmission curve of the filter at frequencies near the middle of the passband closely approximates a theoretical third-order Butterworth filter profile, as shown in Figure 2.

Savchenkov, Anatoliy↗

DOE OSTI2023

Spin–phonon interactions and magnetoelectric coupling in Co 4 B 2 O 9 ( B = Nb, Ta)

In order to explore the consequences of spin–orbit coupling on spin–phonon interactions in a set of chemically similar mixed metal oxides, we measured the infrared vibrational properties of Co 4 B 2 O 9 (B = Nb, Ta) as a function of temperature and compared our findings with lattice dynamics calculations and several different models of spin–phonon coupling. Frequency vs temperature trends for the Co 2+ shearing mode near 150 cm –1 reveal significant shifts across the magnetic ordering temperature that are especially large in relative terms. Bringing these results together and accounting for noncollinearity, we obtain spin–phonon coupling constants of –3.4 and –4.3 cm –1 for Co 4 Nb 2 O 9 and the Ta analog, respectively. Analysis reveals that these coupling constants are derived from interlayer (rather than intralayer) exchange interactions and that the interlayer interactions contain competing antiferromagnetic and ferromagnetic contributions. At the same time, beyond-Heisenberg terms are minimized due to fortuitous symmetry considerations, different from most other 4d- and 5d-containing oxides. Comparison with other contemporary oxides shows that spin–phonon coupling in this family of materials is among the strongest ever reported, suggesting an origin for magnetoelectric coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Unified Approach to Polarons and Phonon-Induced Band Structure Renormalization

Ab initio calculations of the phonon-induced band structure renormalization are currently based on the perturbative Allen-Heine theory and its many-body generalizations. These approaches are unsuitable to describe materials where electrons form localized polarons. Here, we develop a self-consistent, many-body Green’s function theory of band structure renormalization that incorporates localization and self-trapping. Here, we show that the present approach reduces to the Allen-Heine theory in the weak-coupling limit, and to total energy calculations of self-trapped polarons in the strong-coupling limit. To demonstrate this methodology, we reproduce the path-integral results of Feynman and diagrammatic Monte Carlo calculations for the Fröhlich model at all couplings, and we calculate the zero point renormalization of the band gap of an ionic insulator including polaronic effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Predicting Mössbauer Parameters of Nonheme Diiron Complexes with Density Functional Theory

Mössbauer spectroscopy provides significant insights into the electronic structure and environment of the metal centers. In this work, we investigate the electronic structures of a set of nonheme diiron complexes by evaluating two key parameters pertaining to Mössbauer spectroscopy, namely, the isomer shift (δ) and quadrupole splitting (|ΔE Q |), using different levels of density functional theory (DFT). The diiron systems investigated here span diverse oxidation states, bridging motifs, and spin coupling patterns, which present a challenging case for theoretical predictions. We demonstrate that the combination of B97-D3/def2-TZVP is an efficient approach in modeling both the δ and |ΔEQ| values with high accuracy for the representative nonheme diiron complexes. We also show that δ is accurately predicted irrespective of the choice of approximate density functional while the |ΔE Q | is sensitive to the level of theory employed. Further investigation shows that the present methodology assessed using synthetic nonheme diiron complexes could be extended to nonheme diiron enzyme active sites, featuring both ferromagnetic and antiferromagnetic coupling between the iron centers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗