Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complexes”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Gaseous cyclodextrin- closo -dodecaborate complexes χCD·B 12 X 12 2– (χ = α, β, and γ; X = F, Cl, Br, and I): electronic structures and intramolecular interactions

A fundamental understanding of cyclodextrin- closo -dodecaborate inclusion complexes is of great interest in supramolecular chemistry. Herein, we report a systematic investigation on the electronic structures and intramolecular interactions of perhalogenated closo -dodecaborate dianions B 12 X 12 2– (X = F, Cl, Br and I) binding to α-, β-, and γ-cyclodextrins (CDs) in the gas phase using combined negative ion photoelectron spectroscopy (NIPES) and density functional theory (DFT) calculations. The vertical detachment energy (VDE) of each complex and electronic stabilization of each dianion due to the CD binding (ΔVDE, relative to the corresponding isolated B 12 X 12 2– ) are determined from the experiments along α-, β- and γ-CD in the form of VDE (ΔVDE): 4.00 (2.10), 4.33 (2.43), and 4.30 (2.40) eV in X = F; 4.09 (1.14), 4.64 (1.69), and 4.69 (1.74) eV in X = Cl; 4.11 (0.91), 4.58 (1.38), and 4.70 (1.50) eV in X = Br; and 3.54 (0.74), 3.88 (1.08), and 4.05 (1.25) eV in X = I, respectively. All complexes have significantly higher VDEs than the corresponding isolated dodecaborate dianions with ΔVDE spanning from 0.74 eV at (α, I) to 2.43 eV at (β, F), sensitive to both host CD size and guest substituent X. DFT-optimized complex structures indicate that all B 12 X 12 2– prefer binding to the wide openings of CDs with the insertion depth and binding motif strongly dependent on the CD size and halogen X. Dodecaborate anions with heavy halogens, i.e. , X = Cl, Br, and I, are found outside of α-CD, while B 12 F 12 2– is completely wrapped by γ-CD. Partial embedment of B 12 X 12 2– into CDs is observed for the other complexes via multipronged B–X … H–O/C interlocking patterns. The simulated spectra based on the density of states agree well with those of the experiments and the calculated VDEs well reproduce the experimental trends. Molecular orbital analyses suggest that the spectral features at low binding energies originated from electrons detached from the dodecaborate dianion, while those at higher binding energies are derived from electron detachment from CDs. Energy decomposition analyses reveal that the electrostatic interaction plays a dominating role in contributing to the host–guest interactions for the X = F series partially due to the formation of a O/C–H … X–B hydrogen bonding network, and the dispersion forces gradually become important with the increase of halogen size.

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

Near-infrared C-term MCD spectroscopy of octahedral uranium (V) complexes

C-term magnetic circular dichroism (MCD) spectroscopy is a powerful method for probing d–d and f–f transitions in paramagnetic metal complexes. However, this technique remains underdeveloped both experimentally and theoretically for studies of U(V) complexes of Oh symmetry, which have been of longstanding interest for probing electronic structure, bonding, and covalency in 5f systems. In this study, C-term NIR MCD of the Laporte forbidden f–f transitions of [UCl 6 ] - and [UF 6 ] - are reported, demonstrating the significant fine structure resolution possible with this technique including for the low energy Γ 7 → Γ 8 transitions in [UF 6 ] - . The experimental NIR MCD studies were further extended to [U(OC 6 F 5 ) 6 ] - , [U(CH 2 SiMe 3 ) 6 ] - , and [U(NC(tBu)(Ph)) 6 ] - to evaluate the effects of ligand-type on the f–f MCD fine structure features. Theoretical calculations were conducted to determine the Laporte forbidden f–f transitions and their MCD intensity experimentally observed in the NIR spectra of the U(V) hexahalide complexes, via the inclusion of vibronic coupling, to better understand the underlying spectral fine structure features for these complexes. These spectra and simulations provide an important platform for the application of MCD spectroscopy to this widely studied class of U(V) complexes and identify areas for continued theoretical development.

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

Actinide arene-metalates: ion pairing effects on the electronic structure of unsupported uranium–arenide sandwich complexes

Addition of [UI 2 (THF) 3 (μ-OMe)] 2 ·THF (2·THF) to THF solutions containing 6 equiv. of K[C 14 H 10 ] generates the heteroleptic dimeric complexes [K(18-crown-6)(THF) 2 ] 2 [U(η 6 -C 14 H 10 )(η 4 -C 14 H 10 )(μ-OMe)] 2 ·4THF (1 18C6 ·4THF) and {[K(THF) 3 ][U(η 6 -C 14 H 10 )(η 4 -C 14 H 10 )(μ-OMe)]} 2 (1 THF ) upon crystallization of the products in THF in the presence or absence of 18-crown-6, respectively. Both 1 18C6 ·4THF and 1 THF are thermally stable in the solid-state at room temperature; however, after crystallization, they become insoluble in THF or DME solutions and instead gradually decompose upon standing. X-ray diffraction analysis reveals 1 18C6 ·4THF and 1 THF to be structurally similar, possessing uranium centres sandwiched between bent anthracenide ligands of mixed tetrahapto and hexahapto ligation modes. Yet, the two complexes are distinguished by the close contact potassium-arenide ion pairing that is seen in 1 THF but absent in 1 18C6 ·4THF, which is observed to have a significant effect on the electronic characteristics of the two complexes. Structural analysis, SQUID magnetometry data, XANES spectral characterization, and computational analyses are generally consistent with U(IV) formal assignments for the metal centres in both 1 18C6 ·4THF and 1 THF , though noticeable differences are detected between the two species. For instance, the effective magnetic moment of 1 THF (3.74 μ B ) is significantly lower than that of 1 18C6 ·4THF (4.40 μ B ) at 300 K. Furthermore, the XANES data shows the U L III -edge absorption energy for 1 THF to be 0.9 eV higher than that of 1 18C6 ·4THF, suggestive of more oxidized metal centres in the former. Of note, CASSCF calculations on the model complex {[U(η 6 -C 14 H 10 )(η 4 -C 14 H 10 )(μ-OMe)] 2 } 2– (1*) shows highly polarized uranium–arenide interactions defined by π-type bonds where the metal contributions are primarily comprised by the 6d-orbitals (7.3 ± 0.6%) with minor participation from the 5f-orbitals (1.5 ± 0.5%). These unique complexes provide new insights into actinide–arenide bonding interactions and show the sensitivity of the electronic structures of the uranium atoms to coordination sphere effects.

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

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗
DOE OSTI2025

Dicopper( i ) complexes of a binucleating, dianionic, naphthyridine bis(amide) ligand

The dinucleating ligand, 1,8-naphthyridine-2,7-bis(2,6-diisopropylphenyl)carboxamide (NBDA), was synthesized by palladium-catalyzed aminocarbonylation. This ligand was treated with two equivalents of mesitylcopper(I) in the presence of [nBu4N]X (X = Cl, N3) to give the anionic complexes [nBu4N][Cu2(NBDA)(μ-Cl)] and [nBu4N][Cu2(NBDA)(μ-N3)]. Treatment of H2NBDA with mesitylcopper(I) and two equivalents of xylyl isocyanide led to the formation of a charge-neutral dicopper(I) complex, [Cu2(NBDA)(CNXyl)2], displaying two isocyanide ligands, each terminally bound to one of the copper atoms. The complexes were characterized by NMR and IR spectroscopy, as well as by single-crystal X-ray diffraction analysis. Electrochemical characterization of the complexes using cyclic voltammetry revealed a reversible ligand-based reduction between -1.65 and -2.0 V vs. Fc/Fc+. DFT calculations suggest a more ionic bonding character and weaker Cu-Cu interactions in the NBDA complexes compared to those with other 1,8-naphthyridine-based ligands. This is congruent with intermetallic separations of over 3 Å induced by relatively strong coordination of the copper atoms to the amide nitrogen donor atoms observed in the solid state molecular structures.

Sévery, Laurent↗
DOE OSTI2025

Structure and ultrafast dynamics of tri-nuclear Ag-/Tl–Pt 2 POP 4 complexes in solution

The energetics and dynamics of ion assembly in solution has broad influence in nanomaterials and inorganic synthesis. To investigate the fundamental processes involved, we present a time-resolved x-ray solution scattering (TR-XSS) study of the trinuclear silver and thallium complexes of the diplatinum ion PtPOP [Pt 2 (H 2 P 2 O 5 )$_4^{4−}$] in aqueous solution. These complexes, their structural properties, and their electronic structure are not well understood and afford a unique opportunity to study the metal–metal bond formation that influences molecular and material assembly in solution. We present model-independent analysis of the observed dynamics as well as an analysis incorporating time-resolved structural refinements of key bond lengths with $<$100 fs time resolution. We find that upon photoexcitation, the Pt atoms contract ∼0.25 Å toward the center of both the Ag- and the Tl-PtPOP complexes, as previously observed for the PtPOP anion. For the AgPtPOP system, an ultrafast Ag-Pt bond expansion of ∼0.2 Å is observed, whereas in contrast, the TlPtPOP system exhibits a Tl-Pt bond contraction of ∼0.3 Å upon photoexcitation. For both complexes, the change in electronic state leads to coherent (“wave-packet”) oscillations along the metal–Pt coordinates. Based on these structural dynamics, we propose an electronic structure model that describes the metal–metal bonding behavior in both the ground and excited state for both complexes.

Lenzen, Philipp [Technical Univ. of Denmark, Lyngb↗
DOE OSTI2021

Coupled cluster benchmarks of large noncovalent complexes: The L7 dataset as well as DNA–ellipticine and buckycatcher–fullerene

In this work, benchmark binding energies for dispersion-bound complexes in the L7 dataset, the DNA–ellipticine intercalation complex, and the buckycatcher–C60 complex with 120 heavy atoms using a focal-point method based on the canonical form of second-order Møller−Plesset theory (MP2) and the domain based local pair natural orbital scheme for the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set (CBS) limit are reported. This work allows for increased confidence given the agreement with respect to values recently obtained using the local natural orbital CCSD(T) for L7 and the canonical CCSD(T)/CBS result for the coronene dimer (C2C2PD). Therefore, these results can be considered pushing the CCSD(T)/CBS binding benchmark to the hundred-atom scale. The disagreements between the two state-of-the-art methods, CCSD(T) and fixed-node diffusion Monte Carlo, are substantial with at least 2.0 (∼10%), 1.9 (∼5%), and 10.3 kcal/mol (∼25%) differences for C2C2PD in L7, DNA–ellipticine, and buckycatcher–C60, respectively. Such sizable discrepancy above “chemical accuracy” for large noncovalent complexes indicates how challenging it is to obtain benchmark binding interactions for systems beyond small molecules, although the three up-to-date density functionals, PBE0+D4, ωB97M-V, and B97M-V, agree better with CCSD(T) for these large systems. In addition to reporting these values, different basis sets and various CBS extrapolation parameters for Hartree–Fock and MP2 correlation energies were tested for the first time in large noncovalent complexes with the goal of providing some indications toward optimal cost effective routes to approach the CBS limit without substantial loss in quality.

Ballesteros, Francisco (ORCID:0000000325732307)↗
DOE OSTI2021

Revealing the nature of electron correlation in transition metal complexes with symmetry breaking and chemical intuition

Here, in this work, we provide a nuanced view of electron correlation in the context of transition metal complexes, reconciling computational characterization via spin and spatial symmetry breaking in single-reference methods with qualitative concepts from ligand-field and molecular orbital theories. These insights provide the tools to reliably diagnose the multi-reference character, and our analysis reveals that while strong (i.e., static) correlation can be found in linear molecules (e.g., diatomics) and weakly bound and antiferromagnetically coupled (monometal-noninnocent ligand or multi-metal) complexes, it is rarely found in the ground-states of mono-transition-metal complexes. This leads to a picture of static correlation that is no more complex for transition metals than it is, e.g., for organic biradicaloids. In contrast, the ability of organometallic species to form more complex interactions, involving both ligand-to-metal σ-donation and metal-to-ligand π-backdonation, places a larger burden on a theory's treatment of dynamic correlation. We hypothesize that chemical bonds in which inter-electron pair correlation is non-negligible cannot be adequately described by theories using MP2 correlation energies and indeed find large errors vs experiment for carbonyl-dissociation energies from double-hybrid density functionals. A theory's description of dynamic correlation (and to a less important extent, delocalization error), which affects relative spin-state energetics and thus spin symmetry breaking, is found to govern the efficacy of its use to diagnose static correlation.

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

A bottom-up perspective on photodynamics and photoprotection in light-harvesting complexes using anti-Brownian trapping

Single-molecule fluorescence spectroscopy allows direct, real-time observation of dynamic photophysical changes in light harvesting complexes. The Anti-Brownian ELectrokinetic (ABEL) trap is one such single-molecule method with useful advantages. This approach is particularly well-suited to make detailed spectroscopic measurements of pigment-protein complexes in solution phase because it enables extended-duration single-molecule observation by counteracting Brownian motion. This Perspective summarizes recent contributions by the authors and others that have utilized the unique capabilities of the ABEL trap to advance our understanding of phycobiliproteins and the phycobilisome complex, the primary light-harvesting apparatus of cyanobacteria. Monitoring the rich spectroscopic data from these measurements, which include brightness, fluorescence lifetime, polarization, and emission spectra, among other measurable parameters, has provided direct characterization of pigments and energy transfer pathways in the phycobilisome, spanning scales from single pigments and monomeric phycobiliproteins to higher order oligomers and protein-protein interactions of the phycobilisome complex. Importantly, new photophysical states and photodynamics were observed to modulate the flow of energy through the phycobilisome and suggest a previously unknown complexity in phycobilisome light harvesting and energy transport with a possible link to photoadaptive or photoprotective functions in cyanobacteria. Here, beyond deepening our collective understanding of natural light-harvesting systems, these and future discoveries may serve as inspiration for engineering improved artificial light-harvesting technologies.

14 SOLAR ENERGY↗
DOE OSTI2026

Complexity analysis of a CT injection experiment on BRB

In this work, we use Jensen–Shannon complexity and permutation entropy to analyze the magnetic field fluctuations of an astrophysically scaled plasma experiment. The experiment was intended to emulate an interplanetary coronal mass ejection event in the lab, recreating the major sections seen in satellite data. We also use a technique called “delay,” in which we use select elements, skipping one or more data points at a time, in our time series data to obtain Jensen–Shannon complexity as a function of frequency and investigate the frequency of maximized complexity. We then compare the delay frequencies to other frequencies in the plasma. We found that the frequencies for maximum complexity do not correspond to the frequencies investigated, implying that other physical mechanisms lead to an increase in complexity at these frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Redox conditions correlated with vibronic coupling modulate quantum beats in photosynthetic pigment–protein complexes

Significance Photosynthetic organisms evolved their light-harvesting antenna complexes to optimize energy transfer. It was recently shown that the redox environment can tune the mixing of electronic and vibrational states to steer energy through different pathways of a pigment–protein complex. Quantum beating signals in the spectra of pigment–protein complexes have been used to probe the excited-state dynamics within the complexes, but the microscopic dynamics that generate these signals and their role in promoting energy transfer are not fully understood. Here, we show that the redox environment that tunes energy transfer similarly tunes the quantum beating signals in the same complex. We find that the beats report on excited-state vibrations that maintain coherence through the vibronically enhanced energy transfer process.

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

Mechanistic roles of metal- and ligand-protonated species in hydrogen evolution with [Cp*Rh] complexes

Protonation reactions involving organometallic complexes are ubiquitous in redox chemistry and often result in the generation of reactive metal hydrides. However, some organometallic species supported by η 5 -pentamethylcyclopentadienyl (Cp*) ligands have recently been shown to undergo ligand-centered protonation by direct proton transfer from acids or tautomerization of metal hydrides, resulting in the generation of complexes bearing the uncommon η 4 -pentamethylcyclopentadiene (Cp*H) ligand. Here, time-resolved pulse radiolysis (PR) and stopped-flow spectroscopic studies have been applied to examine the kinetics and atomistic details involved in the elementary electron- and proton-transfer steps leading to complexes ligated by Cp*H, using Cp*Rh(bpy) as a molecular model (where bpy is 2,2′-bipyridyl). Stopped-flow measurements coupled with infrared and UV-visible detection reveal that the sole product of initial protonation of Cp*Rh(bpy) is [Cp*Rh(H)(bpy)] + , an elusive hydride complex that has been spectroscopically and kinetically characterized here. Tautomerization of the hydride leads to the clean formation of [(Cp*H)Rh(bpy)] + . Variable-temperature and isotopic labeling experiments further confirm this assignment, providing experimental activation parameters and mechanistic insight into metal-mediated hydride-to-proton tautomerism. Spectroscopic monitoring of the second proton transfer event reveals that both the hydride and related Cp*H complex can be involved in further reactivity, showing that [(Cp*H)Rh] is not necessarily an off-cycle intermediate, but, instead, depending on the strength of the acid used to drive catalysis, an active participant in hydrogen evolution. Identification of the mechanistic roles of the protonated intermediates in the catalysis studied here could inform design of optimized catalytic systems supported by noninnocent cyclopentadienyl-type ligands.

08 HYDROGEN↗
DOE OSTI2024

Path integrals, complex probabilities and the discrete Weyl representation

Abstract A discrete formulation of the real-time path integral as the expectation value of a functional of paths with respect to a complex probability on a sample space of discrete valued paths is explored. The formulation in terms of complex probabilities is motivated by a recent reinterpretation of the real-time path integral as the expectation value of a potential functional with respect to a complex probability distribution on cylinder sets of paths. The discrete formulation in this work is based on a discrete version of the Weyl algebra that can be applied to any observable with a finite number of outcomes. The origin of the complex probability in this work is the completeness relation. In the discrete formulation the complex probability exactly factors into products of conditional probabilities and exact unitarity is maintained at each level of approximation. The approximation of infinite dimensional quantum systems by discrete systems is discussed. The method is illustrated by applying it to scattering theory and quantum field theory. The implications of these applications for quantum computing is discussed.

Physics↗
DOE OSTI2023

The structure of the high-affinity nickel-binding site in the Ni,Zn-HypA•UreE2 complex

Abstract The maturation pathway for the nickel-dependent enzyme urease utilizes the protein UreE as a metallochaperone to supply Ni(II) ions. In Helicobacter pylori urease maturation also requires HypA and HypB, accessory proteins that are commonly associated with hydrogenase maturation. Herein we report on the characterization of a protein complex formed between HypA and the UreE2 dimer. Nuclear magnetic resonance (NMR) coupled with molecular modelling show that the protein complex apo, Zn-HypA•UreE2, forms between the rigorously conserved Met-His-Glu (MHE motif) Ni-binding N-terminal sequence of HypA and the two conserved His102A and His102B located at the dimer interface of UreE2. This complex forms in the absence of Ni(II) and is supported by extensive protein contacts that include the use of the C-terminal sequences of UreE2 to form additional strands of β-sheet with the Ni-binding domain of HypA. The Ni-binding properties of apo, Zn-HypA•UreE2 and the component proteins were investigated by isothermal titration calorimetry using a global fitting strategy that included all of the relevant equilibria, and show that the Ni,Zn-HypA•UreE2 complex contains a single Ni(II)-binding site with a sub-nanomolar KD. The structural features of this novel Ni(II) site were elucidated using proteins produced with specifically deuterated amino acids, protein point mutations, and the analyses of X-ray absorption spectroscopy, hyperfine shifted NMR features, as well as molecular modeling coupled with quantum-mechanical calculations. The results show that the complex contains a six-coordinate, high-spin Ni(II) site with ligands provided by both component proteins.

Zambelli, Barbara (ORCID:0000000238760051)↗
DOE OSTI2021

Carbon complexes in highly C-doped GaN

We investigate the properties of heavily C-doped GaN grown by hydride vapor phase epitaxy using both optical experiments and hybrid density functional theory calculations. Previous work has established that carbon acceptors (C N ) give rise to a yellow luminescence band near 2.2 eV along with a blue luminescence band near 2.9 eV. Photoluminescence measurements show the yellow band shifting as a function of carbon concentration, suggesting a change in the behavior of carbon species as carbon content increases. With hybrid density functional theory we calculate the electrical and optical behavior of carbon centers containing multiple carbon impurities, which may arise in heavily doped material. We compare the behavior of these complexes to the isolated centers, and find that the dicarbon donor-acceptor (C Ga –C N ) complex is a candidate to explain the shift in the yellow luminescence peak. Tricarbon complexes have high formation energies and modest binding energies, and also give rise to optical transitions that are inconsistent with the observed spectra. We also identify the split dicarbon interstitial on the gallium site as a low-energy species with a large binding energy that may act to compensate carbon acceptors. Here, local vibrational modes are calculated for carbon impurity centers, and we compare these results to recent experiments. Dicarbon and tricarbon complexes involving C Ga and C N exhibit modes that are only slightly higher than the isolated species, while carbon interstitials and related complexes give rise to vibrational modes much higher than C Ga and C N .

36 MATERIALS SCIENCE↗
DOE OSTI2026

Carbon-nitrogen complex as an alternative to the 𝑇 center in Si

A first-principles study of a carbon-nitrogen (CN) impurity complex in silicon is presented as an isoelectronic alternative to the T center [(CCH) Si ]. The latter has been pursued for applications in quantum information science, yet its sensitivity to the presence of hydrogen is still problematic. Our proposed complex has no hydrogen, thereby eliminating this issue. First, we show that the CN complex is stable against decomposition into substitutional and interstitial defects. Next, we show that due to being isoelectronic to the T center, the CN complex has a similar electronic structure, and therefore could be used in similar applications. We assess several low-energy configurations of the CN complex, finding (CN) Si to be stable and have the largest Debye-Waller factor. We predict a zero-phonon line (ZPL) of 828 meV (in the telecom S-band) and a radiative lifetime of 4.2 μs, comparable to the T center. Due to the presence of a bound exciton, choice of the exchange-correlation functional and also supercell-size scaling of the ZPL and transition dipole moment require special scrutiny; we rigorously justify our extrapolation schemes that allow computing values in the dilute limit.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular Structure and Electron Affinity of Metal-Solvent Complexes: Insights from Density Functional Theory Simulations

A molecular level understanding of the structure and energetics of the monovalent and divalent metal ion complexes is of great importance for development of next-generation batteries. Here, Density Functional Theory (DFT) simulations at the ωb97xD/6-31 + G(d,p) level of theory are performed to investigate the interaction of metal ions (Li + , Na + , K + , Mg 2+ , Ca 2+ , Zn 2+ ) with 26 organic solvent molecules. The reduction energetics (electron affinity and reduction potential) and structural responses of the solvent molecules and the molecular complexes are discussed. The DFT calculations are carried out to investigate the structure, energetics and electron affinities of chelated complexes of water (H 2 O), tetrahydrofuran (THF) and di-methoxy ethane (DME) solvent molecules. Additionally, ab initio dynamic simulations (AIMD) at 298 K using atom centered density matrix propagation (ADMP) formalism are performed to understand the spontaneous structure formation upon electron attachment of the metal ion-solvent complexes. The ADMP simulations indicate the decomposition of Mg + -(DME) 3 complex via cleavage of C–O bond of one of the three DME molecules indicating irreversible decomposition of DME in the presence of the Mg + radical. We believe that the data collected as part of this investigation serves as a library of fundamental knowledge towards a deeper understanding of the electrode-electrolyte interfacial reactions.

25 ENERGY STORAGE↗
DOE OSTI2024

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Large-scale multiphysics simulations are computationally challenging due to the coupling of multiple processes with widely disparate time scales. The advent of exascale computing systems exacerbates these challenges since these systems enable ever-increasing size and complexity. In recent years, there has been renewed interest in developing multirate methods as a means to handle the large range of time scales, as these methods may afford greater accuracy and efficiency than more traditional approaches of using implicit-explicit (IMEX) and low-order operator splitting schemes. However, to date there have been few performance studies that compare different classes of multirate integrators on complex application problems. In this work, we study the performance of several newly developed multirate infinitesimal (MRI) methods, implemented in the SUNDIALS solver package, on two reacting flow model problems built on structured mesh frameworks. The first model revisits prior work on a compressible reacting flow problem with complex chemistry that is implemented using BoxLib but where we now include comparisons between a new explicit MRI scheme with the multirate spectral deferred correction (SDC) methods in the original paper. The second problem uses the same complex chemistry as the first problem, combined with a simplified flow model, but runs at a large spatial scale where explicit methods become infeasible due to stability constraints. Two recently developed IMEX MRI multirate methods are tested. These methods rely on advanced features of the AMReX framework on which the model is built, such as multilevel grids and multilevel preconditioners. The results from these two problems show that MRI multirate methods can offer significant performance benefits on complex multiphysics application problems and that these methods may be combined with advanced spatial discretization to compound the advantages of both.

97 MATHEMATICS AND COMPUTING↗